Hi Moet and dennis

 

Thanks alot. Script is now working ok and  making new file with variable on regular lat lon grid. One another thing that i have seen in original and regrided file is that the values are slightly changed. Is this due to the interpolation of the data i.e rcm2grid fuction .

 

Siraj

 

On 4/3/07, Moetasim <mashfaq@purdue.edu> wrote:

Siraj, I think you got everything right after Dennis reply. I just forgot to
remove that extra False from there while modifying the 3d data script to 2d
data.



Quoting Siraj Ul Islam <sirajkhan78@gmail.com>:

> 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