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 |