pkg/seaice/seaice_itd_pickup.F

C $Header: /u/gcmpack/MITgcm/pkg/seaice/seaice_itd_pickup.F,v 1.2 2012/10/23 13:20:49 jmc Exp $
C $Name:  $

#include "SEAICE_OPTIONS.h"

C !ROUTINE: SEAICE_ITD_PICKUP

C !INTERFACE: ==========================================================
      SUBROUTINE SEAICE_ITD_PICKUP(
     I     myIter, myThid )

C !DESCRIPTION: \bv
C     *===========================================================*
C     | SUBROUTINE SEAICE_ITD_PICKUP
C     | o called in case pickup file does not contain
C     |   ITD variables but mean ice thickness and concentration
C     |
C     | o choose between two schemes:
C     |
C     |   a) a simple scheme where the mean values are just put
C     |      into the first ITD category and then redustributed
C     |      into the correct category by SEAICE_ITD_REDIST
C     |      -> simpleSchemeFlag = .TRUE.
C     |
C     |   b) a scheme that assumes a log-normal distribution based
C     |      on the mean ice thickness and a standard decviation
C     |      of LND_sigma=0.25
C     |      -> simpleSchemeFlag = .FALSE.
C     |
C     | Torge Martin, Mai 2012, torge@mit.edu
C     *===========================================================*
C \ev

C !USES: ===============================================================
      IMPLICIT NONE

C     === Global variables needed ===
C     AREA      :: total sea ice area fraction
C     HEFF      :: mean in-situ sea ice thickness
C     HSNOW     :: mean in-situ snow layer depth
C
C     === Global variables to be changed ===
C     AREAITD   :: sea ice area      by category
C     HEFFITD   :: sea ice thickness by category
C     HSNOWITD  :: snow thickness    by category
C
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "GRID.h"
#include "SEAICE_SIZE.h"
#include "SEAICE_PARAMS.h"
#include "SEAICE.h"

#ifdef ALLOW_AUTODIFF_TAMC
# include "tamc.h"
#endif

C !INPUT PARAMETERS: ===================================================
C     === Routine arguments ===
C     myIter    :: iteration number
C     myThid    :: Thread no. that called this routine.
      INTEGER myIter
      INTEGER myThid
CEndOfInterface

#ifdef SEAICE_ITD

C !LOCAL VARIABLES: ====================================================
C     === Local variables ===
C     i,j,bi,bj,k :: Loop counters
C     nITD        :: number of sea ice thickness categories
C
      INTEGER i, j, bi, bj, k
#ifdef ALLOW_AUTODIFF_TAMC
      INTEGER itmpkey
#endif /* ALLOW_AUTODIFF_TAMC */
      _RL dummyTime

C     local variables for picking up ITD from single category pickup file
      INTEGER LND_i, LND_iend
C     parameters for log-normal distribution (LND)
      _RL LND_sigma, LND_mu
      PARAMETER(LND_sigma=0.25)
      _RL LND_dx
      _RL LND_tmp
C     bin width of distribution
      PARAMETER( LND_iend = 1000 )
      PARAMETER( LND_dx = 100.D0 / LND_iend )
c     PARAMETER(LND_dx=0.1)
c     PARAMETER(LND_iend=INT(100./LND_dx))
      _RL LND_x  (LND_iend)
      _RL LND_pdf(LND_iend)
C     flag for pickup scheme
      LOGICAL simpleSchemeFlag

      simpleSchemeFlag = .TRUE.
      dummyTime = 1.0

C---+-|--1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
      IF (simpleSchemeFlag) THEN
C--      Put all ice into one bin:
C
         DO bj=myByLo(myThid),myByHi(myThid)
          DO bi=myBxLo(myThid),myBxHi(myThid)
           DO j=1-OLy,sNy+OLy
            DO i=1-OLx,sNx+OLx
             AREAITD(i,j,1,bi,bj)  = AREA(i,j,bi,bj)
             HEFFITD(i,j,1,bi,bj)  = HEFF(i,j,bi,bj)
             HSNOWITD(i,j,1,bi,bj) = HSNOW(i,j,bi,bj)
            ENDDO
           ENDDO
          ENDDO
         ENDDO

C---+-|--1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
      ELSE
C--      Assume log-normal ITD:

         DO bj=myByLo(myThid),myByHi(myThid)
          DO bi=myBxLo(myThid),myBxHi(myThid)
           DO j=1-OLy,sNy+OLy
            DO i=1-OLx,sNx+OLx
C
C            initialize log-normal distribution
             LND_mu = log(HEFF(i,j,bi,bj)/AREA(i,j,bi,bj))
     &                - 0.5*LND_sigma*LND_sigma
             LND_x(1) = 0.+LND_dx/2.
C            make thickness bins
             DO LND_i=2,LND_iend
              LND_x(LND_i)=LND_x(LND_i-1)+LND_dx
             ENDDO
C            log-normal distribution:
             DO LND_i=2,LND_iend
              LND_tmp = log(LND_x(LND_i))-LND_mu
              LND_pdf(LND_i)= 1.
     &             / (LND_x(LND_i)*LND_sigma*sqrt(2*3.1416))
     &             * exp( -(LND_tmp*LND_tmp)
     &             /       (2*LND_sigma*LND_sigma) )
     &             * AREA(i,j,bi,bj)
             ENDDO
C            assign bins to ice thickness categories
             k=1
             DO LND_i=1,LND_iend
              IF ( LND_x(LND_i).GT.Hlimit(k) ) k=k+1
              AREAITD(i,j,k,bi,bj) = AREAITD(i,j,k,bi,bj)
     &                             + LND_pdf(LND_i)*LND_dx
              HEFFITD(i,j,k,bi,bj) = HEFFITD(i,j,k,bi,bj)
     &                             + LND_pdf(LND_i)*LND_x(LND_i)*LND_dx
             ENDDO
C
            ENDDO
           ENDDO
          ENDDO
         ENDDO

      ENDIF

C---+-|--1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
C     finally sort into correct ice thickness category
C      and compute bulk variables
C      (needed for dynamic solver at beginning of seaice_model.F)
      DO bj=myByLo(myThid),myByHi(myThid)
       DO bi=myBxLo(myThid),myBxHi(myThid)
        CALL SEAICE_ITD_REDIST( bi, bj, dummyTime, myIter, myThid)
        CALL SEAICE_ITD_SUM( bi, bj, dummyTime, myIter, myThid)
       ENDDO
      ENDDO

C---+-|--1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
#endif /* SEAICE_ITD */
      RETURN
      END
1	jmc	1.3	C $Header: /u/gcmpack/MITgcm/pkg/seaice/seaice_itd_pickup.F,v 1.2 2012/10/23 13:20:49 jmc Exp $
2	jmc	1.2	C $Name: $
3	heimbach	1.1
4			#include "SEAICE_OPTIONS.h"
5
6			C !ROUTINE: SEAICE_ITD_PICKUP
7
8			C !INTERFACE: ==========================================================
9			SUBROUTINE SEAICE_ITD_PICKUP(
10			I myIter, myThid )
11
12			C !DESCRIPTION: \bv
13			C ===========================================================
14			C \| SUBROUTINE SEAICE_ITD_PICKUP
15			C \| o called in case pickup file does not contain
16			C \| ITD variables but mean ice thickness and concentration
17			C \|
18			C \| o choose between two schemes:
19			C \|
20			C \| a) a simple scheme where the mean values are just put
21			C \| into the first ITD category and then redustributed
22			C \| into the correct category by SEAICE_ITD_REDIST
23			C \| -> simpleSchemeFlag = .TRUE.
24			C \|
25			C \| b) a scheme that assumes a log-normal distribution based
26	jmc	1.2	C \| on the mean ice thickness and a standard decviation
27	heimbach	1.1	C \| of LND_sigma=0.25
28			C \| -> simpleSchemeFlag = .FALSE.
29			C \|
30			C \| Torge Martin, Mai 2012, torge@mit.edu
31			C ===========================================================
32			C \ev
33
34			C !USES: ===============================================================
35			IMPLICIT NONE
36
37			C === Global variables needed ===
38			C AREA :: total sea ice area fraction
39			C HEFF :: mean in-situ sea ice thickness
40			C HSNOW :: mean in-situ snow layer depth
41			C
42			C === Global variables to be changed ===
43			C AREAITD :: sea ice area by category
44			C HEFFITD :: sea ice thickness by category
45			C HSNOWITD :: snow thickness by category
46			C
47			#include "SIZE.h"
48			#include "EEPARAMS.h"
49			#include "PARAMS.h"
50			#include "GRID.h"
51			#include "SEAICE_SIZE.h"
52			#include "SEAICE_PARAMS.h"
53			#include "SEAICE.h"
54
55			#ifdef ALLOW_AUTODIFF_TAMC
56			# include "tamc.h"
57			#endif
58
59			C !INPUT PARAMETERS: ===================================================
60			C === Routine arguments ===
61			C myIter :: iteration number
62			C myThid :: Thread no. that called this routine.
63			INTEGER myIter
64			INTEGER myThid
65			CEndOfInterface
66
67			#ifdef SEAICE_ITD
68
69			C !LOCAL VARIABLES: ====================================================
70			C === Local variables ===
71			C i,j,bi,bj,k :: Loop counters
72			C nITD :: number of sea ice thickness categories
73			C
74			INTEGER i, j, bi, bj, k
75			#ifdef ALLOW_AUTODIFF_TAMC
76			INTEGER itmpkey
77			#endif /* ALLOW_AUTODIFF_TAMC */
78			_RL dummyTime
79
80	jmc	1.2	C local variables for picking up ITD from single category pickup file
81	heimbach	1.1	INTEGER LND_i, LND_iend
82			C parameters for log-normal distribution (LND)
83			_RL LND_sigma, LND_mu
84			PARAMETER(LND_sigma=0.25)
85			_RL LND_dx
86			_RL LND_tmp
87			C bin width of distribution
88	jmc	1.3	PARAMETER( LND_iend = 1000 )
89			PARAMETER( LND_dx = 100.D0 / LND_iend )
90			c PARAMETER(LND_dx=0.1)
91			c PARAMETER(LND_iend=INT(100./LND_dx))
92			_RL LND_x (LND_iend)
93			_RL LND_pdf(LND_iend)
94	heimbach	1.1	C flag for pickup scheme
95			LOGICAL simpleSchemeFlag
96
97			simpleSchemeFlag = .TRUE.
98			dummyTime = 1.0
99
100			C---+-\|--1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
101			IF (simpleSchemeFlag) THEN
102			C-- Put all ice into one bin:
103			C
104			DO bj=myByLo(myThid),myByHi(myThid)
105			DO bi=myBxLo(myThid),myBxHi(myThid)
106			DO j=1-OLy,sNy+OLy
107			DO i=1-OLx,sNx+OLx
108			AREAITD(i,j,1,bi,bj) = AREA(i,j,bi,bj)
109			HEFFITD(i,j,1,bi,bj) = HEFF(i,j,bi,bj)
110			HSNOWITD(i,j,1,bi,bj) = HSNOW(i,j,bi,bj)
111			ENDDO
112			ENDDO
113			ENDDO
114			ENDDO
115
116			C---+-\|--1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
117			ELSE
118			C-- Assume log-normal ITD:
119
120			DO bj=myByLo(myThid),myByHi(myThid)
121			DO bi=myBxLo(myThid),myBxHi(myThid)
122			DO j=1-OLy,sNy+OLy
123			DO i=1-OLx,sNx+OLx
124			C
125			C initialize log-normal distribution
126	jmc	1.2	LND_mu = log(HEFF(i,j,bi,bj)/AREA(i,j,bi,bj))
127	heimbach	1.1	& - 0.5LND_sigmaLND_sigma
128	jmc	1.2	LND_x(1) = 0.+LND_dx/2.
129	heimbach	1.1	C make thickness bins
130	jmc	1.2	DO LND_i=2,LND_iend
131			LND_x(LND_i)=LND_x(LND_i-1)+LND_dx
132			ENDDO
133	heimbach	1.1	C log-normal distribution:
134	jmc	1.2	DO LND_i=2,LND_iend
135			LND_tmp = log(LND_x(LND_i))-LND_mu
136			LND_pdf(LND_i)= 1.
137			& / (LND_x(LND_i)LND_sigmasqrt(2*3.1416))
138	heimbach	1.1	& * exp( -(LND_tmp*LND_tmp)
139			& / (2LND_sigmaLND_sigma) )
140			& * AREA(i,j,bi,bj)
141	jmc	1.2	ENDDO
142	heimbach	1.1	C assign bins to ice thickness categories
143			k=1
144	jmc	1.2	DO LND_i=1,LND_iend
145			IF ( LND_x(LND_i).GT.Hlimit(k) ) k=k+1
146	heimbach	1.1	AREAITD(i,j,k,bi,bj) = AREAITD(i,j,k,bi,bj)
147			& + LND_pdf(LND_i)*LND_dx
148			HEFFITD(i,j,k,bi,bj) = HEFFITD(i,j,k,bi,bj)
149			& + LND_pdf(LND_i)LND_x(LND_i)LND_dx
150	jmc	1.2	ENDDO
151	heimbach	1.1	C
152			ENDDO
153	jmc	1.2	ENDDO
154	heimbach	1.1	ENDDO
155	jmc	1.2	ENDDO
156	heimbach	1.1
157			ENDIF
158
159			C---+-\|--1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
160	jmc	1.2	C finally sort into correct ice thickness category
161			C and compute bulk variables
162	heimbach	1.1	C (needed for dynamic solver at beginning of seaice_model.F)
163			DO bj=myByLo(myThid),myByHi(myThid)
164			DO bi=myBxLo(myThid),myBxHi(myThid)
165			CALL SEAICE_ITD_REDIST( bi, bj, dummyTime, myIter, myThid)
166			CALL SEAICE_ITD_SUM( bi, bj, dummyTime, myIter, myThid)
167			ENDDO
168			ENDDO
169
170			C---+-\|--1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
171			#endif /* SEAICE_ITD */
172			RETURN
173	jmc	1.2	END