1 |
! This is a script for Ferret. |
2 |
|
3 |
! Description: repeat average plot |
4 |
|
5 |
! This script creates a set of GIF files for a 3D (longitude, |
6 |
! latitude, time) variable, averaged on longitude. |
7 |
! There is one GIF image for each time value. |
8 |
! The script takes 4 input arguments: |
9 |
! - variable name, |
10 |
! - directory of data file, |
11 |
! - lower limit of vertical axis, |
12 |
! - upper limit of vertical axis. |
13 |
|
14 |
!!set mode verify |
15 |
!!set mode journal |
16 |
|
17 |
! Preserve Ferret state: |
18 |
define region/default save |
19 |
set data/save |
20 |
set grid/save |
21 |
|
22 |
say Variable name: $1 |
23 |
say Directory: $2 |
24 |
|
25 |
use histday.nc |
26 |
use/title="Run $2, beginning at day `int(t_ave_00086400[l=1]/86400)`" histday.nc |
27 |
|
28 |
let var_ave = $1[x=@ave] |
29 |
|
30 |
say minimum = `var_ave[y=@min, t=@min]` |
31 |
say maximum = `var_ave[y=@max, t=@max]` |
32 |
|
33 |
repeat/l=1:90 \ |
34 |
(\ |
35 |
plot/vlimits=$3:$4 var_ave;\ |
36 |
frame/file=`l`.gif;\ |
37 |
) |
38 |
|
39 |
! Restore Ferret state: |
40 |
set region save |
41 |
set data/restore |
42 |
set grid/restore |