[ncl-talk] List subscripting used on non-list variable, can't continue

[Borja Sas González]

I don't see anything weird. The message prints that there is an error on
line 51:

lat1 = a[0]->XLAT(0,nlt,nln) ; nearest grid point

But I also see nothing wrong.

El mié., 7 ago. 2019 a las 13:01, Barry Lynn (<barry.h.lynn at gmail.com>)
escribió:

> Hi:
>
> I am not sure which variable is being printed in error, but you are trying
> to subscript (add indices) to a variable that doesn't have them.  Try to
> solve your problem by finding the variable.  You can printVarSummary to see
> the dimensions of each.
>
> Barry
>
> On Wed, Aug 7, 2019 at 2:40 PM Borja Sas González via ncl-talk <
> ncl-talk at ucar.edu> wrote:
>
>> I am trying to write a csv to represent several variables. To make it
>> easy, copy the example from the NCL page (write_csv_5.ncl) that represents
>> the temperature at 2 m.
>>
>> The only changes I made in the script are the output file of the WRF that
>> you are going to read and the map points that I want.
>>
>> When executing it, the following error is printed:
>>
>> fatal:List subscripting used on non-list variable, can't continue
>> fatal:["Execute.c":8637]:Execute: Error occurred at or near line 51 in
>> file csv.ncl
>>
>> The script is:
>>
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
>>
>> begin
>>
>>   a =
>> addfile("/home/bsas/wrf_wind/WRF/test/em_real/wrfout_d03_2019-08-07_00:00:
>> 00.nc","r")
>>
>>   times = wrf_user_list_times(a) ; "2019-08-01_12:00:00", etc
>>   tk2   = wrf_user_getvar(a,"T2",-1) ; T2 in Kelvin
>>
>>   times  = wrf_user_list_times(a)
>>   ntimes = dimsizes(times)
>>   print("ntimes = " + ntimes)
>>
>> ;---Calculate i,j locations of data closest to set of lat/lon points
>>   lats = (/28.290147,28.379555/)
>>   lons = (/-16.377643,-16.585276/)
>>   nlatlon = dimsizes(lats)
>>   loc = wrf_user_ll_to_xy(a, lons, lats, True)   ; 2 x nlatnlon
>>
>> ;  loc = wrf_user_ll_to_ij(a, lons, lats, True)   ; 2 x nlatnlon
>> ;  loc = loc - 1                                  ; wrf_user_ll_to_ij
>> returns 1-based indexes
>>
>> ;---Set up CSV file and header information for the file
>>   csv_filename = "wrf_2m_temperature.csv"
>>   system("rm -f " + csv_filename)                ; Remove file in case it
>> exists.
>>   fields = (/"TIME", "LAT", "LON", "TEMPERATURE (degC)"/)
>>
>> ;---Create a header line for CSV file
>>   dq     = str_get_dq()
>>   fields = dq + fields + dq                      ; Pre/append quotes to
>> field names
>>   header = [/str_join(fields,",")/]              ; Header is field names
>> separated
>>                                                  ; by commas.
>> ;
>> ; Format to use for writing each variable to CSV file.
>> ; If you don't want spaces in CSV file, use the following
>> ; format string:
>> ;     format = "%s,%g,%g,%g"
>> ;
>> format = "%s,%6.2f,%7.2f,%6.2f"
>>
>> ;
>> ; Loop through each time step and desired list of lat/lon values,
>> ; and write a single line of data to CSV file.
>> ;
>>   write_table(csv_filename, "w", header, "%s")   ; Write header to CSV
>> file.
>>   do it = 0,ntimes-1
>>     do nl = 0,nlatlon-1
>>       nln   = loc(0,nl)
>>       nlt   = loc(1,nl)
>>       lat1  = a[0]->XLAT(0,nlt,nln)    ; nearest grid point
>>       lon1  = a[0]->XLONG(0,nlt,nln)
>>       alist = [/times(it),lat1,lon1,tk2(it,nlt,nln)/]   ; Store data to
>> be written in a list.
>>       write_table(csv_filename, "a", alist, format)     ; Write list to
>> CSV file.
>>     end do
>>   end do
>> end
>> --
>> Borja Sas González
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
>
> --
> Barry H. Lynn, Ph.D
> Senior Associate Scientist, Lecturer,
> The Institute of the Earth Science,
> The Hebrew University of Jerusalem,
> Givat Ram, Jerusalem 91904, Israel
> Tel: 972 547 231 170
> Fax: (972)-25662581
>
> C.E.O, Weather It Is, LTD
> Weather and Climate Focus
> http://weather-it-is.com
> Jerusalem, Israel
> Local: 02 930 9525
> Cell: 054 7 231 170
> Int-IS: x972 2 930 9525
>
>

-- 
Borja Sas González
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190807/781143a8/attachment.html>