44 |
C == Local variables == |
C == Local variables == |
45 |
C Loop counters |
C Loop counters |
46 |
INTEGER I, J |
INTEGER I, J |
47 |
|
C number of surface interface layer |
48 |
|
INTEGER kSurface |
49 |
CEOP |
CEOP |
50 |
|
|
51 |
|
if ( buoyancyRelation .eq. 'OCEANICP' ) then |
52 |
|
kSurface = Nr |
53 |
|
else |
54 |
|
kSurface = 1 |
55 |
|
endif |
56 |
|
|
57 |
C-- Forcing term |
C-- Forcing term |
58 |
C Add windstress momentum impulse into the top-layer |
C Add windstress momentum impulse into the top-layer |
59 |
IF ( kLev .EQ. 1 ) THEN |
IF ( kLev .EQ. kSurface ) THEN |
60 |
DO j=jMin,jMax |
DO j=jMin,jMax |
61 |
DO i=iMin,iMax |
DO i=iMin,iMax |
62 |
gU(i,j,kLev,bi,bj) = gU(i,j,kLev,bi,bj) |
gU(i,j,kLev,bi,bj) = gU(i,j,kLev,bi,bj) |
117 |
C == Local variables == |
C == Local variables == |
118 |
C Loop counters |
C Loop counters |
119 |
INTEGER I, J |
INTEGER I, J |
120 |
|
C number of surface interface layer |
121 |
|
INTEGER kSurface |
122 |
CEOP |
CEOP |
123 |
|
|
124 |
|
if ( buoyancyRelation .eq. 'OCEANICP' ) then |
125 |
|
kSurface = Nr |
126 |
|
else |
127 |
|
kSurface = 1 |
128 |
|
endif |
129 |
|
|
130 |
C-- Forcing term |
C-- Forcing term |
131 |
C Add windstress momentum impulse into the top-layer |
C Add windstress momentum impulse into the top-layer |
132 |
IF ( kLev .EQ. 1 ) THEN |
IF ( kLev .EQ. kSurface ) THEN |
133 |
DO j=jMin,jMax |
DO j=jMin,jMax |
134 |
DO i=iMin,iMax |
DO i=iMin,iMax |
135 |
gV(i,j,kLev,bi,bj) = gV(i,j,kLev,bi,bj) |
gV(i,j,kLev,bi,bj) = gV(i,j,kLev,bi,bj) |
197 |
C == Local variables == |
C == Local variables == |
198 |
C Loop counters |
C Loop counters |
199 |
INTEGER I, J |
INTEGER I, J |
200 |
|
C number of surface interface layer |
201 |
|
INTEGER kSurface |
202 |
CEOP |
CEOP |
203 |
|
|
204 |
|
if ( buoyancyRelation .eq. 'OCEANICP' ) then |
205 |
|
kSurface = Nr |
206 |
|
else |
207 |
|
kSurface = 1 |
208 |
|
endif |
209 |
|
|
210 |
C-- Forcing term |
C-- Forcing term |
211 |
C Add heat in top-layer |
C Add heat in top-layer |
212 |
IF ( kLev .EQ. 1 ) THEN |
IF ( kLev .EQ. kSurface ) THEN |
213 |
DO j=jMin,jMax |
DO j=jMin,jMax |
214 |
DO i=iMin,iMax |
DO i=iMin,iMax |
215 |
gT(i,j,kLev,bi,bj)=gT(i,j,kLev,bi,bj) |
gT(i,j,kLev,bi,bj)=gT(i,j,kLev,bi,bj) |
230 |
DO i=iMin,iMax |
DO i=iMin,iMax |
231 |
gT(i,j,klev,bi,bj) = gT(i,j,klev,bi,bj) |
gT(i,j,klev,bi,bj) = gT(i,j,klev,bi,bj) |
232 |
& -maskC(i,j,klev,bi,bj)*Qsw(i,j,bi,bj)*(swfracb(1)-swfracb(2)) |
& -maskC(i,j,klev,bi,bj)*Qsw(i,j,bi,bj)*(swfracb(1)-swfracb(2)) |
233 |
& *recip_Cp*recip_rhoNil*recip_drF(klev) |
& *recip_Cp*recip_rhoConst*recip_drF(klev) |
234 |
ENDDO |
ENDDO |
235 |
ENDDO |
ENDDO |
236 |
#endif |
#endif |
287 |
C == Local variables == |
C == Local variables == |
288 |
C Loop counters |
C Loop counters |
289 |
INTEGER I, J |
INTEGER I, J |
290 |
|
C number of surface interface layer |
291 |
|
INTEGER kSurface |
292 |
CEOP |
CEOP |
293 |
|
|
294 |
|
if ( buoyancyRelation .eq. 'OCEANICP' ) then |
295 |
|
kSurface = Nr |
296 |
|
else |
297 |
|
kSurface = 1 |
298 |
|
endif |
299 |
|
|
300 |
|
|
301 |
C-- Forcing term |
C-- Forcing term |
302 |
C Add fresh-water in top-layer |
C Add fresh-water in top-layer |
303 |
IF ( kLev .EQ. 1 ) THEN |
IF ( kLev .EQ. kSurface ) THEN |
304 |
DO j=jMin,jMax |
DO j=jMin,jMax |
305 |
DO i=iMin,iMax |
DO i=iMin,iMax |
306 |
gS(i,j,kLev,bi,bj)=gS(i,j,kLev,bi,bj) |
gS(i,j,kLev,bi,bj)=gS(i,j,kLev,bi,bj) |