Re: [arts-users] On the changes around iy-methods

2018-12-19 Thread Patrick Eriksson 

Hi Victoria,

(I wrote this yesterday, but seems that the sending failed! Anyhow, here 
is an answer complementing Oliver's reply.)


Not sure what you mean with retrieve. If you just want to export the 
data to files, easiest is to do the storing inside the agenda. For example:


AgendaSet( iy_main_agenda ){
  ppathCalc
  iyEmissionStandard
  WriteXML( "binary", ppvar_t, "ppvar_t.xml" )
}

This works for a single beam calculation. If you calculate several 
pencil beams in one ARTS call, there is nowadays a way to give each file 
an unique number. You then write to file like this (still inside the 
agenda):


WriteXMLIndexed("binary",iy_id,ppvar_t,"ppvar_t.xml",9)

Do

arts -d iy_id

for information on the numbering scheme.

Bye,

Patrick


On 2018-12-18 18:30, Victoria Sol Galligani wrote:
Hi! I'm updating old codes to use the latest version of arts and I am 
having some trouble. I'm trying to retrieve pressure, temperature, and 
other variables along the ppath. I am trying to use iyEmissionStandard 
the following way:


Copy(iy_main_agenda, iy_main_agenda__Emission)

IndexSet(iy_agenda_call1, 1)
Tensor3SetConstant( iy_transmission, 0, 0, 0, 0.0 )

ppathCalc

iyEmissionStandard( iy, iy_aux, diy_dx, ppvar_p, ppvar_t, ppvar_nlte, 
ppvar_vmr,   ppvar_wind,

     ppvar_mag, ppvar_f, ppvar_iy,
     ppvar_trans_cumulat, ppvar_trans_partial, iy_id,
     stokes_dim, f_grid, atmosphere_dim, p_grid,
     z_field, t_field, nlte_field, vmr_field,
     abs_species, wind_u_field, wind_v_field,
     wind_w_field, mag_u_field, mag_v_field,
     mag_w_field, cloudbox_on, iy_unit, iy_aux_vars,
     jacobian_do, jacobian_quantities, ppath,
     rte_pos2, propmat_clearsky_agenda,
     water_p_eq_agenda, iy_main_agenda,
     iy_space_agenda, iy_surface_agenda,
     iy_cloudbox_agenda, iy_agenda_call1,
     iy_transmission, rte_alonglos_v,
     surface_props_data )


and I can't seem to be able to ignore diy_dx

What is the correct way to retrieve for example ppvar_p, ppvar_t and 
ppvar_vmr


Thank you in advance!!!

Victoria






___
arts_users.mi mailing list
arts_users.mi@lists.uni-hamburg.de
https://mailman.rrz.uni-hamburg.de/mailman/listinfo/arts_users.mi


___
arts_users.mi mailing list
arts_users.mi@lists.uni-hamburg.de
https://mailman.rrz.uni-hamburg.de/mailman/listinfo/arts_users.mi

Re: [arts-users] On the changes around iy-methods

2018-12-18 Thread Oliver Lemke 
Hi Victoria,

Based on the small extract of your controlfile, it seems you're calling 
iyEmissionStandard directly? Not sure that's gonna work because 
iyEmissionStandard is supposed to be called only through the iy_main_agenda 
which is executed by yCalc. yCalc then takes care of the initialization of 
diy_dx.

However, you can try to do a fake initialization of diy_dx by putting 
'Touch(diy_dx)' before your iyEmissionStandard call.

Cheers,
Oliver


> On 18 Dec 2018, at 18:30, Victoria Sol Galligani 
>  wrote:
> 
> Hi! I'm updating old codes to use the latest version of arts and I am having 
> some trouble. I'm trying to retrieve pressure, temperature, and other 
> variables along the ppath. I am trying to use iyEmissionStandard the 
> following way: 
> 
> Copy(iy_main_agenda, iy_main_agenda__Emission)
> 
> IndexSet(iy_agenda_call1, 1)
> Tensor3SetConstant( iy_transmission, 0, 0, 0, 0.0 )
> 
> ppathCalc
> 
> iyEmissionStandard( iy, iy_aux, diy_dx, ppvar_p, ppvar_t, ppvar_nlte, 
> ppvar_vmr,   ppvar_wind,
> ppvar_mag, ppvar_f, ppvar_iy,
> ppvar_trans_cumulat, ppvar_trans_partial, iy_id,
> stokes_dim, f_grid, atmosphere_dim, p_grid,
> z_field, t_field, nlte_field, vmr_field,
> abs_species, wind_u_field, wind_v_field,
> wind_w_field, mag_u_field, mag_v_field,
> mag_w_field, cloudbox_on, iy_unit, iy_aux_vars,
> jacobian_do, jacobian_quantities, ppath,
> rte_pos2, propmat_clearsky_agenda,
> water_p_eq_agenda, iy_main_agenda,
> iy_space_agenda, iy_surface_agenda,
> iy_cloudbox_agenda, iy_agenda_call1,
> iy_transmission, rte_alonglos_v,
> surface_props_data )
> 
> 
> and I can't seem to be able to ignore diy_dx
> 
> What is the correct way to retrieve for example ppvar_p, ppvar_t and ppvar_vmr
> 
> Thank you in advance!!! 
> 
> Victoria
> 
> 
> 
> 
> 
> ___
> arts_users.mi mailing list
> arts_users.mi@lists.uni-hamburg.de
> https://mailman.rrz.uni-hamburg.de/mailman/listinfo/arts_users.mi

___
arts_users.mi mailing list
arts_users.mi@lists.uni-hamburg.de
https://mailman.rrz.uni-hamburg.de/mailman/listinfo/arts_users.mi

[arts-users] On the changes around iy-methods

2018-12-18 Thread Victoria Sol Galligani 
Hi! I'm updating old codes to use the latest version of arts and I am
having some trouble. I'm trying to retrieve pressure, temperature, and
other variables along the ppath. I am trying to use iyEmissionStandard the
following way:

Copy(iy_main_agenda, iy_main_agenda__Emission)

IndexSet(iy_agenda_call1, 1)
Tensor3SetConstant( iy_transmission, 0, 0, 0, 0.0 )

ppathCalc

iyEmissionStandard( iy, iy_aux, diy_dx, ppvar_p, ppvar_t, ppvar_nlte,
ppvar_vmr,   ppvar_wind,
ppvar_mag, ppvar_f, ppvar_iy,
ppvar_trans_cumulat, ppvar_trans_partial, iy_id,
stokes_dim, f_grid, atmosphere_dim, p_grid,
z_field, t_field, nlte_field, vmr_field,
abs_species, wind_u_field, wind_v_field,
wind_w_field, mag_u_field, mag_v_field,
mag_w_field, cloudbox_on, iy_unit, iy_aux_vars,
jacobian_do, jacobian_quantities, ppath,
rte_pos2, propmat_clearsky_agenda,
water_p_eq_agenda, iy_main_agenda,
iy_space_agenda, iy_surface_agenda,
iy_cloudbox_agenda, iy_agenda_call1,
iy_transmission, rte_alonglos_v,
surface_props_data )


and I can't seem to be able to ignore diy_dx

What is the correct way to retrieve for example ppvar_p, ppvar_t and
ppvar_vmr

Thank you in advance!!!

Victoria
___
arts_users.mi mailing list
arts_users.mi@lists.uni-hamburg.de
https://mailman.rrz.uni-hamburg.de/mailman/listinfo/arts_users.mi