check_for_y_lat_coord

From: 朱学明 <557xiaozhu_at_nyahnyahspammersnyahnyah>
Date: Thu Sep 05 2013 - 12:39:17 MDT

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