Thanks for your answer Richard,
I added abs_xsec_per_speciesAddLines into the cross-section agenda and now it
seems to work. But I do not know why lookupCalc method is working on a single
processor, even if I launched ARTS with the option -n 15.
Mattia
Da: Richard Larsson
Inviato: martedì 8 febbraio 2022 11:00
A: Mattia Sabatini
Cc: arts_users.mi@lists.uni-hamburg.de
Oggetto: 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
:
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
