[ncl-talk] DSSAT input preparation from WRF output
Mary Haley
haley at ucar.edu
Sat Apr 18 12:51:19 MDT 2015
I don't know anything about the DSSAT format, but to see how to read WRF
variables, visit the NCL/WRF pages:
http://www.ncl.ucar.edu/Applications/wrf.shtml
Any one of these links shows various ways for reading data from a WRF file.
As a quick example, you can either read a variable directly:
a = addfile("wrfout_xxxxx.nc","r")
rainc = a->RAINC
"rainc" will contain all the data and the metadata.
or you can use the "wrf_user_getvar" function, which calculates specific
diagnostics by reading required WRF variables and doing the specified
calculation.
http://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_getvar.shtml
If you want to stay in the Python world, then you can use PyNIO to read a
WRF output variable. The syntax is different than NCL:
import Nio as nio
a = nio.open_file("wrfout_xxxx.nc")
rainc = a.variables['RAINC'] # this is a NioVariable with metadataT[:]
to access values
# rainc = a.variables['RAINC'][:] # this is a numpy array with the actual
values.
You can learn more about PyNIO at http://www.pyngl.ucar.edu/Nio.shtml
There may be some WRF PyNIO/PyNGL examples at:
http://www.pyngl.ucar.edu/gallery.shtml
<http://www.pyngl.ucar.edu/Nio.shtml>
I can't check at the moment because the NCAR websites are down due to
semi-annual maintenance.
--Mary
On Mon, Apr 13, 2015 at 5:15 AM, Raja Sivaranjan <rajasivaranjan92 at gmail.com
> wrote:
> Dear all,
>
> I need to prepare weather file for DSSAT model from wrf output. Rainfall,
> Min and Max temperature, solar radiation and weed speed are variable i need
> to extract. I'm using python for handling netcdf file which very difficult
> to work with big file. I did some coding to convert wrf output to DSSAT
> weather file format. Now i want to learn NCL, since it is more powerful to
> work with nc files.
>
> Below are the things i need to do in NCL:
>
> 1. Extract variables such as Rainfall, Min and Max Temperature, solar
> radiation and weed speed for required time
> 2. Write the extracted variable in DSSAT format.
>
> Things i have done in python:
>
> 1. convert wrfnc to excel
> 2. excel to Dssat format.
>
> Can anyone help me out in this.
>
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20150418/3d0c884e/attachment.html
More information about the ncl-talk
mailing list