flex_extract WRF data

[anphi]

Message by ahilboll:

I noticed that when downloading data with GAUSS=1 and ETA=0 (which is, if I understand the comments in flexpartTools.py correctly, the recommended setting for ERA-Interim), the U and V parameters are downloaded in SH grid, not on GG - which means that it probably cannot be ingested by WRF's ungrib.exe. Or am I missing something here?

Furthermore, D is contained in the WRF files twice, once in GG and once in SH. That's because it explicitly gets added to the GG list in the

if c.wrf=='1':

block even though it's already part of SH.

So maybe, the way would be to do

                if c.wrf=='1':
                    self.params['OG__ML'][0]+='/Z/VO'
                    for var_ in ['/D', '/U', '/V']:
                        if var_ not in self.params['OG__ML'][0]:
                            self.params['OG__ML'][0]+=var_

to add the GG U and V fields to the download list, and then (~line 1188)

                    if c.wrf=='1' and gridtype=='regular_ll':                   
                        if levtype=='hybrid':                                   
                            if paramId in [129,130,131,132,133,138,155]:        
                                grib_write(gid,fwrf)                            
                        else:                                                   
                            if paramId in wrfpars:                              
                                grib_write(gid,fwrf)                            

[anphi]

I will consider this in version 7.1!

[anphi]

This issue is NOT considered in version 7.1 and is still to be done in a future version.