pkg/thsice/thsice_main.F

C $Header: /u/gcmpack/MITgcm/pkg/thsice/thsice_main.F,v 1.1 2003/11/23 01:20:13 jmc Exp $
C $Name:  $

#include "THSICE_OPTIONS.h"
 
CBOP
C     !ROUTINE: THSICE_MAIN
C     !INTERFACE:
      SUBROUTINE THSICE_MAIN( 
     I                        myTime, myIter, myThid )
C     !DESCRIPTION: \bv
C     *==========================================================*
C     | S/R  THSICE_MAIN             
C     | o Therm_SeaIce main routine. 
C     |   step forward Thermodynamic_SeaIce variables and modify
C     |    ocean surface forcing accordingly.
C     *==========================================================*

C     !USES:
      IMPLICIT NONE

C     === Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "GRID.h"
#include "DYNVARS.h"
#include "FFIELDS.h"
#include "THSICE_PARAMS.h"
#include "THSICE_VARS.h"
#ifdef ALLOW_BULK_FORCE
#include "BULKF.h"
#endif
 
C     !INPUT/OUTPUT PARAMETERS:
C     === Routine arguments ===
C     myIter :: iteration counter for this thread
C     myTime :: time counter for this thread
C     myThid :: thread number for this instance of the routine.
      _RL  myTime
      INTEGER myIter
      INTEGER myThid
CEOP

#ifdef ALLOW_THSICE
C     !LOCAL VARIABLES:
C     === Local variables ===
      INTEGER i,j
      INTEGER bi,bj
      INTEGER iMin, iMax
      INTEGER jMin, jMax
      _RL prcAtm(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
      _RL evpAtm(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
      _RL flxSW (1-OLx:sNx+OLx,1-OLy:sNy+OLy)

      _RL tauFac

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|

      IF ( stressReduction.GT. 0. _d 0 ) THEN
       iMin = 1-Olx
       iMax = sNx+Olx-1
       jMin = 1-Oly
       jMax = sNy+Oly-1
      ELSE
       iMin = 1
       iMax = sNx
       jMin = 1
       jMax = sNy
      ENDIF

      DO bj=myByLo(myThid),myByHi(myThid)
       DO bi=myBxLo(myThid),myBxHi(myThid)

         DO j = jMin, jMax
          DO i = iMin, iMax
           hOceMxL(i,j,bi,bj) = hfacC(i,j,1,bi,bj)*drF(1)
           tOceMxL(i,j,bi,bj) = theta(i,j,1,bi,bj)
           sOceMxL(i,j,bi,bj) = salt (i,j,1,bi,bj)
           v2ocMxL(i,j,bi,bj) = 
     &              ( uvel(i,j,1,bi,bj)*uvel(i,j,1,bi,bj)
     &              + uvel(i+1,j,1,bi,bj)*uvel(i+1,j,1,bi,bj)
     &              + vvel(i,j+1,1,bi,bj)*vvel(i,j+1,1,bi,bj)
     &              + vvel(i,j,1,bi,bj)*vvel(i,j,1,bi,bj)
     &              )*0.5 _d 0
           prcAtm(i,j) = 0.
           evpAtm(i,j) = 0.
           flxSW (i,j) = 0.
           snowPrc(i,j,bi,bj) = 0. _d 0
#ifdef ALLOW_BULK_FORCE
           prcAtm(i,j) = ( rain(i,j,bi,bj)+runoff(i,j,bi,bj) )*rhofw
           flxSW (i,j) = solar(i,j,bi,bj)
           IF ( iceMask(i,j,bi,bj).GT.0. _d 0
     &       .AND. Tair(i,j,bi,bj).LE.Tf0kel )  THEN
             snowPrc(i,j,bi,bj) = rain(i,j,bi,bj)*rhofw
           ENDIF
#endif
          ENDDO
         ENDDO

         CALL THSICE_STEP_FWD( 
     I                     bi, bj, iMin, iMax, jMin, jMax, 
     I                     prcAtm, 
     U                     evpAtm, flxSW,
     I                     myTime, myIter, myThid )

         CALL THSICE_AVE( 
     I                   evpAtm, flxSW, 
     I                   bi,bj, myTime, myIter, myThid )

c      ENDDO
c     ENDDO

c       IF ( .FALSE. ) THEN
        IF ( stressReduction.GT. 0. _d 0 ) THEN
         DO j = jMin, jMax
          DO i = iMin+1,iMax
            tauFac = stressReduction
     &             *(iceMask(i-1,j,bi,bj)+iceMask(i,j,bi,bj))*0.5 _d 0
            fu(i,j,bi,bj) = (1. _d 0 - tauFac)*fu(i,j,bi,bj)
          ENDDO
         ENDDO
         DO j = jMin+1, jMax
          DO i = iMin, iMax
            tauFac = stressReduction
     &             *(iceMask(i,j-1,bi,bj)+iceMask(i,j,bi,bj))*0.5 _d 0
            fv(i,j,bi,bj) = (1. _d 0 - tauFac)*fv(i,j,bi,bj)
          ENDDO
         ENDDO
        ENDIF

C--  end bi,bj loop
       ENDDO
      ENDDO

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
#endif  /*ALLOW_THSICE*/

      RETURN
      END
1	C $Header: /u/gcmpack/MITgcm/pkg/thsice/thsice_main.F,v 1.1 2003/11/23 01:20:13 jmc Exp $
2	C $Name: $
3
4	#include "THSICE_OPTIONS.h"
5
6	CBOP
7	C !ROUTINE: THSICE_MAIN
8	C !INTERFACE:
9	SUBROUTINE THSICE_MAIN(
10	I myTime, myIter, myThid )
11	C !DESCRIPTION: \bv
12	C ==========================================================
13	C \| S/R THSICE_MAIN
14	C \| o Therm_SeaIce main routine.
15	C \| step forward Thermodynamic_SeaIce variables and modify
16	C \| ocean surface forcing accordingly.
17	C ==========================================================
18
19	C !USES:
20	IMPLICIT NONE
21
22	C === Global variables ===
23	#include "SIZE.h"
24	#include "EEPARAMS.h"
25	#include "PARAMS.h"
26	#include "GRID.h"
27	#include "DYNVARS.h"
28	#include "FFIELDS.h"
29	#include "THSICE_PARAMS.h"
30	#include "THSICE_VARS.h"
31	#ifdef ALLOW_BULK_FORCE
32	#include "BULKF.h"
33	#endif
34
35	C !INPUT/OUTPUT PARAMETERS:
36	C === Routine arguments ===
37	C myIter :: iteration counter for this thread
38	C myTime :: time counter for this thread
39	C myThid :: thread number for this instance of the routine.
40	_RL myTime
41	INTEGER myIter
42	INTEGER myThid
43	CEOP
44
45	#ifdef ALLOW_THSICE
46	C !LOCAL VARIABLES:
47	C === Local variables ===
48	INTEGER i,j
49	INTEGER bi,bj
50	INTEGER iMin, iMax
51	INTEGER jMin, jMax
52	_RL prcAtm(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
53	_RL evpAtm(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
54	_RL flxSW (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
55
56	_RL tauFac
57
58	C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
59
60	IF ( stressReduction.GT. 0. _d 0 ) THEN
61	iMin = 1-Olx
62	iMax = sNx+Olx-1
63	jMin = 1-Oly
64	jMax = sNy+Oly-1
65	ELSE
66	iMin = 1
67	iMax = sNx
68	jMin = 1
69	jMax = sNy
70	ENDIF
71
72	DO bj=myByLo(myThid),myByHi(myThid)
73	DO bi=myBxLo(myThid),myBxHi(myThid)
74
75	DO j = jMin, jMax
76	DO i = iMin, iMax
77	hOceMxL(i,j,bi,bj) = hfacC(i,j,1,bi,bj)*drF(1)
78	tOceMxL(i,j,bi,bj) = theta(i,j,1,bi,bj)
79	sOceMxL(i,j,bi,bj) = salt (i,j,1,bi,bj)
80	v2ocMxL(i,j,bi,bj) =
81	& ( uvel(i,j,1,bi,bj)*uvel(i,j,1,bi,bj)
82	& + uvel(i+1,j,1,bi,bj)*uvel(i+1,j,1,bi,bj)
83	& + vvel(i,j+1,1,bi,bj)*vvel(i,j+1,1,bi,bj)
84	& + vvel(i,j,1,bi,bj)*vvel(i,j,1,bi,bj)
85	& )*0.5 _d 0
86	prcAtm(i,j) = 0.
87	evpAtm(i,j) = 0.
88	flxSW (i,j) = 0.
89	snowPrc(i,j,bi,bj) = 0. _d 0
90	#ifdef ALLOW_BULK_FORCE
91	prcAtm(i,j) = ( rain(i,j,bi,bj)+runoff(i,j,bi,bj) )*rhofw
92	flxSW (i,j) = solar(i,j,bi,bj)
93	IF ( iceMask(i,j,bi,bj).GT.0. _d 0
94	& .AND. Tair(i,j,bi,bj).LE.Tf0kel ) THEN
95	snowPrc(i,j,bi,bj) = rain(i,j,bi,bj)*rhofw
96	ENDIF
97	#endif
98	ENDDO
99	ENDDO
100
101	CALL THSICE_STEP_FWD(
102	I bi, bj, iMin, iMax, jMin, jMax,
103	I prcAtm,
104	U evpAtm, flxSW,
105	I myTime, myIter, myThid )
106
107	CALL THSICE_AVE(
108	I evpAtm, flxSW,
109	I bi,bj, myTime, myIter, myThid )
110
111	c ENDDO
112	c ENDDO
113
114	c IF ( .FALSE. ) THEN
115	IF ( stressReduction.GT. 0. _d 0 ) THEN
116	DO j = jMin, jMax
117	DO i = iMin+1,iMax
118	tauFac = stressReduction
119	& (iceMask(i-1,j,bi,bj)+iceMask(i,j,bi,bj))0.5 _d 0
120	fu(i,j,bi,bj) = (1. _d 0 - tauFac)*fu(i,j,bi,bj)
121	ENDDO
122	ENDDO
123	DO j = jMin+1, jMax
124	DO i = iMin, iMax
125	tauFac = stressReduction
126	& (iceMask(i,j-1,bi,bj)+iceMask(i,j,bi,bj))0.5 _d 0
127	fv(i,j,bi,bj) = (1. _d 0 - tauFac)*fv(i,j,bi,bj)
128	ENDDO
129	ENDDO
130	ENDIF
131
132	C-- end bi,bj loop
133	ENDDO
134	ENDDO
135
136	C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
137	#endif /ALLOW_THSICE/
138
139	RETURN
140	END