<html><head></head><body><div style="font-family:times new roman, new york, times, serif;font-size:16px;"><div style="font-family:times new roman, new york, times, serif;font-size:16px;">Hi</div><div style="font-family:times new roman, new york, times, serif;font-size:16px;"><span><span id="ydp431a2f03result_box" class="ydp431a2f03short_text" lang="en"><span>I recently have made the ncl update to 6.4.0</span></span><span id="ydp431a2f03result_box" class="ydp431a2f03short_text" lang="en"><span></span></span></span> , i had to remove this line to run ncl scripts <br></div><div style="font-family:times new roman, new york, times, serif;font-size:16px;"><span><span>load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl</span></span><br></div><div style="font-family:times new roman, new york, times, serif;font-size:16px;">now it runs but there is this error message <br></div><div style="font-family:times new roman, new york, times, serif;font-size:16px;"><br></div><div style="font-family:times new roman, new york, times, serif;font-size:16px;"><br></div><div style="font-family:times new roman, new york, times, serif;font-size:16px;"><span>Undefined identifier: (wrf_times_c) is undefined, can't continue<br><br></span><div>i want extract data by day and hour</div><div>thanks for help<br></div><div><br></div><div>this is ncl scripts</div><div><br></div><div><span>begin<br> a = addfile("199803/wrfout_d01_1998-03-26_12:00:00.nc","r")<br><br> times = wrf_user_list_times(a) ; "2008-09-29_18:30:00", etc<br> <br> tk2 = wrf_user_getvar(a,"T2",-1) ; T2 in Kelvin<br> rain1 = wrf_user_getvar(a,"RAINC",-1) ; get t2 for all times<br> rain2 = wrf_user_getvar(a,"RAINNC",-1) ; get t2 for all times<br><br> <br> Times = wrf_user_getvar(a,"Times",-1)<br> Times = wrf_times_c(a->Times, 3)<br> <br> times = wrf_user_list_times(a)<br> ntimes = dimsizes(times)<br>; print("ntimes = " + ntimes)<br> <br><br>;---Calculate i,j locations of data closest to set of lat/lon points<br> <br> ip_locs = (/ "Borojurd", "Alighodarz" ,"Khoramabad" , "Dezful" , "Kohrang" /)<br> ip_lats = (/ 33.55, 33.24 ,33.26 , 32.24 ,32.26/)<br> ip_lons = (/ 48.45, 49.42 , 48.17 ,48.23 , 50.07/)<br><br> do ip = 0, 4 ; LOOP through above 20 station <br> opt = True<br> loc = wrf_user_ll_to_ij( a, ip_lons(ip), ip_lats(ip), opt)<br> loc = loc-1<br> locs = ip_locs(ip) <br> ; <br> t2t = tk2(:,loc(1),loc(0)) <br> raint1 = rain1(:,loc(1),loc(0))<br> raint2 = rain2(:,loc(1),loc(0))<br> rain3=raint1+raint2<br><br> slist = [/ " station " + " longitude " + " latitude " + "Times" + " temprature" + " rain1 " + "rain2 " + "rain3" /]<br> write_table("filename", "a", slist, "%s") <br> <br> alist=[/ip_locs(ip), ip_lons(ip), ip_lats(ip), times , t2t-273 , raint1 , raint2 , rain3 /]<br> write_table("filename","a" ,alist, "%s %8.2f %8.2f %s %8.2f %8.2f %8.2f %8.2f " )<br> <br> end do<br>end</span><br></div></div><div><br></div></div></body></html>