61 |
CHARACTER*(MAX_LEN_MBUF) msgBuf |
CHARACTER*(MAX_LEN_MBUF) msgBuf |
62 |
CHARACTER*10 HlimitMsgFormat |
CHARACTER*10 HlimitMsgFormat |
63 |
#endif |
#endif |
64 |
C ToM a random number to generate divergence and convergence randomly for the 1-D case |
#ifdef SEAICE_ITD |
|
_RL rand_num |
|
|
_RL divergence |
|
65 |
INTEGER IT |
INTEGER IT |
66 |
|
CToM a constant divergence rate for an otherwise purely thermodynamic 1-D test case (not a permanent addition) |
67 |
|
_RL divergence |
68 |
|
#endif |
69 |
C#if defined(SEAICE_GROWTH_LEGACY) || defined(ALLOW_AUTODIFF_TAMC) |
C#if defined(SEAICE_GROWTH_LEGACY) || defined(ALLOW_AUTODIFF_TAMC) |
70 |
#if defined(SEAICE_GROWTH_LEGACY) || defined(ALLOW_AUTODIFF_TAMC) || defined(SEAICE_ITD) |
#if defined(SEAICE_GROWTH_LEGACY) || defined(ALLOW_AUTODIFF_TAMC) || defined(SEAICE_ITD) |
71 |
C>>>ToM |
C>>>ToM |
317 |
divergence = 1.0/(5*365.*86400./SEAICE_deltaTtherm) |
divergence = 1.0/(5*365.*86400./SEAICE_deltaTtherm) |
318 |
DO bj=myByLo(myThid),myByHi(myThid) |
DO bj=myByLo(myThid),myByHi(myThid) |
319 |
DO bi=myBxLo(myThid),myBxHi(myThid) |
DO bi=myBxLo(myThid),myBxHi(myThid) |
320 |
DO IT=1,nITD |
DO IT=1,nITD |
321 |
DO J=1-OLy,sNy+OLy |
DO J=1-OLy,sNy+OLy |
322 |
DO I=1-OLx,sNx+OLx |
DO I=1-OLx,sNx+OLx |
323 |
AREAITD(I,J,IT,bi,bj)=MAX(0.0,AREAITD(I,J,IT,bi,bj) |
AREAITD(I,J,IT,bi,bj)=MAX(0.0,AREAITD(I,J,IT,bi,bj) |