/[MITgcm]/MITgcm/pkg/thsice/thsice_slab_ocean.F
ViewVC logotype

Diff of /MITgcm/pkg/thsice/thsice_slab_ocean.F

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

revision 1.2 by jmc, Thu Jun 3 16:34:36 2004 UTC revision 1.3 by jmc, Thu Jun 24 23:46:26 2004 UTC
# Line 7  CBOP Line 7  CBOP
7  C     !ROUTINE: THSICE_SLAB_OCEAN  C     !ROUTINE: THSICE_SLAB_OCEAN
8  C     !INTERFACE:  C     !INTERFACE:
9        SUBROUTINE THSICE_SLAB_OCEAN(        SUBROUTINE THSICE_SLAB_OCEAN(
10         O                      dTsurf,
11       I                      bi, bj, myThid )       I                      bi, bj, myThid )
12  C     !DESCRIPTION: \bv  C     !DESCRIPTION: \bv
13  C     *==========================================================*  C     *==========================================================*
# Line 39  C-- Physics package Line 40  C-- Physics package
40    
41  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
42  C     == Routine Arguments ==  C     == Routine Arguments ==
43          _RL dTsurf(sNx,sNy)
44        INTEGER bi,bj        INTEGER bi,bj
45        INTEGER myThid        INTEGER myThid
46  CEOP  CEOP
# Line 60  C--    add heat flux and fresh-water + s Line 62  C--    add heat flux and fresh-water + s
62         DO j=1,sNy         DO j=1,sNy
63          DO i=1,sNx          DO i=1,sNx
64           IF ( hOceMxL(i,j,bi,bj).NE.0. _d 0 ) THEN           IF ( hOceMxL(i,j,bi,bj).NE.0. _d 0 ) THEN
65              dTsurf(i,j) = tOceMxL(i,j,bi,bj)
66            tOceMxL(i,j,bi,bj) = tOceMxL(i,j,bi,bj)            tOceMxL(i,j,bi,bj) = tOceMxL(i,j,bi,bj)
67       &       - heatFac*Qnet(i,j,bi,bj) / hOceMxL(i,j,bi,bj)       &       - heatFac*Qnet(i,j,bi,bj) / hOceMxL(i,j,bi,bj)
68            sOceMxL(i,j,bi,bj) = sOceMxL(i,j,bi,bj)            sOceMxL(i,j,bi,bj) = sOceMxL(i,j,bi,bj)
# Line 93  C--    add restoring (backward) toward c Line 96  C--    add restoring (backward) toward c
96        ENDIF        ENDIF
97  #endif /* ALLOW_AIM */  #endif /* ALLOW_AIM */
98    
99    C-    Diagnose surf. temp. change
100           DO j=1,sNy
101            DO i=1,sNx
102             IF ( hOceMxL(i,j,bi,bj).NE.0. _d 0 ) THEN
103              dTsurf(i,j) = tOceMxL(i,j,bi,bj) - dTsurf(i,j)
104             ENDIF
105            ENDDO
106           ENDDO
107    
108  #endif  /* ALLOW_THSICE */  #endif  /* ALLOW_THSICE */
109    
110        RETURN        RETURN
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22