<div dir="ltr">; WRF_input_worstcase.csv<br>; 1st 3 lines<br>; cell ID,Longtiude,Latitude,H1,H2,H3,H4,H5,H6,H7,H8,H9,H10,H11,H12,H13,H14,H15,H16,H17,H18,H19,H20,H21,H22,H23,H24<br>; 0,103.551796,1.187599182,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0<br>; 1,103.551796,1.190299988,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0<br><br> dir_csv = "./"<br> fil_csv = "WRF_input_worstcase.csv"<br> pth_csv = dir_csv+fil_csv<br><br> ncol = 28<br> data = <a href="http://www.ncl.ucar.edu/Document/Functions/Contributed/readAsciiTable.shtml"><b>readAsciiTable</b></a>(pth_csv, ncol, "float", 1) ; <== maybe double wood be better<br> printVarSummary(data)<br><br> cell = <b>toint</b>( data(:,0) )<br> LON = data(:,1)<br> LAT = data(:,2)<br> H1 = data(:,3)<br>; :<br> H24 = data(:,27)<br><br><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, May 5, 2021 at 2:11 AM ali mughal via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu">ncl-talk@mailman.ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><br></div><div dir="ltr"><div style="margin:0px;padding:0px 0px 20px;width:1119.2px;font-family:Roboto,RobotoDraft,Helvetica,Arial,sans-serif;font-size:medium"><div id="gmail-m_-3938406879847162614gmail-:123"><div></div><div></div><div></div></div><div><div></div><div id="gmail-m_-3938406879847162614gmail-:155"></div><div id="gmail-m_-3938406879847162614gmail-:121" style="margin:8px 0px 0px;padding:0px"><div id="gmail-m_-3938406879847162614gmail-:122"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr">Dear All <div><br></div><div>I am using the following script to convert the data contained in this (<a href="https://we.tl/t-lS4dhtyDs4" target="_blank">https://we.tl/t-lS4dhtyDs4</a>) csv file to add as a 3 dimensional variable to <a href="http://geo_em.d01.nc/" target="_blank">geo_em.d01.nc</a>.</div><div><br></div><div><div>filename = "WRF_input_worstcase.csv"</div><div>a=addfile("../wrfout_d05_2016-03-31_00:00:<a href="http://00.nc/" target="_blank">00.nc</a>","w") ; I am adding this file only to import structure of variable with W/m^2</div><div>ac_heat=a->SF_AC_URB3D<br></div><div>lines = asciiread(filename,-1,"string")<br></div><div>nlines = dimsizes(lines)-1</div><div>ncols = dimsizes(str_split(lines(0),","))<br></div><div>nrows = dimsizes(lines)</div><div><br></div><div>print("This file has " + nrows + " rows and " + ncols + " columns.")<br></div><div><br></div><div> This file has 27091 rows and 27 columns.<br></div><div><br></div><div> lat = str_get_field(lines , 2, ",")</div><div> lon = str_get_field(lines , 1, ",")</div><div> </div><div> lat := stringtofloat(lat)</div><div> lon := stringtofloat(lon)</div><div><br></div><div> sh= traffic=readAsciiTable(filename, 24, "float", 1) <b>I want to read 24 columns here from H1-H24 but instead it selects the first 24 columns. The structure of the file should be latxlonx24</b></div><div> f = addfile("<a href="http://geo_em.d01.nc/" target="_blank">geo_em.d01.nc</a>","w")</div><div><br></div><div> sh@long_name = "sensible heat flux"</div><div> sh@short_name = "heat flux"</div><div> sh@units = "W/m^2"</div><div> filevardef(a, "SH_FLUX", typeof(ac_heat), getvardims(ac_heat))</div><div> f->SH_FLUX = sh</div><div><br></div></div><div>I would be really grateful for help.</div></div></div></div></div></div></div></div></div></div></div></div></div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="https://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">https://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div>