Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Blame
Revision Log

sedmbc.F90 @ 2281

Last change on this file since 2281 was 2281, checked in by smasson, 14 years ago
set proper svn properties to all files...
Property svn:keywords set to `Id`
File size: 14.0 KB

Line
1	MODULE sedmbc
2	#if defined key_sed
3	!!======================================================================
4	!! * MODULE sedmbc *
5	!! Sediment : mass balance calculation
6	!!=====================================================================
7
8	!!----------------------------------------------------------------------
9	!! sed_mbc :
10	!!----------------------------------------------------------------------
11	!! * Modules used
12	USE sed ! sediment global variable
13	USE seddsr
14
15	IMPLICIT NONE
16	PRIVATE
17
18	!! * Routine accessibility
19	PUBLIC sed_mbc
20
21	!! * Module variables
22	REAL(wp), DIMENSION(jpsol) :: &
23	rain_tot , & !
24	fromsed_tot , & !
25	tosed_tot , & !
26	rloss_tot !
27
28	REAL(wp), DIMENSION(jpwat) :: &
29	diss_in_tot , & !
30	diss_out_tot !
31
32	REAL(wp) :: &
33	cons_tot_o2 , & !
34	sour_tot_no3 , & !
35	cons_tot_no3 , & !
36	sour_tot_c13 !
37
38	REAL(wp) :: &
39	src13p , &
40	src13ca
41
42	CONTAINS
43
44
45	SUBROUTINE sed_mbc( kt )
46	!!----------------------------------------------------------------------
47	!! * ROUTINE sed_mbc *
48	!!
49	!! ** Purpose : computation of total tracer inventories for checking
50	!! mass conservation.
51	!!
52	!!
53	!! ** Method : tracer inventories of each reservoir are computed and added
54	!! subsequently.
55	!!
56	!! History :
57	!! ! 04-10 (N. Emprin, M. Gehlen ) Original code
58	!! ! 06-07 (C. Ethe) Re-organization
59	!!----------------------------------------------------------------------
60
61	!! Arguments
62	INTEGER, INTENT(in) :: &
63	kt
64
65	!! local declarations
66	INTEGER :: &
67	ji, js, jw, jk
68
69	REAL(wp) :: &
70	zinit, zfinal, &
71	zinput, zoutput
72
73	REAL(wp) :: &
74	zdsw, zvol
75
76	REAL, DIMENSION(jpsol) :: &
77	zsolcp_inv_i, zsolcp_inv_f
78
79	REAL, DIMENSION(jpwat) :: &
80	zpwcp_inv_i, zpwcp_inv_f
81
82	REAL(wp) :: &
83	zdelta_sil, zdelta_clay, &
84	zdelta_co2, zdelta_oxy, &
85	zdelta_po4, zdelta_no3, &
86	zdelta_c13, zdelta_c13b
87
88	!!----------------------------------------------------------------------
89	! Initilization
90	!---------------
91	IF( kt == nitsed000 ) THEN
92	cons_tot_o2 = 0.
93	sour_tot_no3 = 0.
94	cons_tot_no3 = 0.
95	sour_tot_c13 = 0.
96
97	DO js = 1, jpsol
98	rain_tot (js) = 0.
99	fromsed_tot(js) = 0.
100	tosed_tot (js) = 0.
101	rloss_tot (js) = 0.
102	ENDDO
103
104	DO jw = 1, jpwat
105	diss_in_tot (jw) = 0.
106	diss_out_tot(jw) = 0.
107	ENDDO
108
109	src13p = rc13P * pdb
110	src13ca = rc13Ca * pdb
111	ENDIF
112
113
114	! Calculation of the cumulativ input and output
115	! for mass balance check
116	!----------------------------------------------
117
118	! cumulativ solid
119	DO js = 1, jpsol
120	DO ji = 1, jpoce
121	! input [mol]
122	rain_tot (js) = rain_tot (js) + dtsed * rainrm_dta(ji,js)
123	fromsed_tot(js) = fromsed_tot(js) + fromsed(ji,js)
124	! output [mol]
125	tosed_tot (js) = tosed_tot (js) + tosed(ji,js)
126	rloss_tot (js) = rloss_tot (js) + rloss(ji,js)
127	ENDDO
128	ENDDO
129
130	! cumulativ dissolved
131	DO jw = 1, jpwat
132	DO ji = 1, jpoce
133	! input [mol]
134	diss_in_tot (jw) = diss_in_tot (jw) + pwcp_dta(ji,jw) * 1.e-3 * dzkbot(ji)
135	! output [mol]
136	diss_out_tot(jw) = diss_out_tot(jw) + tokbot(ji,jw)
137	ENDDO
138	ENDDO
139
140	! cumulativ o2 and no3 consomation
141	DO ji = 1, jpoce
142	cons_tot_o2 = cons_tot_o2 + cons_o2 (ji)
143	sour_tot_no3 = sour_tot_no3 + sour_no3(ji)
144	cons_tot_no3 = cons_tot_no3 + cons_no3(ji)
145	sour_tot_c13 = sour_tot_c13 + sour_c13(ji)
146	ENDDO
147
148
149	! Mass balance check
150	!---------------------
151	IF( kt == nitsedend ) THEN
152	! initial and final inventories for solid component (mole/dx.dy) in sediment
153	zsolcp_inv_i(:) = 0.
154	zsolcp_inv_f(:) = 0.
155	zpwcp_inv_i (:) = 0.
156	zpwcp_inv_f (:) = 0.
157	DO js = 1, jpsol
158	zdsw = dens / mol_wgt(js)
159	DO jk = 2, jpksed
160	DO ji = 1, jpoce
161	zvol = vols3d(ji,jk) * zdsw
162	zsolcp_inv_i(js) = zsolcp_inv_i(js) + solcp0(ji,jk,js) * zvol
163	zsolcp_inv_f(js) = zsolcp_inv_f(js) + solcp (ji,jk,js) * zvol
164	ENDDO
165	END DO
166	ENDDO
167
168	! initial and final inventories for dissolved component (mole/dx.dy) in sediment
169	DO jw = 1, jpwat
170	DO jk = 2, jpksed
171	DO ji = 1, jpoce
172	zvol = volw3d(ji,jk) * 1.e-3
173	zpwcp_inv_i(jw) = zpwcp_inv_i(jw) + pwcp0(ji,jk,jw) * zvol
174	zpwcp_inv_f(jw) = zpwcp_inv_f(jw) + pwcp (ji,jk,jw) * zvol
175	ENDDO
176	END DO
177	ENDDO
178
179	! mass balance for Silica/opal
180	zinit = zsolcp_inv_i(jsopal) + zpwcp_inv_i(jwsil)
181	zfinal = zsolcp_inv_f(jsopal) + zpwcp_inv_f(jwsil)
182	zinput = rain_tot (jsopal) + diss_in_tot (jwsil)
183	zoutput = tosed_tot (jsopal) + rloss_tot (jsopal) + diss_out_tot(jwsil)
184	zdelta_sil = ( zfinal + zoutput ) - ( zinit + zinput )
185
186
187	! mass balance for Clay
188	zinit = zsolcp_inv_i(jsclay)
189	zfinal = zsolcp_inv_f(jsclay)
190	zinput = rain_tot (jsclay) + fromsed_tot(jsclay)
191	zoutput = tosed_tot (jsclay) + rloss_tot (jsclay)
192	zdelta_clay= ( zfinal + zoutput ) - ( zinit + zinput )
193
194	! mass balance for carbon ( carbon in POC, CaCo3, DIC )
195	zinit = zsolcp_inv_i(jspoc) + zsolcp_inv_i(jscal) + zpwcp_inv_i(jwdic)
196	zfinal = zsolcp_inv_f(jspoc) + zsolcp_inv_f(jscal) + zpwcp_inv_f(jwdic)
197	zinput = rain_tot (jspoc) + rain_tot (jscal) + diss_in_tot(jwdic)
198	zoutput = tosed_tot (jspoc) + tosed_tot (jscal) + diss_out_tot(jwdic) &
199	& + rloss_tot (jspoc) + rloss_tot (jscal)
200	zdelta_co2 = ( zfinal + zoutput ) - ( zinit + zinput )
201
202	! mass balance for oxygen : O2 is in POC remineralization
203	zinit = zpwcp_inv_i(jwoxy)
204	zfinal = zpwcp_inv_f(jwoxy)
205	zinput = diss_in_tot(jwoxy)
206	zoutput = diss_out_tot(jwoxy) + cons_tot_o2
207	zdelta_oxy = ( zfinal + zoutput ) - ( zinit + zinput )
208
209	! mass balance for phosphate ( PO4 in POC and dissolved phosphates )
210	zinit = zsolcp_inv_i(jspoc) * spo4r + zpwcp_inv_i(jwpo4)
211	zfinal = zsolcp_inv_f(jspoc) * spo4r + zpwcp_inv_f(jwpo4)
212	zinput = rain_tot (jspoc) * spo4r + diss_in_tot(jwpo4)
213	zoutput = tosed_tot (jspoc) * spo4r + diss_out_tot(jwpo4) &
214	& + rloss_tot (jspoc) * spo4r
215	zdelta_po4 = ( zfinal + zoutput ) - ( zinit + zinput )
216
217
218	! mass balance for Nitrate
219	zinit = zpwcp_inv_i (jwno3)
220	zfinal = zpwcp_inv_f (jwno3)
221	zinput = diss_in_tot (jwno3) + sour_tot_no3
222	zoutput = diss_out_tot(jwno3) + cons_tot_no3
223	zdelta_no3 = ( zfinal + zoutput ) - ( zinit + zinput )
224
225	! mass balance for DIC13
226	zinit = zpwcp_inv_i(jwc13) &
227	& + src13p * zsolcp_inv_i(jspoc) + src13Ca * zsolcp_inv_i(jscal)
228	zfinal = zpwcp_inv_f(jwc13) &
229	& + src13p * zsolcp_inv_f(jspoc) + src13Ca * zsolcp_inv_f(jscal)
230	zinput = diss_in_tot (jwc13) &
231	& + src13p * rain_tot(jspoc) + src13Ca * rain_tot(jscal)
232	zoutput = diss_out_tot(jwc13) &
233	& + src13p * tosed_tot(jspoc) + src13Ca * tosed_tot(jscal) &
234	& + src13p * rloss_tot(jspoc) + src13Ca * rloss_tot(jscal)
235	zdelta_c13 = ( zfinal + zoutput ) - ( zinit + zinput )
236
237	! other mass balance for DIC13
238	zinit = zpwcp_inv_i (jwc13)
239	zfinal = zpwcp_inv_f (jwc13)
240	zinput = diss_in_tot (jwc13) + sour_tot_c13
241	zoutput = diss_out_tot(jwc13)
242	zdelta_c13b= ( zfinal + zoutput ) - ( zinit + zinput )
243
244	END IF
245
246	IF( kt == nitsedend) THEN
247
248	WRITE(numsed,*)
249	WRITE(numsed,*)'================== General mass balance ================== '
250	WRITE(numsed,*)' '
251	WRITE(numsed,*)' '
252	WRITE(numsed,*)' Initial total solid Masses (mole/dx.dy) (k=2-11) '
253	WRITE(numsed,*)' ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~'
254	WRITE(numsed,*)' Opale, Clay, POC, CaCO3, C13'
255	WRITE(numsed,'(4x,5(1PE10.3,2X))')zsolcp_inv_i(jsopal),zsolcp_inv_i(jsclay),zsolcp_inv_i(jspoc), &
256	& zsolcp_inv_i(jscal),( src13P * zsolcp_inv_i(jspoc) + src13Ca * zsolcp_inv_i(jscal) )
257	WRITE(numsed,*)' '
258	WRITE(numsed,*)' Initial total dissolved Masses (mole/dx.dy) (k=2-11) '
259	WRITE(numsed,*)' ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~'
260	WRITE(numsed,*)' Si, O2, DIC, Nit Phos, DIC13'
261	WRITE(numsed,'(4x,6(1PE10.3,2X))') zpwcp_inv_i(jwsil), zpwcp_inv_i(jwoxy), &
262	& zpwcp_inv_i(jwdic), zpwcp_inv_i(jwno3), zpwcp_inv_i(jwpo4), zpwcp_inv_i(jwc13)
263	WRITE(numsed,*)' '
264	WRITE(numsed,*)' Solid inputs : Opale, Clay, POC, CaCO3, C13'
265	WRITE(numsed,'(A4,10X,5(1PE10.3,2X))')'Rain : ',rain_tot(jsopal),rain_tot(jsclay),rain_tot(jspoc),&
266	& rain_tot(jscal),( src13P * rain_tot(jspoc) + src13Ca * rain_tot(jscal) )
267	WRITE(numsed,'(A12,6x,4(1PE10.3,2X))')' From Sed : ',fromsed_tot(jsopal), fromsed_tot(jsclay), &
268	& fromsed_tot(jspoc), fromsed_tot(jscal)
269	WRITE(numsed,*)'Diss. inputs : Si, O2, DIC, Nit, Phos, DIC13'
270	WRITE(numsed,'(A9,1x,6(1PE10.3,2X))')' From Pisc : ', diss_in_tot(jwsil), &
271	& diss_in_tot(jwoxy), diss_in_tot(jwdic), diss_in_tot(jwno3), diss_in_tot(jwpo4), &
272	& diss_in_tot(jwc13)
273	WRITE(numsed,*)' '
274	WRITE(numsed,*)'Solid output : Opale, Clay, POC, CaCO3, C13'
275	WRITE(numsed,'(A6,8x,5(1PE10.3,2X))')'To sed', tosed_tot(jsopal),tosed_tot(jsclay),tosed_tot(jspoc),&
276	& tosed_tot(jscal),( src13P * tosed_tot(jspoc) + src13Ca * tosed_tot(jscal) )
277	WRITE(numsed,'(A5,9x,5(1PE10.3,2X))')'Perdu', rloss_tot(jsopal),rloss_tot(jsclay),rloss_tot(jspoc),&
278	& rloss_tot(jscal),( src13P * rloss_tot(jspoc) + src13Ca * rloss_tot(jscal) )
279	WRITE(numsed,*)'Diss. output : Si, O2, DIC, Nit, Phos, DIC13 '
280	WRITE(numsed,'(A7,2x,6(1PE10.3,2X))')'To kbot', diss_out_tot(jwsil), &
281	& diss_out_tot(jwoxy), diss_out_tot(jwdic), diss_out_tot(jwno3), diss_out_tot(jwpo4), &
282	& diss_out_tot(jwc13)
283	WRITE(numsed,*)' '
284	WRITE(numsed,*)' Total consomation in POC remineralization [mol]: O2, NO3'
285	WRITE(numsed,'(51x,2(1PE10.3,2X))') cons_tot_o2,cons_tot_no3
286	WRITE(numsed,*)' '
287	WRITE(numsed,*)'Final solid Masses (mole/dx.dy) (k=2-11)'
288	WRITE(numsed,*)' Opale, Clay, POC, CaCO3, C13'
289	WRITE(numsed,'(4x,5(1PE10.3,2X))')zsolcp_inv_f(jsopal),zsolcp_inv_f(jsclay),zsolcp_inv_f(jspoc), &
290	& zsolcp_inv_f(jscal),( src13P * zsolcp_inv_f(jspoc) + src13Ca * zsolcp_inv_f(jscal) )
291	WRITE(numsed,*)' '
292	WRITE(numsed,*)'Final dissolved Masses (mole/dx.dy) (k=2-11)'
293	WRITE(numsed,*)' Si, O2, DIC, Nit, Phos, DIC13'
294	WRITE(numsed,'(4x,6(1PE10.3,2X))') zpwcp_inv_f(jwsil), zpwcp_inv_f(jwoxy), &
295	& zpwcp_inv_f(jwdic), zpwcp_inv_f(jwno3), zpwcp_inv_f(jwpo4), zpwcp_inv_f(jwc13)
296	WRITE(numsed,*)' '
297	WRITE(numsed,*)'Delta : Opale, Clay, C, O, N, P, C13'
298	WRITE(numsed,'(7x,7(1PE11.3,1X))') zdelta_sil, zdelta_clay, zdelta_co2, zdelta_oxy, zdelta_no3,&
299	& zdelta_po4, zdelta_c13
300	WRITE(numsed,*)' '
301	WRITE(numsed,*)'deltaC13bis : ',zdelta_c13b
302
303	WRITE(numsed,*)'=========================================================================='
304	WRITE(numsed,*)' Composition of final sediment for point jpoce'
305	WRITE(numsed,*)' ========================================='
306	WRITE(numsed,*)'Opale, Clay, POC, CaCo3, hipor, pH, co3por'
307	DO jk = 1,jpksed
308	WRITE(numsed,'(4(F8.4,4X),3(1PE10.3,2X))') solcp(jpoce,jk,jsopal)100.,solcp(jpoce,jk,jsclay)100.,&
309	& solcp(jpoce,jk,jspoc)100.,solcp(jpoce,jk,jscal)100.,&
310	& hipor(jpoce,jk),-LOG10(hipor(jpoce,jk)/densSW(jpoce)),co3por(jpoce,jk)
311	ENDDO
312	WRITE(numsed,*)'Silicic A., Oxygen, DIC, Nitrats, Phosphats, Alkal., DIC13'
313	DO jk = 1, jpksed
314	WRITE(numsed,'(8(1PE10.3,2X))')pwcp(jpoce,jk,jwsil),pwcp(jpoce,jk,jwoxy),&
315	& pwcp(jpoce,jk,jwdic),pwcp(jpoce,jk,jwno3),pwcp(jpoce,jk,jwpo4),pwcp(jpoce,jk,jwalk),pwcp(jpoce,jk,jwc13)
316	ENDDO
317	WRITE(numsed,*)'densSW at the end : ',densSW(jpoce)
318	WRITE(numsed,*)'=========================================================================='
319
320	ENDIF
321
322	END SUBROUTINE sed_mbc
323
324
325	#else
326	!!======================================================================
327	!! MODULE sedmbc : Dummy module
328	!!======================================================================
329	CONTAINS
330	SUBROUTINE sed_mbc( kt ) ! Empty routine
331	INTEGER, INTENT(in) :: kt
332	WRITE(,) 'sed_mbc: You should not have seen this print! error?', kt
333	END SUBROUTINE sed_mbc
334	#endif
335	END MODULE sedmbc

Note: See TracBrowser for help on using the repository browser.

New URL for NEMO forge! http://forge.nemo-ocean.eu

Context Navigation

source: branches/nemo_v3_3_beta/NEMOGCM/NEMO/TOP_SRC/SED/sedmbc.F90 @ 2281

Download in other formats: