dgoldberg/streamice/streamice_upd_ffrac_uncoupled.F

C $Header: /u/gcmpack/MITgcm_contrib/dgoldberg/streamice/streamice_upd_ffrac_uncoupled.F,v 1.1 2012/03/29 15:59:21 heimbach Exp $
C $Name:  $


C this needs changes

#include "STREAMICE_OPTIONS.h"

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
CBOP 0
      SUBROUTINE STREAMICE_UPD_FFRAC_UNCOUPLED ( myThid )

C     !DESCRIPTION:
C     Initialize STREAMICE variables and constants.

C     !USES:
      IMPLICIT NONE
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "STREAMICE.h"
#include "GRID.h"

C     !INPUT PARAMETERS:
      INTEGER myThid
CEOP

#ifdef ALLOW_STREAMICE

      INTEGER bi, bj, i, j
      _RL OD, rhoi, rhow

      rhoi = streamice_density
      rhow = streamice_density_ocean_avg

      DO bj = myByLo(myThid), myByHi(myThid)
       DO bi = myBxLo(myThid), myBxHi(myThid)
        DO j=1-Oly,sNy+Oly
         DO i=1-Olx,sNx+Olx
          OD = -1.0 * R_low (i,j,bi,bj) - 
     &     H_streamice(i,j,bi,bj) * rhoi/rhow

!           if (i.eq.1) then
!             print *, R_low(1,1,1,1)
!           endif


          IF (OD .ge. 0. _d 0) THEN          

c         ice thickness does not take up whole ocean column -> floating
           float_frac_streamice(i,j,bi,bj) = 0.0
           base_el_streamice(i,j,bi,bj) = R_low(i,j,bi,bj)+OD
           surf_el_streamice(i,j,bi,bj) = 
     &      (1-rhoi/rhow)*H_streamice(i,j,bi,bj)
          ELSE


           float_frac_streamice(i,j,bi,bj) = 1.0
           base_el_streamice(i,j,bi,bj) = R_low(i,j,bi,bj)
           surf_el_streamice(i,j,bi,bj) = R_low(i,j,bi,bj) 
     &      + H_streamice(i,j,bi,bj)
          ENDIF
         ENDDO
        ENDDO
       ENDDO
      ENDDO

      _EXCH_XY_RL(float_frac_streamice, myThid )
      _EXCH_XY_RL(base_el_streamice, myThid )
      _EXCH_XY_RL(surf_el_streamice, myThid )

      
#endif /* ALLOW_STREAMICE */

      RETURN
      END
1	dgoldberg	1.2	C $Header: /u/gcmpack/MITgcm_contrib/dgoldberg/streamice/streamice_upd_ffrac_uncoupled.F,v 1.1 2012/03/29 15:59:21 heimbach Exp $
2	heimbach	1.1	C $Name: $
3
4
5			C this needs changes
6
7			#include "STREAMICE_OPTIONS.h"
8
9			C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
10			CBOP 0
11			SUBROUTINE STREAMICE_UPD_FFRAC_UNCOUPLED ( myThid )
12
13			C !DESCRIPTION:
14			C Initialize STREAMICE variables and constants.
15
16			C !USES:
17			IMPLICIT NONE
18			#include "SIZE.h"
19			#include "EEPARAMS.h"
20			#include "PARAMS.h"
21			#include "STREAMICE.h"
22			#include "GRID.h"
23
24			C !INPUT PARAMETERS:
25			INTEGER myThid
26			CEOP
27
28			#ifdef ALLOW_STREAMICE
29
30			INTEGER bi, bj, i, j
31			_RL OD, rhoi, rhow
32
33			rhoi = streamice_density
34			rhow = streamice_density_ocean_avg
35
36			DO bj = myByLo(myThid), myByHi(myThid)
37			DO bi = myBxLo(myThid), myBxHi(myThid)
38			DO j=1-Oly,sNy+Oly
39			DO i=1-Olx,sNx+Olx
40			OD = -1.0 * R_low (i,j,bi,bj) -
41			& H_streamice(i,j,bi,bj) * rhoi/rhow
42
43			! if (i.eq.1) then
44			! print *, R_low(1,1,1,1)
45			! endif
46
47
48			IF (OD .ge. 0. _d 0) THEN
49
50			c ice thickness does not take up whole ocean column -> floating
51			float_frac_streamice(i,j,bi,bj) = 0.0
52			base_el_streamice(i,j,bi,bj) = R_low(i,j,bi,bj)+OD
53			surf_el_streamice(i,j,bi,bj) =
54			& (1-rhoi/rhow)*H_streamice(i,j,bi,bj)
55			ELSE
56
57
58			float_frac_streamice(i,j,bi,bj) = 1.0
59			base_el_streamice(i,j,bi,bj) = R_low(i,j,bi,bj)
60			surf_el_streamice(i,j,bi,bj) = R_low(i,j,bi,bj)
61			& + H_streamice(i,j,bi,bj)
62			ENDIF
63			ENDDO
64			ENDDO
65			ENDDO
66			ENDDO
67
68			_EXCH_XY_RL(float_frac_streamice, myThid )
69			_EXCH_XY_RL(base_el_streamice, myThid )
70			_EXCH_XY_RL(surf_el_streamice, myThid )
71
72
73
74
75			#endif /* ALLOW_STREAMICE */
76
77			RETURN
78	dgoldberg	1.2	END