/[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.20 by edhill, Thu Oct 9 04:19:18 2003 UTC revision 1.21 by jmc, Sat Dec 6 00:08:35 2003 UTC
# Line 52  C     number of surface interface layer Line 52  C     number of surface interface layer
52        INTEGER kSurface        INTEGER kSurface
53  CEOP  CEOP
54    
55        if ( buoyancyRelation .eq. 'OCEANICP' ) then        if ( buoyancyRelation .eq. 'AMTMOSPHERIC' ) then
56           kSurface = 0
57          elseif ( buoyancyRelation .eq. 'OCEANICP' ) then
58         kSurface = Nr         kSurface = Nr
59        else        else
60         kSurface = 1         kSurface = 1
61        endif        endif
62    
63  C--   Forcing term  C--   Forcing term
64    #ifdef ALLOW_AIM
65          IF ( useAIM ) CALL AIM_TENDENCY_APPLY_U(
66         &                      iMin,iMax, jMin,jMax, bi,bj, kLev,
67         &                      myCurrentTime, myThid )
68    #endif /* ALLOW_AIM */
69    
70  C     Add windstress momentum impulse into the top-layer  C     Add windstress momentum impulse into the top-layer
71        IF ( kLev .EQ. kSurface ) THEN        IF ( kLev .EQ. kSurface ) THEN
72         DO j=jMin,jMax         DO j=jMin,jMax
# Line 125  C     number of surface interface layer Line 133  C     number of surface interface layer
133        INTEGER kSurface        INTEGER kSurface
134  CEOP  CEOP
135    
136        if ( buoyancyRelation .eq. 'OCEANICP' ) then        if ( buoyancyRelation .eq. 'AMTMOSPHERIC' ) then
137           kSurface = 0
138          elseif ( buoyancyRelation .eq. 'OCEANICP' ) then
139         kSurface = Nr         kSurface = Nr
140        else        else
141         kSurface = 1         kSurface = 1
142        endif        endif
143    
144  C--   Forcing term  C--   Forcing term
145    #ifdef ALLOW_AIM
146          IF ( useAIM ) CALL AIM_TENDENCY_APPLY_V(
147         &                      iMin,iMax, jMin,jMax, bi,bj, kLev,
148         &                      myCurrentTime, myThid )
149    #endif /* ALLOW_AIM */
150    
151  C     Add windstress momentum impulse into the top-layer  C     Add windstress momentum impulse into the top-layer
152        IF ( kLev .EQ. kSurface ) THEN        IF ( kLev .EQ. kSurface ) THEN
153         DO j=jMin,jMax         DO j=jMin,jMax
# Line 205  C     number of surface interface layer Line 221  C     number of surface interface layer
221        INTEGER kSurface        INTEGER kSurface
222  CEOP  CEOP
223    
224        if ( buoyancyRelation .eq. 'OCEANICP' ) then        if ( buoyancyRelation .eq. 'AMTMOSPHERIC' ) then
225           kSurface = 0
226          elseif ( buoyancyRelation .eq. 'OCEANICP' ) then
227         kSurface = Nr         kSurface = Nr
228        else        else
229         kSurface = 1         kSurface = 1
230        endif        endif
231    
232  C--   Forcing term  C--   Forcing term
233    #ifdef ALLOW_AIM
234          IF ( useAIM ) CALL AIM_TENDENCY_APPLY_T(
235         &                      iMin,iMax, jMin,jMax, bi,bj, kLev,
236         &                      myCurrentTime, myThid )
237    #endif /* ALLOW_AIM */
238    
239  C     Add heat in top-layer  C     Add heat in top-layer
240        IF ( kLev .EQ. kSurface ) THEN        IF ( kLev .EQ. kSurface ) THEN
241         DO j=jMin,jMax         DO j=jMin,jMax
# Line 295  C     number of surface interface layer Line 319  C     number of surface interface layer
319        INTEGER kSurface        INTEGER kSurface
320  CEOP  CEOP
321    
322        if ( buoyancyRelation .eq. 'OCEANICP' ) then        if ( buoyancyRelation .eq. 'AMTMOSPHERIC' ) then
323           kSurface = 0
324          elseif ( buoyancyRelation .eq. 'OCEANICP' ) then
325         kSurface = Nr         kSurface = Nr
326        else        else
327         kSurface = 1         kSurface = 1
# Line 303  CEOP Line 329  CEOP
329    
330    
331  C--   Forcing term  C--   Forcing term
332    #ifdef ALLOW_AIM
333          IF ( useAIM ) CALL AIM_TENDENCY_APPLY_S(
334         &                      iMin,iMax, jMin,jMax, bi,bj, kLev,
335         &                      myCurrentTime, myThid )
336    #endif /* ALLOW_AIM */
337    
338  C     Add fresh-water in top-layer  C     Add fresh-water in top-layer
339        IF ( kLev .EQ. kSurface ) THEN        IF ( kLev .EQ. kSurface ) THEN
340         DO j=jMin,jMax         DO j=jMin,jMax
Legend:
Removed from v.1.20  
changed lines
  Added in v.1.21
  ViewVC Help
Powered by ViewVC 1.1.22