Matlab Scripts for constructing Lagrangian Trajectories in HOPS
Below are links to Matlab scripts for computing trajectories from the Harvard Ocean Model data
for the California Current including Monterey Bay.
To start, download the files which include:
- HOPS.mat Matlab save file containing variables in arrays u and v, and the matlab structure hops
holding grid point and discretization and coast line information.
- mb_example.m Matlab script which runs an example computing 8 trajectories in and around Monterey Bay.
- trajectories.m Matlab function wrapper for ODE45 or ODE15 calling archive_path.m which is the actual matlab ODE function
to compute the RHS of odes from gridded velocity data.
- archive_path.m See above.
To begin, download the 4 files into a directory and launch MatLab. To run the example script, first load the velocity and
grid data and then simply execute the script mb_example:
>> load HOPS
>> mb_example
The result should be a sequence of images showing the current trajectories (8) and the underlying, time-dependent vector
field.
Please give this a try and let me know if there are problems/questions.
SOME MOVIES TO LOOK AT:
- Exit_Time_1.mov
- FTLE_1.mov
- Exit_Time_Hourly.mov
- Exit_Time_Hourly_wTraj.mov
- Etimes_Control_Nominal.jpg Comparison of Exitimes along Nominal and Control Trajectory.
MOVIES at Day 6 start:
- Exit_Day6.mov
- FTLE_Day6.mov
MOVIES at Day 1 with Control Trajectory:
- Exit_Day1_wtraj_HiRes.mov
- Exit_Day1_ZOOM.mov
MOVIES at Day 5 (t = 121) with Control Trajectories:
- Exit_Day5_wtrajs.mov
- Traj3_Input_Etime.jpg
- Exit_Day5_New.mov With new Blane data for trajectory 3 ...
- CompP1P3_Input_Etime.jpg With new trajectory 3 ... note trajectory 1 leaves control region.
File translated from
TEX
by
TTH,
version 3.44.
On 8 Sep 2008, 13:44.