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	! 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