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

Reply via email to