[Thread Prev][Thread Next][Index]

[ferret_users] REGRIDDING netcdf



Greetings all 


I need help in understanding a weird error that occurs after I regrid data: 

ERROR: yes? fill /x=40E:100E/y=10S:40S SST1[l=1]
              **ERROR: illegal limits: "SST1[L=1]" does not exist at X=40E:100E
              Axis extremes are X=179.8E:179.8W(540.3)

if I did fill where X=W:E and Y=N:S, this will work but if X values are in the East/West only or Y in North/South only , then the error occurs.


Here is the ncdump of the original data :

netcdf ROMS {
dimensions:
LON_AXIS = 230 ;
LAT_AXIS = 194 ;
bnds = 2 ;
TDIM = UNLIMITED ; // (120 currently)
variables:
double LON_AXIS(LON_AXIS) ;
LON_AXIS:long_name = "longitude" ;
LON_AXIS:units = "degrees_east" ;
LON_AXIS:axis = "X" ;
LON_AXIS:modulo = 360. ;
LON_AXIS:point_spacing = "even" ;
LON_AXIS:standard_name = "longitude" ;
double LAT_AXIS(LAT_AXIS) ;
LAT_AXIS:long_name = "latitude" ;
LAT_AXIS:units = "degrees_north" ;
LAT_AXIS:axis = "Y" ;
LAT_AXIS:point_spacing = "uneven" ;
LAT_AXIS:standard_name = "latitude" ;
LAT_AXIS:bounds = "LAT_AXIS_bnds" ;
double LAT_AXIS_bnds(LAT_AXIS, bnds) ;
double TDIM(TDIM) ;
TDIM:units = "months since 1999-01-16 00:00:00" ;
TDIM:calendar = "365_day" ;
TDIM:axis = "T" ;
TDIM:standard_name = "time" ;
double SST(TDIM, LAT_AXIS, LON_AXIS) ;
SST:missing_value = -1.e+34 ;
SST:_FillValue = -1.e+34 ;
SST:long_name = "SST_CPL" ;
SST:history = "From ROMSData" ;

// global attributes:
:history = "FERRET V6.82   14-Aug-14" ;
:Conventions = "CF-1.0" ;
}


Here is the regridding I performed which regridded the file successfully: 

SET MEMORY/SIZE=600
SET DATA ROMS.nc
let SST1 =SST[X=180E:180W:0.5@ave,Y=90N:90S:0.5@ave]
SAVE/FILE =0.5_ROMS.nc SST1


Regards

Molulaqhooa




--
"Try and Fail than to Fail to Try" - Israel Maoyi

[Thread Prev][Thread Next][Index]
Contact Us
Dept of Commerce / NOAA / OAR / PMEL / Ferret

Privacy Policy | Disclaimer | Accessibility Statement