<div dir="ltr"><div>Can you use the 1d interp routine instead of the 3d routine.  Interpolating the whole array for one grid point seems unnecessary. </div><div>wrf_interp_1d(rh(:,60,41), p(:,60,41), (/500,850,1000/) )</div><div><br></div><div><div>Assuming its the interpolation lines slowing you down.  </div><div><br></div></div><div><br></div><div>Good luck, </div><div>    ~Alan</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Oct 13, 2014 at 2:25 PM, Saurabh Singh <span dir="ltr">&lt;<a href="mailto:saurabhsingh123op@gmail.com" target="_blank">saurabhsingh123op@gmail.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">hello all<div><br></div><div>i have a script to extract variables from WRF file at a certain pressure level</div><div>there are around 20 variables that i need to extract , but it takes really long to run the scrip,i am sending the code for 1 variable please if any one can help me to optimise it</div><div><br></div><div>this scripy tries to extract point data for a wrf file generated every 3 hours and total duration is a year.</div><div><br></div><div><div>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl&quot;</div><div>    load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl&quot;</div><div><br></div><div>  files = systemfunc(&quot;ls -1  /home/cee1/Music/Build_WRF_and_WPS.V351/WRF_and_WPS_src/WRFV3/test/ em_real/2k/wrfout_d01_2000*&quot;) + &quot;.nc&quot;</div><div>   a = addfiles(files,&quot;r&quot;)</div><div><div> ;extraction of variables</div><div>   time = -1</div><div>   p  = wrf_user_getvar(a,&quot;pressure&quot;,time) ; total pressure<br></div><div><br></div><div>  rh = wrf_user_getvar(a,&quot;rh&quot;,time)   ; calculate RH</div></div><div>  rh850 =wrf_user_intrp3d( rh,p,&quot;h&quot;,850,0.,False)</div><div>  rh500 =wrf_user_intrp3d( rh,p,&quot;h&quot;,500,0.,False)</div><div>  rh1000 =wrf_user_intrp3d( rh,p,&quot;h&quot;,1000,0.,False)</div><div><br></div><div><div>  rhpoint500 = rh500(:,60,41)             </div><div>  rhpoint850 = rh850(:,60,41)           </div><div>  rhpoint1000 = rh1000(:,60,41)  </div><div> </div></div><div><div>   asciiwrite(&quot;rh0.txt&quot;,rhpoint500)</div><div>   asciiwrite(&quot;rh8.txt&quot;,rhpoint850)</div><div>   asciiwrite(&quot;rh13.txt&quot;,rhpoint1000)</div><div>   </div></div><div>i have the similar code inside for 19 other variables,please help me to optimize this script and i just need my output at a single station and 3 pressure levels i.e 500.850 and 1000</div><span class="HOEnZb"><font color="#888888"><div><br></div>-- <br><div dir="ltr"><font style="color:rgb(0,0,0);font-family:arial,helvetica,sans-serif"><span>With regards</span><span></span></font><font style="font-family:arial,helvetica,sans-serif"><span style="color:rgb(0,0,0)"><span><br></span></span></font><font style="font-family:arial,helvetica,sans-serif"><span style="color:rgb(0,0,0)"><span>Saurabh Kumar singh</span><span></span> </span></font><br style="color:rgb(51,204,0)"><br><font style="color:rgb(102,102,102)" size="1"> <b><i><span style="font-family:Webdings" lang="EN-GB">P</span></i></b><b><i><span style="font-family:Tahoma,sans-serif" lang="EN-GB"> : </span></i></b><b><i><span style="font-family:Arial,sans-serif" lang="EN-GB"> </span></i></b><b><i><span style="font-family:Tahoma,sans-serif" lang="EN-GB">Please consider the environment before printing this e-mail</span></i></b> </font><br style="color:rgb(102,102,102)"><br><br><br></div>
</font></span></div></div>
<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>
<br></blockquote></div><br></div>