/[MITgcm]/MITgcm/pkg/gmredi/gmredi_diagnostics_fill.F
ViewVC logotype

Annotation of /MITgcm/pkg/gmredi/gmredi_diagnostics_fill.F

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

Revision 1.3 - (hide annotations) (download)
Fri May 30 02:50:16 2008 UTC (15 years, 11 months ago) by gforget
Branch: MAIN
CVS Tags: checkpoint64i, checkpoint64h, checkpoint64a, checkpoint64c, checkpoint64b, checkpoint64e, checkpoint64d, checkpoint64g, checkpoint64f, checkpoint63p, checkpoint63q, checkpoint63r, checkpoint63s, checkpoint63l, checkpoint63m, checkpoint63n, checkpoint63o, checkpoint63h, checkpoint63i, checkpoint63j, checkpoint63k, checkpoint63d, checkpoint63e, checkpoint63f, checkpoint63g, checkpoint63a, checkpoint63b, checkpoint63c, checkpoint64, 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.2: +4 -4 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 gforget 1.3 C $Header: /u/gcmpack/MITgcm/pkg/gmredi/gmredi_diagnostics_fill.F,v 1.2 2007/08/20 17:58:57 gforget Exp $
2 jmc 1.1 C $Name: $
3    
4     #include "GMREDI_OPTIONS.h"
5    
6     CStartOfInterface
7     SUBROUTINE GMREDI_DIAGNOSTICS_FILL(
8     I bi, bj, myThid )
9     C *==========================================================*
10     C | SUBROUTINE GMREDI_DIAGNOSTICS_FILL
11     C | o fill GM-Redi diagnostics
12     C *==========================================================*
13     C | note: formerly was part of S/R GMREDI_CALC_TENSOR
14     C | and was isolated in this S/R for TAF reasons
15     C *==========================================================*
16     IMPLICIT NONE
17    
18     C == Global variables ==
19     #include "SIZE.h"
20     #include "EEPARAMS.h"
21     #include "PARAMS.h"
22     #include "GRID.h"
23     #include "GMREDI.h"
24    
25     C == Routine arguments ==
26     C
27     INTEGER bi,bj
28     INTEGER myThid
29 gforget 1.2 #ifdef ALLOW_DIAGNOSTICS
30     LOGICAL DIAGNOSTICS_IS_ON
31     EXTERNAL DIAGNOSTICS_IS_ON
32     #endif /* ALLOW_DIAGNOSTICS */
33 jmc 1.1 CEndOfInterface
34    
35     #ifdef ALLOW_GMREDI
36    
37     C == Local variables ==
38 gforget 1.3 #ifdef ALLOW_EDDYPSI
39 jmc 1.1 INTEGER i,j,k
40     _RL tmpfld3dloc (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
41     #endif
42    
43     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
44    
45    
46     #ifdef ALLOW_DIAGNOSTICS
47     IF ( useDiagnostics ) THEN
48    
49     #ifdef GM_VISBECK_VARIABLE_K
50     IF ( GM_Visbeck_alpha.NE.0. ) THEN
51     CALL DIAGNOSTICS_FILL(VisbeckK,'GM_VisbK',0,1,1,bi,bj,myThid)
52     ENDIF
53     #endif
54     #ifdef GM_NON_UNITY_DIAGONAL
55     CALL DIAGNOSTICS_FILL(Kux,'GM_Kux ',0,Nr,1,bi,bj,myThid)
56     CALL DIAGNOSTICS_FILL(Kvy,'GM_Kvy ',0,Nr,1,bi,bj,myThid)
57     #endif
58     #ifdef GM_EXTRA_DIAGONAL
59     IF ( GM_ExtraDiag ) THEN
60     CALL DIAGNOSTICS_FILL(Kuz,'GM_Kuz ',0,Nr,1,bi,bj,myThid)
61     CALL DIAGNOSTICS_FILL(Kvz,'GM_Kvz ',0,Nr,1,bi,bj,myThid)
62     ENDIF
63     #endif
64     CALL DIAGNOSTICS_FILL(Kwx,'GM_Kwx ',0,Nr,1,bi,bj,myThid)
65     CALL DIAGNOSTICS_FILL(Kwy,'GM_Kwy ',0,Nr,1,bi,bj,myThid)
66     CALL DIAGNOSTICS_FILL(Kwz,'GM_Kwz ',0,Nr,1,bi,bj,myThid)
67     #ifdef GM_BOLUS_ADVEC
68     IF ( GM_AdvForm ) THEN
69     CALL DIAGNOSTICS_FILL(GM_PsiX,'GM_PsiX ',0,Nr,1,bi,bj,myThid)
70     CALL DIAGNOSTICS_FILL(GM_PsiY,'GM_PsiY ',0,Nr,1,bi,bj,myThid)
71     ENDIF
72     #endif
73    
74 gforget 1.3 #ifdef ALLOW_EDDYPSI
75 jmc 1.1 IF ( DIAGNOSTICS_IS_ON('GMEdTauX',myThid) ) THEN
76     DO k = 1, Nr
77     DO j = 1, sny
78     DO i = 1, snx
79     tmpfld3dloc(i,j,k,bi,bj) =
80     & 0.5*rhoConst*fCori(i,j,bi,bj)*
81     & Kwy(i,j,k,bi,bj)
82     ENDDO
83     ENDDO
84     ENDDO
85     CALL DIAGNOSTICS_FILL(tmpfld3dloc,'GMEdTauX',
86     & 0,Nr,1,bi,bj,myThid)
87     ENDIF
88     c
89     IF ( DIAGNOSTICS_IS_ON('GMEdTauY',myThid) ) THEN
90     DO k = 1, Nr
91     DO j = 1, sny
92     DO i = 1, snx
93     tmpfld3dloc(i,j,k,bi,bj) =
94     & -0.5*rhoConst*fCori(i,j,bi,bj)*
95     & Kwx(i,j,k,bi,bj)
96     ENDDO
97     ENDDO
98     ENDDO
99     CALL DIAGNOSTICS_FILL(tmpfld3dloc,'GMEdTauY',
100     & 0,Nr,1,bi,bj,myThid)
101     ENDIF
102 gforget 1.3 #endif /* ALLOW_EDDYPSI */
103 jmc 1.1
104     ENDIF
105     #endif /* ALLOW_DIAGNOSTICS */
106    
107     #endif /* ALLOW_GMREDI */
108    
109     RETURN
110     END
  ViewVC Help
Powered by ViewVC 1.1.22