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
