longitudes array

> 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
Received on Wed Jan 18 12:32:30 2012