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