/[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.13.6.8 by heimbach, Tue Jun 24 23:05:28 2003 UTC revision 1.23 by molod, Mon Dec 15 23:00:44 2003 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3    
4    #include "PACKAGES_CONFIG.h"
5  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
6  #ifdef ALLOW_OBCS  #ifdef ALLOW_OBCS
7  # include "OBCS_OPTIONS.h"  # include "OBCS_OPTIONS.h"
# Line 51  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. 'ATMOSPHERIC' ) 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    C AMM
70    #ifdef ALLOW_FIZHI
71          IF ( useFIZHI ) CALL FIZHI_TENDENCY_APPLY_U(
72         &                      iMin,iMax, jMin,jMax, bi,bj, kLev,
73         &                      myCurrentTime, myThid )
74    #endif /* ALLOW_FIZHI */
75    C AMM
76    
77  C     Add windstress momentum impulse into the top-layer  C     Add windstress momentum impulse into the top-layer
78        IF ( kLev .EQ. kSurface ) THEN        IF ( kLev .EQ. kSurface ) THEN
79         DO j=jMin,jMax         DO j=jMin,jMax
# Line 124  C     number of surface interface layer Line 140  C     number of surface interface layer
140        INTEGER kSurface        INTEGER kSurface
141  CEOP  CEOP
142    
143        if ( buoyancyRelation .eq. 'OCEANICP' ) then        if ( buoyancyRelation .eq. 'ATMOSPHERIC' ) then
144           kSurface = 0
145          elseif ( buoyancyRelation .eq. 'OCEANICP' ) then
146         kSurface = Nr         kSurface = Nr
147        else        else
148         kSurface = 1         kSurface = 1
149        endif        endif
150    
151  C--   Forcing term  C--   Forcing term
152    #ifdef ALLOW_AIM
153          IF ( useAIM ) CALL AIM_TENDENCY_APPLY_V(
154         &                      iMin,iMax, jMin,jMax, bi,bj, kLev,
155         &                      myCurrentTime, myThid )
156    #endif /* ALLOW_AIM */
157    
158    C AMM
159    #ifdef ALLOW_FIZHI
160          IF ( useFIZHI ) CALL FIZHI_TENDENCY_APPLY_V(
161         &                      iMin,iMax, jMin,jMax, bi,bj, kLev,
162         &                      myCurrentTime, myThid )
163    #endif /* ALLOW_FIZHI */
164    C AMM
165  C     Add windstress momentum impulse into the top-layer  C     Add windstress momentum impulse into the top-layer
166        IF ( kLev .EQ. kSurface ) THEN        IF ( kLev .EQ. kSurface ) THEN
167         DO j=jMin,jMax         DO j=jMin,jMax
# Line 204  C     number of surface interface layer Line 235  C     number of surface interface layer
235        INTEGER kSurface        INTEGER kSurface
236  CEOP  CEOP
237    
238        if ( buoyancyRelation .eq. 'OCEANICP' ) then        if ( buoyancyRelation .eq. 'ATMOSPHERIC' ) then
239           kSurface = 0
240          elseif ( buoyancyRelation .eq. 'OCEANICP' ) then
241         kSurface = Nr         kSurface = Nr
242        else        else
243         kSurface = 1         kSurface = 1
244        endif        endif
245    
246  C--   Forcing term  C--   Forcing term
247    #ifdef ALLOW_AIM
248          IF ( useAIM ) CALL AIM_TENDENCY_APPLY_T(
249         &                      iMin,iMax, jMin,jMax, bi,bj, kLev,
250         &                      myCurrentTime, myThid )
251    #endif /* ALLOW_AIM */
252    
253    C AMM
254    #ifdef ALLOW_FIZHI
255          IF ( useFIZHI ) CALL FIZHI_TENDENCY_APPLY_T(
256         &                      iMin,iMax, jMin,jMax, bi,bj, kLev,
257         &                      myCurrentTime, myThid )
258    #endif /* ALLOW_FIZHI */
259    C AMM
260  C     Add heat in top-layer  C     Add heat in top-layer
261        IF ( kLev .EQ. kSurface ) THEN        IF ( kLev .EQ. kSurface ) THEN
262         DO j=jMin,jMax         DO j=jMin,jMax
# Line 294  C     number of surface interface layer Line 340  C     number of surface interface layer
340        INTEGER kSurface        INTEGER kSurface
341  CEOP  CEOP
342    
343        if ( buoyancyRelation .eq. 'OCEANICP' ) then        if ( buoyancyRelation .eq. 'ATMOSPHERIC' ) then
344           kSurface = 0
345          elseif ( buoyancyRelation .eq. 'OCEANICP' ) then
346         kSurface = Nr         kSurface = Nr
347        else        else
348         kSurface = 1         kSurface = 1
# Line 302  CEOP Line 350  CEOP
350    
351    
352  C--   Forcing term  C--   Forcing term
353    #ifdef ALLOW_AIM
354          IF ( useAIM ) CALL AIM_TENDENCY_APPLY_S(
355         &                      iMin,iMax, jMin,jMax, bi,bj, kLev,
356         &                      myCurrentTime, myThid )
357    #endif /* ALLOW_AIM */
358    
359    C AMM
360    #ifdef ALLOW_FIZHI
361          IF ( useFIZHI ) CALL FIZHI_TENDENCY_APPLY_S(
362         &                      iMin,iMax, jMin,jMax, bi,bj, kLev,
363         &                      myCurrentTime, myThid )
364    #endif /* ALLOW_FIZHI */
365    C AMM
366  C     Add fresh-water in top-layer  C     Add fresh-water in top-layer
367        IF ( kLev .EQ. kSurface ) THEN        IF ( kLev .EQ. kSurface ) THEN
368         DO j=jMin,jMax         DO j=jMin,jMax
Legend:
Removed from v.1.13.6.8  
changed lines
  Added in v.1.23
  ViewVC Help
Powered by ViewVC 1.1.22