Re: rcm2rgrid problem

The rcm2rgrid assumes the latitudes are order south-to-north
like in the REGcm grid.

The MM5 is ordered north-to-south. Try

    precip = f->precip ; (time,lat,lon)
    lat2d = f->xlat ; lat order is S->N
    lon2d = f->xlon

    lat2d = (/ lat2d(::-1,:) /) ; rearrange values only
    precip = (/ precip(:,::-1,:) /)

The documentation has been updated to reflect this.o
Good luck

IMRAN NADEEM wrote:
> Hi everybody,
>
> I have precipitation fields from RegCM3 and MM5 models. I have
> written a script to interpolate from model grid to regular grid. The
> script works fine for RegCM3 data
> but produces only missing values for mm5 data. Both datasets covers the
> same region and have almost same number of grid points. The script is
> given below and few time
> steps of both datasets are here (Size 600KB).
>
> http://homepage.boku.ac.at/imrannad/Precip_regcm.nc
> http://homepage.boku.ac.at/imrannad/Precip_mm5.nc
>
>
> Any Help/Suggestion is highly appreciated!
>
> Regards
> Imran Nadeem
> BOKU, Vienna
>
>
>
> ------------------------------
> ------------------------Rcm2rgrid.ncl----------------------------------------------
>
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
>
>
> begin
>
> ; ---------Read input file--------------------
> filein="Precip_regcm.nc" ; input file
> data_file= addfile(filein,"r")
>
> precip = data_file->precip
> lat2d = data_file->xlat
> lon2d = data_file->xlon
> time =data_file->time
>
> lat_g = latGlobeF(901, "lat", "latitude",
> "degrees_north") ;0.20 degree spacing
> lon_g = lonGlobeF(1800, "lon", "longitude", "degrees_east")
> ;0.20 degree spacing
>
> newlon = lon_g({2:19})
> newlat = lat_g({42:50})
>
> precip_rgrid = rcm2rgrid(lat2d,lon2d,precip,newlat,newlon,0) ;
> Regridding
>
> ;---------------- assign named dimensions----------------
> precip_rgrid!0 = "time"
> precip_rgrid!1 = "lat"
> precip_rgrid!2 = "lon"
>
> ;---------------- assign attributes----------------
> precip_rgrid_at_long_name="Precipitation on Regular Grid"
> precip_rgrid_at_units="mm/day"
> precip_rgrid@_FillValue=-1.e+30
>
> ;---------------- assign coordinate variables----------------
> precip_rgrid&lat = newlat
> precip_rgrid&lon = newlon
> precip_rgrid&time =time
>
> ;-----------create output file---------------
> outfile_name = "Precip_regcm_rgrid.nc"
> system(" 'rm' -f " + outfile_name)
> outfile = addfile(outfile_name, "c")
> filedimdef(outfile,"time",-1,True) ; make time and
> record dimension.
> outfile_at_title = "regcm precip on regular lat lon"
>
> ;------------write to output file---------------
> outfile->precip = precip_rgrid
> outfile->lon=newlon
> outfile->lat=newlat
>
> end
> <http://homepage.boku.ac.at/imrannad/Precip_mm5.nc>
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Oct 23 2008 - 08:38:35 MDT