Dear David,
You are reading two dimensional latitude and longitude: lat2d(nlat,mlon) and
lon2d(nlat,mlon)
and dT(zlev,nlat,mlon)
First try to make your latitude and longitude one-dimensional
lat1d(:) = lat2d(:,1)
lon1d(:) = lon2d(1,:)
Then try,
dT!0 = "lev"
dT!1 = "lat1d"
dT!2 = "lon1d"
dT_at_lat1d = lat1d
dT_at_lon1d = lon1d
dT&lat1d_at_units="degrees_north"
dT&lon1d_at_units="degrees_east"
Thanks
Debasish
_____
From: ncl-talk-bounces_at_ucar.edu [mailto:ncl-talk-bounces_at_ucar.edu] On Behalf
Of David Small
Sent: Tuesday, February 10, 2009 7:41 PM
To: ncl-talk_at_ucar.edu
Subject: NARR mapping
Hello,
Can anyone tell me what I'm doing wrong? Thank you very much in advance.
I am having trouble plotting NARR data on a map. When I do my calculation,
the coordinate information in the data gets lost. I try to define the
dimensions, add coordinate variables but there is a problem. In the
following, I read in the 2d latitude and longitude from the NARR temperature
file, calculate the temperature difference from climatology (dT) and then
try to define the dimensions and coordinate variable. I get the following
error:
fatal:No coordinate variable exists for dimension (lat2d) in variable (dT)
fatal:(lat2d) is not coordinate variable in variable(dT).
fatal:Execute: Error occurred at or near line 58 in file TUK_T.ncl
NCL Script:
lat2d=e->lat
lon2d=e->lon
dT=Tw-Tavg
dT!0 = "lev"
dT!1 = "lat2d"
dT!2 = "lon2d"
dT_at_lat2d = lat2d
dT_at_lon2d = lon2d
dT&lat2d_at_units="degrees_north"
dT&lon2d_at_units="degrees_east"
-- David Small Graduate Research Assistant Department of Atmospheric and Oceanic Sciences McGill University Montreal, Quebec, Canada
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Feb 10 2009 - 22:36:30 MST
This archive was generated by hypermail 2.2.0 : Thu Feb 12 2009 - 10:16:44 MST