Re: Fixed grid from ncep map plotting

I see what you did [at least, I think I see what you did].

It speculate you used 'ncl_convert2nc' to
create the "flx*nc" files. You selected "PRATE_98_SFC"
but you did not specify the the corresponding lat_98/lon_98
coordinate information should also be copied.

No problem.

http://www.ncl.ucar.edu/Document/Functions/Contributed/latGau.shtml
http://www.ncl.ucar.edu/Document/Functions/Contributed/lonGlobeF.shtml

can be used to generate the correct lat/lon coordinates

    f = addfile ("test.nc", "r")
    p = f->PRATE_98_SFC_ave ; (Time, lat_98,lon_98)

    dimp = dimsizes( p )
    ntim = dimp(0)
    nlat = dimp(1)
    mlon = dimp(2)

    lon_98 = lonGlobeF(mlon, "lon_98", "longitude", "degrees_east")
    lat_98 = latGau(nlat, "lat_98", "latitude", "degrees_north") ; S=>N

If data are N=>S add the line
    lat_98 = lat_98(::-1) ; make N=>S

    p&lon_98 = lon98
    p&lat_98 = lat98

    printVarSummary( p )

Correia, James wrote:
> By directly, you just use gsn_csm_contour_map_ce ?
> The data doesn't seem to line up very well with the map background.
> Do I need to set gsnAddCyclic = False and also set tfDoNDCOverlay = False?
> Thanks for any help.
> jimmyc
>
> Below is the netcdf header of the data:
> netcdf test {
> dimensions:
> Time = UNLIMITED ; // (62 currently)
> lat_98 = 94 ;
> lon_98 = 192 ;
> variables:
> float PRATE_98_SFC_ave(Time, lat_98, lon_98) ;
> PRATE_98_SFC_ave:forecast_time_units = "hours" ;
> PRATE_98_SFC_ave:forecast_time = 0 ;
> PRATE_98_SFC_ave:parameter_number = 59 ;
> PRATE_98_SFC_ave:parameter_table_version = 2 ;
> PRATE_98_SFC_ave:grid_number = 98 ;
> PRATE_98_SFC_ave:level_indicator = 1 ;
> PRATE_98_SFC_ave:_FillValue = 1.e+20f ;
> PRATE_98_SFC_ave:units = "kg/m^2/s" ;
> PRATE_98_SFC_ave:long_name = "Precipitation rate" ;
> PRATE_98_SFC_ave:center = "US National Weather Service - NCEP (WMC)"
> ;
>
> // global attributes:
> :creation_date = "Thu Jul 30 15:04:18 PDT 2009" ;
> :NCL_Version = "5.1.0" ;
> :system = "Linux nwiceb 2.6.18-92.1.13.el5lustre #1 SMP Fri Sep 26
> 09:45:16 PDT 2008 x86_64 x86_64 x86_64 GNU/Linux" ;
> :conventions = "None" ;
> :grib_source = "flxf00.grb" ;
> :title = "NCL: convert-GRIB-to-netCDF" ;
> :history = "Thu Jul 30 15:09:41 2009: ncrcat flxf00.nc flxf06.nc
> flxf12.nc flxf18.nc flxf24.nc flxf30.nc flxf36.nc flxf42.nc flxf48.nc
> flxf54.nc flxf60.nc flxf66.nc flxf72.nc flxf78.nc flxf84.nc flxf90.nc
> flxf96.nc flxf102.nc flxf108.nc flxf114.nc flxf120.nc flxf126.nc flxf132.nc
> flxf138.nc flxf144.nc flxf150.nc flxf156.nc flxf162.nc flxf168.nc flxf174.nc
> flxf180.nc flxf186.nc flxf192.nc flxf198.nc flxf204.nc flxf210.nc flxf216.nc
> flxf222.nc flxf228.nc flxf234.nc flxf240.nc flxf246.nc flxf252.nc flxf258.nc
> flxf264.nc flxf270.nc flxf276.nc flxf282.nc flxf288.nc flxf294.nc flxf288.nc
> flxf300.nc flxf306.nc flxf312.nc flxf318.nc flxf324.nc flxf330.nc flxf336.nc
> flxf342.nc flxf348.nc flxf354.nc flxf360.nc test.nc" ;
> :nco_openmp_thread_number = 1 ;
>
> On 7/30/09 4:03 PM, "Dennis Shea" <shea_at_ucar.edu> wrote:
>
>> It is my feeling that you should be able to plot the fields directly.
>>
>> On import, the variables have the coordinate information
>> appropriate for the variable. The 94x192 is a standard gaussian grid.
>> Sould be trivial.
>>
>> D
>>
>> Correia, James wrote:
>>> Hi All-
>>> I am using ncl_convert2nc on a number of grib files from NCEPs CFS model. I
>>> am not sure how to associate the data with the correct mapping. The data
>>> appear to be of a similar form of the reanalysis data (2.5 by 2.5 degrees)
>>> with 73 lat and 144 lon points for the pgb files, however the flx files have
>>> a different lat lon grid (grid number = 98) with 94 lat and 192 lon.
>>>
>>> Anyone out there have any experience with this dataset to help me plot it
>>> correctly on a map?
>>> Thanks
>>>
>>>
>>>
>>> James Correia Jr., PhD
>>> Climate Physics Group
>>> Post. Doc.
>>> Pacific Northwest National Lab
>>>
>>> "Wisdom. Strength. Courage. Generosity. Each of us are born with one of
>>> these. It is up to us to find the other three inside of us."
>>> -Into the West
>>>
>>>
>>> _______________________________________________
>>> ncl-talk mailing list
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Jul 31 2009 - 09:45:48 MDT