/[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.2 by jmc, Fri May 7 21:33:34 2004 UTC revision 1.3 by jmc, Mon Mar 13 03:53:40 2006 UTC
# Line 2  C $Header$ Line 2  C $Header$
2  C $Name$  C $Name$
3    
4  #include "THSICE_OPTIONS.h"  #include "THSICE_OPTIONS.h"
5    
6  C     !ROUTINE: THSICE_IMPL_TEMP  C     !ROUTINE: THSICE_IMPL_TEMP
7  C     !INTERFACE:  C     !INTERFACE:
8        SUBROUTINE THSICE_IMPL_TEMP(        SUBROUTINE THSICE_IMPL_TEMP(
9       I                netSW, sFlx,       I                netSW, sFlx,
10       O                dTsurf,       O                dTsurf,
11       I                bi, bj, myTime, myIter, myThid)       I                bi, bj, myTime, myIter, myThid)
12  C     *==========================================================*  C     *==========================================================*
13  C     | S/R  THSICE_IMPL_TEMP              C     | S/R  THSICE_IMPL_TEMP
14  C     | o Calculate sea-ice and surface temp. implicitly  C     | o Calculate sea-ice and surface temp. implicitly
15  C     *==========================================================*  C     *==========================================================*
16  C     | o return surface fluxes for atmosphere boundary layer  C     | o return surface fluxes for atmosphere boundary layer
17  C     |  physics (and therefore called within atmospheric physics)  C     |  physics (and therefore called within atmospheric physics)
18  C     *==========================================================*  C     *==========================================================*
19    
# Line 27  C     === Global variables === Line 27  C     === Global variables ===
27  #include "THSICE_SIZE.h"  #include "THSICE_SIZE.h"
28  #include "THSICE_PARAMS.h"  #include "THSICE_PARAMS.h"
29  #include "THSICE_VARS.h"  #include "THSICE_VARS.h"
30    
31  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
32  C     === Routine arguments ===  C     === Routine arguments ===
33  C     netSW   :: net Short Wave surf. flux (+=down) [W/m2]  C     netSW   :: net Short Wave surf. flux (+=down) [W/m2]
# Line 72  C     evptmp   :: evaporation to the atm Line 72  C     evptmp   :: evaporation to the atm
72        jMin = 1        jMin = 1
73        jMax = sNy        jMax = sNy
74        dBug = .FALSE.        dBug = .FALSE.
75   1010 FORMAT(A,1P4E11.3)   1010 FORMAT(A,1P4E14.6)
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.15 .AND. j.EQ.11 )  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 ;
83  C     Solve for surface and ice temperature (implicitly) ; compute surf. fluxes  C     Solve for surface and ice temperature (implicitly) ; compute surf. fluxes
84  C-------  C-------
85          IF (iceMask(i,j,bi,bj).GT.0. _d 0) THEN          IF (iceMask(i,j,bi,bj).GT.0. _d 0) THEN
86            Tf      = -mu_Tf*sOceMxL(i,j,bi,bj)            Tf      = -mu_Tf*sOceMxL(i,j,bi,bj)
87            hIce    = iceHeight(i,j,bi,bj)            hIce    = iceHeight(i,j,bi,bj)
88            hSnow   = snowHeight(i,j,bi,bj)            hSnow   = snowHeight(i,j,bi,bj)
# Line 97  C------- Line 97  C-------
97             WRITE(6,'(A,2I4,2I2)') 'ThSI_IMPL_T: i,j=',i,j,bi,bj             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=',             WRITE(6,1010) 'ThSI_IMPL_T:-0- iceMask,hIc,hSn,Tsf=',
99       &                   iceMask(i,j,bi,bj), hIce, hSnow, Tsf       &                   iceMask(i,j,bi,bj), hIce, hSnow, Tsf
100               WRITE(6,1010) 'ThSI_IMPL_T:-0- Tice(1,2),Qice(1,2)=',
101         &           Tice1(i,j,bi,bj), Tice2(i,j,bi,bj), qicen
102            ENDIF            ENDIF
103    
104            CALL THSICE_SOLVE4TEMP(            CALL THSICE_SOLVE4TEMP(
# Line 119  C--    Update Sea-Ice state : Line 121  C--    Update Sea-Ice state :
121       &                              Tsf, Tice, dTsurf(i,j)       &                              Tsf, Tice, dTsurf(i,j)
122             WRITE(6,1010) 'ThSI_IMPL_T: sHeat, flxCndBt, Qice =',             WRITE(6,1010) 'ThSI_IMPL_T: sHeat, flxCndBt, Qice =',
123       &           sHeating(i,j,bi,bj), flxCndBt(i,j,bi,bj), qicen       &           sHeating(i,j,bi,bj), flxCndBt(i,j,bi,bj), qicen
124               WRITE(6,1010) 'ThSI_IMPL_T: flxA, fxSW_bf,af=',
125         &           flxtmp, netSW(i,j), flxSW
126            ENDIF            ENDIF
127          ELSE          ELSE
128            dTsurf(i,j) = 0. _d 0            dTsurf(i,j) = 0. _d 0
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22