trunk/Movie/rep_ave_f.jnl

! This is a script for Ferret.

! Description: repeat average fill

! This script uses the data file "histday.nc". The data file should be
! in the current directory. The script creates a set of GIF files for
! a 4D variable. It averages on longitude. There is one GIF image
! for each time value.

! The script takes 5 input arguments:
! - variable name,
! - Identification of run,
! - minimun level of color palette,
! - maximum level of color palette,
! - interval between levels.

!!set mode verify
!!set mode journal

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

say Variable name: $1
say Run: $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, z=@min, t=@min]`
say maximum = `var_ave[y=@max, z=@max, t=@max]`

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