2 |
C $Name$ |
C $Name$ |
3 |
|
|
4 |
#include "SEAICE_OPTIONS.h" |
#include "SEAICE_OPTIONS.h" |
5 |
|
#ifdef ALLOW_AUTODIFF |
6 |
|
# include "AUTODIFF_OPTIONS.h" |
7 |
|
#endif |
8 |
|
|
9 |
C !ROUTINE: SEAICE_ITD_PICKUP |
C !ROUTINE: SEAICE_ITD_PICKUP |
10 |
|
|
88 |
_RL LND_dx |
_RL LND_dx |
89 |
_RL LND_tmp |
_RL LND_tmp |
90 |
C bin width of distribution |
C bin width of distribution |
91 |
PARAMETER(LND_dx=0.1) |
PARAMETER( LND_iend = 1000 ) |
92 |
PARAMETER(LND_iend=INT(100./LND_dx)) |
PARAMETER( LND_dx = 100.D0 / LND_iend ) |
93 |
_RL LND_x (INT(LND_iend)) |
c PARAMETER(LND_dx=0.1) |
94 |
_RL LND_pdf(INT(LND_iend)) |
c PARAMETER(LND_iend=INT(100./LND_dx)) |
95 |
|
_RL LND_x (LND_iend) |
96 |
|
_RL LND_pdf(LND_iend) |
97 |
C flag for pickup scheme |
C flag for pickup scheme |
98 |
LOGICAL simpleSchemeFlag |
LOGICAL simpleSchemeFlag |
99 |
|
|
101 |
dummyTime = 1.0 |
dummyTime = 1.0 |
102 |
|
|
103 |
C---+-|--1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
C---+-|--1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
104 |
|
C reset ITD variables to zero for safety |
105 |
|
DO k = 1, nITD |
106 |
|
DO bj=myByLo(myThid),myByHi(myThid) |
107 |
|
DO bi=myBxLo(myThid),myBxHi(myThid) |
108 |
|
DO j=1-OLy,sNy+OLy |
109 |
|
DO i=1-OLx,sNx+OLx |
110 |
|
AREAITD(i,j,k,bi,bj) = 0. _d 0 |
111 |
|
HEFFITD(i,j,k,bi,bj) = 0. _d 0 |
112 |
|
HSNOWITD(i,j,k,bi,bj) = 0. _d 0 |
113 |
|
ENDDO |
114 |
|
ENDDO |
115 |
|
ENDDO |
116 |
|
ENDDO |
117 |
|
ENDDO |
118 |
IF (simpleSchemeFlag) THEN |
IF (simpleSchemeFlag) THEN |
119 |
C-- Put all ice into one bin: |
C-- Put all ice into one bin: |
120 |
C |
C |
121 |
DO bj=myByLo(myThid),myByHi(myThid) |
DO bj=myByLo(myThid),myByHi(myThid) |
122 |
DO bi=myBxLo(myThid),myBxHi(myThid) |
DO bi=myBxLo(myThid),myBxHi(myThid) |
123 |
DO j=1-OLy,sNy+OLy |
DO j=1-OLy,sNy+OLy |
124 |
DO i=1-OLx,sNx+OLx |
DO i=1-OLx,sNx+OLx |
125 |
AREAITD(i,j,1,bi,bj) = AREA(i,j,bi,bj) |
AREAITD(i,j,1,bi,bj) = AREA(i,j,bi,bj) |
126 |
HEFFITD(i,j,1,bi,bj) = HEFF(i,j,bi,bj) |
HEFFITD(i,j,1,bi,bj) = HEFF(i,j,bi,bj) |
127 |
HSNOWITD(i,j,1,bi,bj) = HSNOW(i,j,bi,bj) |
HSNOWITD(i,j,1,bi,bj) = HSNOW(i,j,bi,bj) |
|
ENDDO |
|
|
ENDDO |
|
128 |
ENDDO |
ENDDO |
129 |
ENDDO |
ENDDO |
130 |
|
ENDDO |
131 |
|
ENDDO |
132 |
|
|
133 |
C---+-|--1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
C---+-|--1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
134 |
ELSE |
ELSE |