pkg/thsice/thsice_slab_ocean.F

C $Header: /u/gcmpack/MITgcm/pkg/thsice/thsice_slab_ocean.F,v 1.7 2007/10/01 13:38:08 jmc Exp $
C $Name:  $

#include "THSICE_OPTIONS.h"

CBOP
C     !ROUTINE: THSICE_SLAB_OCEAN
C     !INTERFACE:
      SUBROUTINE THSICE_SLAB_OCEAN(
     I                      aim_sWght0, aim_sWght1,
     O                      dTsurf,
     I                      bi, bj, myTime, myIter, myThid )
C     !DESCRIPTION: \bv
C     *==========================================================*
C     | S/R  THSICE_SLAB_OCEAN
C     | o Slab ocean for atmosphere (and sea-ice) model
C     *==========================================================*
C     | o add ocean-surface fluxes + restoring term
C     |   and step forward ocean mixed-layer Temp. & Salinity
C     *==========================================================*
C     \ev

C     !USES:
      IMPLICIT NONE

C     == Global variables ==
C-- MITgcm
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "FFIELDS.h"

C-- Sea-Ice package
#include "THSICE_PARAMS.h"
#include "THSICE_VARS.h"

C-- Physics package
#ifdef ALLOW_AIM
#include "AIM_FFIELDS.h"
#endif

C     !INPUT/OUTPUT PARAMETERS:
C     == Routine Arguments ==
C     aim_sWght0 :: weight for time interpolation of surface BC
C     aim_sWght1 :: 0/1 = time period before/after the current time
C     dTsurf     :: diagnostics of slab-ocean temperature change [K/iter]
C     bi,bj      :: tile indices
C     myTime     :: Current time of simulation ( s )
C     myIter     :: Current iteration number in simulation
C     myThid     :: my Thread number Id.
      _RL     aim_sWght0, aim_sWght1
      _RL     dTsurf(sNx,sNy)
      _RL     myTime
      INTEGER bi,bj
      INTEGER myIter, myThid
CEOP

#ifdef ALLOW_THSICE

C     == Local variables ==
C     i,j          :: Loop counters
      _RL dtFac, fwFac, heatFac
#ifdef ALLOW_AIM
      _RL oceTfreez, locTemp, dtFacR
#endif
      INTEGER i,j

cph the following structure is not supported by TAF
cph      IF ( .NOT.stepFwd_oceMxL ) RETURN
      IF ( stepFwd_oceMxL ) THEN

C--    add heat flux and fresh-water + salt flux :
       dtFac   = ocean_deltaT/rhosw
       fwFac   = ocean_deltaT*sMxL_default/rhosw
       heatFac = ocean_deltaT/(cpwater*rhosw)
       DO j=1,sNy
        DO i=1,sNx
         IF ( hOceMxL(i,j,bi,bj).NE.0. _d 0 ) THEN
          dTsurf(i,j) = tOceMxL(i,j,bi,bj)
          tOceMxL(i,j,bi,bj) = tOceMxL(i,j,bi,bj)
     &       - heatFac*Qnet(i,j,bi,bj) / hOceMxL(i,j,bi,bj)
          sOceMxL(i,j,bi,bj) = sOceMxL(i,j,bi,bj)
     &       + (fwFac*EmPmR(i,j,bi,bj) - dtFac*saltFlux(i,j,bi,bj))
     &                                 / hOceMxL(i,j,bi,bj)
         ENDIF
        ENDDO
       ENDDO

#ifdef ALLOW_AIM
       IF ( tauRelax_MxL .GT. 0. _d 0 ) THEN
C--    add restoring (backward) toward climatological Temp.
        dtFac   = ocean_deltaT/tauRelax_MxL
        dtFacR  = 1. _d 0 /(1. _d 0 + dtFac)
        oceTfreez = - 1.9 _d 0
        DO j=1,sNy
         DO i=1,sNx
          IF ( hOceMxL(i,j,bi,bj).NE.0. _d 0 ) THEN
           oceTfreez = -mu_Tf*sOceMxL(i,j,bi,bj)
           locTemp = ( aim_sWght0*aim_sst0(i,j,bi,bj)
     &               + aim_sWght1*aim_sst1(i,j,bi,bj)
     &               ) - celsius2K
           locTemp = MAX( locTemp , oceTfreez )
           tOceMxL(i,j,bi,bj) =
     &         (tOceMxL(i,j,bi,bj) + dtFac*locTemp)*dtFacR
          ENDIF
         ENDDO
        ENDDO
       ENDIF
       IF ( tauRelax_MxL_salt .GT. 0. _d 0 ) THEN
C--    add restoring (backward) toward climatological fixed Salinity
        dtFac   = ocean_deltaT/tauRelax_MxL_salt
        dtFacR  = 1. _d 0 /(1. _d 0 + dtFac)
        DO j=1,sNy
         DO i=1,sNx
          IF ( hOceMxL(i,j,bi,bj).NE.0. _d 0 ) THEN
           sOceMxL(i,j,bi,bj) =
     &         (sOceMxL(i,j,bi,bj) + dtFac*sMxL_default)*dtFacR
          ENDIF
         ENDDO
        ENDDO
       ENDIF
#endif /* ALLOW_AIM */

C-    Diagnose surf. temp. change
       DO j=1,sNy
        DO i=1,sNx
         IF ( hOceMxL(i,j,bi,bj).NE.0. _d 0 ) THEN
          dTsurf(i,j) = tOceMxL(i,j,bi,bj) - dTsurf(i,j)
         ENDIF
        ENDDO
       ENDDO

c-- End of IF ( stepFwd_oceMxL ) THEN
      ENDIF

#endif  /* ALLOW_THSICE */

      RETURN
      END
1	C $Header: /u/gcmpack/MITgcm/pkg/thsice/thsice_slab_ocean.F,v 1.7 2007/10/01 13:38:08 jmc Exp $
2	C $Name: $
3
4	#include "THSICE_OPTIONS.h"
5
6	CBOP
7	C !ROUTINE: THSICE_SLAB_OCEAN
8	C !INTERFACE:
9	SUBROUTINE THSICE_SLAB_OCEAN(
10	I aim_sWght0, aim_sWght1,
11	O dTsurf,
12	I bi, bj, myTime, myIter, myThid )
13	C !DESCRIPTION: \bv
14	C ==========================================================
15	C \| S/R THSICE_SLAB_OCEAN
16	C \| o Slab ocean for atmosphere (and sea-ice) model
17	C ==========================================================
18	C \| o add ocean-surface fluxes + restoring term
19	C \| and step forward ocean mixed-layer Temp. & Salinity
20	C ==========================================================
21	C \ev
22
23	C !USES:
24	IMPLICIT NONE
25
26	C == Global variables ==
27	C-- MITgcm
28	#include "SIZE.h"
29	#include "EEPARAMS.h"
30	#include "PARAMS.h"
31	#include "FFIELDS.h"
32
33	C-- Sea-Ice package
34	#include "THSICE_PARAMS.h"
35	#include "THSICE_VARS.h"
36
37	C-- Physics package
38	#ifdef ALLOW_AIM
39	#include "AIM_FFIELDS.h"
40	#endif
41
42	C !INPUT/OUTPUT PARAMETERS:
43	C == Routine Arguments ==
44	C aim_sWght0 :: weight for time interpolation of surface BC
45	C aim_sWght1 :: 0/1 = time period before/after the current time
46	C dTsurf :: diagnostics of slab-ocean temperature change [K/iter]
47	C bi,bj :: tile indices
48	C myTime :: Current time of simulation ( s )
49	C myIter :: Current iteration number in simulation
50	C myThid :: my Thread number Id.
51	_RL aim_sWght0, aim_sWght1
52	_RL dTsurf(sNx,sNy)
53	_RL myTime
54	INTEGER bi,bj
55	INTEGER myIter, myThid
56	CEOP
57
58	#ifdef ALLOW_THSICE
59
60	C == Local variables ==
61	C i,j :: Loop counters
62	_RL dtFac, fwFac, heatFac
63	#ifdef ALLOW_AIM
64	_RL oceTfreez, locTemp, dtFacR
65	#endif
66	INTEGER i,j
67
68	cph the following structure is not supported by TAF
69	cph IF ( .NOT.stepFwd_oceMxL ) RETURN
70	IF ( stepFwd_oceMxL ) THEN
71
72	C-- add heat flux and fresh-water + salt flux :
73	dtFac = ocean_deltaT/rhosw
74	fwFac = ocean_deltaT*sMxL_default/rhosw
75	heatFac = ocean_deltaT/(cpwater*rhosw)
76	DO j=1,sNy
77	DO i=1,sNx
78	IF ( hOceMxL(i,j,bi,bj).NE.0. _d 0 ) THEN
79	dTsurf(i,j) = tOceMxL(i,j,bi,bj)
80	tOceMxL(i,j,bi,bj) = tOceMxL(i,j,bi,bj)
81	& - heatFac*Qnet(i,j,bi,bj) / hOceMxL(i,j,bi,bj)
82	sOceMxL(i,j,bi,bj) = sOceMxL(i,j,bi,bj)
83	& + (fwFacEmPmR(i,j,bi,bj) - dtFacsaltFlux(i,j,bi,bj))
84	& / hOceMxL(i,j,bi,bj)
85	ENDIF
86	ENDDO
87	ENDDO
88
89	#ifdef ALLOW_AIM
90	IF ( tauRelax_MxL .GT. 0. _d 0 ) THEN
91	C-- add restoring (backward) toward climatological Temp.
92	dtFac = ocean_deltaT/tauRelax_MxL
93	dtFacR = 1. _d 0 /(1. _d 0 + dtFac)
94	oceTfreez = - 1.9 _d 0
95	DO j=1,sNy
96	DO i=1,sNx
97	IF ( hOceMxL(i,j,bi,bj).NE.0. _d 0 ) THEN
98	oceTfreez = -mu_Tf*sOceMxL(i,j,bi,bj)
99	locTemp = ( aim_sWght0*aim_sst0(i,j,bi,bj)
100	& + aim_sWght1*aim_sst1(i,j,bi,bj)
101	& ) - celsius2K
102	locTemp = MAX( locTemp , oceTfreez )
103	tOceMxL(i,j,bi,bj) =
104	& (tOceMxL(i,j,bi,bj) + dtFaclocTemp)dtFacR
105	ENDIF
106	ENDDO
107	ENDDO
108	ENDIF
109	IF ( tauRelax_MxL_salt .GT. 0. _d 0 ) THEN
110	C-- add restoring (backward) toward climatological fixed Salinity
111	dtFac = ocean_deltaT/tauRelax_MxL_salt
112	dtFacR = 1. _d 0 /(1. _d 0 + dtFac)
113	DO j=1,sNy
114	DO i=1,sNx
115	IF ( hOceMxL(i,j,bi,bj).NE.0. _d 0 ) THEN
116	sOceMxL(i,j,bi,bj) =
117	& (sOceMxL(i,j,bi,bj) + dtFacsMxL_default)dtFacR
118	ENDIF
119	ENDDO
120	ENDDO
121	ENDIF
122	#endif /* ALLOW_AIM */
123
124	C- Diagnose surf. temp. change
125	DO j=1,sNy
126	DO i=1,sNx
127	IF ( hOceMxL(i,j,bi,bj).NE.0. _d 0 ) THEN
128	dTsurf(i,j) = tOceMxL(i,j,bi,bj) - dTsurf(i,j)
129	ENDIF
130	ENDDO
131	ENDDO
132
133	c-- End of IF ( stepFwd_oceMxL ) THEN
134	ENDIF
135
136	#endif /* ALLOW_THSICE */
137
138	RETURN
139	END