/[MITgcm]/MITgcm/pkg/seaice/seaice_itd_pickup.F
ViewVC logotype

Diff of /MITgcm/pkg/seaice/seaice_itd_pickup.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.2 by jmc, Tue Oct 23 13:20:49 2012 UTC revision 1.5 by mlosch, Fri Apr 29 14:28:38 2016 UTC
# Line 2  C $Header$ Line 2  C $Header$
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    
# Line 85  C     parameters for log-normal distribu Line 88  C     parameters for log-normal distribu
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    
# Line 96  C     flag for pickup scheme Line 101  C     flag for pickup scheme
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

Legend:
Removed from v.1.2  
changed lines
  Added in v.1.5

  ViewVC Help
Powered by ViewVC 1.1.22