[ncl-talk] Plotting Meteogram From WRF Data?
Ronald Stenz
rds238 at cornell.edu
Mon Aug 6 20:12:17 MDT 2018
Great, thanks for the help! I'll give that a try!
On Mon, Aug 6, 2018 at 5:05 PM, Dennis Shea <shea at ucar.edu> wrote:
> Actually, I created some good' netCDF files and a script
>
> %> gzip -d meteo*nc.gz
> %> ncl meteo_3.ncl
>
> WRF data has different dimensionality
>
> A3(Time, south_north, west_east)
> A4(Time, bottom_top, south_north, west_east)
>
> The meteogram script wants 'Time' to be the 'x-axis' and 'Pressure' to be
> the y-axis.
>
> a4 = A4( bottom_top|:, Time|:, south_north|j , west_east|i)
> printVarSummary(a4) ; a4(bottom_top, Time) ===> a4(Y,X)
>
> You would have to interpolate to pressure.
>
> Good luck
>
> On Sun, Aug 5, 2018 at 10:01 PM, Dennis Shea <shea at ucar.edu> wrote:
>
>> *re: "how is this file created"*
>>
>> The file and script were donated 15 years ago.
>>
>> [snp]
>> // global attributes:
>> :title = "*Meteogram data from John Ertl*" ;
>> :date = "Tue Sep 14 11:22:51 MDT *2004*" ;
>>
>> The file is a poorly created netCDF file.
>>
>> ====
>> re: "*I would like to make those meteograms from my simulations, the
>> format of meteo_data.nc <http://meteo_data.nc> seems different?*
>>
>> Not sure what you mean by "different". Here is a translation of the file.
>>
>> netcdf meteo_data.new {
>> dimensions:
>> levels = 8 ;
>> time = 25 ;
>> variables:
>> float rh(levels, time) ;
>> rh:_FillValue = -999.f ;
>> rh:long_name = "relative humidity" ;
>> rh:units = "%" ;
>> float levels(levels) ;
>> levels:_FillValue = -999.f ;
>> levels:units = "hPa" ;
>> levels:long_name = "pressure" ;
>> float time(time) ;
>> time:_FillValue = -999.f ;
>> time:units = "hours" ;
>> float ugrid(levels, time) ;
>> ugrid:_FillValue = -999.f ;
>> ugrid:long_name = "zonal wind" ;
>> ugrid:units = "???" ;
>> float vgrid(levels, time) ;
>> vgrid:_FillValue = -999.f ;
>> vgrid:long_name = "meridional wind" ;
>> vgrid:units = "???" ;
>> float tempisobar(levels, time) ;
>> tempisobar:_FillValue = -999.f ;
>> tempisobar:long_name = "temperature" ;
>> tempisobar:units = "degF" ;
>> float rain03(time) ;
>> rain03:_FillValue = -999.f ;
>> rain03:long_name = "*3hr rain total*" ;
>> rain03:units = "inches" ;
>> float tempht(time) ;
>> tempht:_FillValue = -999.f ;
>> tempht:long_name = "temperature at 2m" ;
>> tempht:units = "degF" ;
>>
>> =====
>> If you have (say)
>>
>> RH(time,levels)
>>
>> your can reorder using named dimension reordering.
>>
>> rh = RH(levels|:,time|:)
>> printVarSummary(rh)
>>
>> =====
>>
>> If you want to basically copy that code you must cahnge the units, etc
>>
>> =======
>>
>> Please respond *only* to nal-talk,
>>
>> Good luck
>>
>> On Sun, Aug 5, 2018 at 8:31 PM, Ronald Stenz <rds238 at cornell.edu> wrote:
>>
>>> Thanks, but how is this file created? I would like to make those
>>> meteograms from my simulations, the format of meteo_data.nc seems
>>> different?
>>>
>>> Thanks
>>>
>>> On Sun, Aug 5, 2018 at 9:24 PM, Dennis Shea <shea at ucar.edu> wrote:
>>>
>>>> The file *meteo_data.nc <http://meteo_data.nc>* is distributed with
>>>> NCL.
>>>>
>>>> %> ncdump -h $NCARG_ROOT/lib/ncarg/data/cdf/meteo_data.nc
>>>>
>>>> or
>>>>
>>>> %> ncl_filedump -h $NCARG_ROOT/lib/ncarg/data/cdf/meteo_data.nc
>>>>
>>>>
>>>> On Sun, Aug 5, 2018 at 5:33 PM, Ronald Stenz <rds238 at cornell.edu>
>>>> wrote:
>>>>
>>>>> I'm looking at meteo_1.ncl and meteo_2.ncl and am wondering where the
>>>>> file meteo_data.nc comes from? Is there a separate program to
>>>>> generate this file from the wrfout files? The wrfout files will not work
>>>>> for these programs, but do in meteo_5.ncl (which uses wrfout files).
>>>>>
>>>>> https://www.ncl.ucar.edu/Applications/Scripts/meteo_2.ncl
>>>>>
>>>>> Here is the code I am referring to for reference.
>>>>>
>>>>> Thanks
>>>>>
>>>>> _______________________________________________
>>>>> ncl-talk mailing list
>>>>> ncl-talk at ucar.edu
>>>>> 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/20180806/ceabf15c/attachment.html>
More information about the ncl-talk
mailing list