<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<STYLE type=text/css> <!--@import url(scrollbar.css); --></STYLE>
<META content="text/html; charset=utf-8" http-equiv=Content-Type>
<STYLE> body{FONT-SIZE:12pt; FONT-FAMILY:宋体,serif;} </STYLE>
<META name=GENERATOR content="MSHTML 8.00.7601.17514"><BASE
target=_blank></HEAD>
<BODY
style="LINE-HEIGHT: 1.3; BORDER-RIGHT-WIDTH: 0px; MARGIN: 12px; BORDER-TOP-WIDTH: 0px; BORDER-BOTTOM-WIDTH: 0px; BORDER-LEFT-WIDTH: 0px"
marginheight="0" marginwidth="0">
<DIV><FONT color=#000000 size=3 face=宋体></FONT> </DIV>
<DIV> </DIV>
<DIV><FONT size=2 face=Verdana>
<HR>
</FONT></DIV>
<DIV><FONT size=2 face=Verdana><STRONG>发件人：</STRONG>Si ZH
<wh.lakes.2015@gmail.com></FONT></DIV>
<DIV><FONT size=2
face=Verdana><STRONG>发送时间：</STRONG>2015-10-25 16:36</FONT></DIV>
<DIV><FONT size=2 face=Verdana><STRONG>主题：</STRONG>about "check_for_y_lat_coord
and valid latitude coordinate array"</FONT></DIV>
<DIV><FONT size=2
face=Verdana><STRONG>收件人：</STRONG>"ncl-talk"<ncl-talk@ucar.edu></FONT></DIV>
<DIV><FONT size=2 face=Verdana><STRONG>抄送：</STRONG></FONT></DIV>
<DIV><FONT size=2 face=Verdana></FONT> </DIV>
<DIV><FONT size=2 face=Verdana>
<DIV dir=ltr>Hi there,
<DIV><BR></DIV>
<DIV>I am a new learner of UCL and was blocked at the gate of .nc file "read in"
step. I use "ncl_filedump" and can see the dimensions and variables in my .nc
file, but I don't know how to read those dimensions and variables correctly in
NCL,</DIV>
<DIV><BR></DIV>
<DIV>NCL version I am using is 6.3.0, the contents of "<A
href="http://myfile.nc">myfile.nc</A>" like below:</DIV>
<DIV><BR></DIV>
<DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>$
ncl_filedump <A href="http://myfile.nc">myfile.nc</A></FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)"
color=#000000> Copyright (C) 1995-2015 - All Rights Reserved</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)"
color=#000000> University Corporation for Atmospheric Research</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000> NCAR
Command Language Version 6.3.0</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000> The
use of this software is governed by a License Agreement.</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000> See <A
href="http://www.ncl.ucar.edu/">http://www.ncl.ucar.edu/</A> for more
details.</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)"
color=#000000><BR></FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>Variable:
f</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>Type:
file</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>filename:
myfile</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>path:
<A href="http://myfile.nc">myfile.nc</A></FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
file global attributes:</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
CDI : Climate Data Interface version 1.6.1 (<A
href="http://code.zmaw.de/projects/cdi">http://code.zmaw.de/projects/cdi</A>)</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
Conventions : CF-1.5</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
history : Thu Jan 22 19:25:20 2015: cdo sub <A
href="http://SM_MAM.zhyaa.pm.03236.206912-209911.nc">SM_MAM.zhyaa.pm.03236.206912-209911.nc</A>
<A
href="http://SM_MAM.zhyaa.03236.198012-201011.nc">SM_MAM.zhyaa.03236.198012-201011.nc</A>
<A
href="http://zhyaa.MAM.mean.diff.2080.03236.nc">zhyaa.MAM.mean.diff.2080.03236.nc</A></FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>Thu Jan 15
16:48:56 2015: cdo -timmean -selmon,3,4,5 <A
href="http://zhyaa.03236.monthly.198012-201011.nc">zhyaa.03236.monthly.198012-201011.nc</A>
<A
href="http://SM_MAM.zhyaa.03236.198012-201011.nc">SM_MAM.zhyaa.03236.198012-201011.nc</A></FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>Converted
from PP by cf-python v0.9.8.3</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
source : UM</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
runid : zhyaa</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
stash_code : 3236</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
lbproc : 0</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
submodel : 1</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
CDO : Climate Data Operators version 1.6.1 (<A
href="http://code.zmaw.de/projects/cdo">http://code.zmaw.de/projects/cdo</A>)</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
dimensions:</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
x = 134</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
y = 102</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
time = 1 // unlimited</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
nb2 = 2</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
variables:</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
float longitude ( y, x )</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
standard_name :
longitude</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
long_name : longitude</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
units :
degrees_east</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
_CoordinateAxisType : Lon</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)"
color=#000000><BR></FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
float latitude ( y, x )</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
standard_name :
latitude</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
long_name : latitude</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
units :
degrees_north</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
_CoordinateAxisType : Lat</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)"
color=#000000><BR></FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
double time ( time )</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
standard_name :
time</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
bounds : time_bnds</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
units : days since
1950-12-01 00:00:00</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
calendar : 360_day</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)"
color=#000000><BR></FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
double time_bnds ( time, nb2 )</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
units : days since
1950-12-01 00:00:00</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
calendar : 360_day</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)"
color=#000000><BR></FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
float air_temperature ( time, y, x )</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
standard_name :
air_temperature</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
long_name : TEMPERATURE AT
1.5M</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
units : K</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
coordinates : longitude latitude</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
_FillValue : -1.073742e+09</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
missing_value :
-1.073742e+09</FONT></DIV>
<DIV><FONT style="BACKGROUND-COLOR: rgb(217,234,211)" color=#000000>
cell_methods : time: mean</FONT></DIV></DIV>
<DIV><BR></DIV>
<DIV><BR></DIV>
<DIV>My NCL scripts I wrote like below:</DIV>
<DIV><BR></DIV>
<DIV><BR></DIV>
<DIV>
<DIV><FONT color=#0000ff>load
"$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"</FONT></DIV>
<DIV><FONT color=#0000ff>load
"$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"</FONT></DIV>
<DIV><FONT color=#0000ff><BR></FONT></DIV>
<DIV><FONT color=#0000ff>begin</FONT></DIV>
<DIV><FONT color=#0000ff><BR></FONT></DIV>
<DIV><FONT color=#0000ff>;;;read in temp data</FONT></DIV>
<DIV><FONT color=#0000ff>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;</FONT></DIV>
<DIV><FONT color=#0000ff><BR></FONT></DIV>
<DIV><FONT color=#0000ff>f = addfile("<A
href="http://myfile.nc">myfile.nc</A>","r")</FONT></DIV>
<DIV><FONT color=#0000ff><BR></FONT></DIV>
<DIV><FONT color=#0000ff>temp = f->air_temperature (0,:,:)</FONT></DIV>
<DIV><FONT color=#0000ff>lon = f->longitude (:,:)</FONT></DIV>
<DIV><FONT color=#0000ff>lat = f->latitude (:,:)</FONT></DIV>
<DIV><FONT color=#0000ff>timebnds = f->time_bnds</FONT></DIV>
<DIV><FONT color=#0000ff><BR></FONT></DIV>
<DIV><FONT color=#0000ff>;temp!0 = "time"</FONT></DIV>
<DIV><FONT color=#0000ff>;temp!1 = "y"</FONT></DIV>
<DIV><FONT color=#0000ff>;temp!2 = "x"</FONT></DIV>
<DIV><FONT color=#0000ff> </FONT></DIV>
<DIV><FONT color=#0000ff><BR></FONT></DIV>
<DIV><FONT color=#0000ff> ; temp&y@units = "degrees_north"</FONT></DIV>
<DIV><FONT color=#0000ff> ;temp&x@units = "degrees_east"</FONT></DIV>
<DIV><FONT color=#0000ff><BR></FONT></DIV>
<DIV><FONT color=#0000ff>;temp = temp - 273.15</FONT></DIV>
<DIV><FONT color=#0000ff>lon = lon-180</FONT></DIV>
<DIV><FONT color=#0000ff><BR></FONT></DIV>
<DIV><FONT color=#0000ff>temp@units = "(C)"</FONT></DIV>
<DIV><FONT color=#0000ff><BR></FONT></DIV>
<DIV><FONT color=#0000ff>wks = gsn_open_wks("x11", "myfile")</FONT></DIV>
<DIV><FONT color=#0000ff><BR></FONT></DIV>
<DIV><FONT color=#0000ff>wks@wkWidth = 6400</FONT></DIV>
<DIV><FONT color=#0000ff>wks@wkHeight = 4000</FONT></DIV>
<DIV><FONT color=#0000ff><BR></FONT></DIV>
<DIV><FONT color=#0000ff>res = True
;Indicate you want to set some
resources.</FONT></DIV>
<DIV><FONT color=#0000ff>res@cnFillOn = True ; Turn
on contour line fill.</FONT></DIV>
<DIV><FONT color=#0000ff>;res@cnLinesOn = False ; no
contour lines. </FONT></DIV>
<DIV><FONT color=#0000ff>res@cnLevelSpacingF = 0.5 ;
contour spacing.</FONT></DIV>
<DIV><FONT color=#0000ff>res@gsnSpreadColors = True
; full colors</FONT></DIV>
<DIV><FONT color=#0000ff>res@cnLineLabelsOn = False ;
turn off line label.</FONT></DIV>
<DIV><FONT color=#0000ff><BR></FONT></DIV>
<DIV><FONT color=#0000ff>res@cnMonoFillColor = False
; Use multiple colors.</FONT></DIV>
<DIV><FONT color=#0000ff><BR></FONT></DIV>
<DIV><FONT color=#0000ff><BR></FONT></DIV>
<DIV><FONT color=#0000ff>;res@mpMinLatF = 40</FONT></DIV>
<DIV><FONT color=#0000ff>;res@mpMaxLatF = 50</FONT></DIV>
<DIV><FONT color=#0000ff>;res@mpMinLonF = -92.5</FONT></DIV>
<DIV><FONT color=#0000ff>;res@mpMaxLonF = -75.5</FONT></DIV>
<DIV><FONT color=#0000ff><BR></FONT></DIV>
<DIV><FONT color=#0000ff><BR></FONT></DIV>
<DIV><FONT color=#0000ff>res@pmLabelBarDisplayMode = "Always"</FONT></DIV>
<DIV><FONT color=#0000ff><BR></FONT></DIV>
<DIV><FONT color=#0000ff>res@pmLabelBarOrthogonalPosF = 0.08</FONT></DIV>
<DIV><FONT color=#0000ff><BR></FONT></DIV>
<DIV><FONT color=#0000ff>;res@tiXAxisString =
lon@long_name</FONT></DIV>
<DIV><FONT color=#0000ff>;res@tiYAxisString =
lat@long_name</FONT></DIV>
<DIV><FONT color=#0000ff><BR></FONT></DIV>
<DIV><FONT color=#0000ff>;res@sfXArray =
lon</FONT></DIV>
<DIV><FONT color=#0000ff>;res@sfYArray =
lat</FONT></DIV>
<DIV><FONT color=#0000ff><BR></FONT></DIV>
<DIV><FONT color=#0000ff><BR></FONT></DIV>
<DIV><FONT color=#0000ff>;res@lbAutoLabelStride = True</FONT></DIV>
<DIV><FONT color=#0000ff><BR></FONT></DIV>
<DIV><FONT color=#0000ff>plot = gsn_csm_contour_map (wks,temp,res)</FONT></DIV>
<DIV><FONT color=#0000ff><BR></FONT></DIV>
<DIV><FONT color=#0000ff><BR></FONT></DIV>
<DIV><FONT color=#0000ff>end</FONT><BR></DIV></DIV>
<DIV><BR></DIV>
<DIV><BR></DIV>
<DIV>and NCL error message like below:</DIV>
<DIV><BR></DIV>
<DIV><BR></DIV>
<DIV>
<DIV><SPAN style="BACKGROUND-COLOR: rgb(0,255,255)">(0)
check_for_y_lat_coord: Warning: Data either does not contain a valid latitude
coordinate array or doesn't contain one at all.</SPAN></DIV>
<DIV><SPAN style="BACKGROUND-COLOR: rgb(0,255,255)">(0) A valid
latitude coordinate array should have a 'units' attribute equal to one of the
following values:</SPAN></DIV>
<DIV><SPAN style="BACKGROUND-COLOR: rgb(0,255,255)">(0)
'degrees_north' 'degrees-north' 'degree_north' 'degrees north'
'degrees_N' 'Degrees_north' 'degree_N' 'degreeN' 'degreesN' 'deg
north'</SPAN></DIV>
<DIV><SPAN style="BACKGROUND-COLOR: rgb(0,255,255)">(0)
check_for_lon_coord: Warning: Data either does not contain a valid longitude
coordinate array or doesn't contain one at all.</SPAN></DIV>
<DIV><SPAN style="BACKGROUND-COLOR: rgb(0,255,255)">(0) A valid
longitude coordinate array should have a 'units' attribute equal to one of the
following values:</SPAN></DIV>
<DIV><SPAN style="BACKGROUND-COLOR: rgb(0,255,255)">(0)
'degrees_east' 'degrees-east' 'degree_east' 'degrees east' 'degrees_E'
'Degrees_east' 'degree_E' 'degreeE' 'degreesE' 'deg east'</SPAN></DIV></DIV>
<DIV><BR></DIV>
<DIV><BR></DIV>
<DIV>The data in "<A href="http://myfile.nc">myfile.nc</A>" are North American
Great Lakes area temperature contour data (the large boundary is middle and East
of US and Canada, and I have narrowed the map latitude and longitude to Great
Lakes area, but the results looks nothing about my data, it looks like all over
the world or East Asian area, it is why I used lon = lon - 180). My object is to
figure out Great lakes Area temperature contour.</DIV>
<DIV><BR></DIV>
<DIV>I think my error is like NCL message said (I didn't read in correct
dimensions and coordinates), but I can't 100% understand "<A
href="http://myfile.nc">myfile.nc</A>" via NCL mini manuals, and don't know how
to write tha correct scripts that can read in "<A
href="http://myfile.nc">myfile.nc</A>" (with dimensions: y, x, time and
time_bnd)</DIV>
<DIV><BR></DIV>
<DIV>I attach "<A href="http://myfile.nc">myfile.nc</A>" inside the attachment,
hope you can help this,</DIV>
<DIV><BR></DIV>
<DIV>Thank you,</DIV>
<DIV><BR></DIV>
<DIV>Liang Zhang</DIV>
<DIV><BR></DIV>
<DIV><BR></DIV>
<DIV> </DIV>
<DIV><BR></DIV></DIV></FONT></DIV></BODY></HTML>