/[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.1 by heimbach, Mon Oct 22 21:14:07 2012 UTC revision 1.4 by gforget, Mon Oct 20 03:20:57 2014 UTC
# Line 1  Line 1 
1    C $Header$
2    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 21  C     |      into the correct category b Line 26  C     |      into the correct category b
26  C     |      -> simpleSchemeFlag = .TRUE.  C     |      -> simpleSchemeFlag = .TRUE.
27  C     |  C     |
28  C     |   b) a scheme that assumes a log-normal distribution based  C     |   b) a scheme that assumes a log-normal distribution based
29  C     |      on the mean ice thickness and a standard decviation  C     |      on the mean ice thickness and a standard decviation
30  C     |      of LND_sigma=0.25  C     |      of LND_sigma=0.25
31  C     |      -> simpleSchemeFlag = .FALSE.  C     |      -> simpleSchemeFlag = .FALSE.
32  C     |  C     |
# Line 75  C Line 80  C
80  #endif /* ALLOW_AUTODIFF_TAMC */  #endif /* ALLOW_AUTODIFF_TAMC */
81        _RL dummyTime        _RL dummyTime
82    
83  C     local variables for picking up ITD from single category pickup file  C     local variables for picking up ITD from single category pickup file
84        INTEGER LND_i, LND_iend        INTEGER LND_i, LND_iend
85  C     parameters for log-normal distribution (LND)  C     parameters for log-normal distribution (LND)
86        _RL LND_sigma, LND_mu        _RL LND_sigma, LND_mu
# Line 83  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 119  C--      Assume log-normal ITD: Line 126  C--      Assume log-normal ITD:
126              DO i=1-OLx,sNx+OLx              DO i=1-OLx,sNx+OLx
127  C  C
128  C            initialize log-normal distribution  C            initialize log-normal distribution
129               LND_mu = log(HEFF(i,j,bi,bj)/AREA(i,j,bi,bj))               LND_mu = log(HEFF(i,j,bi,bj)/AREA(i,j,bi,bj))
130       &                - 0.5*LND_sigma*LND_sigma       &                - 0.5*LND_sigma*LND_sigma
131               LND_x(1) = 0.+LND_dx/2.               LND_x(1) = 0.+LND_dx/2.
132  C            make thickness bins  C            make thickness bins
133               DO LND_i=2,LND_iend               DO LND_i=2,LND_iend
134                LND_x(LND_i)=LND_x(LND_i-1)+LND_dx                LND_x(LND_i)=LND_x(LND_i-1)+LND_dx
135               ENDDO               ENDDO
136  C            log-normal distribution:  C            log-normal distribution:
137               DO LND_i=2,LND_iend               DO LND_i=2,LND_iend
138                LND_tmp = log(LND_x(LND_i))-LND_mu                LND_tmp = log(LND_x(LND_i))-LND_mu
139                LND_pdf(LND_i)= 1.                LND_pdf(LND_i)= 1.
140       &             / (LND_x(LND_i)*LND_sigma*sqrt(2*3.1416))       &             / (LND_x(LND_i)*LND_sigma*sqrt(2*3.1416))
141       &             * exp( -(LND_tmp*LND_tmp)       &             * exp( -(LND_tmp*LND_tmp)
142       &             /       (2*LND_sigma*LND_sigma) )       &             /       (2*LND_sigma*LND_sigma) )
143       &             * AREA(i,j,bi,bj)       &             * AREA(i,j,bi,bj)
144               ENDDO               ENDDO
145  C            assign bins to ice thickness categories  C            assign bins to ice thickness categories
146               k=1               k=1
147               DO LND_i=1,LND_iend               DO LND_i=1,LND_iend
148                IF ( LND_x(LND_i).GT.Hlimit(k) ) k=k+1                IF ( LND_x(LND_i).GT.Hlimit(k) ) k=k+1
149                AREAITD(i,j,k,bi,bj) = AREAITD(i,j,k,bi,bj)                AREAITD(i,j,k,bi,bj) = AREAITD(i,j,k,bi,bj)
150       &                             + LND_pdf(LND_i)*LND_dx       &                             + LND_pdf(LND_i)*LND_dx
151                HEFFITD(i,j,k,bi,bj) = HEFFITD(i,j,k,bi,bj)                HEFFITD(i,j,k,bi,bj) = HEFFITD(i,j,k,bi,bj)
152       &                             + LND_pdf(LND_i)*LND_x(LND_i)*LND_dx       &                             + LND_pdf(LND_i)*LND_x(LND_i)*LND_dx
153               ENDDO               ENDDO
154  C  C
155              ENDDO              ENDDO
156             ENDDO             ENDDO
157            ENDDO            ENDDO
158           ENDDO           ENDDO
159    
160        ENDIF        ENDIF
161    
162  C---+-|--1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+-|--1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
163  C     finally sort into correct ice thickness category  C     finally sort into correct ice thickness category
164  C      and compute bulk variables  C      and compute bulk variables
165  C      (needed for dynamic solver at beginning of seaice_model.F)  C      (needed for dynamic solver at beginning of seaice_model.F)
166        DO bj=myByLo(myThid),myByHi(myThid)        DO bj=myByLo(myThid),myByHi(myThid)
167         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
# Line 166  C      (needed for dynamic solver at beg Line 173  C      (needed for dynamic solver at beg
173  C---+-|--1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+-|--1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
174  #endif /* SEAICE_ITD */  #endif /* SEAICE_ITD */
175        RETURN        RETURN
176        END SUBROUTINE SEAICE_ITD_PICKUP        END
   
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22