Hi ALL,
I kindly request for your help if you have some knowledge about the NCL. I am plotting zenith hydrostatic delay, zenith wet delay, total zenith path delay and integrated water vapor over the entire domain from the attached script.

; 2D sigma obtained by taking on the data on a single level (same in all levels), in this case, level index 0.
pg_sum = (/pg_levels_sum(0,0:193,0:241)/) ;copy only the data
pgt_sum = (/pgt_levels_sum(0,0:193,0:241)/) ;copy only the data
pgr_sum = (/pgr_levels_sum(0,0:193,0:241)/) ;copy only the data

; plot terrain height over the entire domain using the geopotential height variable at full-levels obtained
; from the last netcdf file opened through the loop above
; full-level geopotential corresponding to level index 0 is equivalent to the terrain height

But when i run it, it fails to plot. Here is the message received from the terminal;

Copyright (C) 1995-2015 - All Rights Reserved
University Corporation for Atmospheric Research
NCAR Command Language Version 6.3.0
The use of this software is governed by a License Agreement.
See http://www.ncl.ucar.edu/ for more details.
fatal:Subscript out of range, error in subscript #1
fatal:An error occurred reading xlat
fatal:["Execute.c":8575]:Execute: Error occurred at or near line 212 in file 7_tropodelay_iwv_and_terrainht_plots.ncl

I kindly request for help since i am stuck a this stage. I also want to have these plots seasonally like a plot DJF, MAM, JJA, SON
I will be grateful to have a solution.
Thank you very much.