[ncl-talk] reshape error
    Sunmin Park 
    mireiyue at gmail.com
       
    Thu Jul 24 11:26:23 MDT 2014
    
    
  
Dear all nil users	
I am trying to get plots from each month and year. One file contains all 7 years and 12 months values. So, I used reshape but I got a problem. I try to follow the reshape example in the website but it is not working. How can I fix it?
Thank you for any help
Sun-
This is what I got 
fatal:Argument type mismatch on argument (0) of (reshape) can not coerce
fatal:["Execute.c":7743]:Execute: Error occurred at or near line 53 in file CMIP5_ec.ncl
; ==============================================================
; Open the file:
; ==============================================================
  fili    = "ec550aer_aero_NorESM1-M_historicalExt_"+version(ii)+"i1p1_200601-201212.nc"
  varname = "ec550aer” 
 filenc       = addfile (diri+fili+".nc", "r")      ; entire file
;  print (varname)
;  print (filenc)
  data   = filenc->$varname$                          ; (time, ver, lat, lon)
;  print (dimsizes(data))
  time =filenc->time
  lev =filenc->lev
  lat =filenc->lat
  lon =filenc->lon
  data at _FillValue= 0
;  printVarSummary(data)
;  print(dimsizes(time))
;  return
  ntim = dimsizes(time)
  nlev = dimsizes(lev)
  nlat = dimsizes(lat)
  nlon = dimsizes(lon)
  nyears = ntim/12
;  printVarSummary(ntim)
;  printVarSummary(nyears)
;  print(ntim)
;  print(nyears)
  ecreshape = reshape(varname,(/nyears,12,nlev,nlat,nlon/))
;  reec550nm = reshape(ecreshape(:,0:2,:,:,:),(/nyears*12,nlev,nlat,nlon/))
  printVarSummary(ecreshape)
;  printVarSummary(ec550nm)
  return
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140724/e1f25c8e/attachment.html 
    
    
More information about the ncl-talk
mailing list