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