[arts-users] ARTS OEM
Dear ARTS users, I am writing to address a question with the use of OEM in the ARTS development version. From the "TestOEM.arts" in "ARTS/controlfiles/artscomponents/oem", this time I changed the codes as follows. 1) From (line 115): retrievalAddAbsSpecies( species = "O3", unit = "vmr", g1 = p_ret_grid, g2 = lat_grid, g3 = lon_grid ) To: retrievalAddTemperature ( g1 = p_ret_grid, g2 = lat_grid, g3 = lon_grid, hse = "off" ) 2) From (line 125): VectorSetConstant(vars, nelem, 1e-12) DiagonalMatrix(sparse_block, vars) covmat_sxAddBlock(block = sparse_block) To: VectorSetConstant(vars, nelem, 0.1) DiagonalMatrix(sparse_block, vars) covmat_sxAddBlock(block = sparse_block) 3) From (line 37): NumericSet( f_start, 110.436e9 ) NumericSet( f_end, 111.236e9 ) IndexSet( nf, 801 ) IndexSet( np, 81 ) To: NumericSet( f_start, 110.436e9 ) NumericSet( f_end, 111.236e9 ) IndexSet( nf, 18 ) IndexSet( np, 18 ) 4) From (line 42): VectorNLinSpace( f_grid, nf, f_start, f_end ) VectorNLogSpace( p_grid,361, 500e2, 0.1 ) VectorNLogSpace( p_ret_grid, np, 500e2, 0.1 ) To: VectorNLinSpace( f_grid, nf, f_start, f_end ) VectorNLogSpace( p_grid,18, 500e2, 0.1 ) VectorNLogSpace( p_ret_grid, np, 500e2, 0.1 ) Overall, I changed the ozone VMR retrieval to temperature retrieval and reduced the number of frequencies and pressure layers. When I run the controlfile multiple times, I irregularly get errors. The error messages are Run-time error in oem computation: Forward Model Evaluation Error Run-time error in agenda: inversion_iterate_agenda Run-time error in method: atmfields_checkedCalc The variable *t_field* contains one or several NaNs. This is not allowed! or All temperatures in *t_field* must be > 0. or seldom the bad partition function warning. Here are my questions: 1) Is this an anticipated error, simply coming from failing to find a converging solution? Is it expected that the OEM sometimes throws an error as above? 2) Are there mechanisms (e.g. workspace methods) within the ARTS that can prevent the t_field or vmr_field from becoming either negative or NaN during the OEM iterations? Thank you for your assistance. Regards, Byungsuk Lee ARA Consulting & Technology 30 Songdomirae-ro D-1510, Yeonsu-gu, Incheon 21990, South Korea Email: cb...@aracnt.com ___ 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] Recommendation
Dear Guigue, yes, ARTS should be perfectly suitable for that task. Just use the HITRAN catalogue and the CKD_MT continua. (You could look at the HIRS example controlfiles for spectroscopic settings.) There are different ways to get out absorption coefficients, or perhaps even optical depth directly. I’m not sure which way is best for you, it depends a bit on what you want to do with it. Best wishes, Stefan > On 13. Sep 2018, at 15:05, Guillermo Gimenez de Castro > wrote: > > Dear everybody, > > I need to compute atmospheric optical depths from submillimeter to wavelengths > near 10 um. Is it arts suitable for the task? What would you recommend > instead? > > Thank you very much for your help, > > Guigue > > -- > > Si la tolérance naît du doute, qu'on enseigne > à douter des modèles et des utopies, > à récuser les prophètes de salut, > les annonciateurs de catastrophes. > > Appelons de nos voeux la venue des sceptiques > s'ils doivent éteindre le fanatisme. > > Raymond Aron > > ~ > Guillermo Giménez de Castro > Centro de Rádio-Astronomia e Astrofísica Mackenzie (CRAAM) > > Rua da Consolação 896Tel: ++55(11) 2114-8717 > Predio M. Carvalhosa - Cobertura Fax: ++55(11) 3255-3123 > 01302-907, S.Paulo, SP http://www.guigue.gcastro.net > Brasil > > > AVISO LEGAL: Esta mensagem corporativa é destinada somente a quem ela é > dirigida e pode conter informação confidencial e legalmente protegida. Caso > você não seja o destinatário, fica notificado de que é ilegal examiná-la, > utilizá-la, divulgá-la, copiá-la ou distribuí-la no todo ou em partes. Se a > recebeu por engano, pedimos que a retorne, apagando-a dos seus registros. > Fica desprovida de validade a mensagem emitida por quem não detenha poderes > de representação. > DISCLAIMER: This corporate message is intended solely for the addressees and > is confidential. It may contain information that is privileged or otherwise > protected from disclosure. If you have received this transmission in error, > please delete it and immediately notify the sender. Any use not in accord > with its purpose, any dissemination or disclosure, either whole or partial, > is prohibited unless formal approval is granted. > ___ > arts_users.mi mailing list > arts_users.mi@lists.uni-hamburg.de > https://mailman.rrz.uni-hamburg.de/mailman/listinfo/arts_users.mi smime.p7s Description: S/MIME cryptographic signature ___ 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] Recommendation
Dear everybody, I need to compute atmospheric optical depths from submillimeter to wavelengths near 10 um. Is it arts suitable for the task? What would you recommend instead? Thank you very much for your help, Guigue -- Si la tolérance naît du doute, qu'on enseigne à douter des modèles et des utopies, à récuser les prophètes de salut, les annonciateurs de catastrophes. Appelons de nos voeux la venue des sceptiques s'ils doivent éteindre le fanatisme. Raymond Aron ~ Guillermo Giménez de Castro Centro de Rádio-Astronomia e Astrofísica Mackenzie (CRAAM) Rua da Consolação 896Tel: ++55(11) 2114-8717 Predio M. Carvalhosa - Cobertura Fax: ++55(11) 3255-3123 01302-907, S.Paulo, SP http://www.guigue.gcastro.net Brasil AVISO LEGAL: Esta mensagem corporativa é destinada somente a quem ela é dirigida e pode conter informação confidencial e legalmente protegida. Caso você não seja o destinatário, fica notificado de que é ilegal examiná-la, utilizá-la, divulgá-la, copiá-la ou distribuí-la no todo ou em partes. Se a recebeu por engano, pedimos que a retorne, apagando-a dos seus registros. Fica desprovida de validade a mensagem emitida por quem não detenha poderes de representação. DISCLAIMER: This corporate message is intended solely for the addressees and is confidential. It may contain information that is privileged or otherwise protected from disclosure. If you have received this transmission in error, please delete it and immediately notify the sender. Any use not in accord with its purpose, any dissemination or disclosure, either whole or partial, is prohibited unless formal approval is granted. ___ arts_users.mi mailing list arts_users.mi@lists.uni-hamburg.de https://mailman.rrz.uni-hamburg.de/mailman/listinfo/arts_users.mi