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