[RegCNET] questions about rcm2rgrid

Siraj Ul Islam sirajkhan78 at gmail.com
Mon Apr 2 13:06:59 CEST 2007


Hi Moet,

thanks for your comments. Well I have modify the script as you have
suggested. Now the rcm2grid funtion is working but the errors comes
from the line 56( ; ####ERROR LINE)  now where the output RT_re.nc
file is written. I am sending you the error as well as the modified
script with NC dump of HEAD_OUT and SRF files.
Hope you will again help me to solve this problem.

cheers

Siraj
----------------ERROR-------------------------------
Variable: newvar
Type: float
Total Size: 6791136 bytes
            1697784 values
Number of Dimensions: 3
Dimensions and sizes:   [132] x [109] x [118]
Coordinates:
Number Of Attributes: 1
  _FillValue :  -1e+30
fatal:Execute: Error occurred at or near line 56 in file rcm2grid.ncl

---------SCRIPT-------------------------------------
load "/usr/local/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "/usr/local/lib/ncarg/nclscripts/csm/shea_util.ncl"
load "/usr/local/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "/usr/local/lib/ncarg/nclscripts/csm/contributed.ncl"
;********************************************
begin
;Set path to HEAD_OUT.nc
dirhead="/root/"
filhead="head_out.nc"
;Set path to RegCM file
dirrcm="/root/"
filrcm="SRF.nc"
;Pointers to input files
fhead=addfile(dirhead+filhead,"r")
f1 = addfile(dirrcm+filrcm,"r")
;Read in Data
var =  f1->RT
lat2d = fhead->XLAT(0,:,:)
lon2d = fhead->XLON(0,:,:)
newlon = new((/118/),typeof(lon2d))
newlat = new((/109/),typeof(lat2d))
do j=0,117
newlon(j)= 50+(0.5)*j
end do
do i=0,108
newlat(i)=5+(0.5)*i
end do
newvar = rcm2rgrid(lat2d,lon2d,var,newlat,newlon,1)
printVarSummary(newvar)
ncdf=addfile("RT_re.nc","c")
dimNames=(/"time","lat","lon"/)
dimSizes=(/-1,109,118/)
dimUnlim=(/True,False,False,False/)
filedimdef(ncdf,dimNames,dimSizes,dimUnlim)   ; ####ERROR LINE
filevardef(ncdf,"RT",typeof(var),(/"time","lat","lon"/))
ncdf->lat=(/newlat/)
ncdf->lon=(/newlon/)
ncdf->RT=(/newvar/)
end

----------------OUT_HEAD----------------

filename:       head_out
path:   head_out.nc
   file global attributes:
      domxmin : 52.37275
      domxmax : 108.0947
      domymin : -2.378861
      domymax : 44.79095
      domzmin : 1050
      domzmax :  0
   dimensions:
      lon = 118
      lat = 109
      time = 1  // unlimited
   variables:
      float lon ( lon )
         long_name :    Longitude
         units :        degrees_east
         actual_range : <ARRAY>

---------------------------------------------------
> filename:       SRF
> path:   SRF.nc
>  file global attributes:
>     domxmin : 52.37275
>     domxmax : 108.0947
>     domymin : -2.378861
>     domymax : 44.79095
>     domzmin : 1050
>     domzmax : 1050
>  dimensions:
>     lon = 118
>     lat = 109
>     time = 132  // unlimited
>  variables:
>     float lon ( lon )
>        long_name :    Longitude
>        units :        degrees_east
>        actual_range : <ARRAY>

>     float lat ( lat )
>        long_name :    Latitude
>        units :        degrees_north
>        actual_range : <ARRAY>
>
--------------------------------------------
Siraj Ul Islam
Research Fellow (Climate Section)
Global Change Impact Studies Centre
First Floor Saudi Pak tower
Blue Area Islamabad Pakistan
Ph :   92-51-9219785
Fax : 92-51-9219787
http://www.gcisc.org.pk



More information about the RegCNET mailing list