Question about wkSpaceTime_3.ncl script.

*Hi there,
   I would like to use the wkSpaceTime_3.ncl script in order to have some
plots of the space-time spectra.
So I downloaded data from NCEP for Omega variable for the period 01 Jan
2001 to 31 Dec 2004 and this is the result I am getting:*

*victor@hermitage:~/Desktop/temporal$ ncl kiladis.ncl
 Copyright (C) 1995-2010 - All Rights Reserved
 University Corporation for Atmospheric Research
 NCAR Command Language Version 5.2.0
 The use of this software is governed by a License Agreement.
 See http://www.ncl.ucar.edu/ for more details.

Variable: x (parameter)
Type: short
Total Size: 21879936 bytes
            10939968 values
Number of Dimensions: 3
Dimensions and sizes: [time | 5844] x [lat | 13] x [lon | 144]
Coordinates:
            time: [17531688.. 0]
            lat: [-15..15]
            lon: [ 0..357.5]
Number Of Attributes: 18
  level : 500
  long_name : 4xDaily Omega
  unpacked_valid_range : ( -3, 7 )
  actual_range : ( -1.759, 1.18 )
  units : Pascal/s
  add_offset : 29.765
  scale_factor : 0.001
  missing_value : 32766
  precision : 3
  least_significant_digit : 3
  GRIB_id : 39
  GRIB_name : VVEL
  var_desc : Omega (dp/dt)
  dataset : NMC Reanalysis
  level_desc : Multiple levels
  statistic : Individual Obs
  parent_stat : Other
  valid_range : ( -32765, -22765 )
(0)
(0) 4xDaily Omega: min=-30719 max=0
fatal:Assignment type mismatch, right hand side can't be coerced to type of
left hand side
fatal:Execute: Error occurred at or near line 626 in file
$NCARG_ROOT/lib/ncarg/nclscripts/csm/diagnostics_cam.ncl

fatal:Execute: Error occurred at or near line 35 in file kiladis.ncl*

*The script I used is the following:*

*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/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/diagnostics_cam.ncl"

  case = "NCEP: 2001-2005:"
  diro = "./" ; output dir: location of plots

  latN = 15
  latS = -latN ; symmetric about the equator

  nDayWin = 96 ; Wheeler-Kiladis [WK] temporal window length (days)
  nDaySkip = 30 ; time (days) between temporal windows [segments]
                     ; negative means there will be overlapping
                     ; temporal segments

  opt = True
  opt@debug= True
  opt@pltType = "eps"
  opt@cnLinesOn = False ; turn off contour lines
                     ; specify the plot intervals

  spd = 4 ; samples per day (6 hour sampling rate)
  var = "OMEGA500" ; user specified

  diri = "/home/victor/Desktop/temporal/"
  fili = "OMEGA.2001-2005.nc"

  f = addfile (diri+fili, "r")
  omega = f->omega(:,0,{latS:latN},:) ; (time,lat,lon)
  ;omega2 = (omega*0.00100000004749745)+29.7649993896484

                     ; calculation and create plots
  wkSpaceTime (omega, diro, case, var \
              ,latN, spd, nDayWin, nDaySkip, opt )*

*At first I thought this error was because I was not using the scale
factors, and that's why I wrote omega2 in the script, however that didn't
work either.
Do you have any idea about this issue.
Please let me know, any information will be valuable
Thanks in advance.

Victor Torres.
*

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Nov 16 14:51:11 2010