/[MITgcm]/MITgcm/pkg/generic_advdiff/gad_diagnostics_state.F
ViewVC logotype

Contents of /MITgcm/pkg/generic_advdiff/gad_diagnostics_state.F

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

Revision 1.2 - (show annotations) (download)
Fri May 9 21:43:16 2008 UTC (15 years, 11 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint64y, checkpoint64x, checkpoint64z, checkpoint64q, checkpoint64p, checkpoint64s, checkpoint64r, checkpoint64u, checkpoint64t, checkpoint64w, checkpoint64v, checkpoint64i, checkpoint64h, checkpoint64k, checkpoint64j, checkpoint64m, checkpoint64l, checkpoint64o, checkpoint64n, 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, checkpoint65, checkpoint60, checkpoint61, checkpoint62, checkpoint63, checkpoint66g, checkpoint66f, checkpoint66e, checkpoint66d, checkpoint66c, checkpoint66b, checkpoint66a, checkpoint66o, checkpoint66n, checkpoint66m, checkpoint66l, checkpoint66k, checkpoint66j, checkpoint66i, checkpoint66h, checkpoint65z, checkpoint65x, checkpoint65y, checkpoint65r, checkpoint65s, checkpoint65p, checkpoint65q, checkpoint65v, checkpoint65w, checkpoint65t, checkpoint65u, checkpoint65j, checkpoint65k, checkpoint65h, checkpoint65i, checkpoint65n, checkpoint65o, checkpoint65l, checkpoint65m, checkpoint65b, checkpoint65c, checkpoint65a, checkpoint65f, checkpoint65g, checkpoint65d, checkpoint65e, checkpoint59r, 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, HEAD
Changes since 1.1: +3 -3 lines 
change option: GAD_ALLOW_SOM_ADVECT to GAD_ALLOW_TS_SOM_ADV which only
applies to files where Temperature & salinity 2nd Order moments are used
1 C $Header: /u/gcmpack/MITgcm/pkg/generic_advdiff/gad_diagnostics_state.F,v 1.1 2008/02/12 20:26:02 jmc Exp $
2 C $Name: $
3
4 #include "GAD_OPTIONS.h"
5
6 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
7 CBOP 1
8 C !ROUTINE: GAD_DIAGNOSTICS_STATE
9
10 C !INTERFACE:
11 SUBROUTINE GAD_DIAGNOSTICS_STATE( myTime, myIter, myThid )
12
13 C !DESCRIPTION:
14 C Fill-in the diagnostics array for GAD state variables
15
16 C !USES:
17 IMPLICIT NONE
18 #include "SIZE.h"
19 #include "EEPARAMS.h"
20 #include "PARAMS.h"
21 #include "GRID.h"
22 #include "GAD.h"
23 #include "GAD_SOM_VARS.h"
24
25 C !INPUT PARAMETERS:
26 C myTime :: Current time in simulation
27 C myIter :: Current iteration number in simulation
28 C myThid :: my Thread Id. number
29 _RL myTime
30 INTEGER myIter
31 INTEGER myThid
32 CEOP
33
34 #ifdef ALLOW_DIAGNOSTICS
35 #ifdef GAD_ALLOW_TS_SOM_ADV
36 C ! FUNCTIONS:
37 CHARACTER*4 GAD_DIAG_SUFX
38 EXTERNAL GAD_DIAG_SUFX
39 LOGICAL DIAGNOSTICS_IS_ON
40 EXTERNAL DIAGNOSTICS_IS_ON
41
42 C !LOCAL VARIABLES:
43 CHARACTER*8 diagName
44 CHARACTER*4 diagSufx
45 INTEGER i,j,k,bi,bj,n
46 _RL locVar(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
47 _RL recipVol
48 _RL oneThird, oneFifth, oneNinth
49
50 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
51
52 oneThird = 1. _d 0 / 3. _d 0
53 oneFifth = 1. _d 0 / 5. _d 0
54 oneNinth = 1. _d 0 / 9. _d 0
55
56 C--- Temperature:
57 diagSufx = GAD_DIAG_SUFX( GAD_TEMPERATURE, myThid )
58
59 DO n=1,nSOM
60 diagName = 'SM'//somSfx(n)//diagSufx
61 CALL DIAGNOSTICS_FILL( som_T(1-OLx,1-OLy,1,1,1,n), diagName,
62 & 0, Nr, 0, 1, 1, myThid )
63 ENDDO
64
65 diagName = 'SM_v'//diagSufx
66 IF ( DIAGNOSTICS_IS_ON( diagName, myThid ) ) THEN
67 DO bj = myByLo(myThid), myByHi(myThid)
68 DO bi = myBxLo(myThid), myBxHi(myThid)
69 C-- Compute Internal (sub-grid) variance :
70 C [Sx^2 + Sy^2 + Sz^2]/3 +[Sxx^2 + Syy^2 + Szz^2]/5 +[Sxy^2 + Sxz^2 + Syz^2]/9
71 DO k=1,Nr
72 DO j = 1-OLy,sNy+OLy
73 DO i = 1-OLx,sNx+OLx
74 recipVol = recip_rA(i,j,bi,bj)*recip_deepFac2C(k)
75 & *recip_drF(k)*_recip_hFacC(i,j,k,bi,bj)
76 & *recip_rhoFacC(k)
77 locVar(i,j,k) = (
78 & ( som_T(i,j,k,bi,bj,1)*som_T(i,j,k,bi,bj,1)
79 & + som_T(i,j,k,bi,bj,2)*som_T(i,j,k,bi,bj,2)
80 & + som_T(i,j,k,bi,bj,3)*som_T(i,j,k,bi,bj,3)
81 & )*oneThird
82 & + ( som_T(i,j,k,bi,bj,4)*som_T(i,j,k,bi,bj,4)
83 & + som_T(i,j,k,bi,bj,5)*som_T(i,j,k,bi,bj,5)
84 & + som_T(i,j,k,bi,bj,6)*som_T(i,j,k,bi,bj,6)
85 & )*oneFifth
86 & + ( som_T(i,j,k,bi,bj,7)*som_T(i,j,k,bi,bj,7)
87 & + som_T(i,j,k,bi,bj,8)*som_T(i,j,k,bi,bj,8)
88 & + som_T(i,j,k,bi,bj,9)*som_T(i,j,k,bi,bj,9)
89 & )*oneNinth
90 & )*recipVol*recipVol
91
92 ENDDO
93 ENDDO
94 ENDDO
95 CALL DIAGNOSTICS_FILL( locVar, diagName,
96 & 0, Nr, 2, bi, bj, myThid )
97 ENDDO
98 ENDDO
99 ENDIF
100
101 C--- Salinity:
102 diagSufx = GAD_DIAG_SUFX( GAD_SALINITY, myThid )
103
104 DO n=1,nSOM
105 diagName = 'SM'//somSfx(n)//diagSufx
106 CALL DIAGNOSTICS_FILL( som_S(1-OLx,1-OLy,1,1,1,n), diagName,
107 & 0, Nr, 0, 1, 1, myThid )
108 ENDDO
109
110 diagName = 'SM_v'//diagSufx
111 IF ( DIAGNOSTICS_IS_ON( diagName, myThid ) ) THEN
112 DO bj = myByLo(myThid), myByHi(myThid)
113 DO bi = myBxLo(myThid), myBxHi(myThid)
114 C-- Compute Internal (sub-grid) variance :
115 C [Sx^2 + Sy^2 + Sz^2]/3 +[Sxx^2 + Syy^2 + Szz^2]/5 +[Sxy^2 + Sxz^2 + Syz^2]/9
116 DO k=1,Nr
117 DO j = 1-OLy,sNy+OLy
118 DO i = 1-OLx,sNx+OLx
119 recipVol = recip_rA(i,j,bi,bj)*recip_deepFac2C(k)
120 & *recip_drF(k)*_recip_hFacC(i,j,k,bi,bj)
121 & *recip_rhoFacC(k)
122 locVar(i,j,k) = (
123 & ( som_S(i,j,k,bi,bj,1)*som_S(i,j,k,bi,bj,1)
124 & + som_S(i,j,k,bi,bj,2)*som_S(i,j,k,bi,bj,2)
125 & + som_S(i,j,k,bi,bj,3)*som_S(i,j,k,bi,bj,3)
126 & )*oneThird
127 & + ( som_S(i,j,k,bi,bj,4)*som_S(i,j,k,bi,bj,4)
128 & + som_S(i,j,k,bi,bj,5)*som_S(i,j,k,bi,bj,5)
129 & + som_S(i,j,k,bi,bj,6)*som_S(i,j,k,bi,bj,6)
130 & )*oneFifth
131 & + ( som_S(i,j,k,bi,bj,7)*som_S(i,j,k,bi,bj,7)
132 & + som_S(i,j,k,bi,bj,8)*som_S(i,j,k,bi,bj,8)
133 & + som_S(i,j,k,bi,bj,9)*som_S(i,j,k,bi,bj,9)
134 & )*oneNinth
135 & )*recipVol*recipVol
136 ENDDO
137 ENDDO
138 ENDDO
139 CALL DIAGNOSTICS_FILL( locVar, diagName,
140 & 0, Nr, 2, bi, bj, myThid )
141 ENDDO
142 ENDDO
143 ENDIF
144
145 #endif /* GAD_ALLOW_TS_SOM_ADV */
146 #endif /* ALLOW_DIAGNOSTICS */
147
148 RETURN
149 END
  ViewVC Help
Powered by ViewVC 1.1.22