/[MITgcm]/MITgcm_contrib/torge/itd/code/seaice_itd_pickup.F
ViewVC logotype

Annotation of /MITgcm_contrib/torge/itd/code/seaice_itd_pickup.F

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

Revision 1.2 - (hide annotations) (download)
Wed Sep 26 17:30:06 2012 UTC (12 years, 10 months ago) by torge
Branch: MAIN
Changes since 1.1: +5 -1 lines 
change seaice_itd_sum and seaice_itd_redist to exclude bi,bj loops;
both routines must now be called within bi,bj loops;
this is required in seaice_growth
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 torge 1.2 DO bj=myByLo(myThid),myByHi(myThid)
113     DO bi=myBxLo(myThid),myBxHi(myThid)
114     CALL SEAICE_ITD_REDIST( bi, bj, dummyTime, myIter, myThid)
115     ENDDO
116     ENDDO
117 torge 1.1
118     C---+-|--1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
119     ELSE
120     C-- Assume log-normal ITD:
121    
122     DO bj=myByLo(myThid),myByHi(myThid)
123     DO bi=myBxLo(myThid),myBxHi(myThid)
124     DO j=1-OLy,sNy+OLy
125     DO i=1-OLx,sNx+OLx
126     C
127     C initialize log-normal distribution
128     LND_mu = log(HEFF(i,j,bi,bj)/AREA(i,j,bi,bj))
129     & - 0.5*LND_sigma*LND_sigma
130     LND_x(1) = 0.+LND_dx/2.
131     C make thickness bins
132     DO LND_i=2,LND_iend
133     LND_x(LND_i)=LND_x(LND_i-1)+LND_dx
134     ENDDO
135     C log-normal distribution:
136     DO LND_i=2,LND_iend
137     LND_tmp = log(LND_x(LND_i))-LND_mu
138     LND_pdf(LND_i)= 1.
139     & / (LND_x(LND_i)*LND_sigma*sqrt(2*3.1416))
140     & * exp( -(LND_tmp*LND_tmp)
141     & / (2*LND_sigma*LND_sigma) )
142     & * AREA(i,j,bi,bj)
143     ENDDO
144     C assign bins to ice thickness categories
145     k=1
146     DO LND_i=1,LND_iend
147     IF ( LND_x(LND_i).GT.Hlimit(k) ) k=k+1
148     AREAITD(i,j,k,bi,bj) = AREAITD(i,j,k,bi,bj)
149     & + LND_pdf(LND_i)*LND_dx
150     HEFFITD(i,j,k,bi,bj) = HEFFITD(i,j,k,bi,bj)
151     & + LND_pdf(LND_i)*LND_x(LND_i)*LND_dx
152     ENDDO
153     C
154     ENDDO
155     ENDDO
156     ENDDO
157     ENDDO
158    
159     ENDIF
160     C---+-|--1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
161     #endif /* SEAICE_ITD */
162     RETURN
163     END SUBROUTINE SEAICE_ITD_PICKUP
164    
  ViewVC Help
Powered by ViewVC 1.1.22