--- MITgcm/pkg/thsice/thsice_impl_temp.F 2006/03/13 03:53:40 1.3 +++ MITgcm/pkg/thsice/thsice_impl_temp.F 2006/05/25 18:03:24 1.4 @@ -1,4 +1,4 @@ -C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/thsice/thsice_impl_temp.F,v 1.3 2006/03/13 03:53:40 jmc Exp $ +C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/thsice/thsice_impl_temp.F,v 1.4 2006/05/25 18:03:24 jmc Exp $ C $Name: $ #include "THSICE_OPTIONS.h" @@ -27,6 +27,9 @@ #include "THSICE_SIZE.h" #include "THSICE_PARAMS.h" #include "THSICE_VARS.h" + INTEGER siLo, siHi, sjLo, sjHi + PARAMETER ( siLo = 1-OLx , siHi = sNx+OLx ) + PARAMETER ( sjLo = 1-OLy , sjHi = sNy+OLy ) C !INPUT/OUTPUT PARAMETERS: C === Routine arguments === @@ -49,89 +52,64 @@ #ifdef ALLOW_THSICE C !LOCAL VARIABLES: C === Local variables === -C flxSW :: net Short-Wave (+=down) at surface [W/m2] -C flxExcSw :: surf. heat flux (+=down) except SW, function of surf. temp Ts: -C 0: Flx(Ts=0) ; 1: Flx(Ts=Ts^n) ; 2: d.Flx/dTs(Ts=Ts^n) -C sHeating :: surf heating left to melt snow or ice (= Atmos-conduction) -C flxCnB :: heat flux conducted through the ice to bottom surface -C flxtmp :: net heat flux from the atmosphere ( >0 downward) -C evptmp :: evaporation to the atmosphere [kg/m2/s] (>0 if evaporate) +C tFrzOce :: sea-water freezing temperature [oC] (function of S) +C dTsrf :: surf. temp adjusment: Ts^n+1 - Ts^n INTEGER i,j INTEGER iMin, iMax INTEGER jMin, jMax - _RL flxSW - _RL Tf - _RL flxtmp, evptmp - _RL flxExcSw(0:2) - _RL hIce, hSnow, Tsf, Tice(nlyr), qicen(nlyr) + _RL tFrzOce(1-OLx:sNx+OLx,1-OLy:sNy+OLy) +c _RL dTsrf (1-OLx:sNx+OLx,1-OLy:sNy+OLy) + LOGICAL dBugFlag + +C- define grid-point location where to print debugging values +#include "THSICE_DEBUG.h" + + 1010 FORMAT(A,1P4E14.6) - LOGICAL dBug +C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| iMin = 1 iMax = sNx jMin = 1 jMax = sNy - dBug = .FALSE. - 1010 FORMAT(A,1P4E14.6) - - DO j = jMin, jMax - DO i = iMin, iMax -c dBug = ( bi.EQ.3 .AND. i.EQ.15 .AND. j.EQ.11 ) + dBugFlag = debugLevel.GE.debLevB C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| C part.1 : ice-covered fraction ; C Solve for surface and ice temperature (implicitly) ; compute surf. fluxes C------- + DO j = jMin, jMax + DO i = iMin, iMax IF (iceMask(i,j,bi,bj).GT.0. _d 0) THEN - Tf = -mu_Tf*sOceMxL(i,j,bi,bj) - hIce = iceHeight(i,j,bi,bj) - hSnow = snowHeight(i,j,bi,bj) - flxExcSw(0) = sFlx(i,j,0) - flxExcSw(1) = sFlx(i,j,1) - flxExcSw(2) = sFlx(i,j,2) - flxSW = netSW(i,j) - Tsf = Tsrf(i,j,bi,bj) - qicen(1)= Qice1(i,j,bi,bj) - qicen(2)= Qice2(i,j,bi,bj) - IF ( dBug ) THEN + tFrzOce(i,j) = -mu_Tf*sOceMxL(i,j,bi,bj) + icFlxSW(i,j,bi,bj) = netSW(i,j) +#ifdef ALLOW_DBUG_THSICE + IF ( dBug(i,j,bi,bj) ) THEN WRITE(6,'(A,2I4,2I2)') 'ThSI_IMPL_T: i,j=',i,j,bi,bj WRITE(6,1010) 'ThSI_IMPL_T:-0- iceMask,hIc,hSn,Tsf=', - & iceMask(i,j,bi,bj), hIce, hSnow, Tsf + & iceMask(i,j,bi,bj), iceHeight(i,j,bi,bj), + & snowHeight(i,j,bi,bj), Tsrf(i,j,bi,bj) WRITE(6,1010) 'ThSI_IMPL_T:-0- Tice(1,2),Qice(1,2)=', - & Tice1(i,j,bi,bj), Tice2(i,j,bi,bj), qicen + & Tice1(i,j,bi,bj), Tice2(i,j,bi,bj), + & Qice1(i,j,bi,bj), Qice2(i,j,bi,bj) ENDIF - - CALL THSICE_SOLVE4TEMP( - I useBulkforce, flxExcSw, Tf, hIce, hSnow, - U flxSW, Tsf, qicen, - O Tice, sHeating(i,j,bi,bj), flxCndBt(i,j,bi,bj), - O dTsurf(i,j), flxtmp, evptmp, - I i,j, bi,bj, myThid) -C-- Update Fluxes : - Qsw(i,j,bi,bj) = -flxSW -C-- Update Sea-Ice state : - Tsrf(i,j,bi,bj) =Tsf - Tice1(i,j,bi,bj)=Tice(1) - Tice2(i,j,bi,bj)=Tice(2) - Qice1(i,j,bi,bj)=qicen(1) - Qice2(i,j,bi,bj)=qicen(2) - - IF ( dBug ) THEN - WRITE(6,1010) 'ThSI_IMPL_T: Tsf, Tice(1,2), dTsurf=', - & Tsf, Tice, dTsurf(i,j) - WRITE(6,1010) 'ThSI_IMPL_T: sHeat, flxCndBt, Qice =', - & sHeating(i,j,bi,bj), flxCndBt(i,j,bi,bj), qicen - WRITE(6,1010) 'ThSI_IMPL_T: flxA, fxSW_bf,af=', - & flxtmp, netSW(i,j), flxSW - ENDIF - ELSE - dTsurf(i,j) = 0. _d 0 - Qsw(i,j,bi,bj) = 0. _d 0 +#endif ENDIF - ENDDO ENDDO + CALL THSICE_SOLVE4TEMP( + I bi, bj, siLo, siHi, sjLo, sjHi, + I iMin,iMax, jMin,jMax, dBugFlag, .FALSE., + I iceMask(siLo,sjLo,bi,bj), iceHeight(siLo,sjLo,bi,bj), + I snowHeight(siLo,sjLo,bi,bj), tFrzOce, sFlx, + U icFlxSW(siLo,sjLo,bi,bj), Tsrf(siLo,sjLo,bi,bj), + U Qice1(siLo,sjLo,bi,bj), Qice2(siLo,sjLo,bi,bj), + O Tice1(siLo,sjLo,bi,bj), Tice2(siLo,sjLo,bi,bj), dTsurf, + O sHeating(siLo,sjLo,bi,bj), flxCndBt(siLo,sjLo,bi,bj), + O icFlxAtm(siLo,sjLo,bi,bj), icFrwAtm(siLo,sjLo,bi,bj), + I myTime, myIter, myThid ) + C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| #endif /* ALLOW_THSICE */