[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