Dear Dennis <div><br></div><div>Thank you for your kind response.</div><div><br></div><div>When replying I thought I have copied ncl_talk but was not in cc.</div><div><br></div><div>Thank you for the code it solves the problem.</div><div><br></div><div>If you allow can I ask a related question in this mail or should I write a separate email with a new subject.</div><div><br></div><div><br></div><div><br>On Monday, May 10, 2021, Dennis Shea <<a href="mailto:shea@ucar.edu">shea@ucar.edu</a>> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>The utility 'dos2unix' must be used on the text file from excel.</div><div>----</div><div>%> ncl AliMughal.WRF_csv.ncl | less<br></div><div>----</div><div><br></div><div>Please ***ALL*** comments should go through ncl_talk.</div><div><br></div><div>Otherwise, I become your personal programmer.</div><div><br></div><div>Good luck<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, May 7, 2021 at 7:40 PM Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@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>[1] "excel file" .... I think there are special characters embedded in the csv file. Possibly 'dos' characters. <br></div><div> ****These must be removed before anything will work.**** This is your responsibility<br></div><div><br></div><div>[2] There are 27091 lines in the file. one header + 27090 lines with data</div><div><br></div><div>[3] data = <a href="http://www.ncl.ucar.edu/Document/Functions/Contributed/readAsciiTable.shtml" target="_blank"><b>readAsciiTable</b></a>(pth_csv, ncol, "double", 1) <br></div><div><br></div><div>[4] This looks like gridded data. If so, do you know the dimensions: (nlat,mlon) ?</div><div><br></div><div>possibly: nlat=210, mlon=129 ?????</div><div><br></div><div> N = 27090</div><div> lat = LAT(0:nlat-1)</div><div> lon = LON(::mlon) ; same as (0:N-1:mlon)<br></div><div><br></div><div> lat!0 = "lat"</div><div> lon!0 = "lon"</div><div><br></div><div> lat@long_name = "latitude"</div><div> lon:@ong_name = "longitude"</div><div><br></div><div> lat@units = "degrees_north" :</div><div> lon@units = "degrees_east"<br></div><div><br></div><div> delete([/LAT, LON/])</div><div><br></div><div> print(lat)</div><div> print(lon)<br></div><div><br></div><div> h1 = <a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/onedtond.shtml" target="_blank"><b>onedtond</b></a>(tofloat(H1),(/nlat,<wbr>mlon/))</div><div> h1@long_name = "..."</div><div> h1@units = "..."</div><div> h1!0 = "lat"</div><div> h1!1 = "lon"</div><div> <br></div><div>=====</div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, May 7, 2021 at 3:18 AM ali mughal <<a href="mailto:mughalali655@gmail.com" target="_blank">mughalali655@gmail.com</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"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr">Dear Dennis<div><br></div><div>Thank you for the reply.</div><div><br></div><div>I am getting two warning when executing the above lines </div><div><br></div><div>One is after the line </div><div><br></div><div>cell = <b>toint</b>( data(:,0) )<br></div><div><br></div><div><div>warning:toint: there are 967 doubles larger than INT_MAX, which has been flagged missing.</div></div><div><br></div><div>and the second one comes after the line H24<span style="white-space:pre-wrap"> </span>=data(:,26)</div><div>warning:asciiread: End of file reached and only (731430) elements were read from the file, filling remaining elements with the default missing value for the req <br></div><div><br></div><div>Apart from these warnings the data structure is different to what I want. Each of the variables shows similar dimensions as below. </div><div><div>Variable: LAT</div><div>Type: double</div><div>Total Size: 216720 bytes</div><div> 27090 values</div><div>Number of Dimensions: 1</div><div>Dimensions and sizes: [27090]</div><div>Coordinates:</div><div>Number Of Attributes: 1</div><div> _FillValue : 9.969209968386869e+36</div></div><div><br></div><div>The new sh variable which I intend to create should have the dimensions (LAT,LON, H-H24) similar to the excel file. </div><div><br></div><div><br></div><div>Still I cannot add the variable sh to netcdf file. I have tried to add H values one by one to the variable sh using the following command but then I get the following error.</div><div><br></div><div><div>sh=(/(/LAT/),(/LON/), (/H1/)/)</div></div><div><br></div><div><div>fatal:["NclFile.c":432]:<wbr>FileAddVar: an error occurred while adding a variable to a file, check to make sure data type is supported by the output format</div><div>fatal:["Execute.c":8578]:<wbr>Execute: Error occurred at or near line 43 in file test.ncl</div></div><div><br></div><div>Could you kindly help. So sorry for the trouble.</div><div><br></div><div><br></div><div><br></div></div></div></div></div></div></div></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, May 5, 2021 at 9:17 PM Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@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">; WRF_input_worstcase.csv<br>; 1st 3 lines<br>; cell ID,Longtiude,Latitude,H1,H2,<wbr>H3,H4,H5,H6,H7,H8,H9,H10,H11,<wbr>H12,H13,H14,H15,H16,H17,H18,<wbr>H19,H20,H21,H22,H23,H24<br>; 0,103.551796,1.187599182,0,0,<wbr>0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,<wbr>0,0,0,0,0,0,0<br>; 1,103.551796,1.190299988,0,0,<wbr>0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,<wbr>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" target="_blank"><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" target="_blank">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><div></div><div></div><div></div></div><div><div></div><div></div><div style="margin:8px 0px 0px;padding:0px"><div><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-<wbr>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"<wbr>)<br></div><div>nlines = dimsizes(lines)-1</div><div>ncols = dimsizes(str_split(lines(0),",<wbr>"))<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(<wbr>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>
______________________________<wbr>_________________<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/<wbr>mailman/listinfo/ncl-talk</a></blockquote></div>
</blockquote></div>
</blockquote></div>
</blockquote></div>
</blockquote></div>