C $Header: /u/gcmpack/MITgcm/pkg/ecco/ecco_cost_summary.F,v 1.12 2014/10/18 18:15:45 gforget Exp $
C $Name: $
#include "ECCO_OPTIONS.h"
subroutine ECCO_COST_SUMMARY( mythid )
c ==================================================================
c SUBROUTINE ecco_cost_summary
c ==================================================================
c
c o Summarize the cost function related parts of the ECCO release.
c
c started: Christian Eckert eckert@mit.edu 30-Jun-1999
c
c changed: Christian Eckert eckert@mit.edu 25-Feb-2000
c
c - Restructured the code in order to create a package
c for the MITgcmUV.
c
c ==================================================================
c SUBROUTINE ecco_cost_summary
c ==================================================================
implicit none
c == global variables ==
#if (defined (ALLOW_ECCO) defined (ECCO_CTRL_DEPRECATED))
#include "EEPARAMS.h"
#include "SIZE.h"
#include "ecco_cost.h"
#ifdef ALLOW_CTRL
# include "CTRL_OBCS.h"
#endif
#endif
c == routine arguments ==
integer mythid
#if (defined (ALLOW_ECCO) defined (ECCO_CTRL_DEPRECATED))
c == local variables ==
integer il
character*(max_len_mbuf) msgbuf
c == external ==
integer ilnblnk
external
c == end of interface ==
write(msgbuf,'(a)')
&' '
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a)')
&'// ======================================================='
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a)')
&'// ECCO cost function configuration >>> START <<<'
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a)')
&'// ======================================================='
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a)')
&' '
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a)')
&' Multipliers for the indivdual cost function contributions:'
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a)')
&' '
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a,e10.3)')
&' Net heat flux: ',mult_hflux
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a,e10.3)')
&' Salt flux: ',mult_sflux
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a,e10.3)')
&' Zonal wind stress: ',mult_tauu
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a,e10.3)')
&' Meridional wind stress: ',mult_tauv
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a,e10.3)')
&' Mean sea surface height: ',mult_hmean
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a,e10.3)')
&' Sea surface height anomalies: ',mult_h
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a,e10.3)')
&' Temperature Lev.: ',mult_temp
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a,e10.3)')
&' Salinity Lev.: ',mult_salt
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a,e10.3)')
&' Temperature ini.: ',mult_temp0
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a,e10.3)')
&' Salinity ini.: ',mult_salt0
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a,e10.3)')
&' Sea level ini.: ',mult_etan0
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a,e10.3)')
&' zonal velocity ini.: ',mult_uvel0
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a,e10.3)')
&' merid velocity ini.: ',mult_vvel0
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a,e10.3)')
&' TMI Sea surface temperature: ',mult_tmi
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a,e10.3)')
&' Sea surface temperature: ',mult_sst
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a,e10.3)')
&' Sea surface salinity: ',mult_sss
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a,e10.3)')
&' CTD temperature: ',mult_ctdt
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a,e10.3)')
&' CTD salinity: ',mult_ctds
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a,e10.3)')
&' CTD clim temperature: ',mult_ctdtclim
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a,e10.3)')
&' CTD clim salinity: ',mult_ctdsclim
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a,e10.3)')
&' XBT Temperature: ',mult_xbt
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a,e10.3)')
&' ARGO Temperature: ',mult_argot
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a,e10.3)')
&' ARGO Salt: ',mult_argos
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a,e10.3)')
&' drifter velocities: ',mult_drift
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a,e10.3)')
&' drift between last and 1st year:',mult_tdrift
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a,e10.3)')
&' drift between last and 1st year:',mult_sdrift
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a,e10.3)')
&' Ageostrophic bdy flow: ',mult_ageos
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
#ifdef ALLOW_CTRL
write(msgbuf,'(a,e10.3)')
&' OB North: ',mult_obcsn
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a,e10.3)')
&' OB South: ',mult_obcss
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a,e10.3)')
&' OB West: ',mult_obcsw
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a,e10.3)')
&' OB East: ',mult_obcse
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
#endif
write(msgbuf,'(a)')
&' '
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a)')
&' '
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
il = ilnblnk(tdatfile)
write(msgbuf,'(a,a)')
&' Temperature data are read from: ',tdatfile(1:il)
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
il = ilnblnk(sdatfile)
write(msgbuf,'(a,a)')
&' Salinity data are read from: ',sdatfile(1:il)
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
#ifdef ALLOW_SSH_MEAN_COST_CONTRIBUTION
write(msgbuf,'(A,L5)')
& ' using_cost_altim[T=mdt,F=no]: ',using_cost_altim
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
if (using_cost_altim) then
il = ilnblnk(mdtdatfile)
write(msgbuf,'(a,a)')
& ' MDT is read from: ',mdtdatfile(1:il)
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a,i10,i8)')
& ' MDT startdate are: ',
& mdtstartdate(1),
& mdtstartdate(2)
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a,i10,i8)')
& ' MDT enddate are: ',
& mdtenddate(1),
& mdtenddate(2)
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
endif
#endif /* ALLOW_SSH_MEAN_COST_CONTRIBUTION */
#ifdef ALLOW_SSH_TPANOM_COST_CONTRIBUTION
if (using_tpj) then
il = ilnblnk(topexfile)
write(msgbuf,'(a,a)')
& ' T/P data are read from: ',topexfile(1:il)
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a,i10,i8)')
& ' T/P start date is: ',
& topexstartdate(1),
& topexstartdate(2)
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a,e18.11)')
& ' T/P sampling period is: ',
& topexperiod
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
endif
#endif /* ALLOW_SSH_TPANOM_COST_CONTRIBUTION */
#ifdef ALLOW_SSH_ERSANOM_COST_CONTRIBUTION
if (using_ers) then
il = ilnblnk(ersfile)
write(msgbuf,'(a,a)')
& ' ERS data are read from: ',ersfile(1:il)
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a,i10,i8)')
& ' ERS start date is: ',
& ersstartdate(1),
& ersstartdate(2)
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a,e18.11)')
& ' ERS sampling period is: ',
& ersperiod
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
endif
#endif /* ALLOW_SSH_ERSANOM_COST_CONTRIBUTION */
#ifdef ALLOW_SSH_GFOANOM_COST_CONTRIBUTION
if (using_gfo) then
il = ilnblnk(gfofile)
write(msgbuf,'(a,a)')
& ' GFO data are read from: ',gfofile(1:il)
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a,i10,i8)')
& ' GFO start date is: ',
& gfostartdate(1),
& gfostartdate(2)
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
write(msgbuf,'(a,e18.11)')
& ' GFO sampling period is: ',
& gfoperiod
call PRINT_MESSAGE( msgbuf, standardmessageunit,
& SQUEEZE_RIGHT , mythid)
endif
#endif /* ALLOW_SSH_GFOANOM_COST_CONTRIBUTION */
#endif /* ALLOW_ECCO and ECCO_CTRL_DEPRECATED */
return
end