Noelia
I suspect that the lines below ,may not actually be doing anything.
> correl&lat@units = "degrees_north"
> correl&lon@units = "degrees_east"
Instead, make sure the lat/lon arrays have the required attributes BEFORE you assign them as coordinate arrays:
> lat@long_name = "latitude"
>lat@units = "degrees_north"
> correl&lat = lat
;do same thing for 'lon'
Will
From: Noelia otero <noeli1680@gmail.com<mailto:noeli1680@gmail.com>>
Date: Friday, 4 July 2014 4:08 PM
To: "ncl-talk@ucar.edu<mailto:ncl-talk@ucar.edu>" <ncl-talk@ucar.edu<mailto:ncl-talk@ucar.edu>>
Subject: message warning check_for_y_lat_coord
Hi,
I am looking for some support to fix the next problem. I have already read some posts concerning this message, but so far, I did not solve this (I am start learning the language).
after running the script, this is the message:
Variable: correl
Type: float
Total Size: 2016 bytes
504 values
Number of Dimensions: 2
Dimensions and sizes: [lon | 18] x [lat | 28]
Coordinates:
lon: [37.67308962904533..-9.767145559195573]
lat: [-28.125..47.8125]
Number Of Attributes: 1
_FillValue : 1e+20
(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'
And this is the code :
;*****************************************************************************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
;*******************************************************************************
begin
a = addfile("/home/nof/pruebas/ensem_means/pr_Amon_bcc-csm1-1_historical_ensmean_t42_westA_meanJAS.nc","r")
b = addfile("/home/nof/pruebas/WAMI_CEOF/WAMI_Amon_bcc-csm1-1_historical_ensmean.nc","r")
pr =a->pr
wami =b->WAMI(:,0,0,0)
lat = a->lat
lon = a->lon
pr!0="time" ; assign attributes
pr!1="lon"
pr!2="lat"
wami!0="time"
;pr&lat=lat
;pr&lon=lon
prT = transpose(pr)
copy_VarCoords(pr,wami) ; copy attributes from temp to pcp
correl=escorc(wami,prT(lon|:,lat|:,time|:))
correl!1="lat"
correl!0="lon"
correl&lat=pr&lat
correl&lon=pr&lon
correl&lat@units = "degrees_north"
correl&lon@units = "degrees_east"
printVarSummary(correl)
wks = gsn_open_wks("x11","example")
plot = gsn_csm_contour_map_ce(wks,correl,False)
end
I don't understand the message very much,because looking into the variables and dimensions, the names seem correct to me, so I don't know where the problem is...what am I missing?
Thanks a lot in advance!
Cheers,
Noelia
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Jul 04 06:14:12 2014
This archive was generated by hypermail 2.1.8 : Wed Jul 23 2014 - 15:33:47 MDT