Back to home page

MITgcm

 
 

    


File indexing completed on 2021-04-08 05:12:34 UTC

view on githubraw file Latest commit ba0b0470 on 2021-04-08 01:06:32 UTC
6d54cf9ca1 Ed H*0001 #include "TIMEAVE_OPTIONS.h"
                0002 
aa7db3783b Jean*0003 #ifdef ALLOW_TIMEAVE
                0004 
                0005 CBOP
                0006 C     !ROUTINE: TIMEAVE_STATV.h
                0007 C     !INTERFACE:
                0008 C     include "TIMEAVE_STATV.h"
                0009 C     !DESCRIPTION: \bw
                0010 C     *================================================================*
                0011 C     | TIMEAVE_STATV.h
9a31010562 Jean*0012 C     | o Time averages of model state-variables
aa7db3783b Jean*0013 C     |   (common block TAVE_STATEVARS)
                0014 C     *================================================================*
                0015 C     | Time average of state variables is (generally) centered on the
9a31010562 Jean*0016 C     |  middle of the time step (time average interval = timeAve_half)
                0017 C     | Time average of intermediate and tendency variables is centered
                0018 C     |  on the time step (time average interval=timeAve_full)
aa7db3783b Jean*0019 C     *================================================================*
                0020 C     \ev
                0021 CEOP
                0022 
                0023 C     TimeAve_*    :: time of temporal integration (s) *** for each thread ***
9a31010562 Jean*0024 C     timeAve_half :: half time_step multiple (used for state variables)
                0025 C     timeAve_full :: full time_step multiple (used for for intermediate var.)
9d36aec500 Dimi*0026 C     uFluxtave    :: zonal surface wind stress (N/m^2,
                0027 C                     >0 for increase in uVel, i=1 held at western face)
                0028 C     vFluxtave    :: meridional surface wind stress (N/m^2,
                0029 C                     >0 for increase in vVel, j=1 held at southern face)
                0030 C     tFluxtave    :: net surface heat flux (W/m^2, >0 for increase in theta)
                0031 C     sFluxtave    :: net surface salt flux (g/m^2/s, >0 for increase in salt)
aa7db3783b Jean*0032 C     etatave      :: surface displacement (r unit, i.e. ocean:z, atmos:p)
                0033 C     uVeltave     :: zonal velocity (m/s, i=1 held at western face)
                0034 C     vVeltave     :: meridional velocity (m/s, j=1 held at southern face)
                0035 C     wVeltave     :: vertical velocity ([r]/s, i.e.: ocean:m/s atmos:Pa/s)
                0036 C     thetatave    :: potential temperature (oC, held at pressure/tracer point)
ba0b047096 Mart*0037 C     salttave     :: salinity (g/kg, held at pressure/tracer point)
9d36aec500 Dimi*0038 C     Eta2tave     :: eta * eta
aa7db3783b Jean*0039 C     TTtave       :: theta * theta
                0040 C     UUtave       :: uVel * uVel (used to compute the averaged KE)
                0041 C     VVtave       :: vVel * vVel (used to compute the averaged KE)
4237a8e6b2 Jean*0042 C     UVtave       :: uVel * vVel (at vorticity point, i.e. grid-corner)
aa7db3783b Jean*0043 C     KEtave       :: Kinetic Energy
                0044 C     UTtave       :: uVel * theta (* hFacW)
                0045 C     VTtave       :: vVel * theta (* hFacS)
1a4a1e8371 Alis*0046 C     WTtave       :: wVel * theta
03afacb44c Alis*0047 C     UStave       :: uVel * salt (* hFacW)
                0048 C     VStave       :: vVel * salt (* hFacS)
1a4a1e8371 Alis*0049 C     WStave       :: wVel * salt
4237a8e6b2 Jean*0050 C     tDiffRtave   :: vertical diffusion flux of Temperature (theta)
9d36aec500 Dimi*0051 C     uZetatave    :: uVel*Relativ_Vorticity_3 (computed at v point)
                0052 C     vZetatave    :: vVel*Relativ_Vorticity_3 (computed at u point)
aa7db3783b Jean*0053 C     phiHydtave   :: Hydrostatic (ocean) pressure / (atmos) geo- Potential
60c223928f Mart*0054 C     phiHydLowtave:: Hydrostatic (ocean) pressure / (atmos) geo- Potential
                0055 C                     at the fixed boundary: (ocean) bottom pressure
                0056 C                     (atmos) geo- Potential
aa7db3783b Jean*0057 C     ConvectCountTave :: Average number of convective adjustment event
cf6eb84240 Jean*0058 
9a31010562 Jean*0059       COMMON /TAVE_TIME/ timeAve_half, timeAve_full
                0060       _RL timeAve_half(nSx,nSy)
                0061       _RL timeAve_full(nSx,nSy)
cf6eb84240 Jean*0062 
9d36aec500 Dimi*0063       COMMON /TAVE_STATEVARS/
e7f0a9dbe1 Dimi*0064      &                  uFluxtave,vFluxtave,tFluxtave,sFluxtave
                0065      &                 ,etatave,uVeltave,vVeltave,wVeltave
1bf6bedaf8 Dimi*0066      &                 ,thetatave,salttave,phiHydLowtave
1a4a1e8371 Alis*0067      &                 ,UTtave,VTtave,WTtave,UStave,VStave,WStave
e7f0a9dbe1 Dimi*0068      &                 ,Eta2tave,TTtave,UUtave,VVtave,UVtave
52cc3612c0 Jean*0069      &                 ,TdiffRtave
                0070 #ifdef ALLOW_MOM_VECINV
                0071      &                 ,uZetatave, vZetatave
                0072 #endif /* ALLOW_MOM_VECINV */
                0073      &                 ,phiHydtave
1bf6bedaf8 Dimi*0074      &                 ,phiHydLow2Tave
e7f0a9dbe1 Dimi*0075      &                 ,ConvectCountTave
9d36aec500 Dimi*0076       _RL  uFluxtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0077       _RL  vFluxtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0078       _RL  tFluxtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0079       _RL  sFluxtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
aa7db3783b Jean*0080       _RL  etatave  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
cf6eb84240 Jean*0081       _RL  uVeltave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
                0082       _RL  vVeltave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
aa7db3783b Jean*0083       _RL  wVeltave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
cf6eb84240 Jean*0084       _RL  thetatave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
                0085       _RL  salttave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
1bf6bedaf8 Dimi*0086       _RL phiHydLowtave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
e7f0a9dbe1 Dimi*0087       _RL  UTtave   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
                0088       _RL  VTtave   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
1a4a1e8371 Alis*0089       _RL  WTtave   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
03afacb44c Alis*0090       _RL  UStave   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
                0091       _RL  VStave   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
1a4a1e8371 Alis*0092       _RL  WStave   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
e7f0a9dbe1 Dimi*0093       _RL  eta2Tave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
aea29c8517 Alis*0094       _RL  TTtave   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
                0095       _RL  UUtave   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
                0096       _RL  VVtave   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
4237a8e6b2 Jean*0097       _RL  UVtave   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
                0098       _RL TdiffRtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
99bc607d7a Ed H*0099 #ifdef ALLOW_MOM_VECINV
6add245ab0 Jean*0100       _RL uZetatave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
                0101       _RL vZetatave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
e7f0a9dbe1 Dimi*0102 #endif /* ALLOW_MOM_VECINV */
cf6eb84240 Jean*0103       _RL phiHydtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
60c223928f Mart*0104       _RL phiHydLow2Tave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
cf6eb84240 Jean*0105       _RL ConvectCountTave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
aa7db3783b Jean*0106 
                0107 #ifdef NONLIN_FRSURF
                0108 C     hUtave       :: average zonal flow (=hFacW*uVel) (still in m/s !)
                0109 C     hVtave       :: average merid.flow (=hFacS*vVel) (still in m/s !)
                0110 C     hFacCtave    :: average thickness fraction of open water, Center
                0111 C     hFacWtave    :: average thickness fraction of open water, West side
                0112 C     hFacStave    :: average thickness fraction of open water, South side
                0113 
9a31010562 Jean*0114       COMMON /TAVE_THICKNESS/
aa7db3783b Jean*0115      &              hUtave, hVtave
9a31010562 Jean*0116 c    &            , hFacCtave, hFacWtave, hFacStave
aa7db3783b Jean*0117       _RL  hUtave   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
                0118       _RL  hVtave   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
                0119 c     _RL  hFacCtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
                0120 c     _RL  hFacWtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
                0121 c     _RL  hFacStave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
                0122 #endif /* NONLIN_FRSURF */
                0123 
9a31010562 Jean*0124 #endif /* ALLOW_TIMEAVE */