Meteoinfolab Script Example: Data projection-flexpart

Flexpart is a hysplit-like diffusion pattern, its output netcdf file refers to the WRF, but the global properties are not written, such as only a projection name (such as Lambert), there is no relevant projection parameters: central longitude, standard latitude and so on. The WRF header file must be consulted to reconstruct the projection (why not tiger the property all over?). ）。 The latitude and longitude coordinates of the data are some, but the coordinates under the Lambert projection are not, the x, y coordinates of the projection can be obtained by the PROJECTXY function, and the LAT is the latitude and longitude of the lower left corner of the data in Lon.

Script Program:

f = AddFile ('D:/temp/nc/header_d01.nc') Topov= f['Topography']topo=topov[:,:]#Set ProjectionFromproj = Projinfo (proj='Longlat') Toproj= Projinfo (proj='LCC', lon_0=117.0092, lat_0=40.11041, lat_1=30.0, lat_2=60.0)#Set x/yLon = 112.14053344726562lat= 35.62180709838867xn= 216yn= 240DX= 4000Dy= 4000x, y=Projectxy (lon, lat, xn, yn, dx, dy, toproj, fromproj)#PlotAxesm (Projinfo=toproj, griddx=2, griddy=2, gridline=True) Mlayer= Shaperead ('d:/temp/map/bou2_4p.shp') geoshow (mlayer) layer= IMSHOWM (x, Y, topo, proj=toproj) title ('Flexpart Topography') Colorbar (layer) show ()

Meteoinfolab Script Example: Data projection-flexpart