[gradsusr] Trajectory Analysis

[Bernat Codina]

I'm afraid it makes no sense computing trajectories with snapshots
taken at 1-month intervals.

Bernat Codina
Dept. d'Astronomia i Meteorologia
Universitat de Barcelona




2010/4/16 James T. Potemra <jimp at hawaii.edu>:
> 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
>>