<div dir="ltr"><div><div><div>You can download the entire *beta* 6.4.0 library which includes: wgt_vertical_n<br><br><a href="http://www.cgd.ucar.edu/~shea/contributed.ncl_beta_640">http://www.cgd.ucar.edu/~shea/contributed.ncl_beta_640</a><br><br>===<br><br></div>In your NCL script:<br><br></div>load "/path/contributed.ncl_beta_640" ; NOTE: this is a *BETA* version<br><br>===<br></div>The wgt_vertical_n' function is attached as a separate entity.<br><br><br><div><br><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Jun 11, 2016 at 1:04 AM, Guido Cioni <span dir="ltr"><<a href="mailto:guidocioni@gmail.com" target="_blank">guidocioni@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word">Hey Erick,<div>if you need to integrate over a certain depth I believe you can write the function yourself. I did something similar in one of my scripts:</div><div><br></div><div><font style="font-size:11px" face="Menlo">;;;; first compute the pressure layers thicknesses ;;;</font></div><div><div><font style="font-size:11px" face="Menlo"> dp=new(dim sizes(height),typeof(p))</font></div><div><font style="font-size:11px" face="Menlo"> do i=0,dimsizes(time)-1</font></div><div><font style="font-size:11px" face="Menlo"> dp(i,:)=dpres_plevel(p(i,:),pres_sfc(i),p(i,0),0)</font></div><div><font style="font-size:11px" face="Menlo"> end do</font></div><div><font style="font-size:11px" face="Menlo"><br></font></div><div><font style="font-size:11px" face="Menlo">; then compute the integral by summing the contribution of every layer</font></div><div><font style="font-size:11px" face="Menlo">; note that NCL does the sum without looping, so everything should be fast</font></div><div><font style="font-size:11px" face="Menlo">; also with big arrays</font></div><div><div><font style="font-size:11px" face="Menlo">var_int=dim_sum_n(var*dp, 1)/g</font></div></div><div><br></div><div>I was using domain-averaged data, so with coordinates (time,z) but you can easily adapt the script for 2D (lat,lon) data. </div><div>Cheers</div><div>
<div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><br>Guido Cioni</div><div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><a href="http://guidocioni.altervista.org" target="_blank">http://guidocioni.altervista.org</a> </div>
</div>
<br><div><blockquote type="cite"><div><div class="h5"><div>On 11 Jun 2016, at 03:36, Erick Gordon <<a href="mailto:erick.gordon24@gmail.com" target="_blank">erick.gordon24@gmail.com</a>> wrote:</div><br></div></div><div><div><div class="h5"><div dir="ltr"><div>Hi</div><div>I need to know how I can calculate the integral of a variable in NCL. In the case of wgt_vertical_n function will only be available from version 6.4.0 and later. For example, I need to integrate the divergence of moisture flux from 1000 to 500 hPa. If anyone can help me I would be grateful</div><div>Thank you very much</div></div></div></div>
_______________________________________________<br>ncl-talk mailing list<br><a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><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></div></blockquote></div><br></div></div><br>_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>