<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 12pt;
font-family:Calibri
}
--></style></head>
<body class='hmmessage'><div dir='ltr'><div>Dear Mary</div><div>thank you very much for yoru response. i run the revised script from you but it is again giving output as a single ascii ouput as 18T1.... how could i get 19 individual stations files like 1<span style="font-family: monospace, monospace; font-size: 12pt;">T2oben_1kmtemp_rainnonc_rainc_raintot_wspd_wdir_snow_snowh.txt", 2</span><span style="font-family: monospace, monospace; font-size: 12pt;">T2oben_1kmtemp_rainnonc_rainc_raintot_wspd_wdir_snow_snowh.txt etc.</span></div><div><span style="font-family: monospace, monospace; font-size: 12pt;">regards</span></div><div><span style="font-family: monospace, monospace; font-size: 12pt;">ram</span></div><div><br></div><div><br></div><div><br><div dir="ltr"><div class="ecxgmail_default" style="font-size:small;">Please respond back to ncl-talk if you have follow-up questions. </div><div class="ecxgmail_default" style="font-size:small;"><br></div><div class="ecxgmail_default" style="font-size:small;">You have two sets of double do loops, where the first set of loops reads in the data, and the second loop tries to write it.</div><div class="ecxgmail_default" style="font-size:small;"><br></div><div class="ecxgmail_default" style="font-size:small;">The problem is that all the data you're reading in the first set of loops gets reset every time, so by the time you get to your second set of do loops, you don't have any of the original data left to write.</div><div class="ecxgmail_default" style="font-size:small;"><br></div><div class="ecxgmail_default" style="font-size:small;">You want to use just one set of do loops, so that you read in your data, and then immediately write it to a file before you then read the next set of data.</div><div class="ecxgmail_default" style="font-size:small;"><br></div><div class="ecxgmail_default" style="font-size:small;">I moved all your 'wrf_user_getvar' calls to the outside do loop, because they only change for every value of "it" and not for every value of "ip". It's time-consuming to have code in a loop that you don't need, so it's important to keep your "do" loops as clutter free as possible.</div><div class="ecxgmail_default" style="font-size:small;"><br></div><div class="ecxgmail_default" style="font-size:small;">This is UNTESTED:</div><div class="ecxgmail_default" style="font-size:small;"><br></div><div class="ecxgmail_default" style=""><div class="ecxgmail_default" style=""><br></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> ; -------------- LOAD FUNCTIONS AND PROCEDURES ----------------</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"><br></font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"><br></font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> ; -------------- BEGINING OF NCL SCRIPT ----------------</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace">begin</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace">DATADir = "/mnt/gpfs1/data/private/fgvy024/Build_WRF/WRFV6/run/"</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace">FILES = systemfunc (" ls -1 " + DATADir + "wrfout_d03* ")</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace">numFILES = dimsizes(FILES)</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace">a = addfiles(FILES+".nc","r")</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace">times = wrf_user_list_times(a) ; get times in the file</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace">ntimes = dimsizes(times) ; number of times in the file or 1584 for domain1 4750</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"><br></font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace">;new parameters naming</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace">temp = new(ntimes,float) </font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace">rain_expp = new(ntimes,float) </font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace">rain_conn = new(ntimes,float) </font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace">rain_tot = new(ntimes,float)</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace">wspd = new(ntimes,float)</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace">wdir = new(ntimes,float)</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace">snoww_eq = new(ntimes,float)</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace">snowh_h = new(ntimes,float)</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace">;rhh = new(ntimes,float)</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"><br></font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace">ip_lats = (/ 27.3081, 27.6331,27.86,27.959167, \</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> 27.696667, 27.8024,27.8953, 27.99,27.9818, \ </font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> 27.90506667, 27.89671667, 27.901, 27.8987, \</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> 27.89261111, 27.87828333, 27.85911111, \</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> 27.60111, 27.74661, 27.85722 /)</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace">ip_lons = (/ 86.5042, 86.2331, 86.46, 86.81278, 86.721389, 86.7144, \</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> 86.8188, 86.83, 86.7649, 86.37536667, 86.37428333, 86.37625, 86.37911667, \</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> 86.37594444,86.43361667, 86.43388889, 86.74028, 86.713, 86.79417 /)</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace">res = True </font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace">res@returnInt = True </font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"><br></font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace">data = new(ntimes, "string") </font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"><br></font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace">;******************************* time loop</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace">do it = 0, ntimes-1 </font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> tc2 = wrf_user_getvar(a,"T2",it) ; temp extraction</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> tc2 = tc2 -273.16</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> rain_exp = wrf_user_getvar(a,"RAINNC",it) </font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> rain_con = wrf_user_getvar(a,"RAINC",it) </font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> u10 = wrf_user_getvar(a,"U10",it) ; u at 10 m, mass point</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> v10 = wrf_user_getvar(a,"V10",it) ; v at 10 m, mass point</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> snow_eq = wrf_user_getvar(a,"SNOW",it)</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> snow_h = wrf_user_getvar(a,"SNOWH",it)</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"><br></font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> ;************************************************</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> do ip = 0, 18 ;loop for number of station locations</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> point = wrf_user_ll_to_ij(a, ip_lons(ip), ip_lats(ip), True)</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> x = point(0) - 1</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> y = point(1) - 1 </font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"><br></font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> temp(it) = tc2(y,x) ; temperature extraction for station location</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> rain_expp(it) = rain_exp(y,x)</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> rain_conn(it) = rain_con(y,x)</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> rain_tot(it) = rain_exp(y,x) + rain_con(y,x)</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> wspd(it) = ndtooned( sqrt(u10(y,x)^2 + v10(y,x)^2) )</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> wdir(it) = ndtooned( atan2(u10(y,x),v10(y,x))/0.01745329 +180. ) </font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> snoww_eq(it) = snow_eq(y,x)</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> snowh_h (it) = snow_h(y,x)</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"><br></font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace">;---Create filename to write to</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> fName = ip+"T2oben_1kmtemp_rainnonc_rainc_raintot_wspd_wdir_snow_snowh.txt"</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"><br></font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> data (it)= (sprintf("%5.0f",it) +" " \</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> +sprintf("%21.2f", temp(it)) +" " \</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> +sprintf("%22.2f", rain_expp(it)) +" " \</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> +sprintf("%23.2f", rain_conn(it)) +" " \</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> +sprintf("%24.2f", rain_tot(it)) +" " \</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> +sprintf("%17.2f", wspd(it)) +" " \ </font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> +sprintf("%25.2f", wdir(it)) +" " \ </font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> +sprintf("%18.2f", snoww_eq(it)) +" " \</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> +sprintf("%19.2f", snowh_h(it)) +" " )</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> end do ; end of "ip" loop</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace"> asciiwrite (fName , data)</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace">end do ; end of time loop</font></div><div class="ecxgmail_default" style=""><font face="monospace, monospace">en</font>d</div><div style="font-size:small;"><br></div></div></div><div class="ecxgmail_extra"><br><div class="ecxgmail_quote">On Mon, May 9, 2016 at 10:23 AM, Ramchandra Karki <span dir="ltr"><<a href="mailto:rammetro@hotmail.com" target="_blank">rammetro@hotmail.com</a>></span> wrote:<br><blockquote class="ecxgmail_quote" style="border-left:1px #ccc solid;padding-left:1ex;">
<div><div dir="ltr">Dear Mary<div>thank you very much. <span style="font-size:12pt;">i am new to these things. the script is not giving errors but i am not getting what i am expecting. i want to get </span></div><div>-individual ascii file outputs for each stations (19 ascii files for 19 stations containing time series for whole period). would you please kindly suggest me where i need to change to get it?<br><div><div><br><i>Regards<br></i><div><i>ram</i></div><br><div class="ecxhm ecxHOEnZb"><br></div><div><div class="ecxhm ecxHOEnZb"><hr>Date: Mon, 9 May 2016 10:13:48 -0600<br>Subject: Re: [ncl-talk] multiple ascii file outputs<br>From: <a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a><br>To: <a href="mailto:rammetro@hotmail.com" target="_blank">rammetro@hotmail.com</a><br>CC: <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a></div><div><div class="h5"><br><br><div dir="ltr"><div>It helps if you provide more information about what the problem is. Are you getting an error of any kind? If so, please indicate which line the error occurs on.</div><div><br></div><div>I think the problem is here:</div><div><br></div><div><div style="line-height:21.3px;color:rgb(68,68,68);"><font face="monospace, monospace"> do ip = 0, 18</font></div><div style="line-height:21.3px;color:rgb(68,68,68);"><font face="monospace, monospace"> do it = 0,ntimes-1 ;ntimes-1</font></div><div style="line-height:21.3px;color:rgb(68,68,68);"><font face="monospace, monospace"> data (it)= (sprintf("%5.0f",it) +" " \</font></div><div style="line-height:21.3px;color:rgb(68,68,68);"><font face="monospace, monospace"> +sprintf("%21.2f", temp(it)) +" " \</font></div><div style="line-height:21.3px;color:rgb(68,68,68);"><font face="monospace, monospace"> +sprintf("%22.2f", rain_expp(it)) +" " \</font></div><div style="line-height:21.3px;color:rgb(68,68,68);"><font face="monospace, monospace"> +sprintf("%23.2f", rain_conn(it)) +" " \</font></div><div style="line-height:21.3px;color:rgb(68,68,68);"><font face="monospace, monospace"> +sprintf("%24.2f", rain_tot(it)) +" " \</font></div><div style="line-height:21.3px;color:rgb(68,68,68);"><font face="monospace, monospace"> +sprintf("%17.2f", wspd(it)) +" " \ </font></div><div style="line-height:21.3px;color:rgb(68,68,68);"><font face="monospace, monospace"> +sprintf("%25.2f", wdir(it)) +" " \ </font></div><div style="line-height:21.3px;color:rgb(68,68,68);"><font face="monospace, monospace"> +sprintf("%18.2f", snoww_eq(it)) +" " \</font></div><div style="line-height:21.3px;color:rgb(68,68,68);"><font face="monospace, monospace"> +sprintf("%19.2f", snowh_h(it)) +" " )</font></div><div style="line-height:21.3px;color:rgb(68,68,68);"><font face="monospace, monospace">end do</font></div><div style="line-height:21.3px;color:rgb(68,68,68);"><font face="monospace, monospace"> end do ; end of time loop</font></div><div style="line-height:21.3px;color:rgb(68,68,68);"><font face="monospace, monospace"><br></font></div><div style="line-height:21.3px;color:rgb(68,68,68);"><font face="monospace, monospace"> asciiwrite (fName , data)</font></div><div style="line-height:21.3px;color:rgb(68,68,68);"><br></div><div style="line-height:21.3px;color:rgb(68,68,68);">You have a double "do" loop, where "data" is getting created inside the inner do loop. For each iteration of "ip", then, the previous strings assigned to "data" are getting clobbered. In fact, I'm not even sure why you are looping across "ip=0,18", because the inside loop is doing the exact same thing every time.</div><div style="line-height:21.3px;color:rgb(68,68,68);"><br></div><div style="line-height:21.3px;color:rgb(68,68,68);">Finally, you call "asciiwrite", but this is only going to contain "data" from the last iteration of "ip" (ip=18).</div><div style="line-height:21.3px;color:rgb(68,68,68);">--Mary</div><div style="line-height:21.3px;color:rgb(68,68,68);"><br></div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;"><br></div></div></div><div><br><div>On Sat, May 7, 2016 at 8:52 AM, Ramchandra Karki <span dir="ltr"><<a href="mailto:rammetro@hotmail.com" target="_blank">rammetro@hotmail.com</a>></span> wrote:<br><blockquote style="border-left:1px #ccc solid;padding-left:1ex;">
<div><div dir="ltr"><div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"> Dear NCL help,</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">I am trying to extract multiple stations time series (ascci output for each stations) from wrfoutput file but i could not succeed in doing so. Would you please kindly guide me to get the solution?</div></div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"><br></div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"> ; -------------- LOAD FUNCTIONS AND PROCEDURES ----------------</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"><br style="line-height:21.3px;"></div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"><br style="line-height:21.3px;"></div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"> ; -------------- BEGINING OF NCL SCRIPT ----------------</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">begin</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">DATADir = "/mnt/gpfs1/data/private/fgvy024/Build_WRF/WRFV6/run/"</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">FILES = systemfunc (" ls -1 " + DATADir + "wrfout_d03* ")</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">numFILES = dimsizes(FILES)</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">a = addfiles(FILES+".nc","r")</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">times = wrf_user_list_times(a) ; get times in the file</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">ntimes = dimsizes(times) ; number of times in the file or 1584 for domain1 4750</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"><br style="line-height:21.3px;"></div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">;new parameters naming</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">temp = new(ntimes,float) </div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">rain_expp = new(ntimes,float) </div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">rain_conn = new(ntimes,float) </div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">rain_tot = new(ntimes,float)</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">wspd = new(ntimes,float)</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">wdir = new(ntimes,float)</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">snoww_eq = new(ntimes,float)</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">snowh_h = new(ntimes,float)</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">;rhh = new(ntimes,float)</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"><br style="line-height:21.3px;"></div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">;******************************* time loop</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">do it = 0, ntimes-1 </div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">time = it</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"><br style="line-height:21.3px;"></div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"> ;************************************************</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"><span style="line-height:22.72px;font-size:12pt;">do ip = 0, 18 ;loop for number of station locations</span></div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">ip_lats = (/ 27.3081, 27.6331,<span style="line-height:21.3px;white-space:pre-wrap;"> </span>27.86,<span style="line-height:21.3px;white-space:pre-wrap;"> </span>27.959167, \</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"> 27.696667, 27.8024,<span style="line-height:21.3px;white-space:pre-wrap;"> </span>27.8953, 27.99,<span style="line-height:21.3px;white-space:pre-wrap;"> </span>27.9818, \ </div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"> 27.90506667, 27.89671667, 27.901, 27.8987, \</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"> 27.89261111, 27.87828333, 27.85911111, \</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"> 27.60111, 27.74661, 27.85722 /)</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">ip_lons = (/ 86.5042, 86.2331, 86.46, 86.81278, 86.721389, 86.7144, \</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"> 86.8188, 86.83, 86.7649, 86.37536667, 86.37428333, 86.37625, 86.37911667, \</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"><span style="line-height:21.3px;white-space:pre-wrap;"> </span> 86.37594444,<span style="line-height:21.3px;white-space:pre-wrap;"> </span>86.43361667, 86.43388889, 86.74028, 86.713, 86.79417 /)</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">res = True </div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">res@returnInt = True </div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">point = wrf_user_ll_to_ij(a, ip_lons(ip), ip_lats(ip), True)</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">x = point(0) - 1</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">y = point(1) - 1 </div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"><br style="line-height:21.3px;"></div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"><span style="line-height:22.72px;font-size:12pt;">tc2 = wrf_user_getvar(a,"T2",it) ; temp extraction</span></div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">tc2 = tc2 -273.16</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">temp(it) = tc2(y,x) ; temperature extraction for station location</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">rain_exp = wrf_user_getvar(a,"RAINNC",it) </div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">rain_con = wrf_user_getvar(a,"RAINC",it) </div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">rain_expp(it) = rain_exp(y,x)</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">rain_conn(it) = rain_con(y,x)</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">rain_tot(it) = rain_exp(y,x) + rain_con(y,x)</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">u10 = wrf_user_getvar(a,"U10",it) ; u at 10 m, mass point</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">v10 = wrf_user_getvar(a,"V10",it) ; v at 10 m, mass point</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">wspd(it) = ndtooned( sqrt(u10(y,x)^2 + v10(y,x)^2) )</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">wdir(it) = ndtooned( atan2(u10(y,x),v10(y,x))/0.01745329 +180. ) </div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"><span style="line-height:22.72px;font-size:12pt;">snow_eq = wrf_user_getvar(a,"SNOW",it)</span></div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">snow_h = wrf_user_getvar(a,"SNOWH",it)</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">snoww_eq(it) = snow_eq(y,x)</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">snowh_h (it) = snow_h(y,x)</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"><br style="line-height:21.3px;"></div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">end do ; end of station loop</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">end do ; end of time loop</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"><br style="line-height:21.3px;"></div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"> npts=ntimes</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"> fName = ip+"T2oben_1kmtemp_rainnonc_rainc_raintot_wspd_wdir_snow_snowh.txt"</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"> data = new( npts, "string") </div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"> ;print(" Time temp ")</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"> do ip = 0, 18</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"> do it = 0,ntimes-1 ;ntimes-1</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"> data (it)= (sprintf("%5.0f",it) +" " \</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"> +sprintf("%21.2f", temp(it)) +" " \</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"> +sprintf("%22.2f", rain_expp(it)) +" " \</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"> +sprintf("%23.2f", rain_conn(it)) +" " \</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"> +sprintf("%24.2f", rain_tot(it)) +" " \</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"> +sprintf("%17.2f", wspd(it)) +" " \ </div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"> +sprintf("%25.2f", wdir(it)) +" " \ </div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"> +sprintf("%18.2f", snoww_eq(it)) +" " \</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"> +sprintf("%19.2f", snowh_h(it)) +" " )</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">end do</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"> end do ; end of time loop</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"><br></div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"> asciiwrite (fName , data)</div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);"><br></div><div style="line-height:21.3px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255);">end</div><div><br></div><div><i><br></i></div> </div></div>
<br>_______________________________________________<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" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div></div></div></div></div></div></div> </div></div>
</blockquote></div><br></div></div> </div></body>
</html>