[ncl-talk] plotting model output where x/y are not lat/lon

[mberdahl at envsci.rutgers.edu]

Hi all,

I'm trying to make a simple plot of cumulative snowfall with data from a
model I did not run.  When I try to make a simple plot, I receive the
following error:

(0)     check_for_y_lat_coord: Warning: Data either does not contain a
valid latitude coordinate array or doesn't contain one at all.
(0)     A valid latitude coordinate array should have a 'units' attribute
equal to one of the following values:
(0)         'degrees_north' 'degrees-north' 'degree_north' 'degrees north'
'degrees_N' 'Degrees_north' 'degree_N' 'degreeN' 'degreesN' 'deg north'
(0)     check_for_lon_coord: Warning: Data either does not contain a valid
longitude coordinate array or doesn't contain one at all.
(0)     A valid longitude coordinate array should have a 'units' attribute
equal to one of the following values:
(0)         'degrees_east' 'degrees-east' 'degree_east' 'degrees east'
'degrees_E' 'Degrees_east' 'degree_E' 'degreeE' 'degreesE' 'deg east'


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


When I dump the file it looks like this:



netcdf ICE.1983.01-12.f77 {
dimensions:
        TIME = UNLIMITED ; // (365 currently)
        X10_69 = 60 ;
        Y18_127 = 110 ;
        SECTOR = 2 ;
        SECTOR1_1 = 1 ;
        ATMLAY14_23 = 10 ;
        bnds = 2 ;
        ATMLAY23_23 = 1 ;
        ATM6H3_3 = 1 ;
        OUTLAY = 20 ;
        OUTLAY1_1 = 1 ;

variables:
        double TIME(TIME) ;
                TIME:units = "HOURS since 1901-01-15 00:00:00" ;
                TIME:long_name = "time" ;
                TIME:time_origin = "15-JAN-1901 00:00:00" ;
                TIME:axis = "T" ;
        float DATE(TIME) ;
                DATE:units = "YYYYMMDDHH" ;
                DATE:long_name = "Date(YYYY MM DD HH)" ;
                DATE:history = "From ICE.f77.1983.01.01-15" ;
        double X10_69(X10_69) ;
                X10_69:units = "km" ;
                X10_69:long_name = "x" ;
                X10_69:point_spacing = "even" ;
                X10_69:axis = "X" ;

 double Y18_127(Y18_127) ;
                Y18_127:units = "km" ;
                Y18_127:long_name = "y" ;
                Y18_127:point_spacing = "even" ;
                Y18_127:axis = "Y" ;
        float LON(Y18_127, X10_69) ;
                LON:units = "degrees" ;
                LON:long_name = "Longitude" ;
                LON:history = "From ICE.f77.1983.01.01-15" ;
        float LAT(Y18_127, X10_69) ;
                LAT:units = "degrees" ;
                LAT:long_name = "Latitude" ;
                LAT:history = "From ICE.f77.1983.01.01-15" ;


Clearly there are lats and lons associated with the x and y spacing, but
I'm not sure how to tell NCL how to use those.  Does anyone have advice on
this?

Thanks in advance,
Mira