Re: map projections

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Fri, 14 Dec 2007 15:13:46 -0700

Remove:
tfDoNDCOverlay = True

Also .... did you unpack the precip data?
 It is of type short?

     prc = short2flt( f->precip )
     printVarSummary(prc)

jimmyc_at_iastate.edu wrote:
> All-
> I am trying to plot the Unified precipitation dataset from CDC. It is
> a netcdf file but has no information on the map projection. The data
> is on a lat lon grid:
> netcdf precip.1993 {
> dimensions:
> lon = 321 ;
> lat = 161 ;
> time = UNLIMITED ; // (365 currently)
> variables:
> float lat(lat) ;
> lat:units = "degrees_north" ;
> lat:long_name = "Latitude" ;
> lat:actual_range = 20.f, 60.f ;
> float lon(lon) ;
> lon:units = "degrees_east" ;
> lon:long_name = "Longitude" ;
> lon:actual_range = 220.f, 300.f ;
> double time(time) ;
> time:units = "hours since 1700-1-1 00:00:0.0" ;
> time:long_name = "Time" ;
> time:actual_range = 2568384., 2577120. ;
> time:delta_t = "0000-00-01 00:00:00" ;
> time:avg_period = "0000-00-01 00:00:00" ;
> short precip(time, lat, lon) ;
> precip:long_name = "Daily Accumulated Precipitation" ;
>
>
> When I attempt to plot it using the
> mpProjection = "LambertConformal"
> mpLimitMode = "Corners"
> mpLeftCornerLatF = lat(0)
> mpLeftCornerLonF = lon(0)
> mpRightCornerLatF = lat(160)
> mpRightCornerLonF = lon(320)
>
> The image that is produced has the map background rotated 90 degrees
> but the displayed data is in the correct orientation.
>
> I do have this:
> tfDoNDCOverlay = True
>
> I have no idea if this is a Lambert grid but other projection types
> didnt seem to work either.
>
> Any ideas?
>

-- 
======================================================
Dennis J. Shea                  tel: 303-497-1361    |
P.O. Box 3000                   fax: 303-497-1333    |
Climate Analysis Section                             |
Climate & Global Dynamics Div.                       |
National Center for Atmospheric Research             |
Boulder, CO  80307                                   |
USA                        email: shea 'at' ucar.edu |
======================================================
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Dec 14 2007 - 15:13:46 MST

This archive was generated by hypermail 2.2.0 : Mon Dec 31 2007 - 09:18:01 MST