itd/code/seaice_itd_pickup.F


#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_dx=0.1)
      PARAMETER(LND_iend=INT(100./LND_dx))
      _RL LND_x  (INT(LND_iend))
      _RL LND_pdf(INT(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         ... then sort into correct ice thickness category 
         CALL SEAICE_ITD_REDIST( dummyTime, myIter, myThid)

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-|--+----|
#endif /* SEAICE_ITD */
      RETURN
      END SUBROUTINE SEAICE_ITD_PICKUP

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