/[MITgcm]/MITgcm/verification/tidal_basin_2d/code/external_forcing.F
ViewVC logotype

Diff of /MITgcm/verification/tidal_basin_2d/code/external_forcing.F

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

revision 1.7 by jmc, Tue May 6 15:52:51 2014 UTC revision 1.8 by jmc, Sat Jul 5 15:11:36 2014 UTC
# Line 61  CEOP Line 61  CEOP
61        ENDIF        ENDIF
62    
63  C--   Forcing term  C--   Forcing term
 #ifdef ALLOW_AIM  
       IF ( useAIM ) CALL AIM_TENDENCY_APPLY_U(  
      &                      iMin,iMax, jMin,jMax, bi,bj, kLev,  
      &                      myTime, myThid )  
 #endif /* ALLOW_AIM */  
   
 #ifdef ALLOW_FIZHI  
       IF ( useFIZHI ) CALL FIZHI_TENDENCY_APPLY_U(  
      &                      iMin,iMax, jMin,jMax, bi,bj, kLev,  
      &                      myTime, myThid )  
 #endif /* ALLOW_FIZHI */  
64    
65  C     Add windstress momentum impulse into the top-layer  C     Add windstress momentum impulse into the top-layer
66        IF ( kLev .EQ. kSurface ) THEN        IF ( kLev .EQ. kSurface ) THEN
# Line 114  c    &   *min( myTime/86400. , 1.) Line 103  c    &   *min( myTime/86400. , 1.)
103         ENDDO         ENDDO
104        ENDDO        ENDDO
105    
 #if (defined (ALLOW_TAU_EDDY))  
        CALL TAUEDDY_EXTERNAL_FORCING_U(  
      I           iMin,iMax, jMin,jMax, bi,bj, kLev,  
      I           myTime, myThid )  
 #endif  
   
 #ifdef ALLOW_OBCS  
       IF (useOBCS) THEN  
        CALL OBCS_SPONGE_U(  
      I           iMin,iMax, jMin,jMax, bi,bj, kLev,  
      I           myTime, myThid )  
       ENDIF  
 #endif  
   
106        RETURN        RETURN
107        END        END
108    
# Line 187  CEOP Line 162  CEOP
162        ENDIF        ENDIF
163    
164  C--   Forcing term  C--   Forcing term
 #ifdef ALLOW_AIM  
       IF ( useAIM ) CALL AIM_TENDENCY_APPLY_V(  
      &                      iMin,iMax, jMin,jMax, bi,bj, kLev,  
      &                      myTime, myThid )  
 #endif /* ALLOW_AIM */  
   
 #ifdef ALLOW_FIZHI  
       IF ( useFIZHI ) CALL FIZHI_TENDENCY_APPLY_V(  
      &                      iMin,iMax, jMin,jMax, bi,bj, kLev,  
      &                      myTime, myThid )  
 #endif /* ALLOW_FIZHI */  
165    
166  C     Add windstress momentum impulse into the top-layer  C     Add windstress momentum impulse into the top-layer
167        IF ( kLev .EQ. kSurface ) THEN        IF ( kLev .EQ. kSurface ) THEN
# Line 212  C-jmc: Without CD-scheme, this is OK ; b Line 176  C-jmc: Without CD-scheme, this is OK ; b
176         ENDDO         ENDDO
177        ENDIF        ENDIF
178    
 #if (defined (ALLOW_TAU_EDDY))  
        CALL TAUEDDY_EXTERNAL_FORCING_V(  
      I           iMin,iMax, jMin,jMax, bi,bj, kLev,  
      I           myTime, myThid )  
 #endif  
   
 #ifdef ALLOW_OBCS  
       IF (useOBCS) THEN  
        CALL OBCS_SPONGE_V(  
      I           iMin,iMax, jMin,jMax, bi,bj, kLev,  
      I           myTime, myThid )  
       ENDIF  
 #endif  
   
179        RETURN        RETURN
180        END        END
181    
# Line 294  CEOP Line 244  CEOP
244        recip_Cp = 1. _d 0 / HeatCapacity_Cp        recip_Cp = 1. _d 0 / HeatCapacity_Cp
245    
246  C--   Forcing term  C--   Forcing term
 #ifdef ALLOW_AIM  
       IF ( useAIM ) CALL AIM_TENDENCY_APPLY_T(  
      &                      iMin,iMax, jMin,jMax, bi,bj, kLev,  
      &                      myTime, myThid )  
 #endif /* ALLOW_AIM */  
   
 #ifdef ALLOW_FIZHI  
       IF ( useFIZHI ) CALL FIZHI_TENDENCY_APPLY_T(  
      &                      iMin,iMax, jMin,jMax, bi,bj, kLev,  
      &                      myTime, myThid )  
 #endif /* ALLOW_FIZHI */  
247    
248  C     Add heat in top-layer  C     Add heat in top-layer
249        IF ( kLev .EQ. kSurface ) THEN        IF ( kLev .EQ. kSurface ) THEN
# Line 343  c     IF ( usePenetratingSW ) THEN Line 282  c     IF ( usePenetratingSW ) THEN
282  c     ENDIF  c     ENDIF
283  #endif  #endif
284    
 #ifdef ALLOW_OBCS  
       IF (useOBCS) THEN  
        CALL OBCS_SPONGE_T(  
      I           iMin,iMax, jMin,jMax, bi,bj, kLev,  
      I           myTime, myThid )  
       ENDIF  
 #endif  
   
285        RETURN        RETURN
286        END        END
287    
# Line 411  CEOP Line 342  CEOP
342        ENDIF        ENDIF
343    
344  C--   Forcing term  C--   Forcing term
 #ifdef ALLOW_AIM  
       IF ( useAIM ) CALL AIM_TENDENCY_APPLY_S(  
      &                      iMin,iMax, jMin,jMax, bi,bj, kLev,  
      &                      myTime, myThid )  
 #endif /* ALLOW_AIM */  
   
 #ifdef ALLOW_FIZHI  
       IF ( useFIZHI ) CALL FIZHI_TENDENCY_APPLY_S(  
      &                      iMin,iMax, jMin,jMax, bi,bj, kLev,  
      &                      myTime, myThid )  
 #endif /* ALLOW_FIZHI */  
345    
346  C     Add fresh-water in top-layer  C     Add fresh-water in top-layer
347        IF ( kLev .EQ. kSurface ) THEN        IF ( kLev .EQ. kSurface ) THEN
# Line 434  C     Add fresh-water in top-layer Line 354  C     Add fresh-water in top-layer
354         ENDDO         ENDDO
355        ENDIF        ENDIF
356    
 #ifdef ALLOW_OBCS  
       IF (useOBCS) THEN  
        CALL OBCS_SPONGE_S(  
      I           iMin,iMax, jMin,jMax, bi,bj, kLev,  
      I           myTime, myThid )  
       ENDIF  
 #endif  
   
357        RETURN        RETURN
358        END        END

Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8

  ViewVC Help
Powered by ViewVC 1.1.22