[RegCNET] questions about rcm2rgrid

Moetasim mashfaq at purdue.edu
Thu Apr 5 18:17:28 CEST 2007


In your case, it's more due to the difference in res of RegCM3 grid and 
regridded grid. You may decrease the difference by improving the newlat,newlon 
declaration. 

The very simplest of doing that is to change multiplier from 0.5 to 0.4545 
(50km ~ 0.4545 degree). Or you may try fspan function instead of ispan. 

Moet

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

> 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 at 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 at 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
> 



More information about the RegCNET mailing list