Re: Regriding Problem

Dear NCl-talk

I am trying to interpolate temperature from an NCEP1 file into new lat and lon but the interpolated temperature shows different min and max with same value of 9.96921e+36. Could someone advise.
The part of the script and output are shown below

SCRIPT
;load variables
air2 = short2flt(f->air(:,:,:,:))
printVarSummary(air2)
air2@_FillValue = 32766
air2@missing_value = air2@_FillValue
print("air2: min="+min(air2)+" max="+max(air2))

;create lat & lon arrays
LAT=fspan(-37.00,-40.00,3)
LON=fspan(6.00,53.00,47)
LAT@units = "degrees_south"
LON@units = "degrees_east"

;interpolate to new lat & lon arrays
regrid=linint2_Wrap(air2&lon,air2&lat,air2,True,LON,LAT,0)
print("regrid: min="+min(regrid)+" max="+max(regrid))
printVarSummary(regrid)

OUTPUT
Variable: air2
Type: float
Total Size: 1043631360 bytes
            260907840 values
Number of Dimensions: 4
Dimensions and sizes: [time | 1460] x [level | 17] x [lat | 73] x [lon | 144]
Coordinates:
            time: [1779456..1788210]
            level: [1000..10]
            lat: [90..-90]
            lon: [ 0..357.5]
Number Of Attributes: 19
  long_name : 6-hourly Air Temperature on Pressure Levels
  unpacked_valid_range : ( 137.5, 362.5 )
  actual_range : ( 179.3, 326.2 )
  units : degK
  missing_value_original : 32766
  _FillValue_original : -32767
  _FillValue : -32767
  precision : 2
  least_significant_digit : 1
  GRIB_id : 11
  GRIB_name : TMP
  var_desc : Air temperature
  dataset : NCEP/DOE AMIP-II Reanalysis (Reanalysis-2)
  level_desc : Pressure Levels
  statistic : Individual Obs
  parent_stat : Other
  standard_name : air_temperature
  valid_range : ( 137.5, 362.55 )
  missing_value : -32767
(0) air2: min=179.3 max=326.2
(0) regrid: min=9.96921e+36 max=9.96921e+36

Variable: regrid
Type: float
Total Size: 13998480 bytes
            3499620 values
Number of Dimensions: 4
Dimensions and sizes: [time | 1460] x [level | 17] x [LAT | 3] x [LON | 47]
Coordinates:
            time: [1779456..1788210]
            level: [1000..10]
            LAT: [-37..-40]
            LON: [ 6..53]
Number Of Attributes: 19
  long_name : 6-hourly Air Temperature on Pressure Levels
  unpacked_valid_range : ( 137.5, 362.5 )
  actual_range : ( 179.3, 326.2 )
  units : degK
  missing_value_original : 32766
  _FillValue_original : -32767
  _FillValue : 32766
  precision : 2
  least_significant_digit : 1
  GRIB_id : 11
  GRIB_name : TMP
  var_desc : Air temperature
  dataset : NCEP/DOE AMIP-II Reanalysis (Reanalysis-2)
  level_desc : Pressure Levels
  statistic : Individual Obs
  parent_stat : Other
  standard_name : air_temperature
  valid_range : ( 137.5, 362.55 )
  missing_value : 32766
warning:int2p_n: 205860 input array(s) contained all missing data. No interpolation performed on these arrays

Thanks
Ditiro

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sat Oct 26 01:10:25 2013