Hi,
Jonas: Yes, I confirm that retrieved winds in qpack2_wind3d_demo are
coming out as zero. I have run this demo several times, but failed to
notice this. So this error has probably been around for quite some time!
Sorry.
Richard: I have confirmed that spectra with zero and non-zero winds
differ, but that wind jacobians seem to be zero. The cfile produced by
qpack attached. Can you have a look? Have I missed to adjust qpack in
any way?
No more time today to look at this. I will try to look at this further
tomorrow. Just wanted to give a quick response.
Bye,
Patrick
On 2016-11-04 11:08, Richard Larsson wrote:
All,
This might be my fault. I have no idea how this particular qpack-case
works but there was an overhaul in the design of the Jacobians in ARTS
last year that made their calculation less forced and more theoretical
(i.e., the jacobian before was from partial perturbations but it is now
from direct derivation as much as possible). For the wind field this
means that the Jacobian is calculated based on the derivative of the
transmission with respect to frequency before being scaled back to wind
field parameters.
There could be edge-cases that I have miss-implemented. It sounds like
this at the least.
Could you explain your test case to me so that I can look into fixing
the problem? What is the atmospheric geometry and what is the a priori
wind field? What kind of results do you expect?
//Richard
2016-11-04 18:01 GMT+09:00 Jonas Hagen <[email protected]
<mailto:[email protected]>>:
Dear all
I'm currently porting my wind retrievals from ARTS/atmlab v2.2 to
v2.3 (developement version). But somehow the wind retrievals just
return zeros (wind_{u|v}_x and wind_{u|v}_A in L2 are exactly 0).
This is also the case for the qpack2_wind3d_demo.m. Can anyone
confirm this? I could need some help with the debugging. There are
no error messages or warnings.
Best regards,
Jonas
--
Jonas Hagen
Institute of Applied Physics
University of Bern
Sidlerstr. 5, CH-3012 Bern, Switzerland
_______________________________________________
qpack.mi mailing list
[email protected] <mailto:[email protected]>
https://mailman.rrz.uni-hamburg.de/mailman/listinfo/qpack.mi
<https://mailman.rrz.uni-hamburg.de/mailman/listinfo/qpack.mi>
_______________________________________________
qpack.mi mailing list
[email protected]
https://mailman.rrz.uni-hamburg.de/mailman/listinfo/qpack.mi
# Control file created by Atmlab function *qarts2cfile*
Arts2{
INCLUDE "/home/patrick/SVN/ARTS/arts/controlfiles/general/general.arts"
INCLUDE "/home/patrick/SVN/ARTS/arts/controlfiles/general/agendas.arts"
INCLUDE "/home/patrick/SVN/ARTS/arts/controlfiles/general/continua.arts"
INCLUDE "/home/patrick/SVN/ARTS/arts/controlfiles/general/planet_earth.arts"
IndexSet(abs_f_interp_order,3)
IndexSet(stokes_dim,1)
ReadXML(f_grid,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/f_grid.xml")
AtmosphereSet3D
abs_speciesSet(species=["O3","H2O"])
ReadXML(p_grid,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/p_grid.xml")
ReadXML(lat_grid,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/lat_grid.xml")
ReadXML(lon_grid,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/lon_grid.xml")
ReadXML(lat_true,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/lat_true.xml")
ReadXML(lon_true,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/lon_true.xml")
ReadXML(t_field,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/t_field.xml")
ReadXML(vmr_field,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/vmr_field.xml")
ReadXML(z_field,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/z_field.xml")
ReadXML(wind_u_field,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/wind_u_field.xml")
ReadXML(wind_v_field,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/wind_v_field.xml")
ReadXML(z_surface,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/z_surface.xml")
atmfields_checkedCalc(bad_partition_functions_ok=1)
atmgeom_checkedCalc
NumericSet(p_hse,10000.000)
NumericSet(z_hse_accuracy,0.100)
z_fieldFromHSE
abs_linesReadFromArts(abs_lines,"/home/patrick/SVN/ARTS/atmlab/demos/data/o3line111ghz",0.000e+00,1.000e+99)
abs_lines_per_speciesCreateFromLines
AgendaSet(abs_xsec_agenda){
abs_xsec_per_speciesInit
abs_xsec_per_speciesAddLines
abs_xsec_per_speciesAddConts
}
abs_xsec_agenda_checkedCalc
AgendaSet(propmat_clearsky_agenda){
Ignore(rtp_mag)
Ignore(rtp_los)
propmat_clearskyInit
propmat_clearskyAddOnTheFly
}
propmat_clearsky_agenda_checkedCalc
nlteOff
Copy(iy_surface_agenda,iy_surface_agenda__UseSurfaceRtprop)
Copy(iy_space_agenda,iy_space_agenda__CosmicBackground)
NumericSet( ppath_lmax, 2.50e+02 )
Copy(ppath_agenda,ppath_agenda__FollowSensorLosPath)
Copy(ppath_step_agenda,ppath_step_agenda__GeometricPath)
ReadXML(sensor_pos,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/sensor_pos.xml")
ReadXML(sensor_los,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/sensor_los.xml")
ReadXML(sensor_time,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/sensor_time.xml")
Copy(iy_main_agenda,iy_main_agenda__Emission)
StringSet(iy_unit,"RJBT")
WriteXML(output_file_format,sensor_pos,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/sensor_pos.xml")
WriteXML(output_file_format,sensor_los,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/sensor_los.xml")
ReadXML(sensor_response,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/sensor_response.xml")
ReadXML(antenna_dim,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/antenna_dim.xml")
ReadXML(mblock_dlos_grid,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/mblock_dlos_grid.xml")
ReadXML(sensor_response_f,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/sensor_response_f.xml")
ReadXML(sensor_response_dlos,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/sensor_response_dlos.xml")
ReadXML(sensor_response_pol,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/sensor_response_pol.xml")
ReadXML(sensor_response_f_grid,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/sensor_response_f_grid.xml")
ReadXML(sensor_response_dlos_grid,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/sensor_response_dlos_grid.xml")
ReadXML(sensor_response_pol_grid,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/sensor_response_pol_grid.xml")
jacobianInit
VectorCreate(gas1_jac_p_grid)
ReadXML(gas1_jac_p_grid,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/gas1_jac_p_grid.xml")
VectorCreate(gas1_jac_lat_grid)
ReadXML(gas1_jac_lat_grid,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/gas1_jac_lat_grid.xml")
VectorCreate(gas1_jac_lon_grid)
ReadXML(gas1_jac_lon_grid,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/gas1_jac_lon_grid.xml")
jacobianAddAbsSpecies(jacobian_quantities,jacobian_agenda,
atmosphere_dim,p_grid,lat_grid,lon_grid,
gas1_jac_p_grid,gas1_jac_lat_grid,gas1_jac_lon_grid,
"O3",
"analytical","rel",1,0)
VectorCreate(wind_u_jac_p_grid)
ReadXML(wind_u_jac_p_grid,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/wind_u_jac_p_grid.xml")
VectorCreate(wind_u_jac_lat_grid)
ReadXML(wind_u_jac_lat_grid,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/wind_u_jac_lat_grid.xml")
VectorCreate(wind_u_jac_lon_grid)
ReadXML(wind_u_jac_lon_grid,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/wind_u_jac_lon_grid.xml")
jacobianAddWind(jacobian_quantities,jacobian_agenda,
atmosphere_dim,p_grid,lat_grid,lon_grid,
wind_u_jac_p_grid,wind_u_jac_lat_grid,wind_u_jac_lon_grid,"u")
VectorCreate(wind_v_jac_p_grid)
ReadXML(wind_v_jac_p_grid,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/wind_v_jac_p_grid.xml")
VectorCreate(wind_v_jac_lat_grid)
ReadXML(wind_v_jac_lat_grid,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/wind_v_jac_lat_grid.xml")
VectorCreate(wind_v_jac_lon_grid)
ReadXML(wind_v_jac_lon_grid,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/wind_v_jac_lon_grid.xml")
jacobianAddWind(jacobian_quantities,jacobian_agenda,
atmosphere_dim,p_grid,lat_grid,lon_grid,
wind_v_jac_p_grid,wind_v_jac_lat_grid,wind_v_jac_lon_grid,"v")
jacobianAddFreqShift(df=2.500000e+04)
jacobianAddPolyfit(jacobian_quantities, jacobian_agenda,
sensor_response_pol_grid,sensor_response_dlos_grid,
sensor_pos,0)
jacobianClose
cloudboxOff
cloudbox_checkedCalc
sensor_checkedCalc
yCalc
WriteXML(output_file_format,y,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/y.xml")
WriteXML(output_file_format,y_f,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/y_f.xml")
WriteXML(output_file_format,y_pol,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/y_pol.xml")
WriteXML(output_file_format,y_pos,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/y_pos.xml")
WriteXML(output_file_format,y_los,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/y_los.xml")
WriteXML(output_file_format,y_geo,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/y_geo.xml")
WriteXML(output_file_format,y_aux,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/y_aux.xml")
WriteXML(output_file_format,jacobian,"/tmp/atmlab-patrick-tpf7474803_8cdf_41c5_9b8c_50381c07e496/jacobian.xml")
}
_______________________________________________
qpack.mi mailing list
[email protected]
https://mailman.rrz.uni-hamburg.de/mailman/listinfo/qpack.mi
