[arts-users] ARTS OEM

2018-09-13 Thread 이병석
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

2018-09-13 Thread Stefan Buehler
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

2018-09-13 Thread Guillermo Gimenez de Castro
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