rcm2points_Wrap

I've a curvilinear grid from the Princeton Ocean Model I'm trying to
interpolate to a rectilinear grid using rcm2points_Wrap.

I'm reading in a model file with the following excertped from
ncl_filedump output:
path: GOMHD-2011.3.100.GOMEX-PPP.nc
    file global attributes:
       file_type : Full_Grid
       Conventions : COARDS
       grid_type : Orthogonal Curv Grid (IM=217,JM=160)
       z_type : Sigma Coordinates (KB= 37)
       model : NOS Gulf of Mexico Forecast System (NGOM
       title : GOMEX-PPP Initial Hindcast Test
       comment : Original MODAS Assimilation
       source : Model ID. GOMHD-2011.3.100
       institution : NOAA/NOS/CSDL
       history : Dynalysis Princeton Ocean Model
       references : rich.patchen@noaa.gov
       creation_date : 2011-01-06 8:34:50 00:00
    dimensions:
       nx = 217
       ny = 160
       sigma = 37
       time = 31 // unlimited
    variables:
       float time ( time )
          long_name : Time
          base_date : ( 2000, 8, 1, 0 )
          units : days since 2000-08-01 0:00:00 00:00
          standard_name : time

       float lon ( ny, nx )
          long_name : Longitude
          units : degrees_east
          standard_name : longitude

       float lat ( ny, nx )
          long_name : Latitude
          units : degrees_north
          standard_name : latitude

       float mask ( ny, nx )
          long_name : Land Mask
          units : nondimensional

       float depth ( ny, nx )
          long_name : Bathymetry
          units : meters
          positive : down
          standard_name : depth

       float sigma ( sigma )
          long_name : Sigma Stretched Vertical Coordinate at Nodes
          units : sigma_level
          positive : down
          standard_name : ocean_sigma_coordinate
          formula_terms : sigma: sigma eta: zeta depth: depth

       float temp ( time, sigma, ny, nx )
          long_name : Temperature
          units : Celsius
          missing_value : -7.762932
          standard_name : sea_water_temperature
          averaging_window : 1
          averaging_offset : 0
          _FillValue : -7.762932

I'm doing the following:
fnameModel6HC =
"/home/gerry/GOMEX-PPP/data/GOMEX-PPP/6_NGOM_Patchen/GOMHD-2011.3.100.GOMEX-PPP.nc"
modelFile6HC=addfile(fnameModel6HC,"r")

model6DepthHC = modelFile6HC->depth
model6LonHC = modelFile6HC->lon
model6LatHC = modelFile6HC->lat
model6WaterTempHC = modelFile6HC->temp(:,:,:,:)
model6TimeHC = modelFile6HC->time
waterTempHC=rcm2points_Wrap(model6LatHC,model6LonHC,model6WaterTempHC,model6LatHC1d,
model6LonHC1d,0)

And I'm getting: fatal:Variable (model6LatHC2d) is undefined

Do I need to declare model6LatHC1d

Thanks, Gerry

-- 
Gerry Creager -- gerry.creager@tamu.edu
Texas Mesonet -- AATLT, Texas A&M University
Cell: 979.229.5301 Office: 979.458.4020 FAX: 979.862.3983
Office: 1700 Research Parkway Ste 160, TAMU, College Station, TX 77843
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk