Thanks Yemi -- that works nicely. I add the illegal print statement into
my script following the variables I want to examine, copy the script to
an interactive ncl session to run, which then crashes, but leaves the
variables accessible in the workspace:
; code snippet ********************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/contrib/calendar_decode2.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
;************************************************
begin
;************************************************
; read in data
;************************************************
f=addfile
("../data/PFEL_winds_and_ocean_transports/erdlasFnWPr_4baa_b7bc_5c31.nc","r")
curl = f->curl(:,0,0)
time = f->time
yearfrac = cd_calendar(time, 4)
print(1/0) ; illegal division
;****************************
************* interactive session in terminal ****************
ncl 36> printVarSummary(yearfrac)
Variable: yearfrac
Type: double
Total Size: 4456 bytes
557 values
Number of Dimensions: 1
Dimensions and sizes: [557]
Coordinates:
Number Of Attributes: 1
calendar : standard
ncl 37> print(time)
Variable: time
Type: double
Total Size: 4456 bytes
557 values
Number of Dimensions: 1
Dimensions and sizes: [time | 557]
Coordinates:
time: [-93355200..1368792000]
Number Of Attributes: 8
_CoordinateAxisType : Time
actual_range : ( -93355200, 1368792000 )
axis : T
ioos_category : Time
long_name : reference time
standard_name : time
time_origin : 01-JAN-1970 00:00:00
units : seconds since 1970-01-01T00:00:00Z
(0) -93355200
(1) -90806400
(2) -88257600
(3) -85622400
(4) -82987200
(5) -80352000
(6) -77716800
ncl 38> print(curl)
Variable: curl
Type: float
Total Size: 2228 bytes
557 values
Number of Dimensions: 1
Dimensions and sizes: [time | 557]
Coordinates:
time: [-93355200..1368792000]
Number Of Attributes: 8
longitude : 242.5
latitude : 32.5
_FillValue : -9999
colorBarMaximum : 3
colorBarMinimum : -3
ioos_category : Wind
long_name : Curl of Surface Wind Stress
units : MPa m-1
(0) 0.06939734
(1) 0.1011225
(2) -0.09163887
(3) -0.2258962
(4) 0.07263987
(5) -0.02644638
ncl 39>
This will make debugging much more efficient.
Best fishes
Sam
-- Sam McClatchie, Supervisory oceanographer, Fisheries oceanography Southwest Fisheries Science Center, NOAA, 8901 La Jolla Shores Dr. La Jolla, CA 92037-1509, U.S.A. email <Sam.McClatchie@noaa.gov> Office: 858 546 7083, Cellular: 858 752 8495 Research home page <www.fishocean.info> /\ ...>><xX(¡> //// \\\\ <¡)Xx><< ///// \\\\\\ ><(((¡> >><(((¡> ...>><xX(¡>O<¡)Xx><< _______________________________________________ ncl-talk mailing list List instructions, subscriber options, unsubscribe: http://mailman.ucar.edu/mailman/listinfo/ncl-talkReceived on Sun Jun 16 10:20:18 2013
This archive was generated by hypermail 2.1.8 : Mon Jun 24 2013 - 11:46:47 MDT