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

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

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

revision 1.17 by jmc, Fri Dec 5 02:40:04 2003 UTC revision 1.18 by dimitri, Wed Mar 3 05:17:44 2004 UTC
# Line 28  C     === Global variables === Line 28  C     === Global variables ===
28  #include "DYNVARS.h"  #include "DYNVARS.h"
29  #include "GRID.h"  #include "GRID.h"
30  #include "SURFACE.h"  #include "SURFACE.h"
31    #ifdef ALLOW_SEAICE
32    #include "SEAICE.h"
33    #endif /* ALLOW_SEAICE */
34    
35  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
36  C     === Routine arguments ===  C     === Routine arguments ===
# Line 88  C--   Surface restoring term : Line 91  C--   Surface restoring term :
91        IF ( doThetaClimRelax .OR. doSaltClimRelax ) THEN        IF ( doThetaClimRelax .OR. doSaltClimRelax ) THEN
92         DO j = jMin, jMax         DO j = jMin, jMax
93          DO i = iMin, iMax          DO i = iMin, iMax
94    #ifdef ALLOW_SEAICE
95    C     Don't restore under sea-ice
96    C     Heat Flux (restoring term) :
97              surfaceTendencyT(i,j,bi,bj) = surfaceTendencyT(i,j,bi,bj)
98         &      -lambdaThetaClimRelax * (1-AREA(i,j,1,bi,bj)
99         &         *(theta(i,j,kSurface,bi,bj)-SST(i,j,bi,bj))
100    C     Salt Flux (restoring term) :
101              surfaceTendencyS(i,j,bi,bj) = surfaceTendencyS(i,j,bi,bj)
102         &      -lambdaSaltClimRelax * (1-AREA(i,j,1,bi,bj)
103         &         *(salt(i,j,kSurface,bi,bj)-SSS(i,j,bi,bj))
104    #else /* ifndef ALLOW_SEAICE */
105  C     Heat Flux (restoring term) :  C     Heat Flux (restoring term) :
106           IF ( abs(yC(i,j,bi,bj)).LE.latBandClimRelax ) THEN           IF ( abs(yC(i,j,bi,bj)).LE.latBandClimRelax ) THEN
107            surfaceTendencyT(i,j,bi,bj) = surfaceTendencyT(i,j,bi,bj)            surfaceTendencyT(i,j,bi,bj) = surfaceTendencyT(i,j,bi,bj)
# Line 98  C     Salt Flux (restoring term) : Line 112  C     Salt Flux (restoring term) :
112       &      -lambdaSaltClimRelax       &      -lambdaSaltClimRelax
113       &         *(salt(i,j,kSurface,bi,bj)-SSS(i,j,bi,bj))       &         *(salt(i,j,kSurface,bi,bj)-SSS(i,j,bi,bj))
114           ENDIF           ENDIF
115    #endif /* ALLOW_SEAICE */
116          ENDDO          ENDDO
117         ENDDO         ENDDO
118        ENDIF        ENDIF
Legend:
Removed from v.1.17  
changed lines
  Added in v.1.18
  ViewVC Help
Powered by ViewVC 1.1.22