/[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.39 by heimbach, Wed Jun 28 21:24:54 2006 UTC revision 1.40 by dfer, Fri Jan 5 01:29:31 2007 UTC
# Line 237  C     == Global data == Line 237  C     == Global data ==
237  #include "GRID.h"  #include "GRID.h"
238  #include "DYNVARS.h"  #include "DYNVARS.h"
239  #include "FFIELDS.h"  #include "FFIELDS.h"
240    #include "SURFACE.h"
241    
242  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
243  C     == Routine arguments ==  C     == Routine arguments ==
# Line 303  C     Add heat in top-layer Line 304  C     Add heat in top-layer
304         ENDDO         ENDDO
305        ENDIF        ENDIF
306    
307          IF (linFSConserveTr) THEN
308           DO j=1,sNy
309            DO i=1,sNx
310              IF (kLev .EQ. ksurfC(i,j,bi,bj)) THEN
311                gT(i,j,kLev,bi,bj)=gT(i,j,kLev,bi,bj)
312         &        +TsurfCor*recip_drF(kLev)*_recip_hFacC(i,j,kLev,bi,bj)
313              ENDIF
314            ENDDO
315           ENDDO
316          ENDIF
317    
318  #ifdef ALLOW_SHELFICE  #ifdef ALLOW_SHELFICE
319        IF ( useShelfIce )        IF ( useShelfIce )
320       &     CALL SHELFICE_FORCING_T(       &     CALL SHELFICE_FORCING_T(
# Line 381  C     == Global data == Line 393  C     == Global data ==
393  #include "GRID.h"  #include "GRID.h"
394  #include "DYNVARS.h"  #include "DYNVARS.h"
395  #include "FFIELDS.h"  #include "FFIELDS.h"
396    #include "SURFACE.h"
397    
398  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
399  C     == Routine arguments ==  C     == Routine arguments ==
# Line 439  C     Add fresh-water in top-layer Line 452  C     Add fresh-water in top-layer
452          ENDDO          ENDDO
453         ENDDO         ENDDO
454        ENDIF        ENDIF
455    
456          IF (linFSConserveTr) THEN
457           DO j=1,sNy
458            DO i=1,sNx
459              IF (kLev .EQ. ksurfC(i,j,bi,bj)) THEN
460                gS(i,j,kLev,bi,bj)=gS(i,j,kLev,bi,bj)
461         &        +SsurfCor*recip_drF(kLev)*_recip_hFacC(i,j,kLev,bi,bj)
462              ENDIF
463            ENDDO
464           ENDDO
465          ENDIF
466    
467  #ifdef ALLOW_SHELFICE  #ifdef ALLOW_SHELFICE
468        IF ( useShelfIce )        IF ( useShelfIce )
Legend:
Removed from v.1.39  
changed lines
  Added in v.1.40
  ViewVC Help
Powered by ViewVC 1.1.22