[RegCNET] questions about rcm2rgrid

[Siraj Ul Islam]

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.ictp.it/pipermail/regcnet/attachments/20070405/abc251bd/attachment-0002.html>