WRFCHEM output python
2014-09-29
###WRFCHEM_output pythonic####
- Related to the last note on getting txt time series data from wrf output here. This note details getting numpy array from wrf output and plotting data in python.
- The variables required from wrfoutput are Temperature, Relative humidity, Wind direction, wind speed, PM2.5 and PM 10. Numpy array from these variables will be used for air earth application, to insert the grid data into Postgresql. Data plots is for a quick view of the wrf output.
- The numpy array creation is based on this, this and this.
- The plotting is primarily based on this and this.
- The numpy array creation and plotting for each wrf output variables is attempted as follows ###Temperature###
- From this, it is required to post-process wrfoutput using
p_interp
utility to get pressure correction for variables such as RH. - For this, the source code was downloaded and compiled using the command
gfortran p_interp.F90 -L/usr/local/lib -lnetcdf -lnetcdff -lm \ -I/usr/local/include -Mfree -o p_interp
, with this p_interp got compiled and an executable of the same name is created.