<html><body><div style="color:#000; background-color:#fff; font-family:HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:12pt"><div class="" style="">Hi all,</div><div class="" style=""><br class="" style=""></div><div class="" style="">I am trying to calculate Vertical Integrated Moisture Transport (VIMT) from surface to 500 hPa in units of kg/m/sec using grads function "vint(psexpr,expr,top)" only on the one time step.</div><div class="" style="">I used the following script to get VIMT:</div><div class="" style=""><br class="" style=""></div><div class="" style=""><span style="font-size: 12pt;">'sdfopen spr-daily.nc' ;* spr is the surface pressure in Pa.</span><br></div><div class="" style="">'set time 20jun1970'</div><div class="" style="">'define surfpres=pres/100' ;* /100 to convert in hPa</div><div class="" style="">'close 1'</div><div class="" style=""><br class=""
style=""></div><div class="" style="">'sdfopen shum-daily.nc'</div><div class="" style="">'sdfopen uwnd-daily.nc'</div><div class="" style="">'sdfopen vwnd-daily.nc'</div><div class="" style=""><br class="" style=""></div><div class="" style="">'set lat 0 45'</div><div class="" style="">'set lon 50 100'</div><div class="" style=""><br class="" style=""></div><div class="" style="">'set time 20jun1970'</div><div class="" style="">'define vintq=vint(surfpres,shum.1,500)' ;*shum is the specific humidity in kg/kg</div><div class="" style="">'define vintu=vint(surfpres,uwnd.2,500)' ;*uwnd is the u component of wind in m/sec</div><div class="" style="">'define vintv=vint(surfpres,vwnd.3,500)' ;*vwnd is the v component of wind in m/sec</div><div class="" style="">'define zon=(vintq*vintu)'</div><div class="" style="">'define mer=(vintq*vintv)'</div><div class="" style=""><br class="" style=""></div><div class=""
style="">'set gxout vector'</div><div class="" style="">'d zon;mer'</div><div class="" style=""><br class="" style=""></div><div class="" style="">Script runs successfully showing VIMT, but the magnitude of vector appears (4e+06) in the vector scale bar.</div><div class="" style="">I guess, this magnitude is not right. The value should lie around 20 or 40 (kg/m/sec). </div><div class="" style=""><br class="" style=""></div><div class="" style="">I'll be very thankful if anyone could help me to sort out the issue mentioned above.</div><div class="" style=""><br class="" style=""></div><div class="" style="">Regards</div><div class="" style=""><br class="" style=""></div><div class="" style="">Muhammad Latif</div><div class="" style="">PhD Scholar</div><div class="" style="">Department of Meteorology</div><div class="" style="">CIIT, Islamabad </div><div class="" style=""><br class="" style=""></div></div></body></html>