[gradsusr] Trajectory Analysis

[James T. Potemra]

Owanda:

The way traj.gs is written is that it scans the ctl file and looks for 
the line with "TDEF" (or "tdef"), and then gets time information from 
that line.  If you are reading the netcdf file without a ctl file, you 
option is to modify the traj.gs script to do something like "q dims" to 
get the time information.  Another option would be to make a ctl file 
and put in a TDEF line.  Of course, you could also "hardwire" the time 
definitions in traj.gs.

Jim
----------------
Here are the relevant lines from traj.gs; the first part reads info from 
the 'query' command, the next part ('fixter = subwrd(rec,2)') reads from 
the ctl file:

  'q file'
  rec=sublin(result,5)
  _Tmax=subwrd(rec,12)
  rec=sublin(result,2)
  fitxer=subwrd(rec,2)
  rc=0
  while(rc=0)
    rec=read(fitxer)
    rc=sublin(rec,1)
    lin=sublin(rec,2)
    if (subwrd(lin,1)='tdef' | subwrd(lin,1)='TDEF')
      interval=subwrd(lin,5)
-------------------

Owanda Otieno wrote:
>
> Hi all,
>
> I am trying to construct moisture track trajectories from RegCM3 
> simulation data, I get the following error when I run the script 
> traj.gs <http://traj.gs>.
>
> TDEF not found in ATMmonera.8788.nc <http://ATMmonera.8788.nc>
>
> Any idea please,
>
> Thank you
> -- 
> Owanda Otieno
> Department of Natural Resources and the Environment
> University of Connecticut
> 1376 Storrs Road, Unit 4087
> Storrs, CT 06269-4087
> Phone: (860) 486-1876
>