1 | !+ |
---|
2 | ! |
---|
3 | ! .. _read.jnl: |
---|
4 | ! |
---|
5 | ! |
---|
6 | ! ======== |
---|
7 | ! read.jnl |
---|
8 | ! ======== |
---|
9 | ! |
---|
10 | ! |
---|
11 | ! DESCRIPTION |
---|
12 | ! =========== |
---|
13 | ! |
---|
14 | ! read TROPFLUX file ${PROJECT_D}/TropFlux_19890101_20091231.nc |
---|
15 | ! |
---|
16 | ! plot the first time step of swr |
---|
17 | ! |
---|
18 | ! EXAMPLES |
---|
19 | ! ======== |
---|
20 | ! |
---|
21 | ! :: |
---|
22 | ! |
---|
23 | ! ferret> GO read |
---|
24 | ! |
---|
25 | ! or:: |
---|
26 | ! |
---|
27 | ! $ ferret -nojnl -script read.jnl |
---|
28 | ! |
---|
29 | ! SEE ALSO |
---|
30 | ! ======== |
---|
31 | ! |
---|
32 | ! matlab/octave similar program : :ref:`read.m` |
---|
33 | ! |
---|
34 | ! IDL+SAXO similar program : :ref:`read.pro` |
---|
35 | ! |
---|
36 | ! TODO |
---|
37 | ! ==== |
---|
38 | ! |
---|
39 | ! usage of $TROFLUX_D |
---|
40 | ! |
---|
41 | ! make it usable for real end-user |
---|
42 | ! |
---|
43 | ! learn ferret |
---|
44 | ! |
---|
45 | ! clean metafiles |
---|
46 | ! |
---|
47 | ! pb with Fprint on zeus probably because of an mtt problem |
---|
48 | ! |
---|
49 | ! pb with ferret on adonis |
---|
50 | ! |
---|
51 | ! pb with underscore in dataset label on zeus |
---|
52 | ! |
---|
53 | ! pas de floatstr function sur zeus ?! |
---|
54 | ! cf http://ferret.pmel.noaa.gov/Ferret/documentation/users-guide/handing-string-data-symbols/STRING-FUNCTIONS?searchterm=floatstr |
---|
55 | ! |
---|
56 | ! pb missing value |
---|
57 | ! EVOLUTIONS |
---|
58 | ! ========== |
---|
59 | ! |
---|
60 | ! - fplod 20110106T164602Z aedon.locean-ipsl.upmc.fr (Darwin) |
---|
61 | ! |
---|
62 | ! * add plot |
---|
63 | ! |
---|
64 | ! - fplod 20110106T090443Z aedon.locean-ipsl.upmc.fr (Darwin) |
---|
65 | ! |
---|
66 | ! * creation thanks to ACP |
---|
67 | ! |
---|
68 | !- |
---|
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 | !++ http://www.pmel.noaa.gov/maillists/tmap/ferret_users/fu_2010/msg00690.html |
---|
79 | !let mm = 1 |
---|
80 | !say `mm,zw=2` |
---|
81 | ! ++ http://www.pmel.noaa.gov/maillists/tmap/ferret_users/fu_2009/msg00771.html |
---|
82 | !repeat/range=1:2/name=m (use "http://dodsserver/sst_`m`.nc"; show dat) |
---|
83 | !repeat/range=1:12/name=mm ( show sym mm; say `mm,zw=2`) |
---|
84 | !repeat/range=1:12/name=mm ( show sym mm; LET CMM B`mm,zw=2`A; LIST CMM) |
---|
85 | |
---|
86 | ! pour les log |
---|
87 | ! http://www.pmel.noaa.gov/maillists/tmap/ferret_users/fu_2009/msg00036.html |
---|
88 | show symbol session* |
---|
89 | let date = "($session_date)" |
---|
90 | let datetime = "($session_date) ($session_time)" |
---|
91 | list date |
---|
92 | list datetime |
---|
93 | |
---|
94 | ! DEFINE VARIABLE yyyymmdd_b = "19890101" |
---|
95 | ! DEFINE VARIABLE yyyymmdd_e = "20091231" |
---|
96 | !++ DEFINE VARIABLE cb = {spawn:printf '%8.8d' `b` "} |
---|
97 | |
---|
98 | DEFINE SYMBOL yyyymmdd_b = 19890101 |
---|
99 | DEFINE SYMBOL yyyymmdd_e = 20091231 |
---|
100 | |
---|
101 | ! ++ pb floatstr si yyyymmdd_b est un entier DEFINE VARIABLE fullfilename = "/usr/zeus/incas/fplod/tropflux_d"+"TropFlux_" + floatstr(yyyymmdd_b,"(I8.8)") + "_" + floatstr(yyyymmdd_e,"(I8.8)") + ".nc" |
---|
102 | !DEFINE VARIABLE filename = "TropFlux_" + yyyymmdd_b + "_" + yyyymmdd_e + ".nc" |
---|
103 | DEFINE SYMBOL filename = TropFlux_($yyyymmdd_b)_($yyyymmdd_e).nc |
---|
104 | ! ++ $PROJECT_D defini ? |
---|
105 | !DEFINE SYMBOL tropflux_d = /usr/zeus/incas/fplod/tropflux_d/ |
---|
106 | !DEFINE SYMBOL fullfilename = ($tropflux_d)($filename) |
---|
107 | !DEFINE VARIABLE tropflux_d = {spawn:"echo $PROJECT_ID"} |
---|
108 | DEFINE SYMBOL tropflux_d = {spawn:"echo $PROJECT_ID"} |
---|
109 | SHOW SYM tropflux_d |
---|
110 | ! DEFINE SYMBOL fullfilename = "$tropflux_d($filename)" |
---|
111 | DEFINE SYMBOL fullfilename = `tropflux_d`($filename) |
---|
112 | CANCEL SYMBOL filename |
---|
113 | ! open |
---|
114 | ! test if the file exists |
---|
115 | !++DEFINE VARIABLE a = {"dummy", SPAWN:"ls `fullfilename`"} |
---|
116 | !++DEFINE VARIABLE nfiles = `a,RETURN=IEND` |
---|
117 | !++CANCEL a |
---|
118 | !++LIST nfiles |
---|
119 | !++IF `nfiles EQ 2` THEN EXIT |
---|
120 | !++CANCEL nfiles |
---|
121 | |
---|
122 | ! ++ use fullfilename |
---|
123 | !use "/usr/zeus/incas/fplod/tropflux_d/TropFlux_19890101_20091231.nc" |
---|
124 | !use "`fullfilename`" |
---|
125 | use "($fullfilename)" |
---|
126 | ! |
---|
127 | ! plot the first timestep of swr |
---|
128 | shade/l=1 swr |
---|
129 | ! |
---|
130 | ! plot world map |
---|
131 | go land |
---|
132 | ! |
---|
133 | ! create a Postscript file |
---|
134 | sp Fprint -l cps -p portrait -o TropFlux_19890101_20091231.ferret.ps metafile.plt |
---|
135 | ! |
---|
136 | ! create a GIF file |
---|
137 | frame/file=TropFlux_19890101_20091231.ferret.gif |
---|
138 | ! |
---|
139 | !++exit |
---|