trunk/Movie/rep_ave_pl.jnl

! This is a script for Ferret.

! Description: repeat average plot

! This script creates a set of GIF files for a 3D (longitude,
! latitude, time) variable, averaged on longitude.
! There is one GIF image for each time value.
! The script takes 4 input arguments:
! - variable name,
! - directory of data file,
! - lower limit of vertical axis,
! - upper limit of vertical axis.

!!set mode verify
!!set mode journal

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

say Variable name: $1
say Directory: $2

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

let var_ave = $1[x=@ave]

say minimum = `var_ave[y=@min, t=@min]`
say maximum = `var_ave[y=@max, t=@max]`

repeat/l=1:90 \
(\
   plot/vlimits=$3:$4 var_ave;\
   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 average plot
4
5	! This script creates a set of GIF files for a 3D (longitude,
6	! latitude, time) variable, averaged on longitude.
7	! There is one GIF image for each time value.
8	! The script takes 4 input arguments:
9	! - variable name,
10	! - directory of data file,
11	! - lower limit of vertical axis,
12	! - upper limit of vertical axis.
13
14	!!set mode verify
15	!!set mode journal
16
17	! Preserve Ferret state:
18	define region/default save
19	set data/save
20	set grid/save
21
22	say Variable name: $1
23	say Directory: $2
24
25	use histday.nc
26	use/title="Run $2, beginning at day `int(t_ave_00086400[l=1]/86400)`" histday.nc
27
28	let var_ave = $1[x=@ave]
29
30	say minimum = `var_ave[y=@min, t=@min]`
31	say maximum = `var_ave[y=@max, t=@max]`
32
33	repeat/l=1:90 \
34	(\
35	plot/vlimits=$3:$4 var_ave;\
36	frame/file=`l`.gif;\
37	)
38
39	! Restore Ferret state:
40	set region save
41	set data/restore
42	set grid/restore