Re: [arts-users] ARTS 2.5.0 - absorption
Dear Mattia, It looks like you are removing line calculations from your lookup calculations. Could you try adding abs_xsec_per_speciesAddLines to your lookup cross-section agenda? I think this will solve the problem. We are currently in a transition in the 2.5-branch of moving away from abs_xsec_agenda entirely. There are some lingering problems when using lookup table calculations at this time. Mainly, since the line calculations should happen in propmat_clearsky_agenda now, they are not part of any of the default cross-section agendas as in the past. The lookup generation interface is currently under active development, so this interface might change in the near future. With hope, //Richard Den fre 21 jan. 2022 kl 14:53 skrev Mattia Sabatini < mattia.sabat...@artov.ismar.cnr.it>: > Dear all, > > > > I have a question concerning the absorption calculation with ARTS 2.5.0. > Here are some details about the controlfile I wrote: > > > > 1D atmosphere, clear sky, 4 IR channels > > > > VectorLinSpace( f_grid, 2.02456915e13, 4.2827494e13, 5e8 ) – i.e. from 7um > to 14.8 um > > > > abs_species = "H2O, H2O-SelfContCKDMT100, H2O-ForeignContCKDMT100", "O3”, > "CO2, CO2-CKDMT100", "N2O", "CH4", "CFC11-HXSEC", "CFC12-HXSEC", > "HCFC22-HXSEC", "CFC113-HXSEC", "CFC114-HXSEC". > > > > abs_linesReadSpeciesSplitCatalog(abs_lines, "cat/") - from > arts-xml-data-2.4/spectroscopy > > ReadXsecData(basename="coefficients_arts/") – from arts-crossfit-main/ > > > > On-the-fly absorption set with > > > > AgendaSet( propmat_clearsky_agenda ){ > > Ignore(rtp_mag) > > Ignore(rtp_los) > > Ignore(rtp_nlte) > > propmat_clearskyInit > > propmat_clearskyAddXsecAgenda > > propmat_clearskyAddLines > > propmat_clearskyAddHitranXsec > > } > > > > AgendaSet( abs_xsec_agenda ){ > > abs_xsec_per_speciesInit > > abs_xsec_per_speciesAddConts > > } > > > > Batch_atm_fields_compact from Eresmaa 137L, gases VMR from FASCOD except > for CFC ones which are set as constant. > > > > Surface set with > > > > AgendaSet( surface_rtprop_agenda ){ > >specular_losCalc > >Extract( surface_skin_t, t_surface_vector, ybatch_index ) > >surfaceFlatScalarReflectivity > > } > > (Surface_vector contains the surface skin temperatures from the Eresmaa > surface data. The correct skin temperature is extracted during ybatchCalc, > accordingly to ybatch_index) > > > > This is my output for the first profile (channel BTs in Kelvin): > > 291.688825481684 > > 288.903990987693 > > 287.623301104291 > > 284.402270055496 > > > > I repeated the calculation by using a lookup table, editing the > controlfile: > > > > AgendaSet( propmat_clearsky_agenda ){ > > Ignore(rtp_mag) > > Ignore(rtp_los) > > Ignore(rtp_nlte) > > propmat_clearskyInit > > propmat_clearskyAddFromLookup > > propmat_clearskyAddHitranXsec > > } > > > > abs_lookupSetupBatch > > abs_xsec_agenda_checkedCalc > > lbl_checkedCalc > > abs_lookupCalc > > WriteXML( "binary", abs_lookup ) > > ReadXML( abs_lookup ) > > abs_lookupAdapt > > > > The lookup calculation is done in few minutes, and the output I get this > time is: > > > > 293.350716564285 > > 291.452571755867 > > 296.771709918824 > > 292.567643470471 > > > > Why are the two outputs for the same profile so different? Also, why is > the lookup table calculation so fast compared to version 2.4 (a couple of > hours with same f_grid and atmospheric profiles)? Where did I go wrong? > > > > Mattia > ___ > 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] read in HITRAN par file
Dear Pengwang, You should be able to read a Hitran file in 2.4 that you could read in 2.3. Now, it sounds like you have downloaded a new version of the Hitran catalog. Arts 2.4 does not directly support Hitran2020, neither did 2.3 to be clear. For water, I still think this should still be Ok but I am not sure. If you want to give it a try, the 'modern' way of dealing with Hitran data in Arts is to download it with the additional tags in the Hitran download interface. On the tab currently called '4. Select or Create Output Format'*, you can create a new output format. To create a new format that works well with Arts, this format should contain at the top the ".par line" tag, as the second option it should contain the " qns' " tag, and as a last option it should have the " qns'' " tag. The field separator and line endings should be the default ([tab] ; (CR LF)). If you manage to download data like this from Hitran, you can give the hitran_type="Online" option in ReadHITRAN. This should give decent results. Again, remember that Arts 2.4 is from before Hitran2020 was released. It is possible this will not work in 2.4. What I describe above works for the development version of Arts (v.2.5). With hope, //Richard * From hitran.org, Data Access -> Line-By-Line -> select H2O and go forward -> go forward (or deselect some isotopologues first) -> go forward (or select some Kayser range first) -> You are now on the right tab to specify output format Den tors 3 feb. 2022 kl 16:06 skrev Pengwang Zhai : > Hello, ARTS Community, > > I used arts 2.3 to generate gas absorption coefficient lookup table for > o2, co2, etc. One key function was: "abs_linesReadFromHitran”. Now I need > to revisit the calculation with HITRAN2020. I thought I would use arts 2.4, > a newer version. After successfully compiling arts2.4, and run the same > .arts file I used before (a modification of TestAbs.arts), arts tells me > that “abs_linesReadFromHitran” is not available any more. I went over the > documentation, and found “READHITRAN”. I tried the following: > > ReadHITRAN( abs_lines, >“./h2o_61fbd3ce.par", > 1.3e+14, > 1.0e+15 ) > > arts reports: > > "Run-time error in method: ReadHITRAN > Error parsing quantum number Sym”. > > Is there any way we can read in HITRAN par file and create a lookup table? > Do you suggest me going back to arts 2.3? > > Thanks for your advice, > > Pengwang > > > Pengwang Zhai (he/his/him) > Graduate Program Director, Atmospheric Physics, > Associate Professor, Physics Department, UMBC > pwz...@umbc.edu > > > > > > > ___ > 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] cannot compile the latest version of ARTS
Hello, I had some issue with arts 2.4 in reading HITRAN par file. I was advised to use the latest version of arts. I did the following: git clone https://github.com/atmtools/arts.git cd arts mkdir build cd build cmake -DENABLE_NETCDF=1 .. make The compilation stoped at 25%, with the following error message: /arts/src/nc_io.h:37:10: fatal error: 'abs_species_tags.h' file not found #include "abs_species_tags.h" ^~~~ 1 error generated. make[2]: *** [src/CMakeFiles/artscore.dir/auto_md.cc.o] Error 1 make[1]: *** [src/CMakeFiles/artscore.dir/all] Error 2 make: *** [all] Error 2 Any suggestion here on how to fix this? Note that my purpose is to create an absorption lookup table from HITRAN 2020. I need the NetCDF support to make the resultant absorption lookup table manageable. Appreciate it very much, Yours, Pengwang ___ 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] cannot compile the latest version of ARTS
Hi Pengwang, Thanks for reporting this bug. I'll merge a PR with a fix soon. See https://github.com/atmtools/arts/pull/394 However, if you don't plan to use or modify the lookup tables outside of ARTS, I recommend that you save them as binary XML instead for best efficiency: WriteXML( "binary", abs_lookup, "lookup.xml" ) Cheers, Oliver > On 8 Feb 2022, at 20:02, Pengwang Zhai wrote: > > Hello, > > I had some issue with arts 2.4 in reading HITRAN par file. I was advised to > use the latest version of arts. I did the following: > > git clone https://github.com/atmtools/arts.git > cd arts > mkdir build > cd build > cmake -DENABLE_NETCDF=1 .. > make > > The compilation stoped at 25%, with the following error message: > > /arts/src/nc_io.h:37:10: fatal error: 'abs_species_tags.h' file not found > #include "abs_species_tags.h" > ^~~~ > 1 error generated. > make[2]: *** [src/CMakeFiles/artscore.dir/auto_md.cc.o] Error 1 > make[1]: *** [src/CMakeFiles/artscore.dir/all] Error 2 > make: *** [all] Error 2 > > Any suggestion here on how to fix this? > > Note that my purpose is to create an absorption lookup table from HITRAN > 2020. I need the NetCDF support to make the resultant absorption lookup table > manageable. > > Appreciate it very much, > > Yours, > > Pengwang > > > > > > > ___ > 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