/[MITgcm]/MITgcm/pkg/ecco/ecco_cost_summary.F
ViewVC logotype

Contents of /MITgcm/pkg/ecco/ecco_cost_summary.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph

Revision 1.3 - (show annotations) (download)
Tue Oct 9 00:02:51 2007 UTC (16 years, 7 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint60, checkpoint61, checkpoint62, checkpoint59q, checkpoint59p, checkpoint59r, checkpoint59m, checkpoint59l, checkpoint59o, checkpoint59n, checkpoint59i, checkpoint59k, checkpoint59j, checkpoint62c, checkpoint62b, checkpoint62a, checkpoint62g, checkpoint62f, checkpoint62e, checkpoint62d, checkpoint62k, checkpoint62j, checkpoint62i, checkpoint62h, checkpoint62o, checkpoint62n, checkpoint62m, checkpoint62l, checkpoint62s, checkpoint62r, checkpoint62q, checkpoint62p, checkpoint62t, checkpoint61f, checkpoint61g, checkpoint61d, checkpoint61e, checkpoint61b, checkpoint61c, checkpoint61a, checkpoint61n, checkpoint61o, checkpoint61l, checkpoint61m, checkpoint61j, checkpoint61k, checkpoint61h, checkpoint61i, checkpoint61v, checkpoint61w, checkpoint61t, checkpoint61u, checkpoint61r, checkpoint61s, checkpoint61p, checkpoint61q, checkpoint61z, checkpoint61x, checkpoint61y
Changes since 1.2: +4 -3 lines 
add missing cvs $Header:$ or $Name:$
1 C $Header: $
2 C $Name: $
3
4 #include "COST_CPPOPTIONS.h"
5
6
7 subroutine ecco_cost_summary( mythid )
8
9 c ==================================================================
10 c SUBROUTINE ecco_cost_summary
11 c ==================================================================
12 c
13 c o Summarize the cost function related parts of the ECCO release.
14 c
15 c started: Christian Eckert eckert@mit.edu 30-Jun-1999
16 c
17 c changed: Christian Eckert eckert@mit.edu 25-Feb-2000
18 c
19 c - Restructured the code in order to create a package
20 c for the MITgcmUV.
21 c
22 c ==================================================================
23 c SUBROUTINE ecco_cost_summary
24 c ==================================================================
25
26 implicit none
27
28 c == global variables ==
29
30 #include "EEPARAMS.h"
31 #include "SIZE.h"
32
33 #include "ecco_cost.h"
34
35 c == routine arguments ==
36
37 integer mythid
38
39 c == local variables ==
40
41 integer il
42
43 character*(max_len_mbuf) msgbuf
44
45 c == external ==
46
47 integer ilnblnk
48 external ilnblnk
49
50 c == end of interface ==
51
52 write(msgbuf,'(a)')
53 &' '
54 call print_message( msgbuf, standardmessageunit,
55 & SQUEEZE_RIGHT , mythid)
56 write(msgbuf,'(a)')
57 &'// ======================================================='
58 call print_message( msgbuf, standardmessageunit,
59 & SQUEEZE_RIGHT , mythid)
60 write(msgbuf,'(a)')
61 &'// ECCO cost function configuration >>> START <<<'
62 call print_message( msgbuf, standardmessageunit,
63 & SQUEEZE_RIGHT , mythid)
64 write(msgbuf,'(a)')
65 &'// ======================================================='
66 call print_message( msgbuf, standardmessageunit,
67 & SQUEEZE_RIGHT , mythid)
68 write(msgbuf,'(a)')
69 &' '
70 call print_message( msgbuf, standardmessageunit,
71 & SQUEEZE_RIGHT , mythid)
72
73 write(msgbuf,'(a)')
74 &' Multipliers for the indivdual cost function contributions:'
75 call print_message( msgbuf, standardmessageunit,
76 & SQUEEZE_RIGHT , mythid)
77 write(msgbuf,'(a)')
78 &' '
79 call print_message( msgbuf, standardmessageunit,
80 & SQUEEZE_RIGHT , mythid)
81 write(msgbuf,'(a,e10.3)')
82 &' Net heat flux: ',mult_hflux
83 call print_message( msgbuf, standardmessageunit,
84 & SQUEEZE_RIGHT , mythid)
85 write(msgbuf,'(a,e10.3)')
86 &' Salt flux: ',mult_sflux
87 call print_message( msgbuf, standardmessageunit,
88 & SQUEEZE_RIGHT , mythid)
89 write(msgbuf,'(a,e10.3)')
90 &' Zonal wind stress: ',mult_tauu
91 call print_message( msgbuf, standardmessageunit,
92 & SQUEEZE_RIGHT , mythid)
93 write(msgbuf,'(a,e10.3)')
94 &' Meridional wind stress: ',mult_tauv
95 call print_message( msgbuf, standardmessageunit,
96 & SQUEEZE_RIGHT , mythid)
97 write(msgbuf,'(a,e10.3)')
98 &' Mean sea surface height: ',mult_hmean
99 call print_message( msgbuf, standardmessageunit,
100 & SQUEEZE_RIGHT , mythid)
101 write(msgbuf,'(a,e10.3)')
102 &' Sea surface height anomalies: ',mult_h
103 call print_message( msgbuf, standardmessageunit,
104 & SQUEEZE_RIGHT , mythid)
105 write(msgbuf,'(a,e10.3)')
106 &' Temperature Lev.: ',mult_temp
107 call print_message( msgbuf, standardmessageunit,
108 & SQUEEZE_RIGHT , mythid)
109 write(msgbuf,'(a,e10.3)')
110 &' Salinity Lev.: ',mult_salt
111 call print_message( msgbuf, standardmessageunit,
112 & SQUEEZE_RIGHT , mythid)
113 write(msgbuf,'(a,e10.3)')
114 &' Temperature ini.: ',mult_temp0
115 call print_message( msgbuf, standardmessageunit,
116 & SQUEEZE_RIGHT , mythid)
117 write(msgbuf,'(a,e10.3)')
118 &' Salinity ini.: ',mult_salt0
119 call print_message( msgbuf, standardmessageunit,
120 & SQUEEZE_RIGHT , mythid)
121 write(msgbuf,'(a,e10.3)')
122 &' TMI Sea surface temperature: ',mult_tmi
123 call print_message( msgbuf, standardmessageunit,
124 & SQUEEZE_RIGHT , mythid)
125 write(msgbuf,'(a,e10.3)')
126 &' Sea surface temperature: ',mult_sst
127 call print_message( msgbuf, standardmessageunit,
128 & SQUEEZE_RIGHT , mythid)
129 write(msgbuf,'(a,e10.3)')
130 &' Sea surface salinity: ',mult_sss
131 call print_message( msgbuf, standardmessageunit,
132 & SQUEEZE_RIGHT , mythid)
133 write(msgbuf,'(a,e10.3)')
134 &' CTD temperature: ',mult_ctdt
135 call print_message( msgbuf, standardmessageunit,
136 & SQUEEZE_RIGHT , mythid)
137 write(msgbuf,'(a,e10.3)')
138 &' CTD salinity: ',mult_ctds
139 call print_message( msgbuf, standardmessageunit,
140 & SQUEEZE_RIGHT , mythid)
141 write(msgbuf,'(a,e10.3)')
142 &' CTD clim temperature: ',mult_ctdtclim
143 call print_message( msgbuf, standardmessageunit,
144 & SQUEEZE_RIGHT , mythid)
145 write(msgbuf,'(a,e10.3)')
146 &' CTD clim salinity: ',mult_ctdsclim
147 call print_message( msgbuf, standardmessageunit,
148 & SQUEEZE_RIGHT , mythid)
149 write(msgbuf,'(a,e10.3)')
150 &' XBT Temperature: ',mult_xbt
151 call print_message( msgbuf, standardmessageunit,
152 & SQUEEZE_RIGHT , mythid)
153 write(msgbuf,'(a,e10.3)')
154 &' ARGO Temperature: ',mult_argot
155 call print_message( msgbuf, standardmessageunit,
156 & SQUEEZE_RIGHT , mythid)
157 write(msgbuf,'(a,e10.3)')
158 &' ARGO Salt: ',mult_argos
159 call print_message( msgbuf, standardmessageunit,
160 & SQUEEZE_RIGHT , mythid)
161 write(msgbuf,'(a,e10.3)')
162 &' drifter velocities: ',mult_drift
163 call print_message( msgbuf, standardmessageunit,
164 & SQUEEZE_RIGHT , mythid)
165 write(msgbuf,'(a,e10.3)')
166 &' drift between last and 1st year:',mult_tdrift
167 call print_message( msgbuf, standardmessageunit,
168 & SQUEEZE_RIGHT , mythid)
169 write(msgbuf,'(a,e10.3)')
170 &' drift between last and 1st year:',mult_sdrift
171 call print_message( msgbuf, standardmessageunit,
172 & SQUEEZE_RIGHT , mythid)
173 write(msgbuf,'(a,e10.3)')
174 &' Ageostrophic bdy flow: ',mult_ageos
175 call print_message( msgbuf, standardmessageunit,
176 & SQUEEZE_RIGHT , mythid)
177 write(msgbuf,'(a,e10.3)')
178 &' OB North: ',mult_obcsn
179 call print_message( msgbuf, standardmessageunit,
180 & SQUEEZE_RIGHT , mythid)
181 write(msgbuf,'(a,e10.3)')
182 &' OB South: ',mult_obcss
183 call print_message( msgbuf, standardmessageunit,
184 & SQUEEZE_RIGHT , mythid)
185 write(msgbuf,'(a,e10.3)')
186 &' OB West: ',mult_obcsw
187 call print_message( msgbuf, standardmessageunit,
188 & SQUEEZE_RIGHT , mythid)
189 write(msgbuf,'(a,e10.3)')
190 &' OB East: ',mult_obcse
191 call print_message( msgbuf, standardmessageunit,
192 & SQUEEZE_RIGHT , mythid)
193
194
195 write(msgbuf,'(a)')
196 &' '
197 call print_message( msgbuf, standardmessageunit,
198 & SQUEEZE_RIGHT , mythid)
199 write(msgbuf,'(a)')
200 &' '
201 call print_message( msgbuf, standardmessageunit,
202 & SQUEEZE_RIGHT , mythid)
203
204 il = ilnblnk(tdatfile)
205 write(msgbuf,'(a,a)')
206 &' Temperature data are read from: ',tdatfile(1:il)
207 call print_message( msgbuf, standardmessageunit,
208 & SQUEEZE_RIGHT , mythid)
209 il = ilnblnk(sdatfile)
210 write(msgbuf,'(a,a)')
211 &' Salinity data are read from: ',sdatfile(1:il)
212 call print_message( msgbuf, standardmessageunit,
213 & SQUEEZE_RIGHT , mythid)
214
215 #ifdef ALLOW_SSH_MEAN_COST_CONTRIBUTION
216 if (using_topex) then
217 il = ilnblnk(topexmeanfile)
218 write(msgbuf,'(a,a)')
219 & ' T/P mean is read from: ',topexmeanfile(1:il)
220 call print_message( msgbuf, standardmessageunit,
221 & SQUEEZE_RIGHT , mythid)
222 endif
223 #endif /* ALLOW_SSH_MEAN_COST_CONTRIBUTION */
224
225 #ifdef ALLOW_SSH_TPANOM_COST_CONTRIBUTION
226 if (using_topex) then
227 il = ilnblnk(topexfile)
228 write(msgbuf,'(a,a)')
229 & ' T/P data are read from: ',topexfile(1:il)
230 call print_message( msgbuf, standardmessageunit,
231 & SQUEEZE_RIGHT , mythid)
232 write(msgbuf,'(a,i10,i8)')
233 & ' T/P start date is: ',
234 & topexstartdate(1),
235 & topexstartdate(2)
236 call print_message( msgbuf, standardmessageunit,
237 & SQUEEZE_RIGHT , mythid)
238 write(msgbuf,'(a,e18.11)')
239 & ' T/P sampling period is: ',
240 & topexperiod
241 call print_message( msgbuf, standardmessageunit,
242 & SQUEEZE_RIGHT , mythid)
243 endif
244 #endif /* ALLOW_SSH_TPANOM_COST_CONTRIBUTION */
245
246 #ifdef ALLOW_SSH_ERSANOM_COST_CONTRIBUTION
247 if (using_ers) then
248 il = ilnblnk(ersfile)
249 write(msgbuf,'(a,a)')
250 & ' ERS data are read from: ',ersfile(1:il)
251 call print_message( msgbuf, standardmessageunit,
252 & SQUEEZE_RIGHT , mythid)
253 write(msgbuf,'(a,i10,i8)')
254 & ' ERS start date is: ',
255 & ersstartdate(1),
256 & ersstartdate(2)
257 call print_message( msgbuf, standardmessageunit,
258 & SQUEEZE_RIGHT , mythid)
259 write(msgbuf,'(a,e18.11)')
260 & ' ERS sampling period is: ',
261 & ersperiod
262 call print_message( msgbuf, standardmessageunit,
263 & SQUEEZE_RIGHT , mythid)
264 endif
265 #endif /* ALLOW_SSH_ERSANOM_COST_CONTRIBUTION */
266
267 #ifdef ALLOW_SSH_GFOANOM_COST_CONTRIBUTION
268 if (using_gfo) then
269 il = ilnblnk(gfofile)
270 write(msgbuf,'(a,a)')
271 & ' GFO data are read from: ',gfofile(1:il)
272 call print_message( msgbuf, standardmessageunit,
273 & SQUEEZE_RIGHT , mythid)
274 write(msgbuf,'(a,i10,i8)')
275 & ' GFO start date is: ',
276 & gfostartdate(1),
277 & gfostartdate(2)
278 call print_message( msgbuf, standardmessageunit,
279 & SQUEEZE_RIGHT , mythid)
280 write(msgbuf,'(a,e18.11)')
281 & ' GFO sampling period is: ',
282 & gfoperiod
283 call print_message( msgbuf, standardmessageunit,
284 & SQUEEZE_RIGHT , mythid)
285 endif
286 #endif /* ALLOW_SSH_GFOANOM_COST_CONTRIBUTION */
287
288 return
289 end
290
  ViewVC Help
Powered by ViewVC 1.1.22