1 | |
---|
2 | CCC $Header$ |
---|
3 | CCC TOP 1.0 , LOCEAN-IPSL (2005) |
---|
4 | C This software is governed by CeCILL licence see modipsl/doc/NEMO_CeCILL.txt |
---|
5 | C --------------------------------------------------------------------------- |
---|
6 | C $Id$ |
---|
7 | CCC--------------------------------------------------------------------- |
---|
8 | CCC |
---|
9 | CCC trcini.lobster1.h |
---|
10 | CCC ****************** |
---|
11 | CCC |
---|
12 | CCC purpose : |
---|
13 | CCC --------- |
---|
14 | CCC specific initialisation for lobster1 model |
---|
15 | CCC |
---|
16 | CCC modifications : |
---|
17 | CC ------------- |
---|
18 | CC original : 99-09 (M. Levy) |
---|
19 | CC additions : 00-12 (0. Aumont, E. Kestenare) |
---|
20 | CC add sediment computations |
---|
21 | CC |
---|
22 | CCC--------------------------------------------------------------------- |
---|
23 | CCC opa8, ipsl (11/96) |
---|
24 | CCC--------------------------------------------------------------------- |
---|
25 | CC local declarations |
---|
26 | CC ================== |
---|
27 | INTEGER ji,jj,jk |
---|
28 | REAL zdm0(jpi,jpj,jpk),zrro(jpi,jpj),zfluo,zfluu |
---|
29 | REAL ztest |
---|
30 | C |
---|
31 | C 1. initialization of fields for optical model |
---|
32 | C -------------------------------------------- |
---|
33 | C |
---|
34 | DO jj=1,jpj |
---|
35 | DO ji=1,jpi |
---|
36 | xze(ji,jj)=5. |
---|
37 | END DO |
---|
38 | END DO |
---|
39 | |
---|
40 | DO jk=1,jpk |
---|
41 | DO jj=1,jpj |
---|
42 | DO ji=1,jpi |
---|
43 | xpar(ji,jj,jk)=0. |
---|
44 | END DO |
---|
45 | END DO |
---|
46 | END DO |
---|
47 | C |
---|
48 | C 2. initialization for passive tracer remineralisation-damping array |
---|
49 | C ------------------------------------------------------------------------- |
---|
50 | C |
---|
51 | DO jn=1,jptra |
---|
52 | DO jk=1,jpk |
---|
53 | remdmp(jk,jn)=tminr |
---|
54 | END DO |
---|
55 | END DO |
---|
56 | C |
---|
57 | IF(lwp) THEN |
---|
58 | WRITE(numout,*) ' ' |
---|
59 | WRITE(numout,*) ' trcini: compute remineralisation-damping ' |
---|
60 | WRITE(numout,*) ' arrays for tracers' |
---|
61 | ENDIF |
---|
62 | C |
---|
63 | C 3. initialization of biological variables |
---|
64 | C ------------------------------------------ |
---|
65 | C |
---|
66 | C Calculate vertical distribution of newly formed biogenic poc |
---|
67 | C in the water column in the case of max. possible bottom depth |
---|
68 | C ------------------------------------------------------------ |
---|
69 | C |
---|
70 | zdm0 = 0. |
---|
71 | zrro = 1. |
---|
72 | DO jk = jpkb,jpkm1 |
---|
73 | C |
---|
74 | DO jj = 1,jpj |
---|
75 | DO ji = 1,jpi |
---|
76 | C |
---|
77 | zfluo = (fsdepw(ji,jj,jk)/fsdepw(ji,jj,jpkb))**xhr |
---|
78 | zfluu = (fsdepw(ji,jj,jk+1)/fsdepw(ji,jj,jpkb))**xhr |
---|
79 | IF (zfluo.gt.1.) zfluo = 1. |
---|
80 | zdm0(ji,jj,jk) = zfluo-zfluu |
---|
81 | IF (jk.le.jpkb-1) zdm0(ji,jj,jk)=0. |
---|
82 | zrro(ji,jj) = zrro(ji,jj)-zdm0(ji,jj,jk) |
---|
83 | C |
---|
84 | ENDDO |
---|
85 | ENDDO |
---|
86 | C |
---|
87 | ENDDO |
---|
88 | C |
---|
89 | DO jj = 1,jpj |
---|
90 | DO ji = 1,jpi |
---|
91 | zdm0(ji,jj,jpk) = zrro(ji,jj) |
---|
92 | ENDDO |
---|
93 | ENDDO |
---|
94 | C |
---|
95 | C Calculate vertical distribution of newly formed biogenic poc |
---|
96 | C in the water column with realistic topography (first "dry" layer |
---|
97 | C contains total fraction, which has passed to the upper layers) |
---|
98 | C ---------------------------------------------------------------------- |
---|
99 | C |
---|
100 | dminl = 0. |
---|
101 | dmin3 = zdm0 |
---|
102 | C |
---|
103 | DO jk = 1,jpk |
---|
104 | DO jj = 1,jpj |
---|
105 | DO ji = 1,jpi |
---|
106 | |
---|
107 | IF(tmask(ji,jj,jk).eq.0.) THEN |
---|
108 | dminl(ji,jj) = dminl(ji,jj)+dmin3(ji,jj,jk) |
---|
109 | dmin3(ji,jj,jk) = 0.0 |
---|
110 | ENDIF |
---|
111 | C |
---|
112 | ENDDO |
---|
113 | ENDDO |
---|
114 | ENDDO |
---|
115 | C |
---|
116 | DO jj = 1,jpj |
---|
117 | DO ji = 1,jpi |
---|
118 | IF (tmask(ji,jj,1).eq.0.) dmin3(ji,jj,1) = 0. |
---|
119 | ENDDO |
---|
120 | ENDDO |
---|
121 | C |
---|
122 | C CALCUL DU MASK DE COTE |
---|
123 | C |
---|
124 | cmask=0. |
---|
125 | do ji=2,jpi-1 |
---|
126 | do jj=2,jpj-1 |
---|
127 | if (tmask(ji,jj,1).eq.1) then |
---|
128 | ztest=tmask(ji+1,jj,1)*tmask(ji-1,jj,1)*tmask(ji,jj+1,1) |
---|
129 | . *tmask(ji,jj-1,1) |
---|
130 | if (ztest.eq.0) cmask(ji,jj)=1. |
---|
131 | endif |
---|
132 | end do |
---|
133 | end do |
---|
134 | |
---|
135 | cmask(1,:)=cmask(jpi-1,:) |
---|
136 | cmask(jpi,:)=cmask(2,:) |
---|
137 | C |
---|
138 | C CALCUL DE LA SURFACE COTIERE |
---|
139 | C |
---|
140 | areacot=0. |
---|
141 | do ji=2,jpi-1 |
---|
142 | do jj=2,jpj-1 |
---|
143 | areacot=areacot+e1t(ji,jj)*e2t(ji,jj)*cmask(ji,jj) |
---|
144 | end do |
---|
145 | end do |
---|
146 | |
---|