<!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> BLOCKQUOTE{margin-Top: 0px; margin-Bottom: 0px; margin-Left: 2em} body{FONT-SIZE:12.1pt; COLOR:#001; FONT-FAMILY:宋体,serif;} </STYLE>
<META name=GENERATOR content="MSHTML 10.00.9200.16660"><BASE
target=_blank></HEAD>
<BODY
style="BORDER-LEFT-WIDTH: 0px; BORDER-RIGHT-WIDTH: 0px; BORDER-BOTTOM-WIDTH: 0px; MARGIN: 12px; LINE-HEIGHT: 1.3; BORDER-TOP-WIDTH: 0px"
marginheight="0" marginwidth="0"><STATIONERY>
<DIV><FONT color=#000000 size=3 face=宋体>Hi All,</FONT></DIV>
<DIV><FONT color=#000000 size=3> I am trying to plot the geostrophic
current fields from aviso. And get some screen output and error messages as
follow,</FONT></DIV>
<DIV><FONT color=#000000 size=3></FONT> </DIV>
<DIV>(0) 19930101</DIV>
<DIV> </DIV>
<DIV>Variable: u<BR>Type: float<BR>Total Size: 664692
bytes<BR>
166173 values<BR>Number of Dimensions: 2<BR>Dimensions and sizes:
[NbLongitudes | 579] x [NbLatitudes |
287]<BR>Coordinates:<BR>
NbLongitudes:
[98.33333333332349..290.9999999999709]<BR>
NbLatitudes: [-17.79766126545592..60.88112837177402]<BR>Number Of Attributes:
5<BR> date : 1993-01-01
00:00:00.000000 UTC<BR> Date_CNES_JD
: 15706<BR> units
: cm/s<BR> long_name :
U<BR> _FillValue : 1.844674e+19</DIV>
<DIV> </DIV>
<DIV>Variable: v<BR>Type: float<BR>Total Size: 664692
bytes<BR>
166173 values<BR>Number of Dimensions: 2<BR>Dimensions and sizes:
[NbLongitudes | 579] x [NbLatitudes |
287]<BR>Coordinates:<BR>
NbLongitudes:
[98.33333333332349..290.9999999999709]<BR>
NbLatitudes: [-17.79766126545592..60.88112837177402]<BR>Number Of Attributes:
5<BR> _FillValue : 1.844674e+19<BR> long_name :
V<BR> units : cm/s<BR>
Date_CNES_JD : 15706<BR> date
: 1993-01-01 00:00:00.000000 UTC</DIV>
<DIV> </DIV>
<DIV>Variable: speed<BR>Type: float<BR>Total Size: 664692
bytes<BR>
166173 values<BR>Number of Dimensions: 2<BR>Dimensions and sizes:
[579] x [287]<BR>Coordinates:<BR>Number Of Attributes: 1<BR> _FillValue
: 1.844674e+19<BR>(0) check_for_y_lat_coord:
Warning: Data either does not contain a valid latitude coordinate array or
doesn't contain one at all.<BR>(0) A valid latitude
coordinate array should have a 'units' attribute equal to one of the following
values:<BR>(0) 'degrees_north'
'degrees-north' 'degree_north' 'degrees north' 'degrees_N' 'Degrees_north'
'degree_N' 'degreeN' 'degreesN' 'deg north'<BR>(0)
check_for_lon_coord: Warning: Data either does not contain a valid longitude
coordinate array or doesn't contain one at all.<BR>(0) A
valid longitude coordinate array should have a 'units' attribute equal to one of
the following values:<BR>(0)
'degrees_east' 'degrees-east' 'degree_east' 'degrees east' 'degrees_E'
'Degrees_east' 'degree_E' 'degreeE' 'degreesE' 'deg
east'<BR>(0) gsn_csm_map_ce: Warning: you set mpMaxLonF
to a value > 180,
but<BR>(0)
didn't set mpCenterLonF. Setting mpCenterLonF to 194.5<BR></DIV>
<DIV>My code is as follow,</DIV>
<DIV> </DIV>
<DIV> <FONT color=#0000ff> </FONT><FONT color=#0000ff
face="Times New Roman">do nt=nb,ne<BR>
in =
addfile(infile(nt-1),"r")<BR> utc_date =
str_get_cols(infile(nt-1),44,51)<BR>
print(""+utc_date)<BR>
<BR>;************************************************<BR>; read in zonal [u] and
meridional [v] winds
<BR>;************************************************<BR>
tmp = in->Grid_0001<BR>
u = tmp({lon_min:lon_max},{lat_min:lat_max})
;/100.0<BR> tmp =
in->Grid_0002<BR> v =
tmp({lon_min:lon_max},{lat_min:lat_max}) ;/100.0<BR>
speed = sqrt(u^2 +
v^2)<BR> <BR> </FONT><A
href="mailto:vcres@tiMainString"><FONT color=#0000ff
face="Times New Roman">vcres@tiMainString</FONT></A><FONT color=#0000ff
face="Times New Roman"> =
"Geostrophic Current Field"<BR> </FONT><A
href="mailto:vcres@gsnLeftString"><FONT color=#0000ff
face="Times New Roman">vcres@gsnLeftString</FONT></A><FONT color=#0000ff
face="Times New Roman"> =
utc_date<BR> </FONT><A
href="mailto:vcres@gsnRightString"><FONT color=#0000ff
face="Times New Roman">vcres@gsnRightString</FONT></A><FONT color=#0000ff
face="Times New Roman"> =
"m/s"</FONT></DIV>
<DIV><FONT color=#0000ff face="Times New Roman"></FONT> </DIV>
<DIV><FONT color=#ffff00 face="Times New Roman"><FONT
color=#0000ff> wks = gsn_open_wks(wks_type
,"uv_"+utc_date)
; open a ps file<BR>
gsn_define_colormap(wks,"rainbow")
;"BlAqGrYeOrRe") ; choose
colormap<BR>
printVarSummary(u)<BR>
printVarSummary(v)<BR>
printVarSummary(speed)<BR>
plot=gsn_csm_vector_scalar_map_ce(wks,u,v,speed,vcres) ; create
plot<BR> end do</FONT><BR></FONT></DIV>
<DIV><FONT face="Times New Roman"><FONT color=#000000>Anyone can help me for
figuring out the error? Thanks.</FONT></FONT></DIV>
<DIV><FONT face="Times New Roman"><FONT color=#000000></FONT> </DIV></FONT>
<DIV align=left><FONT color=#c0c0c0 size=2
face=Verdana>2013-09-06</FONT></DIV><FONT size=2 face=Verdana>
<HR id=SignNameHR style="HEIGHT: 2px; WIDTH: 122px" align=left SIZE=2>
</FONT>
<DIV><FONT color=#c0c0c0 size=2 face=Verdana><SPAN
id=_FlashSignName>朱学明</SPAN></FONT></DIV></STATIONERY></BODY></HTML>