Re: longitudes array

The error message is telling you that you are doing some calculation on the right-hand side that is return an array that is different than the array that you have on the left-hand side.

Your "v" array is dimensioned lat x lon x time. Your "avg_sahel_t" array is dimensioned 21 x 33 x 2 12.

You have to make sure, then, that the leftmost two dimensions of avg_sahel_t (21 x 33) are the same as the leftmost two dimensions as "v". Otherwise, you can't do the assignment:

> avg_sahel_t(:,:,iy,im) = dim_avg_Wrap(v)

--Mary

On Jan 18, 2012, at 12:32 PM, elodyfluck@orange.fr wrote:

>
>
>
>
>
> > Message du 18/01/12 16:39
> > De : elodyfluck@orange.fr
> > A : "ncl-talk [ncl-talk]"
> > Copie à :
> > Objet : longitudes array
> >
> > Hi all NCL users,
>
> I'm a beginner in NCL language.I would like to have an average of t2m in multiple NetCDF files from ERA_INTERIM for a specific array: longitude{-21:27} and latitude{0:30}.
> Here t2m has dimensions: (time,latitude,longitude). Longitudes are going from 0 to 360 (so I convert it in {-180:180} to specify my area), and latitudes goes 90 to -90.
> With the script below, I cannot plot my desired area for the specific array but I have all longitudes "compressed" in region (0 to 30N; -21W to 27E).
>
> Version of NCL is 4.3.0.
>
> Any suggestion?
> Thanks in advance
> Elody
>
>
>
> 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"
>
> begin
>
>
> avg_sahel_t = new((/21,33,2,12/),float)
> yrs=(/"1996","1997"/)
> yr=(/"96","97"/)
> mns=(/"01","02","03","04","05","06","07","08","09","10","11","12"/)
>
> do iy = 0,1
> do im = 0,11
> f = addfile("/archive/crct/shared/ERA_INTERIM/SURFACE"+"/erai_sfc_"+yrs(iy)+"_" +mns(im)+".nc","r")
> print(iy+" "+im)
>
> var = f->t2m(:,{0:30},:)
> ;printVarSummary(var)
> v = var(latitude|:,longitude|:,time|:)
> ;printVarSummary(v)
>
> lat2d = f->latitude
> lon2d = f->longitude
> mlon = 240
> lon = lonGlobeF(mlon,"lon2d","longitude","degrees_east")
> lon = (/lon-180./)
> lon&lon2d = lon
> ;print(lon)
>
> minlon = -21
> maxlon = 27
> olon = 33
>
> lon1 = fspan(minlon,maxlon,olon)
> lon1@units = "degrees_east"
> ;print(lon1)
>
> v!0 = "lat2d"
> v!1 = "lon1"
> v@lon1 = lon1
> printVarSummary(v)
> avg_sahel_t!0 = "lat2d"
> avg_sahel_t!1 = "lon1"
> avg_sahel_t(:,:,iy,im) = dim_avg_Wrap(v) : fatal:Dimension size mismatch on subscript #1, left-hand and right-hand side dimensions do not match
> printVarSummary(avg_sahel_t)
> delete(f)
> delete(var)
> delete(v)
> end do
> end do
>
> vv = dim_avg_Wrap(dim_avg_Wrap(avg_sahel_t))*0.001554884
>
> printVarSummary(vv)
> print(min(vv)+" "+max(vv))
>
>
> wks = gsn_open_wks("ps","sahel2")
> setvalues wks
> "wkColorMap" : "rainbow"
> end setvalues
>
> res = True
> res@cnLinesOn = False
> res@cnFillOn = True
> res@cnLineLabelsOn = False
>
> res@cnRasterModeOn = True
> res@gsnAddCyclic = False
> ;Zoom in Sahel
> res@mpLimitMode = "LatLon"
> res@mpMinLatF = -90
> res@mpMaxLatF = 90
> res@mpMinLonF = -180
> res@mpMaxLonF = 180
>
> res@mpCenterLonF = (res@mpMaxLonF-res@mpMinLonF)/2
>
> res@tfDoNDCOverlay = True
>
> res@mpGeophysicalLineColor = 2
> res@mpPerimOn = True
> res@mpGridLineDashPattern = 2
> res@mpOutlineBoundarySets = "AllBoundaries"
> res@mpUSStateLineColor = "black"
> res@pmTickMarkDisplayMode = "Always"
> res@pmLabelBarWidthF = 0.6
> res@pmLabelBarHeightF = 0.1
> res@lbLabelFontHeightF = 0.013
>
> res@cnLevelSelectionMode = "ExplicitLevels"
> res@cnLevels =(/5,7,9,11,13,15,17,19,21,23,25,27,29,31,33,35,37/)
> res@cnFillColors=(/12,23,45,52,57,62,73,89,98,128,136,147,152,157,160,170,175,186/)
> ;res@cnFillColors=(/12,23,45,52,57,62,73,89,98,128,136,147,152,157,160,170,175,186/)
>
> res@tiMainString = "Annual mean temperature 1996-1997"
>
> plot = gsn_csm_contour_map(wks,vv,res)
>
>
> end
>
> PrintVar Summary of v gives:
>
> Variable: v
> Type: short
> Total Size: 1249920 bytes
> 624960 values
> Number of Dimensions: 3
> Dimensions and sizes: [lat2d | 21] x [lon1 | 240] x [time | 124]
> Coordinates:
> lat2d: [ 0..30]
> lon1: [ 0..358.5]
> time: [841512..842250]
> Number Of Attributes: 7
> lon1 :
> long_name : 2 metre temperature
> units : K
> missing_value : -32767
> _FillValue : -32767
> add_offset : 262.5528775112522
> scale_factor : 0.001703177224306923
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Jan 18 14:19:35 2012