trunk/Movie/rep_fill.jnl

! This is a script for Ferret.

! Description: repeat fill

! This script creates a set of GIF files for a 3D (longitude,
! latitude, time) variable.
! It does not average on longitude.
! There is one GIF image for each time value.
! The script takes 5 input arguments:
! - variable name,
! - directory of data file,
! - 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 Directory: $2

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

say minimum = `$1[x=@min, y=@min, t=@min]`
say maximum = `$1[x=@max, y=@max, t=@max]`

repeat/l=1:90 \
(\
   fill/levels=($3,$4,$5) $1; go land;\
   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 fill
4
5	! This script creates a set of GIF files for a 3D (longitude,
6	! latitude, time) variable.
7	! It does not average on longitude.
8	! There is one GIF image for each time value.
9	! The script takes 5 input arguments:
10	! - variable name,
11	! - directory of data file,
12	! - minimun level of color palette,
13	! - maximum level of color palette,
14	! - interval between levels.
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
25	say Directory: $2
26
27	use histday.nc
28	use/title="Run $2, beginning at day `int(t_ave_00086400[l=1]/86400)`" histday.nc
29
30	say minimum = `$1[x=@min, y=@min, t=@min]`
31	say maximum = `$1[x=@max, y=@max, t=@max]`
32
33	repeat/l=1:90 \
34	(\
35	fill/levels=($3,$4,$5) $1; go land;\
36	frame/file=`l`.gif;\
37	)
38
39	! Restore Ferret state:
40	set region save
41	set data/restore
42	set grid/restore