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 |