Dear Patrick and Oliver,
Thank you for your advice! Indeed, I had used the 2.2.64 stable version
assuming it was a safe choice for a beginner. Now I tried Patrick's piece of
code for FASTEM call with the current ARTS version from Oliver's GitHub
(2.3.1234, I assume), and everything worked flawlessly.
I appreciate your help a lot.
Best regards,
Natalia
Natalia Bliankinshtein
Research Council Officer/Agent du Conseil de recherches
NRC Aerospace/CNRC Aérospatiale
National Research Council Canada/Conceil National de recherches Canada
1200 Montreal Rd, Ottawa (Ontario) K1A 0R6
Bldg/Edifice U-61 | Rm/Ch. 233A
1920 Research Private
Tel.: (613) 998-5349 | Cellphone: (343) 549-4262
natalia.bliankinsht...@nrc-cnrc.gc.ca
________________________________________
From: Patrick Eriksson <patrick.eriks...@chalmers.se>
Sent: March 7, 2019 2:21 PM
To: Bliankinshtein, Natalia; arts_users.mi@lists.uni-hamburg.de
Subject: Re: [arts-users] FASTEM in ARTS
Dear Natalia,
Try this:
VectorCreate( trv )
NumericCreate( wind_speed )
NumericCreate( wind_dir )
AgendaSet( iy_surface_agenda ){
SurfaceDummy
specular_losCalc
ppathCalc( rte_pos=rtp_pos, rte_los=specular_los )
ArrayOfStringSet( iy_aux_vars, ["Optical depth"] )
iyEmissionStandard
transmittanceFromIy_aux( transmittance=trv )
surfaceFastem( wind_speed=wind_speed, wind_direction=wind_dir,
transmittance=trv )
iySurfaceRtpropCalc
}
This should work. I have not actually made a test run (due to lack of
time) but hopefully I have not missed any detail. (If you use this, you
can ignore surface_rtprop_agenda).
FASTEM does not fit into the pattern of ARTS, as it wants the
transmission as input. This makes the agenda definition a bit messy.
FASTEM 6 is default.
My suggestion uses the default for salinity, adopt if you want to use
another value.
You need a quite recent version of ARTS 2.3 for this. I added the method
transmittanceFromIy_aux not a long time ago.
Bye,
Patrick
On 2019-03-07 19:22, Bliankinshtein, Natalia wrote:
> Dear Patrick,
>
> Thanks for your prompt reply! I would like to use FASTEM in a simple 1D
> forward model configuration, so no retrievals of surface properties will be
> needed. I am particularly interested in FASTEM6 - not sure how its versions
> align with those of ARTS. Also, I think a simple prescribing of surface wind,
> temperature and salinity as Numerics would be sufficient for my purposes.
>
> I have started building my controlfile upon the classroom exercise #2
> (rtcalc) included in the release package, so if you could show me how to
> include FASTEM into that example, that would be just great.
> Please let me know if you need more info about the planned usage to give me a
> hint on this.
>
> Thank you for your help,
> Natalia
>
>
>
> Natalia Bliankinshtein
>
> Research Council Officer/Agent du Conseil de recherches
> NRC Aerospace/CNRC Aérospatiale
> National Research Council Canada/Conceil National de recherches Canada
> 1200 Montreal Rd, Ottawa (Ontario) K1A 0R6
> Bldg/Edifice U-61 | Rm/Ch. 233A
> 1920 Research Private
> Tel.: (613) 998-5349 | Cellphone: (343) 549-4262
> natalia.bliankinsht...@nrc-cnrc.gc.ca
>
> ________________________________________
> From: Patrick Eriksson <patrick.eriks...@chalmers.se>
> Sent: March 7, 2019 2:43 AM
> To: Bliankinshtein, Natalia; arts_users.mi@lists.uni-hamburg.de
> Subject: Re: [arts-users] FASTEM in ARTS
>
> Dear Natalia,
>
> Yes, FASTEM is at hand, but it is a bit tricky to give a general answer.
> It depends on e.g.
>
> * Do you want to retrieve any surface properties?
>
> * Do you want to set specific wind speeds etc, or do you want to obtain
> wind speed etc. by interpolation from some data?
>
> Can you give us some more information?
>
> Bye,
>
> Patrick
>
>
>
> On 2019-03-06 20:51, Bliankinshtein, Natalia wrote:
>> Dear ARTS users and developers,
>>
>>
>> I have recently started learning ARTS and appreciate its extreme
>> flexibility.
>>
>> According to the documentation, ARTS does not include an option of using
>> FASTEM6 water surface model, however I see some evidence people have
>> used it, for example here in Mr. Oliver Lemke's repository
>> https://github.com/olemke/arts/tree/master/3rdparty/fastem and in Mr.
>> Stuart Fox's email to the list in November 2018.
>>
>>
>> I am struggling though to come up with a controlfile that would call
>> FASTEM6. So I am kindly asking if anyone could please share a minimum
>> working example of how to do that?
>>
>>
>> Also, a workaround that could work for my purposes would be reading
>> surface reflectivities from a file, which I tried. It seems to me,
>> however, that surface emissivities are always computed as
>> (1-reflectivity), which is not the case for FASTEM. Thus my question: is
>> it possible to override surface emissivity and to read it from a file
>> instead?
>>
>>
>> Looking forward to hear your valuable advice.
>>
>>
>> Best regards,
>>
>> Natalia
>>
>>
>>
>> *Natalia Bliankinshtein*
>>
>> Research Council Officer/Agent du Conseil de recherches
>> NRC Aerospace/CNRC Aérospatiale
>> National Research Council Canada/Conceil National de recherches Canada
>> 1200 Montreal Rd, Ottawa (Ontario) K1A 0R6
>> Bldg/Edifice U-61 | Rm/Ch. 233A
>> 1920 Research Private
>> Tel.: (613) 998-5349 | Cellphone: (343) 549-4262
>> natalia.bliankinsht...@nrc-cnrc.gc.ca
>>
>>
>>
>> _______________________________________________
>> 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
