/[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.3 - (hide annotations) (download)
Tue Oct 9 00:02:51 2007 UTC (16 years, 9 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 jmc 1.3 C $Header: $
2     C $Name: $
3 heimbach 1.1
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 jmc 1.3 integer mythid
38 heimbach 1.1
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 jmc 1.3 & SQUEEZE_RIGHT , mythid)
199 heimbach 1.1 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 heimbach 1.2 #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 heimbach 1.1 return
289     end
290    
  ViewVC Help
Powered by ViewVC 1.1.22