<html><head><meta http-equiv="Content-Type" content="text/html charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">Hi Mary <div><br></div><div>I think the interpolation still generates some errors. I sent you the values and the minimum still shows minus.</div><div>I plotted original data and interpolated data with same legend scale. It showed quite different plotting. I modified the script as you suggested.</div><div><div><div>The original data has following min max => (0) Ambient Aerosol Extinction at 550 nm: min=0 max=0.00242234</div><div>After interpolation => (0)<span class="Apple-tab-span" style="white-space: pre;"> </span>min=-0.000818642 max=0.00353843</div><div>I think it should not be minus value. How can I fix this?</div></div><div><br></div><div>Attached plots. the first plot is from original data and second is after interpolation. </div><div><br></div><div>Any suggestions will help.</div><div><br></div><div>Thank you</div><div>Sun- </div></div><div><br></div><div><div><font color="#e32400">hyam = filenc->a(::-1)</font></div><div><font color="#e32400">hybm = filenc->b(::-1)</font></div><div><br></div><div>;hyam = filenc->a</div><div>;hybm = filenc->b</div><div>plev0 = (/1000,925,850,700,600,500,400,300,250,200,150,100,70,50,30,20,10/)</div><div>psfc = filenc->ps</div><div>p0=1000.</div><div><br></div><div>;printVarSummary(psfc)</div><div>;printVarSummary(hyam)</div><div><br></div><div><font color="#b51a00">ec550aer = data(:,::-1,:,:)</font></div><div>printVarSummary(ec550aer)</div><div>;ec550aer = data ; read ec550aer from all files</div><div>ec550aerNEW=vinth2p(ec550aer,hyam,hybm,plev0,psfc,1,p0,1,True)</div><div>delete(ec550aer)</div><div>ec550aer=ec550aerNEW</div><div>delete(ec550aerNEW)</div><div>printMinMax(ec550aer,0)</div><div>return</div></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div></div></body></html>