CDO

Hi there,

I am trying to use remapcon to interpolate data on a regular lat-lon grid to a target grid that is regular in NAD_1983_HARN_StatePlane_Florida_East_FIPS_0901 (resolution of 2 mi x 2 mi so lat and lon vary with lon and lat). Is this possible with remapcon or remapnn?
I have lat and lon as arrays for the 2 mi x 2 mi target grid. And they are function of dimensions row and column, respectively.
My main variable, rainfall, is a function of dimensions col, row, and time.

When I try to run remapcon on it, I get:
cdo remapcon (Abort): Unsupported target grid type (generic)!

Does this mean I have to make my rainfall variable have lon, lat, and time dimensions instead? How, if my lon and lat are arrays?

This is the output of R's RNetCDF print.nc on my target grid:

netcdf classic {
dimensions:
timestamps = 37621 ;
cells = 4830 ;
col = 70 ;
row = 119 ;
pair = 2 ;
variables:
NC_FLOAT timestamps(timestamps) ;
NC_CHAR timestamps:long_name = "Time in days" ;
NC_CHAR timestamps:units = "minutes since 1914-01-01 00:00:00" ;
NC_INT cells(cells) ;
NC_CHAR cells:long_name = "ID of GRIDIO cells" ;
NC_CHAR cells:units = "Index number" ;
NC_INT col(col) ;
NC_CHAR col:long_name = "GRIDIO column number" ;
NC_CHAR col:units = "Index number" ;
NC_INT row(row) ;
NC_CHAR row:long_name = "GRIDIO row number" ;
NC_CHAR row:units = "Index number" ;
NC_INT roco(pair, cells) ;
NC_CHAR roco:long_name = "Row/Column of GRIDIO cells" ;
NC_CHAR roco:units = "Index number" ;
NC_FLOAT coords(pair, cells) ;
NC_CHAR coords:long_name = "XY cordinates of GRIDIO cell centroid" ;
NC_CHAR coords:units = "State Plane Coordinates (ft)" ;
NC_CHAR coords:Project Coordinate System = "NAD_1983_HARN_StatePlane_Florida_East_FIPS_0901" ;
NC_CHAR coords:Projection = "Transverse Mercator" ;
NC_CHAR coords:False_Easting = "656166.66666700" ;
NC_CHAR coords:False_Northing = "0.00000000" ;
NC_CHAR coords:Central Meridian = "-81.00000000" ;
NC_CHAR coords:Scale_Factor = "0.99994100" ;
NC_CHAR coords:Latitude_Of_Origin = "24.33333300" ;
NC_CHAR coords:Linear Unit = "US Foot" ;
NC_FLOAT dxdy(pair, cells) ;
NC_CHAR dxdy:long_name = "Column/row width of GRIDIO cell" ;
NC_CHAR dxdy:units = "Linear ft" ;
NC_FLOAT area(cells) ;
NC_CHAR area:long_name = "Area of GRIDIO cell" ;
NC_CHAR area:units = "ft^2" ;
NC_FLOAT rainfall(col, row, timestamps) ;
NC_CHAR rainfall:long_name = "Rainfall" ;
NC_CHAR rainfall:units = "inches" ;
NC_FLOAT lat(col, row) ;
NC_CHAR lat:units = "degreeN" ;
NC_FLOAT lon(col, row) ;
NC_CHAR lon:units = "degreeE" ;

// global attributes:
NC_CHAR :Version = "Rainfall v4.7 1914-2016" ;
NC_CHAR :Description = "QA/QC rainfall gauge data using TIN-10 method from 1914 to May 2002. June 2002 to December 2016 use NEXRAD." ;

Thanks for any assistance,

Michelle