[ncl-talk] Interpolation error
Sunmin Park
mireiyue at gmail.com
Wed Aug 6 15:16:44 MDT 2014
Dear NCL users,
I have a problem to do interpolation using NorESM1 model results. I am using extinction coefficient values and when I read the data from the files, the MIN and MAX values are
"(0) Ambient Aerosol Extinction at 550 nm: min=0 max=0.00242234 “
But after interpolation the values are
"(0) min=-2.58212e+22 max=0.0127834”
I do not understand why the values are changing and the MIN value shows minus. How can I fix it?
Thank you for any help!!
Sun-
; ==============================================================
; Open the file:
; ==============================================================
fili = "ec550aer_aero_NorESM1-M_historicalExt_r1i1p1_200601-201212.nc"
; fili = "ec550aer_aero_NorESM1-M_historicalExt_"+version(ii)+"i1p1_200601-201212.nc"
varname = "ec550aer"
varname1 = "a"
varname2 = "b"
varname3 = "ps"
varname4 = "p0"
filenc = addfile (diri+fili+".nc", "r") ; entire file
; printMinMax(varname,0)
; print(varname1)
; print(varname2)
; print(varname3)
; print(varname4)
; print (filenc)
; return
data = filenc->$varname$ ; (time, ver, lat, lon)
; print (dimsizes(data))
time =filenc->time
lev =filenc->lev
lat =filenc->lat
printMinMax(data,0)
printVarSummary(data)
;return
lon =filenc->lon
; print(time)
; print(lev)
; print(lat)
; print(lon)
; data at _FillValue=1e+20
; newlev =lev*1000
; print(newlev)
;return
;############Converting from hybrid levels to standard pressure levels######################################
hyam = filenc->a
hybm = filenc->b
plev0 = (/1000,925,850,700,600,500,400,300,250,200,150,100,70,50,30,20,10/)
psfc = filenc->ps
p0=1000.
;printVarSummary(psfc)
;printVarSummary(hyam)
ec550aer = data;filenc->ec550aer ; read ec550aer from all files
ec550aerNEW=vinth2p(ec550aer,hyam,hybm,plev0,psfc,1,p0,1,True)
delete(ec550aer)
ec550aer=ec550aerNEW
delete(ec550aerNEW)
printVarSummary(ec550aer)
printMinMax(ec550aer,0)
return
;###################################################################################
RESULTS
(0) Ambient Aerosol Extinction at 550 nm: min=0 max=0.00242234
Variable: data
Type: float
Total Size: 120766464 bytes
30191616 values
Number of Dimensions: 4
Dimensions and sizes: [time | 84] x [lev | 26] x [lat | 96] x [lon | 144]
Coordinates:
time: [15.5..2539.5]
lev: [0.9925560999999998..0.00354463800000001]
lat: [ -90.. 90]
lon: [ 0..357.5]
Number Of Attributes: 12
standard_name : volume_extinction_coefficient_in_air_due_to_ambient_aerosol
long_name : Ambient Aerosol Extinction at 550 nm
comment : """ambient"" means ""wetted"". "
units : m-1
original_name : BAK
original_units : 1/km
Variable: ec550aer
Type: float
Total Size: 78962688 bytes
19740672 values
Number of Dimensions: 4
Dimensions and sizes: [time | 84] x [lev_p | 17] x [lat | 96] x [lon | 144]
Coordinates:
time: [15.5..2539.5]
lev_p: [1000..10]
lat: [ -90.. 90]
lon: [ 0..357.5]
Number Of Attributes: 1
_FillValue : 1e+20
(0) min=-2.58212e+22 max=0.0127834
More information about the ncl-talk
mailing list