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

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

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

Revision 1.4 - (show annotations) (download)
Fri May 30 02:49:04 2008 UTC (16 years, 1 month ago) by gforget
Branch: MAIN
CVS Tags: checkpoint63p, checkpoint63q, checkpoint63l, checkpoint63m, checkpoint63n, checkpoint63o, checkpoint63h, checkpoint63i, checkpoint63j, checkpoint63k, checkpoint63d, checkpoint63e, checkpoint63f, checkpoint63g, checkpoint63a, checkpoint63b, checkpoint63c, checkpoint60, checkpoint61, checkpoint62, checkpoint63, checkpoint62c, checkpoint62b, checkpoint62a, checkpoint62g, checkpoint62f, checkpoint62e, checkpoint62d, checkpoint62k, checkpoint62j, checkpoint62i, checkpoint62h, checkpoint62o, checkpoint62n, checkpoint62m, checkpoint62l, checkpoint62s, checkpoint62r, checkpoint62q, checkpoint62p, checkpoint62w, checkpoint62v, checkpoint62u, checkpoint62t, checkpoint62z, checkpoint62y, checkpoint62x, 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.3: +2 -2 lines 
o bridging the gap between eddy stress and GM.
  -> eddyTau is replaced with eddyPsi (eddyTau = f x rho0 x eddyPsi)
      along with a change in CPP option (now ALLOW_EDDYPSI).
  -> when using GM w/ GM_AdvForm:
      The total eddy streamfunction (Psi = eddyPsi + K x Slope)
      is applied either in the tracer Eq. or in momentum Eq.
      depending on data.gmredi (intro. GM_InMomAsStress).
  -> ALLOW_EDDYPSI_CONTROL for estimation purpose.
  The key modifications are in model/src/taueddy_external_forcing.F
  pkg/gmredi/gmredi_calc_*F pkg/gmredi/gmredi_*transport.F
1 C $Header: /u/gcmpack/MITgcm/pkg/ecco/cost_internal_params.F,v 1.3 2008/02/02 02:34:14 gforget Exp $
2 C $Name: $
3
4 #include "COST_CPPOPTIONS.h"
5
6
7 subroutine cost_internal_params( myiter, mytime, mythid )
8
9 c ==================================================================
10 c SUBROUTINE cost_internal_params
11 c ==================================================================
12 c
13 c o Evaluate cost function contributions of internal parameters adjustments.
14 c
15 c started: Gael Forget 17-Aug-2007
16 c
17 c ==================================================================
18 c SUBROUTINE cost_internal_params
19 c ==================================================================
20
21 implicit none
22
23 c == global variables ==
24 #include "EEPARAMS.h"
25 #include "SIZE.h"
26 #include "PARAMS.h"
27 #include "GRID.h"
28
29 #include "cal.h"
30 #include "ctrl.h"
31 #include "ctrl_dummy.h"
32 #include "ecco_cost.h"
33
34 c == routine arguments ==
35
36 integer myiter
37 _RL mytime
38 integer mythid
39
40 c == end of interface ==
41
42
43 c-- Eddy stress penalty term
44 #ifdef ALLOW_EDDYPSI_COST_CONTRIBUTION
45 call timer_start('cost_tau_eddy [ECCO SPIN-DOWN]', mythid)
46 call cost_tau_eddy( myiter, mytime, mythid )
47 call timer_stop ('cost_tau_eddy [ECCO SPIN-DOWN]', mythid)
48 #endif
49 _BARRIER
50
51 c-- GM coeffs penalty term
52 #ifdef ALLOW_KAPGM_COST_CONTRIBUTION
53 call timer_start('cost_kapgm [ECCO SPIN-DOWN]', mythid)
54 call cost_kapgm( myiter, mytime, mythid )
55 call timer_stop ('cost_kapgm [ECCO SPIN-DOWN]', mythid)
56 #endif
57 _BARRIER
58
59 c-- REDI coeffs penalty term
60 #ifdef ALLOW_KAPREDI_COST_CONTRIBUTION
61 call timer_start('cost_kapredi [ECCO SPIN-DOWN]', mythid)
62 call cost_kapredi( myiter, mytime, mythid )
63 call timer_stop ('cost_kapredi [ECCO SPIN-DOWN]', mythid)
64 #endif
65 _BARRIER
66
67 c-- DIFFKR coeffs penalty term
68 #ifdef ALLOW_DIFFKR_COST_CONTRIBUTION
69 call timer_start('cost_diffkr [ECCO SPIN-DOWN]', mythid)
70 call cost_diffkr( myiter, mytime, mythid )
71 call timer_stop ('cost_diffkr [ECCO SPIN-DOWN]', mythid)
72 #endif
73 _BARRIER
74
75 c-- Bottom Drag penalty term
76 #ifdef ALLOW_BOTTOMDRAG_COST_CONTRIBUTION
77 call timer_start('cost_bottomdrag [ECCO SPIN-DOWN]', mythid)
78 call cost_bottomdrag( myiter, mytime, mythid )
79 call timer_stop ('cost_bottomdrag [ECCO SPIN-DOWN]', mythid)
80 #endif
81 _BARRIER
82
83 end
  ViewVC Help
Powered by ViewVC 1.1.22