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

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

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

revision 1.1 by jmc, Wed Apr 7 23:40:34 2004 UTC revision 1.2 by jmc, Fri May 7 21:33:34 2004 UTC
# Line 76  C     evptmp   :: evaporation to the atm Line 76  C     evptmp   :: evaporation to the atm
76    
77        DO j = jMin, jMax        DO j = jMin, jMax
78         DO i = iMin, iMax         DO i = iMin, iMax
79  c       dBug = ( bi.EQ.3 .AND. i.EQ.13 .AND. j.EQ.13 )  c       dBug = ( bi.EQ.3 .AND. i.EQ.15 .AND. j.EQ.11 )
80    
81  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
82  C    part.1 : ice-covered fraction ;  C    part.1 : ice-covered fraction ;
# Line 93  C------- Line 93  C-------
93            Tsf     = Tsrf(i,j,bi,bj)            Tsf     = Tsrf(i,j,bi,bj)
94            qicen(1)= Qice1(i,j,bi,bj)            qicen(1)= Qice1(i,j,bi,bj)
95            qicen(2)= Qice2(i,j,bi,bj)            qicen(2)= Qice2(i,j,bi,bj)
96              IF ( dBug ) THEN
97               WRITE(6,'(A,2I4,2I2)') 'ThSI_IMPL_T: i,j=',i,j,bi,bj
98               WRITE(6,1010) 'ThSI_IMPL_T:-0- iceMask,hIc,hSn,Tsf=',
99         &                   iceMask(i,j,bi,bj), hIce, hSnow, Tsf
100              ENDIF
101    
102            CALL THSICE_SOLVE4TEMP(            CALL THSICE_SOLVE4TEMP(
103       I          useBulkforce, flxExcSw, Tf, hIce, hSnow,       I          useBulkforce, flxExcSw, Tf, hIce, hSnow,
104       U          flxSW, Tsf, qicen,       U          flxSW, Tsf, qicen,
# Line 107  C--    Update Sea-Ice state : Line 113  C--    Update Sea-Ice state :
113            Tice2(i,j,bi,bj)=Tice(2)            Tice2(i,j,bi,bj)=Tice(2)
114            Qice1(i,j,bi,bj)=qicen(1)            Qice1(i,j,bi,bj)=qicen(1)
115            Qice2(i,j,bi,bj)=qicen(2)            Qice2(i,j,bi,bj)=qicen(2)
116    
117              IF ( dBug ) THEN
118               WRITE(6,1010) 'ThSI_IMPL_T: Tsf, Tice(1,2), dTsurf=',
119         &                              Tsf, Tice, dTsurf(i,j)
120               WRITE(6,1010) 'ThSI_IMPL_T: sHeat, flxCndBt, Qice =',
121         &           sHeating(i,j,bi,bj), flxCndBt(i,j,bi,bj), qicen
122              ENDIF
123          ELSE          ELSE
124            dTsurf(i,j) = 0. _d 0            dTsurf(i,j) = 0. _d 0
125            Qsw(i,j,bi,bj) = 0. _d 0            Qsw(i,j,bi,bj) = 0. _d 0
126          ENDIF          ENDIF
127    
128         ENDDO         ENDDO
129        ENDDO        ENDDO
130    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22