17 |
C *================================================================* |
C *================================================================* |
18 |
C | Time average of state variables is (generally) centered on the |
C | Time average of state variables is (generally) centered on the |
19 |
C | middle of the time step (time average interval = TimeAve_half) |
C | middle of the time step (time average interval = TimeAve_half) |
20 |
C | Time average of intermediate and tandancy variables is centered |
C | Time average of intermediate and tendency variables is centered |
21 |
C | on the time step (time average interval=TimeAve_full) |
C | on the time step (time average interval=TimeAve_full) |
22 |
C *================================================================* |
C *================================================================* |
23 |
C \ev |
C \ev |
26 |
C TimeAve_* :: time of temporal integration (s) *** for each thread *** |
C TimeAve_* :: time of temporal integration (s) *** for each thread *** |
27 |
C TimeAve_half :: half time_step multiple (used for state variables) |
C TimeAve_half :: half time_step multiple (used for state variables) |
28 |
C TimeAve_full :: full time_step multiple (used for for intermediate var.) |
C TimeAve_full :: full time_step multiple (used for for intermediate var.) |
29 |
|
C uFluxtave :: zonal surface wind stress (N/m^2, |
30 |
|
C >0 for increase in uVel, i=1 held at western face) |
31 |
|
C vFluxtave :: meridional surface wind stress (N/m^2, |
32 |
|
C >0 for increase in vVel, j=1 held at southern face) |
33 |
|
C tFluxtave :: net surface heat flux (W/m^2, >0 for increase in theta) |
34 |
|
C sFluxtave :: net surface salt flux (g/m^2/s, >0 for increase in salt) |
35 |
C etatave :: surface displacement (r unit, i.e. ocean:z, atmos:p) |
C etatave :: surface displacement (r unit, i.e. ocean:z, atmos:p) |
36 |
C uVeltave :: zonal velocity (m/s, i=1 held at western face) |
C uVeltave :: zonal velocity (m/s, i=1 held at western face) |
37 |
C vVeltave :: meridional velocity (m/s, j=1 held at southern face) |
C vVeltave :: meridional velocity (m/s, j=1 held at southern face) |
38 |
C wVeltave :: vertical velocity ([r]/s, i.e.: ocean:m/s atmos:Pa/s) |
C wVeltave :: vertical velocity ([r]/s, i.e.: ocean:m/s atmos:Pa/s) |
39 |
C thetatave :: potential temperature (oC, held at pressure/tracer point) |
C thetatave :: potential temperature (oC, held at pressure/tracer point) |
40 |
C salttave :: salinity (ppt, held at pressure/tracer point) |
C salttave :: salinity (ppt, held at pressure/tracer point) |
41 |
C Eta2tave :: eta * eta |
C Eta2tave :: eta * eta |
42 |
C TTtave :: theta * theta |
C TTtave :: theta * theta |
43 |
C UUtave :: uVel * uVel (used to compute the averaged KE) |
C UUtave :: uVel * uVel (used to compute the averaged KE) |
44 |
C VVtave :: vVel * vVel (used to compute the averaged KE) |
C VVtave :: vVel * vVel (used to compute the averaged KE) |
48 |
C VTtave :: vVel * theta (* hFacS) |
C VTtave :: vVel * theta (* hFacS) |
49 |
C WTtave :: wVel * theta |
C WTtave :: wVel * theta |
50 |
C tDiffRtave :: vertical diffusion flux of Temperature (theta) |
C tDiffRtave :: vertical diffusion flux of Temperature (theta) |
51 |
C uZetatave :: uVel*Relativ_Vorticity_3 (computed at v point) |
C uZetatave :: uVel*Relativ_Vorticity_3 (computed at v point) |
52 |
C vZetatave :: vVel*Relativ_Vorticity_3 (computed at u point) |
C vZetatave :: vVel*Relativ_Vorticity_3 (computed at u point) |
53 |
C phiHydtave :: Hydrostatic (ocean) pressure / (atmos) geo- Potential |
C phiHydtave :: Hydrostatic (ocean) pressure / (atmos) geo- Potential |
54 |
C phiHydLowtave:: Hydrostatic (ocean) pressure / (atmos) geo- Potential |
C phiHydLowtave:: Hydrostatic (ocean) pressure / (atmos) geo- Potential |
55 |
C at the fixed boundary: (ocean) bottom pressure |
C at the fixed boundary: (ocean) bottom pressure |
60 |
_RL TimeAve_half(Nr,nSx,nSy) |
_RL TimeAve_half(Nr,nSx,nSy) |
61 |
_RL TimeAve_full(Nr,nSx,nSy) |
_RL TimeAve_full(Nr,nSx,nSy) |
62 |
|
|
63 |
COMMON /TAVE_STATEVARS/ |
COMMON /TAVE_STATEVARS/ |
64 |
|
& uFluxtave,vFluxtave,tFluxtave,sFluxtave, |
65 |
& etatave,Eta2tave, |
& etatave,Eta2tave, |
66 |
& uVeltave,vVeltave,wVeltave, |
& uVeltave,vVeltave,wVeltave, |
67 |
& thetatave,salttave, |
& thetatave,salttave, |
71 |
& phiHydLowtave,phiHydLow2Tave, |
& phiHydLowtave,phiHydLow2Tave, |
72 |
& ConvectCountTave |
& ConvectCountTave |
73 |
c & ,KEtave |
c & ,KEtave |
74 |
|
_RL uFluxtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
75 |
|
_RL vFluxtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
76 |
|
_RL tFluxtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
77 |
|
_RL sFluxtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
78 |
_RL etatave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL etatave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
79 |
_RL eta2Tave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL eta2Tave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
80 |
_RL uVeltave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
_RL uVeltave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |