1 | !+ |
---|
2 | ! |
---|
3 | ! .. _average_pb.jnl: |
---|
4 | ! |
---|
5 | ! ============== |
---|
6 | ! average_pb.jnl |
---|
7 | ! ============== |
---|
8 | ! |
---|
9 | ! DESCRIPTION |
---|
10 | ! =========== |
---|
11 | ! |
---|
12 | ! EXAMPLES |
---|
13 | ! ======== |
---|
14 | ! |
---|
15 | ! :: |
---|
16 | ! |
---|
17 | ! ferret> GO average_pb |
---|
18 | ! |
---|
19 | ! or:: |
---|
20 | ! |
---|
21 | ! $ ferret -nojnl -script average_pb.jnl |
---|
22 | ! |
---|
23 | ! SEE ALSO |
---|
24 | ! ======== |
---|
25 | ! |
---|
26 | ! TODO |
---|
27 | ! ==== |
---|
28 | ! |
---|
29 | ! find a solution (with /bad or missing function) |
---|
30 | ! |
---|
31 | ! log |
---|
32 | ! |
---|
33 | ! save figure |
---|
34 | ! |
---|
35 | ! save data to be compare to equivalent nco or cdo or .. implementation |
---|
36 | ! |
---|
37 | ! do the same with levitus, coads (ie a well known dataset) |
---|
38 | ! see http://www.oc.nps.edu/~bird/oc3030/ferret/ferret_notes.html |
---|
39 | ! ocean_atlas_subset.nc seems to be the closest dataset provided by ferret |
---|
40 | ! to tropflux one because of xt,x,y,z axis of temp variable |
---|
41 | ! |
---|
42 | ! see in http://lpas.epfl.ch/MOD/manuals/ferret/chapter3.html : default missing value is 1.E+34 |
---|
43 | ! |
---|
44 | ! get PROJECT_ID |
---|
45 | ! |
---|
46 | ! EVOLUTIONS |
---|
47 | ! ========== |
---|
48 | ! |
---|
49 | ! $Id$ |
---|
50 | ! |
---|
51 | ! $URL$ |
---|
52 | ! |
---|
53 | ! - fplod 20111107T113100Z aedon.locean-ipsl.upmc.fr (Darwin) |
---|
54 | ! |
---|
55 | ! * no more crash |
---|
56 | ! * confirm pb |
---|
57 | ! |
---|
58 | ! - fplod 20111104T130249Z aedon.locean-ipsl.upmc.fr (Darwin) |
---|
59 | ! |
---|
60 | ! * creation to try to answer to [TropFlux:7] mail:: |
---|
61 | ! |
---|
62 | ! I have downloaded TropFlux SST daily data from the site |
---|
63 | ! http://www.locean-ipsl.upmc.fr/~tropflux/data/daily. While using the |
---|
64 | ! data in ferret i encountered some problems. In ferret |
---|
65 | ! list/t=01-may-2001:30-jun-2001 sst[x=66:78@ave,y=6:18@ave] gives no |
---|
66 | ! result whereas list/t=01-may-2001:30-jun-2001 |
---|
67 | ! sst[x=66:70@ave,y=6:18@ave] gives the desired result... |
---|
68 | ! (the only change being the longitude changing from 78 to 70 and hence |
---|
69 | ! less influence of land coming into play). Can you please suggest any |
---|
70 | ! solution/workaround? |
---|
71 | ! |
---|
72 | ! inspired by |
---|
73 | ! /opt/ferret-6.6.5/examples/tutorial.jnl |
---|
74 | ! for plotting |
---|
75 | ! |
---|
76 | !- |
---|
77 | CANCEL DATA/ALL |
---|
78 | CANCEL SYMBOL/ALL |
---|
79 | ! |
---|
80 | ! pour creer des metafile convertible en ps |
---|
81 | set mode metafile |
---|
82 | ! pour avoir des traits plus gros |
---|
83 | GO bold |
---|
84 | ! |
---|
85 | ! pour les log |
---|
86 | ! http://www.pmel.noaa.gov/maillists/tmap/ferret_users/fu_2009/msg00036.html |
---|
87 | let date = "($session_date)" |
---|
88 | let datetime = "($session_date) ($session_time)" |
---|
89 | ! |
---|
90 | DEFINE SYMBOL tropflux_d = /usr/lodyc/incas/fplod/tropflux_d/daily |
---|
91 | DEFINE SYMBOL filename = sst_tropflux_1d_1989_2010.nc |
---|
92 | DEFINE SYMBOL fullfilename = ($tropflux_d)/sst_tropflux_1d_1989_2010.nc |
---|
93 | show symbol fullfilename |
---|
94 | use "($fullfilename)" |
---|
95 | ! |
---|
96 | ! have a look to data structure |
---|
97 | SHOW DATA/VARIABLES |
---|
98 | SHOW GRID sst |
---|
99 | ! have a look to memory |
---|
100 | SHOW MEMORY |
---|
101 | SET MODE DESPERATE ! ++ en attendant de trouver le pb |
---|
102 | ! |
---|
103 | ! original statement in [tropflux:7] |
---|
104 | ! list/t=01-may-2001:30-jun-2001 sst[x=66:78@ave,y=6:18@ave] |
---|
105 | !++ pb format E pour voir les vrais valeurs au lieu de ... à la place de vm |
---|
106 | !++list/format=(E)/t=01-may-2001:30-jun-2001 sst[x=66:78@ave,y=6:18@ave] |
---|
107 | list/t=01-may-2001:30-jun-2001 sst[x=66:78@ave,y=6:18@ave] |
---|
108 | LET sst_ave1 = sst[x=66:78@ave,y=6:18@ave,t=01-may-2001:30-jun-2001] |
---|
109 | PLOT sst_ave1 |
---|
110 | LET/TITLE="SST anomaly x=66:78" sst_anom1 = sst - sst_ave1 |
---|
111 | !++PLOT sst_anom1 |
---|
112 | ! |
---|
113 | ! original statement in [tropflux:7] |
---|
114 | ! list/t=01-may-2001:30-jun-2001 sst[x=66:70@ave,y=6:18@ave] |
---|
115 | list/t=01-may-2001:30-jun-2001 sst[x=66:70@ave,y=6:18@ave] |
---|
116 | LET sst_ave2 = sst[x=66:70@ave,y=6:18@ave,t=01-may-2001:30-jun-2001] |
---|
117 | PLOT sst_ave2 |
---|
118 | ! ++LET/TITLE="SST anomaly x=66:70" sst_anom2 = sst - sst_ave2 |
---|
119 | ! ++PLOT sst_anom2 |
---|
120 | ! |
---|
121 | ! save data |
---|
122 | ! ++ LIST/file=average_pb.output/format=(20E11.3)/order=xy/L=7 sst_ave1 sst_ave2 |
---|
123 | ! |
---|
124 | ! save figure++ |
---|
125 | ! ++exit |
---|