From: Debasish Pai Mazumder <debasish_at_nyahnyahspammersnyahnyah>

Date: Tue, 10 Feb 2009 20:36:30 -0900

Date: Tue, 10 Feb 2009 20:36:30 -0900

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
*