Problem with regriddig 1km curvlinear data set to 8km curvlinear data set

Hi,

I have been trying to re-grid DAYMET 1km dataset into 8km using the
example NCL script prepared to convert from curvlinear to rectilinear
grd (http://www.ncl.ucar.edu/Applications/Scripts/daymet_3.ncl). My
problem is that I want to convert 1km curvlinear DAYMET data into 8km
curvlinear dataset. Therefore, I set the

Opt@DstGridType = "curvilinear"

I have not changed any other option in the example script. The script
works fine when used to convert from curvlinear to rectilinear grd but
failed when trying to use to convert 1km curvlinear DAYMET data into
8km curvlinear data set. It gives the following error message:

gpc-f101n084-$ ncl test.ncl
  Copyright (C) 1995-2012 - All Rights Reserved
  University Corporation for Atmospheric Research
  NCAR Command Language Version 6.1.0
  The use of this software is governed by a License Agreement.
  See http://www.ncl.ucar.edu/ for more details.

Variable: lat
Type: double
Total Size: 168 bytes
             21 values
Number of Dimensions: 1
Dimensions and sizes: [21]
Coordinates:

Variable: lon
Type: double
Total Size: 168 bytes
             21 values
Number of Dimensions: 1
Dimensions and sizes: [21]
Coordinates:
fatal:Number of dimensions in parameter (1) of (curvilinear_to_SCRIP)
is (1), (2) dimensions were expected
fatal:["Execute.c":7743]:Execute: Error occurred at or near line 2978
in file $NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl

fatal:["Execute.c":7743]:Execute: Error occurred at or near line 3098
in file $NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl

fatal:["Execute.c":7743]:Execute: Error occurred at or near line 100
in file test.ncl

==================================================
Here is the source data description:
==================================================
gpc-f101n084-$
/scinet/gpc/Libraries/netcdf-4.1.3/serial_gcc/bin/ncdump -h
prcp.2010_12296.nc
netcdf prcp.2010_12296 {
dimensions:
         x = 212 ;
         y = 250 ;
         time = 365 ;
         nv = 2 ;
variables:
         short lambert_conformal_conic ;
                 lambert_conformal_conic:grid_mapping_name =
"lambert_conformal_conic" ;
                 lambert_conformal_conic:longitude_of_central_meridian
= -100. ;
                 lambert_conformal_conic:latitude_of_projection_origin = 42.5 ;
                 lambert_conformal_conic:false_easting = 0. ;
                 lambert_conformal_conic:false_northing = 0. ;
                 lambert_conformal_conic:standard_parallel = 25., 60. ;
         double x(x) ;
                 x:units = "m" ;
                 x:long_name = "x coordinate of projection" ;
                 x:standard_name = "projection_x_coordinate" ;
         double y(y) ;
                 y:units = "m" ;
                 y:long_name = "y coordinate of projection" ;
                 y:standard_name = "projection_y_coordinate" ;
         double time(time) ;
                 time:long_name = "time" ;
                 time:calendar = "standard" ;
                 time:units = "days since 1980-01-01 00:00:00 UTC" ;
                 time:bounds = "time_bnds" ;
         double lat(y, x) ;
                 lat:units = "degrees_north" ;
                 lat:long_name = "latitude coordinate" ;
                 lat:standard_name = "latitude" ;
         double lon(y, x) ;
                 lon:units = "degrees_east" ;
                 lon:long_name = "longitude coordinate" ;
                 lon:standard_name = "longitude" ;
         short yearday(time) ;
                 yearday:long_name = "yearday" ;
                 yearday:valid_range = 1s, 365s ;
         float prcp(time, y, x) ;
                 prcp:long_name = "daily total precipitation" ;
                 prcp:units = "mm/day" ;
                 prcp:missing_value = -9999.f ;
                 prcp:_FillValue = -9999.f ;
                 prcp:valid_range = 0.f, 200.f ;
                 prcp:coordinates = "lat lon" ;
                 prcp:grid_mapping = "lambert_conformal_conic" ;
                 prcp:cell_methods = "area: sum time: sum" ;
         double time_bnds(time, nv) ;

// global attributes:
                 :tileid = 12296s ;
                 :start_year = 2010s ;
                 :source = "Daymet Software Version 2.0" ;
                 :Version_software = "Daymet Software Version 2.0" ;
                 :Version_data = "Daymet Data Version 2.1" ;
                 :Conventions = "CF-1.4" ;
                 :citation = "Please see http://daymet.ornl.gov/ for
current Daymet data citation information" ;
                 :references = "Please see http://daymet.ornl.gov/ for
current information on Daymet references" ;
}

===========================================================================

I look forward to hearing from you.

Thanks
Kamal

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Feb 20 14:30:05 2013