/[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.1 - (show annotations) (download)
Thu Nov 6 22:10:08 2003 UTC (20 years, 6 months ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint57t_post, checkpoint57o_post, checkpoint52l_pre, checkpoint52e_pre, hrcube4, checkpoint57v_post, checkpoint52n_post, checkpoint52j_post, checkpoint53d_post, checkpoint54a_pre, checkpoint57m_post, checkpoint55c_post, checkpoint54e_post, checkpoint52e_post, checkpoint57s_post, checkpoint54a_post, checkpoint53c_post, checkpoint57k_post, checkpoint55d_pre, checkpoint57d_post, checkpoint57g_post, checkpoint57b_post, checkpoint57c_pre, checkpoint55j_post, checkpoint56b_post, checkpoint57i_post, checkpoint57y_post, hrcube_1, checkpoint57e_post, branch-netcdf, checkpoint52d_pre, checkpoint52l_post, checkpoint55h_post, checkpoint52k_post, checkpoint52b_pre, checkpoint57g_pre, checkpoint54b_post, checkpoint53b_pre, checkpoint55b_post, checkpoint54d_post, checkpoint56c_post, checkpoint52m_post, checkpoint57y_pre, checkpoint55, checkpoint53a_post, checkpoint57f_pre, checkpoint57a_post, checkpoint54, checkpoint54f_post, checkpoint53b_post, checkpoint55g_post, checkpoint52a_pre, checkpoint55f_post, checkpoint57r_post, checkpoint58, checkpoint57a_pre, checkpoint55i_post, checkpoint57, checkpoint56, checkpoint53, checkpoint52, checkpoint52d_post, eckpoint57e_pre, checkpoint52a_post, checkpoint57h_done, checkpoint52b_post, checkpoint53g_post, checkpoint52f_post, checkpoint57x_post, checkpoint57n_post, checkpoint52c_post, checkpoint57w_post, checkpoint57p_post, checkpint57u_post, checkpoint57f_post, checkpoint58a_post, checkpoint57q_post, ecco_c52_e35, hrcube5, checkpoint57z_post, checkpoint57c_post, checkpoint55e_post, checkpoint52i_post, checkpoint52j_pre, checkpoint53f_post, checkpoint55a_post, checkpoint53d_pre, checkpoint54c_post, checkpoint57j_post, checkpoint57h_pre, checkpoint57l_post, checkpoint52i_pre, checkpoint51u_post, checkpoint52h_pre, checkpoint52f_pre, checkpoint57h_post, hrcube_2, hrcube_3, checkpoint56a_post, checkpoint55d_post
Branch point for: netcdf-sm0 
o merging from ecco-branch
o pkg/ecco now containes ecco-specific part of cost function
o top level routines the_main_loop, forward_step
  supersede those in model/src/
  previous input data.cost now in data.ecco
  (new namelist ecco_cost_nml)
1 C $Header: /u/gcmpack/MITgcm/pkg/cost/cost_summary.F,v 1.2.6.3 2003/06/19 15:21:16 heimbach Exp $
2
3 #include "COST_CPPOPTIONS.h"
4
5
6 subroutine ecco_cost_summary( mythid )
7
8 c ==================================================================
9 c SUBROUTINE ecco_cost_summary
10 c ==================================================================
11 c
12 c o Summarize the cost function related parts of the ECCO release.
13 c
14 c started: Christian Eckert eckert@mit.edu 30-Jun-1999
15 c
16 c changed: Christian Eckert eckert@mit.edu 25-Feb-2000
17 c
18 c - Restructured the code in order to create a package
19 c for the MITgcmUV.
20 c
21 c ==================================================================
22 c SUBROUTINE ecco_cost_summary
23 c ==================================================================
24
25 implicit none
26
27 c == global variables ==
28
29 #include "EEPARAMS.h"
30 #include "SIZE.h"
31
32 #include "ecco_cost.h"
33
34 c == routine arguments ==
35
36 integer mythid
37
38 c == local variables ==
39
40 integer il
41
42 character*(max_len_mbuf) msgbuf
43
44 c == external ==
45
46 integer ilnblnk
47 external ilnblnk
48
49 c == end of interface ==
50
51 write(msgbuf,'(a)')
52 &' '
53 call print_message( msgbuf, standardmessageunit,
54 & SQUEEZE_RIGHT , mythid)
55 write(msgbuf,'(a)')
56 &'// ======================================================='
57 call print_message( msgbuf, standardmessageunit,
58 & SQUEEZE_RIGHT , mythid)
59 write(msgbuf,'(a)')
60 &'// ECCO cost function configuration >>> START <<<'
61 call print_message( msgbuf, standardmessageunit,
62 & SQUEEZE_RIGHT , mythid)
63 write(msgbuf,'(a)')
64 &'// ======================================================='
65 call print_message( msgbuf, standardmessageunit,
66 & SQUEEZE_RIGHT , mythid)
67 write(msgbuf,'(a)')
68 &' '
69 call print_message( msgbuf, standardmessageunit,
70 & SQUEEZE_RIGHT , mythid)
71
72 write(msgbuf,'(a)')
73 &' Multipliers for the indivdual cost function contributions:'
74 call print_message( msgbuf, standardmessageunit,
75 & SQUEEZE_RIGHT , mythid)
76 write(msgbuf,'(a)')
77 &' '
78 call print_message( msgbuf, standardmessageunit,
79 & SQUEEZE_RIGHT , mythid)
80 write(msgbuf,'(a,e10.3)')
81 &' Net heat flux: ',mult_hflux
82 call print_message( msgbuf, standardmessageunit,
83 & SQUEEZE_RIGHT , mythid)
84 write(msgbuf,'(a,e10.3)')
85 &' Salt flux: ',mult_sflux
86 call print_message( msgbuf, standardmessageunit,
87 & SQUEEZE_RIGHT , mythid)
88 write(msgbuf,'(a,e10.3)')
89 &' Zonal wind stress: ',mult_tauu
90 call print_message( msgbuf, standardmessageunit,
91 & SQUEEZE_RIGHT , mythid)
92 write(msgbuf,'(a,e10.3)')
93 &' Meridional wind stress: ',mult_tauv
94 call print_message( msgbuf, standardmessageunit,
95 & SQUEEZE_RIGHT , mythid)
96 write(msgbuf,'(a,e10.3)')
97 &' Mean sea surface height: ',mult_hmean
98 call print_message( msgbuf, standardmessageunit,
99 & SQUEEZE_RIGHT , mythid)
100 write(msgbuf,'(a,e10.3)')
101 &' Sea surface height anomalies: ',mult_h
102 call print_message( msgbuf, standardmessageunit,
103 & SQUEEZE_RIGHT , mythid)
104 write(msgbuf,'(a,e10.3)')
105 &' Temperature Lev.: ',mult_temp
106 call print_message( msgbuf, standardmessageunit,
107 & SQUEEZE_RIGHT , mythid)
108 write(msgbuf,'(a,e10.3)')
109 &' Salinity Lev.: ',mult_salt
110 call print_message( msgbuf, standardmessageunit,
111 & SQUEEZE_RIGHT , mythid)
112 write(msgbuf,'(a,e10.3)')
113 &' Temperature ini.: ',mult_temp0
114 call print_message( msgbuf, standardmessageunit,
115 & SQUEEZE_RIGHT , mythid)
116 write(msgbuf,'(a,e10.3)')
117 &' Salinity ini.: ',mult_salt0
118 call print_message( msgbuf, standardmessageunit,
119 & SQUEEZE_RIGHT , mythid)
120 write(msgbuf,'(a,e10.3)')
121 &' TMI Sea surface temperature: ',mult_tmi
122 call print_message( msgbuf, standardmessageunit,
123 & SQUEEZE_RIGHT , mythid)
124 write(msgbuf,'(a,e10.3)')
125 &' Sea surface temperature: ',mult_sst
126 call print_message( msgbuf, standardmessageunit,
127 & SQUEEZE_RIGHT , mythid)
128 write(msgbuf,'(a,e10.3)')
129 &' Sea surface salinity: ',mult_sss
130 call print_message( msgbuf, standardmessageunit,
131 & SQUEEZE_RIGHT , mythid)
132 write(msgbuf,'(a,e10.3)')
133 &' CTD temperature: ',mult_ctdt
134 call print_message( msgbuf, standardmessageunit,
135 & SQUEEZE_RIGHT , mythid)
136 write(msgbuf,'(a,e10.3)')
137 &' CTD salinity: ',mult_ctds
138 call print_message( msgbuf, standardmessageunit,
139 & SQUEEZE_RIGHT , mythid)
140 write(msgbuf,'(a,e10.3)')
141 &' CTD clim temperature: ',mult_ctdtclim
142 call print_message( msgbuf, standardmessageunit,
143 & SQUEEZE_RIGHT , mythid)
144 write(msgbuf,'(a,e10.3)')
145 &' CTD clim salinity: ',mult_ctdsclim
146 call print_message( msgbuf, standardmessageunit,
147 & SQUEEZE_RIGHT , mythid)
148 write(msgbuf,'(a,e10.3)')
149 &' XBT Temperature: ',mult_xbt
150 call print_message( msgbuf, standardmessageunit,
151 & SQUEEZE_RIGHT , mythid)
152 write(msgbuf,'(a,e10.3)')
153 &' ARGO Temperature: ',mult_argot
154 call print_message( msgbuf, standardmessageunit,
155 & SQUEEZE_RIGHT , mythid)
156 write(msgbuf,'(a,e10.3)')
157 &' ARGO Salt: ',mult_argos
158 call print_message( msgbuf, standardmessageunit,
159 & SQUEEZE_RIGHT , mythid)
160 write(msgbuf,'(a,e10.3)')
161 &' drifter velocities: ',mult_drift
162 call print_message( msgbuf, standardmessageunit,
163 & SQUEEZE_RIGHT , mythid)
164 write(msgbuf,'(a,e10.3)')
165 &' drift between last and 1st year:',mult_tdrift
166 call print_message( msgbuf, standardmessageunit,
167 & SQUEEZE_RIGHT , mythid)
168 write(msgbuf,'(a,e10.3)')
169 &' drift between last and 1st year:',mult_sdrift
170 call print_message( msgbuf, standardmessageunit,
171 & SQUEEZE_RIGHT , mythid)
172 write(msgbuf,'(a,e10.3)')
173 &' Ageostrophic bdy flow: ',mult_ageos
174 call print_message( msgbuf, standardmessageunit,
175 & SQUEEZE_RIGHT , mythid)
176 write(msgbuf,'(a,e10.3)')
177 &' OB North: ',mult_obcsn
178 call print_message( msgbuf, standardmessageunit,
179 & SQUEEZE_RIGHT , mythid)
180 write(msgbuf,'(a,e10.3)')
181 &' OB South: ',mult_obcss
182 call print_message( msgbuf, standardmessageunit,
183 & SQUEEZE_RIGHT , mythid)
184 write(msgbuf,'(a,e10.3)')
185 &' OB West: ',mult_obcsw
186 call print_message( msgbuf, standardmessageunit,
187 & SQUEEZE_RIGHT , mythid)
188 write(msgbuf,'(a,e10.3)')
189 &' OB East: ',mult_obcse
190 call print_message( msgbuf, standardmessageunit,
191 & SQUEEZE_RIGHT , mythid)
192
193
194 write(msgbuf,'(a)')
195 &' '
196 call print_message( msgbuf, standardmessageunit,
197 & SQUEEZE_RIGHT , mythid)
198 write(msgbuf,'(a)')
199 &' '
200 call print_message( msgbuf, standardmessageunit,
201 & SQUEEZE_RIGHT , mythid)
202
203 il = ilnblnk(tdatfile)
204 write(msgbuf,'(a,a)')
205 &' Temperature data are read from: ',tdatfile(1:il)
206 call print_message( msgbuf, standardmessageunit,
207 & SQUEEZE_RIGHT , mythid)
208 il = ilnblnk(sdatfile)
209 write(msgbuf,'(a,a)')
210 &' Salinity data are read from: ',sdatfile(1:il)
211 call print_message( msgbuf, standardmessageunit,
212 & SQUEEZE_RIGHT , mythid)
213
214 #ifdef ALLOW_SSH_MEAN_COST_CONTRIBUTION
215 if (using_topex) then
216 il = ilnblnk(topexmeanfile)
217 write(msgbuf,'(a,a)')
218 & ' T/P mean is read from: ',topexmeanfile(1:il)
219 call print_message( msgbuf, standardmessageunit,
220 & SQUEEZE_RIGHT , mythid)
221 endif
222 #endif /* ALLOW_SSH_MEAN_COST_CONTRIBUTION */
223
224 #ifdef ALLOW_SSH_TPANOM_COST_CONTRIBUTION
225 if (using_topex) then
226 il = ilnblnk(topexfile)
227 write(msgbuf,'(a,a)')
228 & ' T/P data are read from: ',topexfile(1:il)
229 call print_message( msgbuf, standardmessageunit,
230 & SQUEEZE_RIGHT , mythid)
231 write(msgbuf,'(a,i10,i8)')
232 & ' T/P start date is: ',
233 & topexstartdate(1),
234 & topexstartdate(2)
235 call print_message( msgbuf, standardmessageunit,
236 & SQUEEZE_RIGHT , mythid)
237 write(msgbuf,'(a,e18.11)')
238 & ' T/P sampling period is: ',
239 & topexperiod
240 call print_message( msgbuf, standardmessageunit,
241 & SQUEEZE_RIGHT , mythid)
242 endif
243 #endif /* ALLOW_SSH_TPANOM_COST_CONTRIBUTION */
244
245 #ifdef ALLOW_SSH_ERSANOM_COST_CONTRIBUTION
246 if (using_ers) then
247 il = ilnblnk(ersfile)
248 write(msgbuf,'(a,a)')
249 & ' ERS data are read from: ',ersfile(1:il)
250 call print_message( msgbuf, standardmessageunit,
251 & SQUEEZE_RIGHT , mythid)
252 write(msgbuf,'(a,i10,i8)')
253 & ' ERS start date is: ',
254 & ersstartdate(1),
255 & ersstartdate(2)
256 call print_message( msgbuf, standardmessageunit,
257 & SQUEEZE_RIGHT , mythid)
258 write(msgbuf,'(a,e18.11)')
259 & ' ERS sampling period is: ',
260 & ersperiod
261 call print_message( msgbuf, standardmessageunit,
262 & SQUEEZE_RIGHT , mythid)
263 endif
264 #endif /* ALLOW_SSH_ERSANOM_COST_CONTRIBUTION */
265
266 return
267 end
268
  ViewVC Help
Powered by ViewVC 1.1.22