--- MITgcm/pkg/thsice/thsice_impl_temp.F	2004/05/07 21:33:34	1.2
+++ MITgcm/pkg/thsice/thsice_impl_temp.F	2006/03/13 03:53:40	1.3
@@ -1,19 +1,19 @@
-C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/thsice/thsice_impl_temp.F,v 1.2 2004/05/07 21:33:34 jmc Exp $
+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 $Name:  $
 
 #include "THSICE_OPTIONS.h"
- 
+
 C     !ROUTINE: THSICE_IMPL_TEMP
 C     !INTERFACE:
-      SUBROUTINE THSICE_IMPL_TEMP( 
+      SUBROUTINE THSICE_IMPL_TEMP(
      I                netSW, sFlx,
      O                dTsurf,
      I                bi, bj, myTime, myIter, myThid)
 C     *==========================================================*
-C     | S/R  THSICE_IMPL_TEMP             
+C     | S/R  THSICE_IMPL_TEMP
 C     | o Calculate sea-ice and surface temp. implicitly
 C     *==========================================================*
-C     | o return surface fluxes for atmosphere boundary layer 
+C     | o return surface fluxes for atmosphere boundary layer
 C     |  physics (and therefore called within atmospheric physics)
 C     *==========================================================*
 
@@ -27,7 +27,7 @@
 #include "THSICE_SIZE.h"
 #include "THSICE_PARAMS.h"
 #include "THSICE_VARS.h"
- 
+
 C     !INPUT/OUTPUT PARAMETERS:
 C     === Routine arguments ===
 C     netSW   :: net Short Wave surf. flux (+=down) [W/m2]
@@ -72,17 +72,17 @@
       jMin = 1
       jMax = sNy
       dBug = .FALSE.
- 1010 FORMAT(A,1P4E11.3)
+ 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 )
 
 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
-C    part.1 : ice-covered fraction ; 
+C    part.1 : ice-covered fraction ;
 C     Solve for surface and ice temperature (implicitly) ; compute surf. fluxes
 C-------
-        IF (iceMask(i,j,bi,bj).GT.0. _d 0) THEN 
+        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)
@@ -97,6 +97,8 @@
            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
+           WRITE(6,1010) 'ThSI_IMPL_T:-0- Tice(1,2),Qice(1,2)=',
+     &           Tice1(i,j,bi,bj), Tice2(i,j,bi,bj), qicen
           ENDIF
 
           CALL THSICE_SOLVE4TEMP(
@@ -119,6 +121,8 @@
      &                              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