Scripts/Movie/rep_vect.jnl

! This is a script for Ferret.

! Description: repeat vector

! This script creates a set of GIF files for a couple of 3D (longitude,
! latitude, time) variables. The  two variables are assembled in a
! vector plot.
! The script does not average on longitude.
! There is one GIF image for each time value.
! The script takes 4 input arguments:
! - variable name 1,
! - variable name 2,
! - directory of data file,
! - standard vector length.

!!set mode verify
!!set mode journal

! Preserve Ferret state:
define region/default save
set data/save
set grid/save

say Variable name 1: $1
say Variable name 2: $2
say Directory: $3

use histday.nc
use/title="Run $3, beginning at day `int(t_ave_00086400[l=1]/86400)`" histday.nc

let norm = ($1^2 + $2^2)^0.5
say maximum of norm = `norm[x=@max, y=@max, t=@max]`

repeat/l=1:90 \
(\
   vector/aspect/length=$4 $1, $2; go land red;\
   frame/file=`l`.gif;\
)

! Restore Ferret state:
set region save
set data/restore
set grid/restore
1	guez	3	! This is a script for Ferret.
2
3			! Description: repeat vector
4
5			! This script creates a set of GIF files for a couple of 3D (longitude,
6			! latitude, time) variables. The two variables are assembled in a
7			! vector plot.
8			! The script does not average on longitude.
9			! There is one GIF image for each time value.
10			! The script takes 4 input arguments:
11			! - variable name 1,
12			! - variable name 2,
13			! - directory of data file,
14			! - standard vector length.
15
16			!!set mode verify
17			!!set mode journal
18
19			! Preserve Ferret state:
20			define region/default save
21			set data/save
22			set grid/save
23
24			say Variable name 1: $1
25			say Variable name 2: $2
26			say Directory: $3
27
28			use histday.nc
29			use/title="Run $3, beginning at day `int(t_ave_00086400[l=1]/86400)`" histday.nc
30
31			let norm = ($1^2 + $2^2)^0.5
32			say maximum of norm = `norm[x=@max, y=@max, t=@max]`
33
34			repeat/l=1:90 \
35			(\
36			vector/aspect/length=$4 $1, $2; go land red;\
37			frame/file=`l`.gif;\
38			)
39
40			! Restore Ferret state:
41			set region save
42			set data/restore
43			set grid/restore