Isaac, I assume that CRTM includes temperatures up to 300K is to cover ice reaching the ground even when the surface temperature is well above 0C. In such cases the ice particles are not in equilibrium with the environment. Internally the temperature will be close to 0C. My suggestion is then to use the data for 0C also for higher temperatures.
Bye, Patrick Isaac Moradi <imor...@umd.edu> skrev: (6 juni 2021 05:25:38 CEST) >Hi Jana, > >I tried to include the za interpolation but it seems that it needs some >touch ups as there are some undefined variables. > >I did manage to do all the interpolations using the code that developed >myself and it seems to be working fine (see the attached file please). > >However, I have another issue which is unrelated to the code. In CRTM, >they have temperatures up to 300 K so I need to do extrapolation as >current SSD database doesn’t provide those temperatures. Those >temperatures may not even make sense in the real world but we don’t >intend to change them at the moment. The extrapolation creates some >negative values (please see the second attached plot). > >I wonder if any of you has any suggestion to avoid this problem? I am >not going to use those temperatures in my own work, but don’t want to >through some junks into CRTM coefficients either as someone may >eventually end up in a situation with T greater than 270 K. > > >The filenames have frequency and temperature included for reference. > >Thank you! > > >Best, >Isaac > > > > > > > > > >> On Jun 4, 2021, at 12:59 PM, Jana Mendrok <jana.mend...@gmail.com> >wrote: >> >> Hi Patrick et al, >> >> >> on the same set of grids. That is, the idea is to use these functions > >> first, before any size interpolation. This is at least indicated in >the >> documentation of assp_interp_size.m. >> >> >> No time for me to dig into the Python code (and anyhow don't have the > >> latest version at hand), and don't remember how closely Jana mimicked >my >> Matlab interface. >> >> Used to follow the Matlab implementation very closely. An >assp_interp_za exists. However, it's obviously nowhere applied (in the >original Zenodo published version; I did not check out newer ones). >> The header of assp_interp_size states the need for common grids. >Though I agree, that's not very prominently, hence easily missed by a >user and definitely missed when not called directly by the user. >Obviously, assp_interp_size misses a proper check that this requirement >is fulfilled... >> >> I guess, back in the days all size instances of one habit had the >same number of angles, so it wasn't necessary to apply that (what >speaks for that is that an IconSnow.rssp file exists; and that had been >created with exactly the code that was published, ie without >az-interpolation). >> >> Isaac/Vasileios - get_assp hence needs to be extended with a call of >assp_interp_za, in between the assp_import_ssdb and >assp.assp_interp_size (not deep enough into it anymore to judge whether >exact location matters...). >> >> In any case, Vasileios, thanks for caring for the code (and cleaning >up after me)! >> Best, >> >> -- >> Jana Mendrok, Ph.D. >> Deutscher Wetterdienst >> Offenbach am Main, Germany >> >> +49 (0)69 8062 3139 >> _______________________________________________ >> 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
