<div dir="ltr">Dear NCL Users,<div>                           I am trying to extract wrf model output parameter from wrfout file corresponding to station latitude longitude from another csv file. If I'll mention latitude and longitude in ncl script as an array then script is writing data to csv files but if I'll try to read latitude and longitude from csv file then it's printing output on terminal but not writing anything to a file. </div><div><br></div><div>Thanks in advance for any help.</div><div><br></div><div>Below is ncl code:</div><div> <br></div><div>begin<br>stn_file = "./sample.csv"<br>f = asciiread(stn_file, -1, "string")      ; -1 is used to read variables as 1 dimensional<br> ; stn  = str_get_field(f4, 4, ",")<br> ; lats = tofloat(str_get_field(f4, 5, ","))<br> ; lons = tofloat(str_get_field(f4, 6, ","))<br> <br> stn    = (/"RRN","SUL","AUR","KAM","MAH","KAP","KOD","KAL","SIR","CHI","CHIT","HUB","CHA","KALG    ","AMI","ALN","KUN"/)<br> lats   = (/12.98,17.39,17.44,17.57,17.53,13.23,17.40,17.35,16.17,17.46,17.69,15.36,15.22,15.18,    15.54,15.42,15.25/)<br> lons   = (/77.51,77.34,76.89,76.98,76.88,74.75,77.19,77.15,77.02,77.42,77.21,75.13,75.15,74.97,    75.05,74.74,75.24/)<br>  <br> ext   = "10_1"                                          ; experiment name<br> lu    = "mod"                                             ; Landuse datafor WRF<br><br></div><div> ; Read the data<br>fo  = addfiles( "wrfout_d03_2017-08-14_00:00:00","r")<br><br>wlat=fo[:]->XLAT(0,:,:)</div><div>wlon=fo[:]->XLONG(0,:,:)<br>rainc     = wrf_user_getvar(fo,"RAINC",-1)<br>time  = wrf_user_getvar(fo,"times",-1)<br>print(time)<br>ntime    = dimsizes(time)<br>nl = dimsizes(lats)<br>ln = nl-1<br><br></div><div>do it=0,nl-1<br>olat = lats(it)<br>olon = lons(it)<br> yl = olat<br> xl = olon<br><br></div><div>; Find the closest point to the site in WRF grids<br> obsij=wrf_latlon_to_ij(wlat,wlon,olat,olon)<br><br> print(obsij)<br> </div><div>; Use the Index to extract data<br> x  = obsij(0)<br> y  = obsij(1)<br>  <br> lm = ntime-1<br> rainCu = rainc(0:lm,x,y)<br> print(stn(it)+"       "+lats(it)+"    "+ lons(it)+"   "+ rainCu)<br><br></div><div> ;;;;Write to csv file<br><br></div><div>csv_filename = "RainC_" +stn(it)+ "_WRF_" +ext+ ".csv"<br>format = "%s  %5.2f   %5.2f   %5.2f"<br>fields =(/"STATION",  "LAT",  "LON",  "RAINC"/)<br>header = [/str_join(fields, " ")/]<br>write_table(csv_filename, "w", header, "%s")<br><br></div><div>alist = [/stn(it), lats(it), lons(it), rainCu/]<br>write_table(csv_filename, "a", alist, format)<br><br></div><div>end do<br>end<br><div><br></div>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div><b>Thanks & Regards,<br></b></div><b>Ajay Bankar.<br></b></div><b><br></b></div></div></div></div>