[ncl-talk] Station obs overlay on map
Debasish Hazra
debasish.hazra5 at gmail.com
Mon Jul 20 11:19:56 MDT 2026
Thanks Dave, attached the revised ncl code with print statements you have
suggested along with print output and error in text message.
On Mon, Jul 20, 2026 at 12:51 PM Dave Allured - NOAA Affiliate <
dave.allured at noaa.gov> wrote:
> Debasish, please check your subscripting. Immediately before the line
> that gets the error, add printVarSummary for lat_new, obs_lat, and
> "indexes". Also include:
>
> printMinMax (indexes, 0)
> print ("num missing = " + num (ismissing (indexes)))
> print (i)
> print (npts_range)
>
> Show us all of this output immediately preceding the error, including the
> error message itself. Do not show repeated outputs, only the last set of
> debug print messages before the error occurs.
>
>
> On Mon, Jul 20, 2026 at 8:05 AM Debasish Hazra via ncl-talk <
> ncl-talk at mailman.ucar.edu> wrote:
>
>> Thanks Karin for your reply. The file is big to attach, here is the link
>> of actual file :
>> https://nomads.ncep.noaa.gov/pub/data/nccf/com/gens/prod/gefs.20260718/00/chem/pgrb2ap25/
>> I am getting error in this line : lat_new(i,0:npts_range-1) =
>> obs_lat(indexes)
>> error says : fatal:Subscript out of range, error in subscript #0
>> fatal:["Execute.c":8637]:Execute: Error occurred at or near line 68 in
>> file Test_read.ncl
>>
>> On Mon, Jul 20, 2026 at 3:17 AM Karin Meier-Fleischer <
>> meier-fleischer at dkrz.de> wrote:
>>
>>> Hi Debasish,
>>>
>>> the file gefs.chem.t00z.a2d_0p25.f018.grib2 is missing.
>>>
>>> Can you provide more Information about the error.
>>>
>>> Regards, Karin
>>>
>>> Am 18.07.26 um 04:26 schrieb Debasish Hazra via ncl-talk:
>>> > Hi,
>>> > I am trying to do a plot with station data in netcdf overlay on model
>>> > output following this example here
>>> > :https://www.ncl.ucar.edu/Applications/Scripts/polyg_8.ncl
>>> > But getting subscript range error. Attached is the ncl code and netcdf
>>> > station data used for this. Any help on this is greatly appreciated.
>>> >
>>> > Thanks
>>> > Debasish
>>
>>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20260720/e37f352c/attachment.htm>
-------------- next part --------------
Copyright (C) 1995-2019 - All Rights Reserved
University Corporation for Atmospheric Research
NCAR Command Language Version 6.6.2
The use of this software is governed by a License Agreement.
See http://www.ncl.ucar.edu/ for more details.
Variable: gsd1
Type: float
Total Size: 4152960 bytes
1038240 values
Number of Dimensions: 2
Dimensions and sizes: [lat_0 | 721] x [lon_0 | 1440]
Coordinates:
lat_0: [90..-90]
lon_0: [ 0..359.75]
Number Of Attributes: 16
center : US National Weather Service - NCEP (WMC)
production_status : Operational products
long_name : Aerosol optical thickness
units : Numeric
_FillValue : 1e+20
grid_type : Latitude/longitude
parameter_discipline_and_category : Meteorological products, Atmospheric chemical constituents
parameter_template_discipline_category_number : ( 48, 0, 20, 102 )
level_type : Entire Atmosphere
level : 0
forecast_time : 18
forecast_time_units : hours
initial_time : 07/16/2026 (00:00)
aerosol_wavelength : >= 5.45e-07 and <= 5.65e-07
aerosol_size : < 2e-05
aerosol_type : Total Aerosol
Variable: gsd1
Type: float
Total Size: 4152960 bytes
1038240 values
Number of Dimensions: 2
Dimensions and sizes: [lat_0 | 721] x [lon_0 | 1440]
Coordinates:
lat_0: [-90..90]
lon_0: [ 0..359.75]
Number Of Attributes: 16
center : US National Weather Service - NCEP (WMC)
production_status : Operational products
long_name : Aerosol optical thickness
units : Numeric
_FillValue : 1e+20
grid_type : Latitude/longitude
parameter_discipline_and_category : Meteorological products, Atmospheric chemical constituents
parameter_template_discipline_category_number : ( 48, 0, 20, 102 )
level_type : Entire Atmosphere
level : 0
forecast_time : 18
forecast_time_units : hours
initial_time : 07/16/2026 (00:00)
aerosol_wavelength : >= 5.45e-07 and <= 5.65e-07
aerosol_size : < 2e-05
aerosol_type : Total Aerosol
Variable: obs_lat
Type: float
Total Size: 8248 bytes
2062 values
Number of Dimensions: 1
Dimensions and sizes: [nhdr | 2062]
Coordinates:
Number Of Attributes: 3
long_name : latitude
_FillValue : -9999
units : degrees_north
Variable: obs_lon
Type: float
Total Size: 8248 bytes
2062 values
Number of Dimensions: 1
Dimensions and sizes: [nhdr | 2062]
Coordinates:
Number Of Attributes: 3
long_name : longitude
_FillValue : -9999
units : degrees_east
Variable: pm25_data
Type: float
Total Size: 15928 bytes
3982 values
Number of Dimensions: 1
Dimensions and sizes: [nobs | 3982]
Coordinates:
Number Of Attributes: 2
long_name : observation value
_FillValue : -9999
Variable: num_distinct_markers
Type: integer
Total Size: 4 bytes
1 values
Number of Dimensions: 1
Dimensions and sizes: [1]
Coordinates:
(0) min=0 max=3981
(0) num missing = 0
Variable: i
Type: integer
Total Size: 4 bytes
1 values
Number of Dimensions: 1
Dimensions and sizes: [1]
Coordinates:
(0) 0
Variable: npts_range
Type: integer
Total Size: 4 bytes
1 values
Number of Dimensions: 1
Dimensions and sizes: [1]
Coordinates:
(0) 3465
Variable: lat_new
Type: float
Total Size: 191136 bytes
47784 values
Number of Dimensions: 2
Dimensions and sizes: [12] x [3982]
Coordinates:
Number Of Attributes: 1
_FillValue : -9999
Variable: obs_lat
Type: float
Total Size: 8248 bytes
2062 values
Number of Dimensions: 1
Dimensions and sizes: [nhdr | 2062]
Coordinates:
Number Of Attributes: 3
long_name : latitude
_FillValue : -9999
units : degrees_north
Variable: indexes
Type: integer
Total Size: 13860 bytes
3465 values
Number of Dimensions: 1
Dimensions and sizes: [3465]
Coordinates:
fatal:Subscript out of range, error in subscript #0
fatal:["Execute.c":8637]:Execute: Error occurred at or near line 73 in file Test_read.ncl
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Test_read.ncl
Type: application/octet-stream
Size: 4514 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20260720/e37f352c/attachment.obj>
More information about the ncl-talk
mailing list