/[MITgcm]/MITgcm/model/src/external_forcing.F
ViewVC logotype

Diff of /MITgcm/model/src/external_forcing.F

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

revision 1.21 by jmc, Sat Dec 6 00:08:35 2003 UTC revision 1.36 by heimbach, Mon Jan 2 21:17:02 2006 UTC
# Line 3  C $Name$ Line 3  C $Name$
3    
4  #include "PACKAGES_CONFIG.h"  #include "PACKAGES_CONFIG.h"
5  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
 #ifdef ALLOW_OBCS  
 # include "OBCS_OPTIONS.h"  
 #endif  
6    
7  CBOP  CBOP
8  C     !ROUTINE: EXTERNAL_FORCING_U  C     !ROUTINE: EXTERNAL_FORCING_U
9  C     !INTERFACE:  C     !INTERFACE:
10        SUBROUTINE EXTERNAL_FORCING_U(        SUBROUTINE EXTERNAL_FORCING_U(
11       I           iMin, iMax, jMin, jMax,bi,bj,kLev,       I           iMin,iMax, jMin,jMax, bi,bj, kLev,
12       I           myCurrentTime,myThid)       I           myTime, myThid )
13  C     !DESCRIPTION: \bv  C     !DESCRIPTION: \bv
14  C     *==========================================================*  C     *==========================================================*
15  C     | S/R EXTERNAL_FORCING_U                                      C     | S/R EXTERNAL_FORCING_U
16  C     | o Contains problem specific forcing for zonal velocity.    C     | o Contains problem specific forcing for zonal velocity.
17  C     *==========================================================*  C     *==========================================================*
18  C     | Adds terms to gU for forcing by external sources            C     | Adds terms to gU for forcing by external sources
19  C     | e.g. wind stress, bottom friction etc..................    C     | e.g. wind stress, bottom friction etc ...
20  C     *==========================================================*  C     *==========================================================*
21  C     \ev  C     \ev
22    
# Line 35  C     == Global data == Line 32  C     == Global data ==
32    
33  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
34  C     == Routine arguments ==  C     == Routine arguments ==
35  C     iMin - Working range of tile for applying forcing.  C     iMin,iMax :: Working range of x-index for applying forcing.
36  C     iMax  C     jMin,jMax :: Working range of y-index for applying forcing.
37  C     jMin  C     bi,bj     :: Current tile indices
38  C     jMax  C     kLev      :: Current vertical level index
39  C     kLev  C     myTime    :: Current time in simulation
40    C     myThid    :: Thread Id number
41        INTEGER iMin, iMax, jMin, jMax, kLev, bi, bj        INTEGER iMin, iMax, jMin, jMax, kLev, bi, bj
42        _RL myCurrentTime        _RL myTime
43        INTEGER myThid        INTEGER myThid
44    
45  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
46  C     == Local variables ==  C     == Local variables ==
47  C     Loop counters  C     i,j       :: Loop counters
48        INTEGER I, J  C     kSurface  :: index of surface layer
49  C     number of surface interface layer        INTEGER i, j
50        INTEGER kSurface        INTEGER kSurface
51  CEOP  CEOP
52    
53        if ( buoyancyRelation .eq. 'AMTMOSPHERIC' ) then        IF ( fluidIsAir ) THEN
54         kSurface = 0         kSurface = 0
55        elseif ( buoyancyRelation .eq. 'OCEANICP' ) then        ELSEIF ( usingPCoords ) THEN
56         kSurface = Nr         kSurface = Nr
57        else        ELSE
58         kSurface = 1         kSurface = 1
59        endif        ENDIF
60    
61  C--   Forcing term  C--   Forcing term
62  #ifdef ALLOW_AIM  #ifdef ALLOW_AIM
63        IF ( useAIM ) CALL AIM_TENDENCY_APPLY_U(        IF ( useAIM ) CALL AIM_TENDENCY_APPLY_U(
64       &                      iMin,iMax, jMin,jMax, bi,bj, kLev,       &                      iMin,iMax, jMin,jMax, bi,bj, kLev,
65       &                      myCurrentTime, myThid )       &                      myTime, myThid )
66  #endif /* ALLOW_AIM */  #endif /* ALLOW_AIM */
67    
68    #ifdef ALLOW_FIZHI
69          IF ( useFIZHI ) CALL FIZHI_TENDENCY_APPLY_U(
70         &                      iMin,iMax, jMin,jMax, bi,bj, kLev,
71         &                      myTime, myThid )
72    #endif /* ALLOW_FIZHI */
73    
74  C     Add windstress momentum impulse into the top-layer  C     Add windstress momentum impulse into the top-layer
75        IF ( kLev .EQ. kSurface ) THEN        IF ( kLev .EQ. kSurface ) THEN
76         DO j=jMin,jMax  c      DO j=1,sNy
77          DO i=iMin,iMax  C-jmc: Without CD-scheme, this is OK ; but with CD-scheme, needs to cover [0:sNy+1]
78           DO j=0,sNy+1
79            DO i=1,sNx+1
80           gU(i,j,kLev,bi,bj) = gU(i,j,kLev,bi,bj)           gU(i,j,kLev,bi,bj) = gU(i,j,kLev,bi,bj)
81       &   +foFacMom*surfaceTendencyU(i,j,bi,bj)       &   +foFacMom*surfaceForcingU(i,j,bi,bj)
82       &   *_maskW(i,j,kLev,bi,bj)       &   *recip_drF(kLev)*recip_hFacW(i,j,kLev,bi,bj)
83          ENDDO          ENDDO
84         ENDDO         ENDDO
85        ENDIF        ENDIF
86    
87  #if (defined (ALLOW_OBCS) && defined (ALLOW_OBCS_SPONGE))  #if (defined (ALLOW_TAU_EDDY))
88           CALL TAUEDDY_EXTERNAL_FORCING_U(
89         I           iMin,iMax, jMin,jMax, bi,bj, kLev,
90         I           myTime, myThid )
91    #endif
92    
93    #ifdef ALLOW_OBCS
94        IF (useOBCS) THEN        IF (useOBCS) THEN
95         CALL OBCS_SPONGE_U(         CALL OBCS_SPONGE_U(
96       I           iMin, iMax, jMin, jMax,bi,bj,kLev,       I           iMin,iMax, jMin,jMax, bi,bj, kLev,
97       I           myCurrentTime,myThid)       I           myTime, myThid )
98        ENDIF        ENDIF
99  #endif  #endif
100    
101        RETURN        RETURN
102        END        END
103    
104    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
105  CBOP  CBOP
106  C     !ROUTINE: EXTERNAL_FORCING_V  C     !ROUTINE: EXTERNAL_FORCING_V
107  C     !INTERFACE:  C     !INTERFACE:
108        SUBROUTINE EXTERNAL_FORCING_V(        SUBROUTINE EXTERNAL_FORCING_V(
109       I           iMin, iMax, jMin, jMax,bi,bj,kLev,       I           iMin,iMax, jMin,jMax, bi,bj, kLev,
110       I           myCurrentTime,myThid)       I           myTime, myThid )
111  C     !DESCRIPTION: \bv  C     !DESCRIPTION: \bv
112  C     *==========================================================*  C     *==========================================================*
113  C     | S/R EXTERNAL_FORCING_V                                      C     | S/R EXTERNAL_FORCING_V
114  C     | o Contains problem specific forcing for merid velocity.    C     | o Contains problem specific forcing for merid velocity.
115  C     *==========================================================*  C     *==========================================================*
116  C     | Adds terms to gV for forcing by external sources            C     | Adds terms to gV for forcing by external sources
117  C     | e.g. wind stress, bottom friction etc..................    C     | e.g. wind stress, bottom friction etc ...
118  C     *==========================================================*  C     *==========================================================*
119  C     \ev  C     \ev
120    
# Line 116  C     == Global data == Line 130  C     == Global data ==
130    
131  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
132  C     == Routine arguments ==  C     == Routine arguments ==
133  C     iMin - Working range of tile for applying forcing.  C     iMin,iMax :: Working range of x-index for applying forcing.
134  C     iMax  C     jMin,jMax :: Working range of y-index for applying forcing.
135  C     jMin  C     bi,bj     :: Current tile indices
136  C     jMax  C     kLev      :: Current vertical level index
137  C     kLev  C     myTime    :: Current time in simulation
138    C     myThid    :: Thread Id number
139        INTEGER iMin, iMax, jMin, jMax, kLev, bi, bj        INTEGER iMin, iMax, jMin, jMax, kLev, bi, bj
140        _RL myCurrentTime        _RL myTime
141        INTEGER myThid        INTEGER myThid
142    
143  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
144  C     == Local variables ==  C     == Local variables ==
145  C     Loop counters  C     i,j       :: Loop counters
146        INTEGER I, J  C     kSurface  :: index of surface layer
147  C     number of surface interface layer        INTEGER i, j
148        INTEGER kSurface        INTEGER kSurface
149  CEOP  CEOP
150    
151        if ( buoyancyRelation .eq. 'AMTMOSPHERIC' ) then        IF ( fluidIsAir ) THEN
152         kSurface = 0         kSurface = 0
153        elseif ( buoyancyRelation .eq. 'OCEANICP' ) then        ELSEIF ( usingPCoords ) THEN
154         kSurface = Nr         kSurface = Nr
155        else        ELSE
156         kSurface = 1         kSurface = 1
157        endif        ENDIF
158    
159  C--   Forcing term  C--   Forcing term
160  #ifdef ALLOW_AIM  #ifdef ALLOW_AIM
161        IF ( useAIM ) CALL AIM_TENDENCY_APPLY_V(        IF ( useAIM ) CALL AIM_TENDENCY_APPLY_V(
162       &                      iMin,iMax, jMin,jMax, bi,bj, kLev,       &                      iMin,iMax, jMin,jMax, bi,bj, kLev,
163       &                      myCurrentTime, myThid )       &                      myTime, myThid )
164  #endif /* ALLOW_AIM */  #endif /* ALLOW_AIM */
165    
166    #ifdef ALLOW_FIZHI
167          IF ( useFIZHI ) CALL FIZHI_TENDENCY_APPLY_V(
168         &                      iMin,iMax, jMin,jMax, bi,bj, kLev,
169         &                      myTime, myThid )
170    #endif /* ALLOW_FIZHI */
171    
172  C     Add windstress momentum impulse into the top-layer  C     Add windstress momentum impulse into the top-layer
173        IF ( kLev .EQ. kSurface ) THEN        IF ( kLev .EQ. kSurface ) THEN
174         DO j=jMin,jMax         DO j=1,sNy+1
175          DO i=iMin,iMax  c       DO i=1,sNx
176    C-jmc: Without CD-scheme, this is OK ; but with CD-scheme, needs to cover [0:sNx+1]
177            DO i=0,sNx+1
178           gV(i,j,kLev,bi,bj) = gV(i,j,kLev,bi,bj)           gV(i,j,kLev,bi,bj) = gV(i,j,kLev,bi,bj)
179       &   +foFacMom*surfaceTendencyV(i,j,bi,bj)       &   +foFacMom*surfaceForcingV(i,j,bi,bj)
180       &   *_maskS(i,j,kLev,bi,bj)       &   *recip_drF(kLev)*recip_hFacS(i,j,kLev,bi,bj)
181          ENDDO          ENDDO
182         ENDDO         ENDDO
183        ENDIF        ENDIF
184    
185  #if (defined (ALLOW_OBCS) && defined (ALLOW_OBCS_SPONGE))  #if (defined (ALLOW_TAU_EDDY))
186           CALL TAUEDDY_EXTERNAL_FORCING_V(
187         I           iMin,iMax, jMin,jMax, bi,bj, kLev,
188         I           myTime, myThid )
189    #endif
190    
191    #ifdef ALLOW_OBCS
192        IF (useOBCS) THEN        IF (useOBCS) THEN
193         CALL OBCS_SPONGE_V(         CALL OBCS_SPONGE_V(
194       I           iMin, iMax, jMin, jMax,bi,bj,kLev,       I           iMin,iMax, jMin,jMax, bi,bj, kLev,
195       I           myCurrentTime,myThid)       I           myTime, myThid )
196        ENDIF        ENDIF
197  #endif  #endif
198    
199        RETURN        RETURN
200        END        END
201    
202    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
203  CBOP  CBOP
204  C     !ROUTINE: EXTERNAL_FORCING_T  C     !ROUTINE: EXTERNAL_FORCING_T
205  C     !INTERFACE:  C     !INTERFACE:
206        SUBROUTINE EXTERNAL_FORCING_T(        SUBROUTINE EXTERNAL_FORCING_T(
207       I           iMin, iMax, jMin, jMax,bi,bj,kLev,       I           iMin,iMax, jMin,jMax, bi,bj, kLev,
208       I           myCurrentTime,myThid)       I           myTime, myThid )
209  C     !DESCRIPTION: \bv  C     !DESCRIPTION: \bv
210  C     *==========================================================*  C     *==========================================================*
211  C     | S/R EXTERNAL_FORCING_T                                      C     | S/R EXTERNAL_FORCING_T
212  C     | o Contains problem specific forcing for temperature.        C     | o Contains problem specific forcing for temperature.
213  C     *==========================================================*  C     *==========================================================*
214  C     | Adds terms to gT for forcing by external sources            C     | Adds terms to gT for forcing by external sources
215  C     | e.g. heat flux, climatalogical relaxation..............    C     | e.g. heat flux, climatalogical relaxation, etc ...
216  C     *==========================================================*  C     *==========================================================*
217  C     \ev  C     \ev
218    
# Line 194  C     == Global data == Line 225  C     == Global data ==
225  #include "GRID.h"  #include "GRID.h"
226  #include "DYNVARS.h"  #include "DYNVARS.h"
227  #include "FFIELDS.h"  #include "FFIELDS.h"
 #ifdef SHORTWAVE_HEATING  
       integer two  
       _RL minusone  
       parameter (two=2,minusone=-1.)  
       _RL swfracb(two)  
 #endif  
228    
229  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
230  C     == Routine arguments ==  C     == Routine arguments ==
231  C     iMin - Working range of tile for applying forcing.  C     iMin,iMax :: Working range of x-index for applying forcing.
232  C     iMax  C     jMin,jMax :: Working range of y-index for applying forcing.
233  C     jMin  C     bi,bj     :: Current tile indices
234  C     jMax  C     kLev      :: Current vertical level index
235  C     kLev  C     myTime    :: Current time in simulation
236    C     myThid    :: Thread Id number
237        INTEGER iMin, iMax, jMin, jMax, kLev, bi, bj        INTEGER iMin, iMax, jMin, jMax, kLev, bi, bj
238        _RL myCurrentTime        _RL myTime
239        INTEGER myThid        INTEGER myThid
 CEndOfInterface  
240    
241  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
242  C     == Local variables ==  C     == Local variables ==
243  C     Loop counters  C     i,j       :: Loop counters
244        INTEGER I, J  C     kSurface  :: index of surface layer
245  C     number of surface interface layer        INTEGER i, j
246        INTEGER kSurface        INTEGER kSurface
247  CEOP  CEOP
248    #ifdef SHORTWAVE_HEATING
249          integer two
250          _RL minusone
251          parameter (two=2,minusone=-1.)
252          _RL swfracb(two)
253          INTEGER kp1
254    #endif
255    
256        if ( buoyancyRelation .eq. 'AMTMOSPHERIC' ) then        IF ( fluidIsAir ) THEN
257         kSurface = 0         kSurface = 0
258        elseif ( buoyancyRelation .eq. 'OCEANICP' ) then        ELSEIF ( usingPCoords ) THEN
259         kSurface = Nr         kSurface = Nr
260        else        ELSE
261         kSurface = 1         kSurface = 1
262        endif        ENDIF
263    
264  C--   Forcing term  C--   Forcing term
265  #ifdef ALLOW_AIM  #ifdef ALLOW_AIM
266        IF ( useAIM ) CALL AIM_TENDENCY_APPLY_T(        IF ( useAIM ) CALL AIM_TENDENCY_APPLY_T(
267       &                      iMin,iMax, jMin,jMax, bi,bj, kLev,       &                      iMin,iMax, jMin,jMax, bi,bj, kLev,
268       &                      myCurrentTime, myThid )       &                      myTime, myThid )
269  #endif /* ALLOW_AIM */  #endif /* ALLOW_AIM */
270    
271    #ifdef ALLOW_FIZHI
272          IF ( useFIZHI ) CALL FIZHI_TENDENCY_APPLY_T(
273         &                      iMin,iMax, jMin,jMax, bi,bj, kLev,
274         &                      myTime, myThid )
275    #endif /* ALLOW_FIZHI */
276    
277  C     Add heat in top-layer  C     Add heat in top-layer
278        IF ( kLev .EQ. kSurface ) THEN        IF ( kLev .EQ. kSurface ) THEN
279         DO j=jMin,jMax         DO j=1,sNy
280          DO i=iMin,iMax          DO i=1,sNx
281           gT(i,j,kLev,bi,bj)=gT(i,j,kLev,bi,bj)           gT(i,j,kLev,bi,bj)=gT(i,j,kLev,bi,bj)
282       &     +maskC(i,j,kLev,bi,bj)*surfaceTendencyT(i,j,bi,bj)       &     +surfaceForcingT(i,j,bi,bj)
283         &     *recip_drF(kLev)*recip_hFacC(i,j,kLev,bi,bj)
284          ENDDO          ENDDO
285         ENDDO         ENDDO
286        ENDIF        ENDIF
287    
288  #ifdef SHORTWAVE_HEATING  #ifdef SHORTWAVE_HEATING
289  C Penetrating SW radiation  C Penetrating SW radiation
290        swfracb(1)=abs(rF(klev))  c     IF ( usePenetratingSW ) THEN
291        swfracb(2)=abs(rF(klev+1))         swfracb(1)=abs(rF(klev))
292        call SWFRAC(         swfracb(2)=abs(rF(klev+1))
293           CALL SWFRAC(
294       I     two,minusone,       I     two,minusone,
295       I     myCurrentTime,myThid,       I     myTime,myThid,
296       U     swfracb)       U     swfracb)
297        DO j=jMin,jMax         kp1 = klev+1
298         DO i=iMin,iMax         IF (klev.EQ.Nr) THEN
299          gT(i,j,klev,bi,bj) = gT(i,j,klev,bi,bj)          kp1 = klev
300       &   -maskC(i,j,klev,bi,bj)*Qsw(i,j,bi,bj)*(swfracb(1)-swfracb(2))          swfracb(2)=0. _d 0
301       &    *recip_Cp*recip_rhoConst*recip_drF(klev)         ENDIF
302           DO j=1,sNy
303            DO i=1,sNx
304             gT(i,j,klev,bi,bj) = gT(i,j,klev,bi,bj)
305         &   -Qsw(i,j,bi,bj)*(swfracb(1)*maskC(i,j,klev,bi,bj)
306         &                   -swfracb(2)*maskC(i,j,kp1, bi,bj))
307         &    *recip_Cp*recip_rhoConst
308         &    *recip_drF(klev)*recip_hFacC(i,j,kLev,bi,bj)
309            ENDDO
310         ENDDO         ENDDO
311        ENDDO  c     ENDIF
312  #endif  #endif
313    
314  #if (defined (ALLOW_OBCS) && defined (ALLOW_OBCS_SPONGE))  #ifdef ALLOW_RBCS
315           if (useRBCS) then
316              call RBCS_ADD_TENDENCY(bi,bj,klev, 1,
317         &                            myTime, myThid )
318           endif
319    #endif
320    
321    #ifdef ALLOW_OBCS
322        IF (useOBCS) THEN        IF (useOBCS) THEN
323         CALL OBCS_SPONGE_T(         CALL OBCS_SPONGE_T(
324       I           iMin, iMax, jMin, jMax,bi,bj,kLev,       I           iMin,iMax, jMin,jMax, bi,bj, kLev,
325       I           myCurrentTime,myThid)       I           myTime, myThid )
326        ENDIF        ENDIF
327  #endif  #endif
328    
329        RETURN        RETURN
330        END        END
331    
332    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
333  CBOP  CBOP
334  C     !ROUTINE: EXTERNAL_FORCING_S  C     !ROUTINE: EXTERNAL_FORCING_S
335  C     !INTERFACE:  C     !INTERFACE:
336        SUBROUTINE EXTERNAL_FORCING_S(        SUBROUTINE EXTERNAL_FORCING_S(
337       I           iMin, iMax, jMin, jMax,bi,bj,kLev,       I           iMin,iMax, jMin,jMax, bi,bj, kLev,
338       I           myCurrentTime,myThid)       I           myTime, myThid )
339    
340  C     !DESCRIPTION: \bv  C     !DESCRIPTION: \bv
341  C     *==========================================================*  C     *==========================================================*
342  C     | S/R EXTERNAL_FORCING_S                                      C     | S/R EXTERNAL_FORCING_S
343  C     | o Contains problem specific forcing for merid velocity.    C     | o Contains problem specific forcing for merid velocity.
344  C     *==========================================================*  C     *==========================================================*
345  C     | Adds terms to gS for forcing by external sources            C     | Adds terms to gS for forcing by external sources
346  C     | e.g. fresh-water flux, climatalogical relaxation.......    C     | e.g. fresh-water flux, climatalogical relaxation, etc ...
347  C     *==========================================================*  C     *==========================================================*
348  C     \ev  C     \ev
349    
# Line 302  C     == Global data == Line 359  C     == Global data ==
359    
360  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
361  C     == Routine arguments ==  C     == Routine arguments ==
362  C     iMin - Working range of tile for applying forcing.  C     iMin,iMax :: Working range of x-index for applying forcing.
363  C     iMax  C     jMin,jMax :: Working range of y-index for applying forcing.
364  C     jMin  C     bi,bj     :: Current tile indices
365  C     jMax  C     kLev      :: Current vertical level index
366  C     kLev  C     myTime    :: Current time in simulation
367    C     myThid    :: Thread Id number
368        INTEGER iMin, iMax, jMin, jMax, kLev, bi, bj        INTEGER iMin, iMax, jMin, jMax, kLev, bi, bj
369        _RL myCurrentTime        _RL myTime
370        INTEGER myThid        INTEGER myThid
371    
372  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
373  C     == Local variables ==  C     == Local variables ==
374  C     Loop counters  C     i,j       :: Loop counters
375        INTEGER I, J  C     kSurface  :: index of surface layer
376  C     number of surface interface layer        INTEGER i, j
377        INTEGER kSurface        INTEGER kSurface
378  CEOP  CEOP
379    
380        if ( buoyancyRelation .eq. 'AMTMOSPHERIC' ) then        IF ( fluidIsAir ) THEN
381         kSurface = 0         kSurface = 0
382        elseif ( buoyancyRelation .eq. 'OCEANICP' ) then        ELSEIF ( usingPCoords ) THEN
383         kSurface = Nr         kSurface = Nr
384        else        ELSE
385         kSurface = 1         kSurface = 1
386        endif        ENDIF
   
387    
388  C--   Forcing term  C--   Forcing term
389  #ifdef ALLOW_AIM  #ifdef ALLOW_AIM
390        IF ( useAIM ) CALL AIM_TENDENCY_APPLY_S(        IF ( useAIM ) CALL AIM_TENDENCY_APPLY_S(
391       &                      iMin,iMax, jMin,jMax, bi,bj, kLev,       &                      iMin,iMax, jMin,jMax, bi,bj, kLev,
392       &                      myCurrentTime, myThid )       &                      myTime, myThid )
393  #endif /* ALLOW_AIM */  #endif /* ALLOW_AIM */
394    
395    #ifdef ALLOW_FIZHI
396          IF ( useFIZHI ) CALL FIZHI_TENDENCY_APPLY_S(
397         &                      iMin,iMax, jMin,jMax, bi,bj, kLev,
398         &                      myTime, myThid )
399    #endif /* ALLOW_FIZHI */
400    
401  C     Add fresh-water in top-layer  C     Add fresh-water in top-layer
402        IF ( kLev .EQ. kSurface ) THEN        IF ( kLev .EQ. kSurface ) THEN
403         DO j=jMin,jMax         DO j=1,sNy
404          DO i=iMin,iMax          DO i=1,sNx
405           gS(i,j,kLev,bi,bj)=gS(i,j,kLev,bi,bj)           gS(i,j,kLev,bi,bj)=gS(i,j,kLev,bi,bj)
406       &   +maskC(i,j,kLev,bi,bj)*surfaceTendencyS(i,j,bi,bj)       &     +surfaceForcingS(i,j,bi,bj)
407         &     *recip_drF(kLev)*recip_hFacC(i,j,kLev,bi,bj)
408          ENDDO          ENDDO
409         ENDDO         ENDDO
410        ENDIF        ENDIF
411    
412  #if (defined (ALLOW_OBCS) && defined (ALLOW_OBCS_SPONGE))  #ifdef ALLOW_RBCS
413           if (useRBCS) then
414              call RBCS_ADD_TENDENCY(bi,bj,klev, 2,
415         &                            myTime, myThid )
416           endif
417    #endif
418    
419    #ifdef ALLOW_OBCS
420        IF (useOBCS) THEN        IF (useOBCS) THEN
421         CALL OBCS_SPONGE_S(         CALL OBCS_SPONGE_S(
422       I           iMin, iMax, jMin, jMax,bi,bj,kLev,       I           iMin,iMax, jMin,jMax, bi,bj, kLev,
423       I           myCurrentTime,myThid)       I           myTime, myThid )
424        ENDIF        ENDIF
425  #endif  #endif
426    
Legend:
Removed from v.1.21  
changed lines
  Added in v.1.36
  ViewVC Help
Powered by ViewVC 1.1.22