<div dir="ltr"><div class="gmail_default" style="font-size:small">My first thought this looked like a _FillValue issue, but I can't tell simply from the output.</div><div class="gmail_default" style="font-size:small"><br>
</div><div class="gmail_default" style="font-size:small">Can you provide the data file for this script? You can use our ftp:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style>
<a href="http://www.ncl.ucar.edu/report_bug.shtml#HowToFTP">http://www.ncl.ucar.edu/report_bug.shtml#HowToFTP</a><br></div><div class="gmail_default" style><br></div><div class="gmail_default" style>Thanks,</div><div class="gmail_default" style>
<br></div><div class="gmail_default" style>--Mary</div><div class="gmail_default" style><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Aug 6, 2014 at 3:16 PM, Sunmin Park <span dir="ltr"><<a href="mailto:mireiyue@gmail.com" target="_blank">mireiyue@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear NCL users,<br>
<br>
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<br>
"(0) Ambient Aerosol Extinction at 550 nm: min=0 max=0.00242234 “<br>
But after interpolation the values are<br>
"(0) min=-2.58212e+22 max=0.0127834”<br>
<br>
I do not understand why the values are changing and the MIN value shows minus. How can I fix it?<br>
<br>
Thank you for any help!!<br>
<br>
Sun-<br>
<br>
; ==============================================================<br>
; Open the file:<br>
; ==============================================================<br>
fili = "ec550aer_aero_NorESM1-M_historicalExt_r1i1p1_200601-201212.nc"<br>
; fili = "ec550aer_aero_NorESM1-M_historicalExt_"+version(ii)+"<a href="http://i1p1_200601-201212.nc" target="_blank">i1p1_200601-201212.nc</a>"<br>
varname = "ec550aer"<br>
varname1 = "a"<br>
varname2 = "b"<br>
varname3 = "ps"<br>
varname4 = "p0"<br>
filenc = addfile (diri+fili+".nc", "r") ; entire file<br>
; printMinMax(varname,0)<br>
; print(varname1)<br>
; print(varname2)<br>
; print(varname3)<br>
; print(varname4)<br>
; print (filenc)<br>
; return<br>
<br>
data = filenc->$varname$ ; (time, ver, lat, lon)<br>
; print (dimsizes(data))<br>
time =filenc->time<br>
lev =filenc->lev<br>
lat =filenc->lat<br>
printMinMax(data,0)<br>
printVarSummary(data)<br>
;return<br>
lon =filenc->lon<br>
; print(time)<br>
; print(lev)<br>
; print(lat)<br>
; print(lon)<br>
; data@_FillValue=1e+20<br>
; newlev =lev*1000<br>
; print(newlev)<br>
<br>
;return<br>
<br>
;############Converting from hybrid levels to standard pressure levels######################################<br>
hyam = filenc->a<br>
hybm = filenc->b<br>
plev0 = (/1000,925,850,700,600,500,400,300,250,200,150,100,70,50,30,20,10/)<br>
psfc = filenc->ps<br>
p0=1000.<br>
<br>
;printVarSummary(psfc)<br>
;printVarSummary(hyam)<br>
<br>
ec550aer = data;filenc->ec550aer ; read ec550aer from all files<br>
ec550aerNEW=vinth2p(ec550aer,hyam,hybm,plev0,psfc,1,p0,1,True)<br>
delete(ec550aer)<br>
ec550aer=ec550aerNEW<br>
delete(ec550aerNEW)<br>
printVarSummary(ec550aer)<br>
printMinMax(ec550aer,0)<br>
return<br>
;###################################################################################<br>
<br>
<br>
RESULTS<br>
<br>
(0) Ambient Aerosol Extinction at 550 nm: min=0 max=0.00242234<br>
<br>
Variable: data<br>
Type: float<br>
Total Size: 120766464 bytes<br>
30191616 values<br>
Number of Dimensions: 4<br>
Dimensions and sizes: [time | 84] x [lev | 26] x [lat | 96] x [lon | 144]<br>
Coordinates:<br>
time: [15.5..2539.5]<br>
lev: [0.9925560999999998..0.00354463800000001]<br>
lat: [ -90.. 90]<br>
lon: [ 0..357.5]<br>
Number Of Attributes: 12<br>
standard_name : volume_extinction_coefficient_in_air_due_to_ambient_aerosol<br>
long_name : Ambient Aerosol Extinction at 550 nm<br>
comment : """ambient"" means ""wetted"". "<br>
units : m-1<br>
original_name : BAK<br>
original_units : 1/km<br>
<br>
<br>
Variable: ec550aer<br>
Type: float<br>
Total Size: 78962688 bytes<br>
19740672 values<br>
Number of Dimensions: 4<br>
Dimensions and sizes: [time | 84] x [lev_p | 17] x [lat | 96] x [lon | 144]<br>
Coordinates:<br>
time: [15.5..2539.5]<br>
lev_p: [1000..10]<br>
lat: [ -90.. 90]<br>
lon: [ 0..357.5]<br>
Number Of Attributes: 1<br>
_FillValue : 1e+20<br>
(0) min=-2.58212e+22 max=0.0127834<br>
<br>
_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
</blockquote></div><br></div>