/[MITgcm]/MITgcm/pkg/thsice/thsice_main.F
ViewVC logotype

Annotation of /MITgcm/pkg/thsice/thsice_main.F

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

Revision 1.3 - (hide annotations) (download)
Thu Jul 8 15:56:36 2004 UTC (19 years, 10 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint55c_post, checkpoint54e_post, checkpoint55d_pre, checkpoint57b_post, checkpoint57c_pre, checkpoint55j_post, checkpoint56b_post, checkpoint55h_post, checkpoint54b_post, checkpoint55b_post, checkpoint54d_post, checkpoint56c_post, checkpoint55, checkpoint57a_post, checkpoint54f_post, checkpoint55g_post, checkpoint55f_post, checkpoint57a_pre, checkpoint55i_post, checkpoint57, checkpoint56, checkpoint57c_post, checkpoint55e_post, checkpoint55a_post, checkpoint54c_post, checkpoint56a_post, checkpoint55d_post
Changes since 1.2: +34 -2 lines 
staggerTimeStep & NLFS: update physical thsice.
1 jmc 1.3 C $Header: /u/gcmpack/MITgcm/pkg/thsice/thsice_main.F,v 1.2 2004/04/07 23:40:34 jmc Exp $
2 jmc 1.1 C $Name: $
3    
4     #include "THSICE_OPTIONS.h"
5    
6 jmc 1.2 CBOP
7 jmc 1.1 C !ROUTINE: THSICE_MAIN
8     C !INTERFACE:
9     SUBROUTINE THSICE_MAIN(
10     I myTime, myIter, myThid )
11 jmc 1.2 C !DESCRIPTION: \bv
12 jmc 1.1 C *==========================================================*
13 jmc 1.2 C | S/R THSICE_MAIN
14 jmc 1.1 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 jmc 1.2
22 jmc 1.1 C === Global variables ===
23     #include "SIZE.h"
24     #include "EEPARAMS.h"
25     #include "PARAMS.h"
26 jmc 1.2 #include "GRID.h"
27 jmc 1.3 #include "SURFACE.h"
28 jmc 1.2 #include "DYNVARS.h"
29 jmc 1.1 #include "FFIELDS.h"
30     #include "THSICE_PARAMS.h"
31 jmc 1.2 #include "THSICE_VARS.h"
32     #ifdef ALLOW_BULK_FORCE
33     #include "BULKF.h"
34     #endif
35 jmc 1.1
36     C !INPUT/OUTPUT PARAMETERS:
37     C === Routine arguments ===
38     C myIter :: iteration counter for this thread
39     C myTime :: time counter for this thread
40     C myThid :: thread number for this instance of the routine.
41     _RL myTime
42     INTEGER myIter
43     INTEGER myThid
44 jmc 1.2 CEOP
45 jmc 1.1
46     #ifdef ALLOW_THSICE
47     C !LOCAL VARIABLES:
48     C === Local variables ===
49     INTEGER i,j
50     INTEGER bi,bj
51     INTEGER iMin, iMax
52     INTEGER jMin, jMax
53 jmc 1.2 _RL prcAtm(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
54     _RL evpAtm(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
55     _RL flxSW (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
56 jmc 1.1
57     _RL tauFac
58    
59     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
60    
61     IF ( stressReduction.GT. 0. _d 0 ) THEN
62     iMin = 1-Olx
63     iMax = sNx+Olx-1
64     jMin = 1-Oly
65     jMax = sNy+Oly-1
66     ELSE
67     iMin = 1
68     iMax = sNx
69     jMin = 1
70     jMax = sNy
71     ENDIF
72    
73     DO bj=myByLo(myThid),myByHi(myThid)
74     DO bi=myBxLo(myThid),myBxHi(myThid)
75 jmc 1.2
76 jmc 1.3 C-- Mixed layer thickness: take the 1rst layer
77     #ifdef NONLIN_FRSURF
78     IF ( staggerTimeStep .AND. nonlinFreeSurf.GT.0 ) THEN
79     IF ( select_rStar.GT.0 ) THEN
80     DO j = jMin, jMax
81     DO i = iMin, iMax
82     hOceMxL(i,j,bi,bj) = drF(1)*h0FacC(i,j,1,bi,bj)
83     & *rStarFacC(i,j,bi,bj)
84     ENDDO
85     ENDDO
86     ELSE
87     DO j = jMin, jMax
88     DO i = iMin, iMax
89     IF ( ksurfC(i,j,bi,bj).EQ.1 ) THEN
90     hOceMxL(i,j,bi,bj) = drF(1)*hFac_surfC(i,j,bi,bj)
91     ELSE
92     hOceMxL(i,j,bi,bj) = drF(1)*hfacC(i,j,1,bi,bj)
93     ENDIF
94     ENDDO
95     ENDDO
96     ENDIF
97     ELSE
98     #else /* ndef NONLIN_FRSURF */
99     IF (.TRUE.) THEN
100     #endif /* NONLIN_FRSURF */
101     DO j = jMin, jMax
102     DO i = iMin, iMax
103     hOceMxL(i,j,bi,bj) = drF(1)*hfacC(i,j,1,bi,bj)
104     ENDDO
105     ENDDO
106     ENDIF
107    
108 jmc 1.2 DO j = jMin, jMax
109     DO i = iMin, iMax
110     tOceMxL(i,j,bi,bj) = theta(i,j,1,bi,bj)
111     sOceMxL(i,j,bi,bj) = salt (i,j,1,bi,bj)
112     v2ocMxL(i,j,bi,bj) =
113     & ( uvel(i,j,1,bi,bj)*uvel(i,j,1,bi,bj)
114     & + uvel(i+1,j,1,bi,bj)*uvel(i+1,j,1,bi,bj)
115     & + vvel(i,j+1,1,bi,bj)*vvel(i,j+1,1,bi,bj)
116     & + vvel(i,j,1,bi,bj)*vvel(i,j,1,bi,bj)
117     & )*0.5 _d 0
118     prcAtm(i,j) = 0.
119     evpAtm(i,j) = 0.
120     flxSW (i,j) = 0.
121     snowPrc(i,j,bi,bj) = 0. _d 0
122     #ifdef ALLOW_BULK_FORCE
123     prcAtm(i,j) = ( rain(i,j,bi,bj)+runoff(i,j,bi,bj) )*rhofw
124     flxSW (i,j) = solar(i,j,bi,bj)
125     IF ( iceMask(i,j,bi,bj).GT.0. _d 0
126     & .AND. Tair(i,j,bi,bj).LE.Tf0kel ) THEN
127     snowPrc(i,j,bi,bj) = rain(i,j,bi,bj)*rhofw
128     ENDIF
129     #endif
130     ENDDO
131     ENDDO
132    
133 jmc 1.1 CALL THSICE_STEP_FWD(
134     I bi, bj, iMin, iMax, jMin, jMax,
135 jmc 1.2 I prcAtm,
136     U evpAtm, flxSW,
137 jmc 1.1 I myTime, myIter, myThid )
138 jmc 1.2
139     CALL THSICE_AVE(
140     I evpAtm, flxSW,
141     I bi,bj, myTime, myIter, myThid )
142    
143 jmc 1.1 c ENDDO
144     c ENDDO
145    
146     c IF ( .FALSE. ) THEN
147     IF ( stressReduction.GT. 0. _d 0 ) THEN
148     DO j = jMin, jMax
149     DO i = iMin+1,iMax
150     tauFac = stressReduction
151     & *(iceMask(i-1,j,bi,bj)+iceMask(i,j,bi,bj))*0.5 _d 0
152     fu(i,j,bi,bj) = (1. _d 0 - tauFac)*fu(i,j,bi,bj)
153     ENDDO
154     ENDDO
155     DO j = jMin+1, jMax
156     DO i = iMin, iMax
157     tauFac = stressReduction
158     & *(iceMask(i,j-1,bi,bj)+iceMask(i,j,bi,bj))*0.5 _d 0
159     fv(i,j,bi,bj) = (1. _d 0 - tauFac)*fv(i,j,bi,bj)
160     ENDDO
161     ENDDO
162     ENDIF
163    
164     C-- end bi,bj loop
165     ENDDO
166     ENDDO
167    
168     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
169     #endif /*ALLOW_THSICE*/
170    
171     RETURN
172     END
  ViewVC Help
Powered by ViewVC 1.1.22