Hi All,
I am trying to plot the geostrophic current fields from aviso. And get some screen output and error messages as follow,
(0) 19930101
Variable: u
Type: float
Total Size: 664692 bytes
166173 values
Number of Dimensions: 2
Dimensions and sizes: [NbLongitudes | 579] x [NbLatitudes | 287]
Coordinates:
NbLongitudes: [98.33333333332349..290.9999999999709]
NbLatitudes: [-17.79766126545592..60.88112837177402]
Number Of Attributes: 5
date : 1993-01-01 00:00:00.000000 UTC
Date_CNES_JD : 15706
units : cm/s
long_name : U
_FillValue : 1.844674e+19
Variable: v
Type: float
Total Size: 664692 bytes
166173 values
Number of Dimensions: 2
Dimensions and sizes: [NbLongitudes | 579] x [NbLatitudes | 287]
Coordinates:
NbLongitudes: [98.33333333332349..290.9999999999709]
NbLatitudes: [-17.79766126545592..60.88112837177402]
Number Of Attributes: 5
_FillValue : 1.844674e+19
long_name : V
units : cm/s
Date_CNES_JD : 15706
date : 1993-01-01 00:00:00.000000 UTC
Variable: speed
Type: float
Total Size: 664692 bytes
166173 values
Number of Dimensions: 2
Dimensions and sizes: [579] x [287]
Coordinates:
Number Of Attributes: 1
_FillValue : 1.844674e+19
(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'
(0) gsn_csm_map_ce: Warning: you set mpMaxLonF to a value > 180, but
(0) didn't set mpCenterLonF. Setting mpCenterLonF to 194.5
My code is as follow,
do nt=nb,ne
in = addfile(infile(nt-1),"r")
utc_date = str_get_cols(infile(nt-1),44,51)
print(""+utc_date)
;************************************************
; read in zonal [u] and meridional [v] winds
;************************************************
tmp = in->Grid_0001
u = tmp({lon_min:lon_max},{lat_min:lat_max}) ;/100.0
tmp = in->Grid_0002
v = tmp({lon_min:lon_max},{lat_min:lat_max}) ;/100.0
speed = sqrt(u^2 + v^2)
vcres@tiMainString = "Geostrophic Current Field"
vcres@gsnLeftString = utc_date
vcres@gsnRightString = "m/s"
wks = gsn_open_wks(wks_type ,"uv_"+utc_date) ; open a ps file
gsn_define_colormap(wks,"rainbow") ;"BlAqGrYeOrRe") ; choose colormap
printVarSummary(u)
printVarSummary(v)
printVarSummary(speed)
plot=gsn_csm_vector_scalar_map_ce(wks,u,v,speed,vcres) ; create plot
end do
Anyone can help me for figuring out the error? Thanks.
2013-09-06
朱学明
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Sep 5 12:39:42 2013
This archive was generated by hypermail 2.1.8 : Mon Sep 16 2013 - 13:43:52 MDT