<div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><br></div><div dir="ltr"><div class="gmail-gs" style="margin:0px;padding:0px 0px 20px;width:1119.2px;font-family:Roboto,RobotoDraft,Helvetica,Arial,sans-serif;font-size:medium"><div id="gmail-:123"><div class="gmail-qQVYZb"></div><div class="gmail-utdU2e"></div><div class="gmail-btm"></div></div><div class="gmail-"><div class="gmail-aHl"></div><div id="gmail-:155" tabindex="-1"></div><div id="gmail-:121" class="gmail-ii gmail-gt" style="margin:8px 0px 0px;padding:0px"><div id="gmail-:122" class="gmail-a3s gmail-aiL"><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">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>