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