<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
So I’m trying to plot idealised Kelvin waves (only) on a WK spectrum, using wkSpaceTime2. This method has worked for me using the ECMWF OA product, but the model output from this analytical KW keeps kicking up this error:
<div class=""><br class="">
</div>
<div class="">
<div style="margin: 0px; font-size: 14px; line-height: normal; font-family: Monaco; background-color: rgb(254, 255, 252);" class="">
<span style="font-variant-ligatures: no-common-ligatures" class="">fatal:["NclFile.c":2100]:Subscript out of range, error in subscript #1</span></div>
<div style="margin: 0px; font-size: 14px; line-height: normal; font-family: Monaco; background-color: rgb(254, 255, 252);" class="">
<span style="font-variant-ligatures: no-common-ligatures" class="">fatal:["Execute.c":8637]:Execute: Error occurred at or near line 25 in file WKST2.ncl</span></div>
<div style="margin: 0px; font-size: 14px; line-height: normal; font-family: Monaco; background-color: rgb(254, 255, 252); min-height: 19px;" class="">
<span style="font-variant-ligatures: no-common-ligatures" class=""></span><br class="">
</div>
<div style="margin: 0px; font-size: 14px; line-height: normal; font-family: Monaco; background-color: rgb(254, 255, 252);" class="">
<span style="font-variant-ligatures: no-common-ligatures" class="">fatal:Variable (U) is undefined</span></div>
<div style="margin: 0px; font-size: 14px; line-height: normal; font-family: Monaco; background-color: rgb(254, 255, 252);" class="">
<span style="font-variant-ligatures: no-common-ligatures" class="">fatal:["Execute.c":8637]:Execute: Error occurred at or near line 35 in file WKST2.ncl</span></div>
</div>
<div class=""><span style="font-variant-ligatures: no-common-ligatures" class=""><br class="">
</span></div>
<div class="">and here is the code:</div>
<div class=""><br class="">
</div>
<div class="">
<div class=""><font face="Monaco" class="">load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"</font></div>
<div class=""><font face="Monaco" class="">load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"</font></div>
<div class=""><font face="Monaco" class="">load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"</font></div>
<div class=""><font face="Monaco" class="">load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/diagnostics_cam.ncl"</font></div>
<div class=""><font face="Monaco" class=""><br class="">
</font></div>
<div class=""><font face="Monaco" class="">;==================================================================</font></div>
<div class=""><font face="Monaco" class="">; USER INTERFACE </font></div>
<div class=""><font face="Monaco" class="">; MAIN DRIVER SCRIPT</font></div>
<div class=""><font face="Monaco" class="">;==================================================================</font></div>
<div class=""><font face="Monaco" class=""> case = "KW_k1-15_3Heqs"</font></div>
<div class=""><font face="Monaco" class=""> diro = "./" ; output dir: location of plots </font></div>
<div class=""><font face="Monaco" class=""><br class="">
</font></div>
<div class=""><font face="Monaco" class=""><br class="">
</font></div>
<div class=""><font face="Monaco" class=""> VAR = "U200"</font></div>
<div class=""><font face="Monaco" class=""> diri = "./"</font></div>
<div class=""><font face="Monaco" class=""> fili = "KW_k1-15_3Heqs.nc"</font></div>
<div class=""><font face="Monaco" class=""><br class="">
</font></div>
<div class=""><font face="Monaco" class=""> nDayWin = 96 ; Wheeler-Kiladis [WK] temporal window length (days)</font></div>
<div class=""><font face="Monaco" class=""> nDaySkip = -65 ; time (days) between temporal windows [segments]</font></div>
<div class=""><font face="Monaco" class=""> ; negative means there will be overlapping temporal segments</font></div>
<div class=""><font face="Monaco" class=""> latBound = 15 ; 15N lat </font></div>
<div class=""><font face="Monaco" class=""><br class="">
</font></div>
<div class=""><font face="Monaco" class=""> ; get data </font></div>
<div class=""><font face="Monaco" class=""> f = addfile (diri+fili, "r")</font></div>
<div class=""><font face="Monaco" class=""> U = f->u(0:1460,1,{latBound:-latBound},:) ; [time | 1461] x [lat | 15] x [lon | 256]</font></div>
<div class=""><font face="Monaco" class=""> spd = 4 ;</font></div>
<div class=""><font face="Monaco" class=""> </font></div>
<div class=""><font face="Monaco" class=""><br class="">
</font></div>
<div class=""><font face="Monaco" class=""> opt = True</font></div>
<div class=""><font face="Monaco" class=""> opt@debug = True </font></div>
<div class=""><font face="Monaco" class=""> opt@pltType = "eps" ; create encapsulated postscript</font></div>
<div class=""><font face="Monaco" class=""><br class="">
</font></div>
<div class=""><font face="Monaco" class=""> opt@pltColorMap = "GMT_haxby" ; change the default color map</font></div>
<div class=""><font face="Monaco" class=""><br class="">
</font></div>
<div class=""><font face="Monaco" class=""> wkSpaceTime (U, diro, case, VAR \</font></div>
<div class=""><font face="Monaco" class=""> ,latBound, spd , nDayWin, nDaySkip, opt )</font></div>
</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
<div class="">and here is the metadata from the model output:</div>
<div class=""><br class="">
</div>
<div class="">
<div style="margin: 0px; font-size: 14px; line-height: normal; font-family: Monaco; background-color: rgb(254, 255, 252);" class="">
<span style="font-variant-ligatures: no-common-ligatures" class="">netcdf KW_k1-15_3Heqs {</span></div>
<div style="margin: 0px; font-size: 14px; line-height: normal; font-family: Monaco; background-color: rgb(254, 255, 252);" class="">
<span style="font-variant-ligatures: no-common-ligatures" class="">dimensions:</span></div>
<div style="margin: 0px; font-size: 14px; line-height: normal; font-family: Monaco; background-color: rgb(254, 255, 252);" class="">
<span style="font-variant-ligatures: no-common-ligatures" class=""><span class="Apple-tab-span" style="white-space:pre"></span>lon = 256 ;</span></div>
<div style="margin: 0px; font-size: 14px; line-height: normal; font-family: Monaco; background-color: rgb(254, 255, 252);" class="">
<span style="font-variant-ligatures: no-common-ligatures" class=""><span class="Apple-tab-span" style="white-space:pre"></span>lat = 128 ;</span></div>
<div style="margin: 0px; font-size: 14px; line-height: normal; font-family: Monaco; background-color: rgb(254, 255, 252);" class="">
<span style="font-variant-ligatures: no-common-ligatures" class=""><span class="Apple-tab-span" style="white-space:pre"></span>lev = 1 ;</span></div>
<div style="margin: 0px; font-size: 14px; line-height: normal; font-family: Monaco; background-color: rgb(254, 255, 252);" class="">
<span style="font-variant-ligatures: no-common-ligatures" class=""><span class="Apple-tab-span" style="white-space:pre"></span>time = 1461 ;</span></div>
<div style="margin: 0px; font-size: 14px; line-height: normal; font-family: Monaco; background-color: rgb(254, 255, 252);" class="">
<span style="font-variant-ligatures: no-common-ligatures" class="">variables:</span></div>
<div style="margin: 0px; font-size: 14px; line-height: normal; font-family: Monaco; background-color: rgb(254, 255, 252);" class="">
<span style="font-variant-ligatures: no-common-ligatures" class=""><span class="Apple-tab-span" style="white-space:pre"></span>float lat(lat) ;</span></div>
<div style="margin: 0px; font-size: 14px; line-height: normal; font-family: Monaco; background-color: rgb(254, 255, 252);" class="">
<span style="font-variant-ligatures: no-common-ligatures" class=""><span class="Apple-tab-span" style="white-space:pre"></span>lat:units = "degrees_north" ;</span></div>
<div style="margin: 0px; font-size: 14px; line-height: normal; font-family: Monaco; background-color: rgb(254, 255, 252);" class="">
<span style="font-variant-ligatures: no-common-ligatures" class=""><span class="Apple-tab-span" style="white-space:pre"></span>float lon(lon) ;</span></div>
<div style="margin: 0px; font-size: 14px; line-height: normal; font-family: Monaco; background-color: rgb(254, 255, 252);" class="">
<span style="font-variant-ligatures: no-common-ligatures" class=""><span class="Apple-tab-span" style="white-space:pre"></span>lon:units = "degrees_east" ;</span></div>
<div style="margin: 0px; font-size: 14px; line-height: normal; font-family: Monaco; background-color: rgb(254, 255, 252);" class="">
<span style="font-variant-ligatures: no-common-ligatures" class=""><span class="Apple-tab-span" style="white-space:pre"></span>float lev(lev) ;</span></div>
<div style="margin: 0px; font-size: 14px; line-height: normal; font-family: Monaco; background-color: rgb(254, 255, 252);" class="">
<span style="font-variant-ligatures: no-common-ligatures" class=""><span class="Apple-tab-span" style="white-space:pre"></span>lev:units = "model level" ;</span></div>
<div style="margin: 0px; font-size: 14px; line-height: normal; font-family: Monaco; background-color: rgb(254, 255, 252);" class="">
<span style="font-variant-ligatures: no-common-ligatures" class=""><span class="Apple-tab-span" style="white-space:pre"></span>float time(time) ;</span></div>
<div style="margin: 0px; font-size: 14px; line-height: normal; font-family: Monaco; background-color: rgb(254, 255, 252);" class="">
<span style="font-variant-ligatures: no-common-ligatures" class=""><span class="Apple-tab-span" style="white-space:pre"></span>time:units = "hours since 2007-01-01 00:00:00" ;</span></div>
<div style="margin: 0px; font-size: 14px; line-height: normal; font-family: Monaco; background-color: rgb(254, 255, 252);" class="">
<span style="font-variant-ligatures: no-common-ligatures" class=""><span class="Apple-tab-span" style="white-space:pre"></span>float gp(time, lev, lat, lon) ;</span></div>
<div style="margin: 0px; font-size: 14px; line-height: normal; font-family: Monaco; background-color: rgb(254, 255, 252);" class="">
<span style="font-variant-ligatures: no-common-ligatures" class=""><span class="Apple-tab-span" style="white-space:pre"></span>gp:units = "m**2 s**-2" ;</span></div>
<div style="margin: 0px; font-size: 14px; line-height: normal; font-family: Monaco; background-color: rgb(254, 255, 252);" class="">
<span style="font-variant-ligatures: no-common-ligatures" class=""><span class="Apple-tab-span" style="white-space:pre"></span>float u(time, lev, lat, lon) ;</span></div>
<div style="margin: 0px; font-size: 14px; line-height: normal; font-family: Monaco; background-color: rgb(254, 255, 252);" class="">
<span style="font-variant-ligatures: no-common-ligatures" class=""><span class="Apple-tab-span" style="white-space:pre"></span>u:units = "m s**-1" ;</span></div>
<div style="margin: 0px; font-size: 14px; line-height: normal; font-family: Monaco; background-color: rgb(254, 255, 252);" class="">
<span style="font-variant-ligatures: no-common-ligatures" class="">}</span></div>
</div>
<div class=""><span style="font-variant-ligatures: no-common-ligatures" class=""><br class="">
</span></div>
<div class=""><span style="font-variant-ligatures: no-common-ligatures" class=""><br class="">
</span></div>
<div class=""><span style="font-variant-ligatures: no-common-ligatures" class="">I have tried so many different attacks to try and get this to work, with varying degrees of success, but I’ve not been able to get any plots using this wkSpaceTime2 skeleton code.
My suspicion is that it may be to do with the units of time in the .nc file, but I don’t think I understand the script well enough to know why, but I imagine there is some NCO or CDO command to change the units of time? Anyway, if this looks like something
familiar to anyone I would be most grateful if you would be able to help me solve this problem.</span></div>
<div class=""><span style="font-variant-ligatures: no-common-ligatures" class=""><br class="">
</span></div>
<div class=""><span style="font-variant-ligatures: no-common-ligatures" class="">Cheers,</span></div>
<div class=""><span style="font-variant-ligatures: no-common-ligatures" class=""><br class="">
</span></div>
<div class=""><span style="font-variant-ligatures: no-common-ligatures" class="">Elliot McKinnon-Gray </span></div>
<div class=""><span style="font-variant-ligatures: no-common-ligatures" class="">—</span></div>
<div class="">Doktorand, International Max Planck Research School for Earth System Modelling IMPRS-ESM</div>
<div class="">Meteorological Institute, Universität Hamburg, Room 4012 GB7</div>
</body>
</html>