<div dir="ltr"><br clear="all"><div>Dear Users,</div><div><br></div><div>I want to plot stream function of wind (v) and I am using zonal_mpsi to calculate the stream function. I keep getting error messages.</div><div><br></div><div>fatal:zonal_mpsi: The pressure array must be between the values of 500 and 100500 (exclusive), and monotonically increasing.<br></div><div><br></div><div>I am not sure how to fix this problem.</div><div><br></div><div>I am attaching my code.</div><div><br></div><div>I read in the V wind and surface pressure</div><div><br></div><div>v-wind has 4 dimensions (time,lev,lat,lon)</div><div><br></div><div>Surface pressure has 3 dimensions (time,lat,lon)</div><div><br></div><div>Here is the code I wrote</div><div><br></div><div>;**************************************************************</div><div><div>;*******************************************</div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"</div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"</div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"</div><div>;*******************************************</div><div><br></div><div>;********************************************</div><div>;This is to read wind speed and wind direction</div><div>;from 1948-2016</div><div>;********************************************</div><div><br></div><div>a=addfile("/Users/ipshita/Documents/grads_NCL/<a href="http://wspd.mon.mean.nc">wspd.mon.mean.nc</a>","r")</div><div>b=addfile("/Users/ipshita/Documents/grads_NCL/<a href="http://uwnd.mon.mean.nc">uwnd.mon.mean.nc</a>","r")</div><div>c=addfile("/Users/ipshita/Documents/grads_NCL/<a href="http://vwnd.mon.mean.nc">vwnd.mon.mean.nc</a>","r")</div><div>d=addfile("/Users/ipshita/Documents/grads_NCL/<a href="http://omega.mon.mean.nc">omega.mon.mean.nc</a>","r")</div><div>e=addfile("/Users/ipshita/Documents/grads_NCL/<a href="http://pres.sfc.mon.mean.nc">pres.sfc.mon.mean.nc</a>","r")</div><div><br></div><div>v=c->vwnd ;Wind in horizontal direction in m/s</div><div><br></div><div>level=c->level ;Pressure level from 1000hpa to 10 hpa</div><div><br></div><div>lat=c->lat</div><div><br></div><div>lon=c->lon</div><div><br></div><div>ps=e->pres</div><div><br></div><div>ps_jan_1948=ps(0,:,:)</div><div><br></div><div>printMinMax(ps_1948_2015,False)</div><div><br></div><div>v_jan_1948=v(0,:,:,:)</div><div>;*********************************************</div><div>;Going from low to high in terms of level</div><div>;*********************************************</div><div><br></div><div>v1=v_jan_1948(::-1,:,:)</div><div>;*********************************************</div><div>;Converting Pressure into Pascals</div><div>;*********************************************</div><div><br></div><div>PsPa=ps_jan_1948*100</div><div><br></div><div>printMinMax(PsPa,False)</div><div><br></div><div>;*********************************************</div><div>zmpsi = zonal_mpsi(v1,v1&lat,v1&level(::-1),PsPa)</div></div><div>;*******************************************************************</div><div><br></div><div>I will be grateful if you can help me resolve this issue.</div><div><br></div><div>Best Regards</div><div>Ipshita</div><div class="gmail_signature"><div dir="ltr"><br></div></div>
</div>