Scripts/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	guez	3	! 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