1 | !+ |
---|
2 | ! |
---|
3 | ! .. _average_demo.jnl: |
---|
4 | ! |
---|
5 | ! ================ |
---|
6 | ! average_demo.jnl |
---|
7 | ! ================ |
---|
8 | ! |
---|
9 | ! DESCRIPTION |
---|
10 | ! =========== |
---|
11 | ! |
---|
12 | ! plot sst anomaly from ferret demo dataset |
---|
13 | ! |
---|
14 | ! no other interest than learning ferret and be able to work on tropflux |
---|
15 | ! dataset (see :ref:`average_pb.jnl`) |
---|
16 | ! |
---|
17 | ! EXAMPLES |
---|
18 | ! ======== |
---|
19 | ! |
---|
20 | ! :: |
---|
21 | ! |
---|
22 | ! ferret> GO average_demo |
---|
23 | ! |
---|
24 | ! or:: |
---|
25 | ! |
---|
26 | ! $ ferret -nojnl -script average_demo.jnl |
---|
27 | ! |
---|
28 | ! SEE ALSO |
---|
29 | ! ======== |
---|
30 | ! |
---|
31 | ! TODO |
---|
32 | ! ==== |
---|
33 | ! |
---|
34 | ! make it work |
---|
35 | ! |
---|
36 | ! save figure |
---|
37 | ! |
---|
38 | ! save data to be compare to equivalent nco or cdo or .. implementation |
---|
39 | ! |
---|
40 | ! do the same with levitus, coads (ie a well known dataset) |
---|
41 | ! see http://www.oc.nps.edu/~bird/oc3030/ferret/ferret_notes.html |
---|
42 | ! ocean_atlas_subset.nc seems to be the closest dataset provided by ferret |
---|
43 | ! to tropflux one because of xt,x,y,z axis of sst variable |
---|
44 | ! |
---|
45 | ! EVOLUTIONS |
---|
46 | ! ========== |
---|
47 | ! |
---|
48 | ! $Id$ |
---|
49 | ! |
---|
50 | ! $URL$ |
---|
51 | ! |
---|
52 | ! - fplod 20111104T150732Z aedon.locean-ipsl.upmc.fr (Darwin) |
---|
53 | ! |
---|
54 | ! * creation to try to answer to [TropFlux:7] see :ref:`average_pb.jnl` |
---|
55 | ! and learning ferret |
---|
56 | ! |
---|
57 | !- |
---|
58 | CANCEL DATA/ALL |
---|
59 | CANCEL SYMBOL/ALL |
---|
60 | ! |
---|
61 | ! pour creer des metafile convertible en ps |
---|
62 | set mode metafile |
---|
63 | ! pour avoir des traits plus gros |
---|
64 | go bold |
---|
65 | ! |
---|
66 | ! pour les log |
---|
67 | ! http://www.pmel.noaa.gov/maillists/tmap/ferret_users/fu_2009/msg00036.html |
---|
68 | let date = "($session_date)" |
---|
69 | let datetime = "($session_date) ($session_time)" |
---|
70 | ! |
---|
71 | SET DATA ocean_atlas_subset |
---|
72 | ! |
---|
73 | ! have a look to data structure |
---|
74 | SHOW DATA/VARIABLES |
---|
75 | SHOW GRID TEMP |
---|
76 | ! have a look to memory |
---|
77 | SHOW MEMORY |
---|
78 | ! |
---|
79 | ! original statement in [tropflux:7] |
---|
80 | ! list/t=01-may-2001:30-jun-2001 TEMP[x=66:78@ave,y=6:18@ave] |
---|
81 | SET REGION/z=1/x=66:78/y=6:18/t=16-jan-0001:17-apr-0001 |
---|
82 | list TEMP[z=1,x=@ave,y=@ave] |
---|
83 | LET sst_ave1 = TEMP[z=1,x=@ave,y=@ave] |
---|
84 | SHOW GRID sst_ave1 |
---|
85 | LIST sst_ave1 |
---|
86 | LET sst_ave1 = TEMP[z=1,x=@ave,y=@ave,t=@ave] |
---|
87 | SHOW GRID sst_ave1 |
---|
88 | LIST sst_ave1 |
---|
89 | PLOT sst_ave1 |
---|
90 | LET/TITLE="TEMP anomaly x=66:78" sst_anom1 = TEMP - sst_ave1 |
---|
91 | SHOW DATA/VARIABLES |
---|
92 | !++ PLOT sst_anom1 |
---|
93 | ! |
---|
94 | ! original statement in [tropflux:7] |
---|
95 | ! list/t=01-may-2001:30-jun-2001 sst[x=66:70@ave,y=6:18@ave] |
---|
96 | list/t=01-may-2001:30-jun-2001 sst[x=66:70@ave,y=6:18@ave] |
---|
97 | LET sst_ave2 = sst[x=66:70@ave,y=6:18@ave,t=01-may-2001:30-jun-2001] |
---|
98 | LET/TITLE="TEMP anomaly x=66:70" sst_anom2 = sst - sst_ave2 |
---|
99 | PLOT sst_anom2 |
---|
100 | ! |
---|
101 | ! save data |
---|
102 | ! ++ LIST/file=average_demo.output/format=(20E11.3)/order=xy/L=7 sst_ave1 sst_ave2 |
---|
103 | ! |
---|
104 | ! save figure++ |
---|
105 | !++exit |
---|