[gradsusr] Trajectory Analysis

[James T. Potemra]

Owanda:

The traj.gs script is written to accept a time-step of either hours or 
minutes (and you are using months).  I think you can change your TDEF 
line to be 720HR, but this of course assumes an even 30 days per month.

Jim

Owanda Otieno wrote:
> Hi all;
> Still have a question on Trajectory plots,
> Following Jim's recomendation, I developed a control (ctl) file for my 
> RegCM3 output. It opens the data and I can plot any variable I want to 
> but when I run _traj.gs <http://traj.gs>_ I get the following error
>  
> HR or MN missing in line TDEF of file ATMmonera.8788.ctl
>  
> my TDEF line is as shown below, which I dont see any problem with
>
> tdef 13 linear 00z15dec1987 1mo
>
> can anyone help figure out this problem
>
> I have also attached the ctl file.
>
> Thanks again.
>
> On Thu, Apr 15, 2010 at 2:46 PM, James T. Potemra <jimp at hawaii.edu 
> <mailto:jimp at hawaii.edu>> wrote:
>
>     Owanda:
>
>     The way traj.gs <http://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 <http://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 <http://traj.gs/>.
>
>     Jim
>     ----------------
>     Here are the relevant lines from traj.gs <http://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/> <http://traj.gs
>         <http://traj.gs/>>.
>
>         TDEF not found in ATMmonera.8788.nc
>         <http://atmmonera.8788.nc/> <http://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
>
>
>
>
> -- 
> 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
>