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

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

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

revision 1.13 by jmc, Mon Mar 13 03:53:40 2006 UTC revision 1.14 by jmc, Tue Mar 14 15:58:27 2006 UTC
# Line 73  C     frw2oc    :: fresh-water flux from Line 73  C     frw2oc    :: fresh-water flux from
73  C     fsalt     :: mass salt flux to the ocean  C     fsalt     :: mass salt flux to the ocean
74  C     frzmltMxL :: ocean mixed-layer freezing/melting potential [W/m2]  C     frzmltMxL :: ocean mixed-layer freezing/melting potential [W/m2]
75  C     TFrzOce   :: sea-water freezing temperature [oC] (function of S)  C     TFrzOce   :: sea-water freezing temperature [oC] (function of S)
76    C     isIceFree :: true for ice-free grid-cell that remains ice-free
77        INTEGER i,j        INTEGER i,j
78        _RL snowPr        _RL snowPr
79        _RL agingTime, ageFac        _RL agingTime, ageFac
# Line 88  C     TFrzOce   :: sea-water freezing te Line 89  C     TFrzOce   :: sea-water freezing te
89        _RL oceV2s, oceTs        _RL oceV2s, oceTs
90        _RL compact, hIce, hSnow, Tsf, Tice(nlyr), qicen(nlyr)        _RL compact, hIce, hSnow, Tsf, Tice(nlyr), qicen(nlyr)
91        _RL tmpflx(0:2), tmpdTs        _RL tmpflx(0:2), tmpdTs
92          LOGICAL isIceFree(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
93  #ifdef ALLOW_DIAGNOSTICS  #ifdef ALLOW_DIAGNOSTICS
94        _RL tmpFac        _RL tmpFac
95  #endif  #endif
# Line 100  C-    Initialise flxAtm Line 102  C-    Initialise flxAtm
102         DO j = 1-Oly, sNy+Oly         DO j = 1-Oly, sNy+Oly
103          DO i = 1-Olx, sNx+Olx          DO i = 1-Olx, sNx+Olx
104            flxAtm(i,j) = 0.            flxAtm(i,j) = 0.
105              isIceFree(i,j) = .FALSE.
106          ENDDO          ENDDO
107         ENDDO         ENDDO
108    
# Line 134  C------- Line 137  C-------
137            siceAlb(i,j,bi,bj) = albedo            siceAlb(i,j,bi,bj) = albedo
138    
139            CALL THSICE_SOLVE4TEMP(            CALL THSICE_SOLVE4TEMP(
140       I          useBulkforce, tmpflx, TFrzOce, hIce, hSnow,       I          useBulkForce, tmpflx, TFrzOce, hIce, hSnow,
141       U          flxSW(i,j), Tsf, qicen,       U          flxSW(i,j), Tsf, qicen,
142       O          Tice, sHeating(i,j,bi,bj), flxCndBt(i,j,bi,bj),       O          Tice, sHeating(i,j,bi,bj), flxCndBt(i,j,bi,bj),
143       O          tmpdTs, flxAtm(i,j), evpAtm(i,j),       O          tmpdTs, flxAtm(i,j), evpAtm(i,j),
# Line 342  C--    Net fluxes : Line 345  C--    Net fluxes :
345  C--   - if esurp > 0 : end  C--   - if esurp > 0 : end
346          ENDIF          ENDIF
347    
348            IF ( hOceMxL(i,j,bi,bj).GT.0. _d 0 )
349         &    isIceFree(i,j) = iceMask(i,j,bi,bj).LE.0. _d 0
350         &                       .AND.   compact .LE.0. _d 0
351          IF ( compact .GT. 0. _d 0 ) THEN          IF ( compact .GT. 0. _d 0 ) THEN
352            iceMask(i,j,bi,bj)=compact            iceMask(i,j,bi,bj)=compact
353            IF ( hSnow .EQ. 0. _d 0 ) snowAge(i,j,bi,bj) = 0. _d 0            IF ( hSnow .EQ. 0. _d 0 ) snowAge(i,j,bi,bj) = 0. _d 0
# Line 371  C--     Compute Sea-Ice Loading (= mass Line 377  C--     Compute Sea-Ice Loading (= mass
377         ENDDO         ENDDO
378        ENDDO        ENDDO
379    
380    #ifdef ALLOW_BULK_FORCE
381          IF ( useBulkForce ) THEN
382            CALL BULKF_FLUX_ADJUST(
383         I                          bi, bj, iMin, iMax, jMin, jMax,
384         I                          isIceFree, myTime, myIter, myThid )
385          ENDIF
386    #endif /* ALLOW_BULK_FORCE */
387    
388  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
389  #endif /* ALLOW_THSICE */  #endif /* ALLOW_THSICE */
390    
Legend:
Removed from v.1.13  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.22