Thank you, I've solved this problem, it is because of the wrong absorption 
models I've choosed. thanks again.



At 2020-02-18 21:33:51, "Lukas Kluft" <lukas.kl...@gmail.com> wrote:

Hello,


Your control file is not running because the `abs_species` workspace variable 
is not set.
Could you please check if you attached the correct control file?


Best wishes,
Lukas



On 18. Feb 2020, at 02:19, Rex <rex...@126.com> wrote:


Dear,
I'm a new user of ARTS, recently I'm using the newest ARTS 2.3 and radiosonde 
data to do a simulation for a groud-based microwave radiometer, the question 
is, the simulated BT have a significant bias compared with observed one at 
water vapor channel, and i don't know how to correct it, is there someone could 
help me?
the following is my script. thanks.


Arts2 {
    INCLUDE "general/general.arts"
    INCLUDE "general/continua.arts"
    INCLUDE "general/agendas.arts"
    INCLUDE "general/planet_earth.arts"


    Copy(abs_xsec_agenda, abs_xsec_agenda__noCIA)
    Copy(iy_main_agenda, iy_main_agenda__Emission)
    Copy(iy_space_agenda, iy_space_agenda__CosmicBackground)
    Copy(iy_surface_agenda, iy_surface_agenda__UseSurfaceRtprop)


    Copy(propmat_clearsky_agenda, propmat_clearsky_agenda__OnTheFly)
    Copy(ppath_agenda, ppath_agenda__FollowSensorLosPath)
    
    Copy(ppath_step_agenda, ppath_step_agenda__GeometricPath)


    IndexSet(stokes_dim, 1)


    jacobianOff
    cloudboxOff
    ReadXML(p_grid, "sonde.p.xml")    
    
    ReadSplitARTSCAT(basename="spectroscopy/Perrin/", fmin=22e9, fmax=58e9)


    abs_linesSetLineShapeType(option="HTP")
    abs_linesSetNormalization(option="None") # this value affect brightness temp
    abs_lines_per_speciesCreateFromLines
    
    AtmRawRead(basename="./AtmRaw/sonde")


    VectorSetConstant(surface_scalar_reflectivity, 1, 0.01)
    Copy(surface_rtprop_agenda, 
surface_rtprop_agenda__Specular_NoPol_ReflFix_SurfTFromt_surface)


    # HTG3 Center Frequency point
    VectorSet(f_grid, [22.24e9, 23.04e9, 23.84e9, 25.44e9, 26.24e9, 27.84e9, 
31.40e9, 51.26e9, 52.28e9, 53.86e9, 54.94e9, 56.66e9, 57.30e9, 58.00e9])


    sensorOff


    StringSet(iy_unit, "PlanckBT")


    AtmosphereSet1D
    AtmFieldsCalc( vmr_zeropadding=1 )
    Extract(z_surface, z_field, 0)
    Extract(t_surface, t_field, 0)


    Copy(sensor_pos, z_surface)
    MatrixSet(sensor_los, [0])


    abs_xsec_agenda_checkedCalc
    lbl_checkedCalc
    propmat_clearsky_agenda_checkedCalc
    atmfields_checkedCalc
    atmgeom_checkedCalc
    cloudbox_checkedCalc
    sensor_checkedCalc
    yCalc


    WriteXML("ascii", y, "results/bt.xml")
    WriteXML("ascii", f_grid, "results/f_grid.xml")
    WriteXML("ascii", abs_species, "results/species.xml")
    WriteXML("ascii", sensor_pos, "results/sensor_pos.xml")
}




<btsim.arts>_______________________________________________
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

Reply via email to