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

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

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

Revision 1.1 - (hide annotations) (download)
Thu Nov 6 22:10:08 2003 UTC (20 years, 7 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 heimbach 1.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