[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