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 | CDIR$ LIST |
---|
7 | SUBROUTINE p4zlim |
---|
8 | #if defined key_passivetrc && defined key_trc_pisces |
---|
9 | CCC--------------------------------------------------------------------- |
---|
10 | CCC |
---|
11 | CCC ROUTINE p4zlim : PISCES MODEL |
---|
12 | CCC ***************************** |
---|
13 | CCC |
---|
14 | CCC PURPOSE : |
---|
15 | CCC --------- |
---|
16 | CCC Compute the co-limitations by the various nutrients |
---|
17 | CCC for the various phytoplankton species |
---|
18 | CCC |
---|
19 | CC INPUT : |
---|
20 | CC ----- |
---|
21 | CC argument |
---|
22 | CC None |
---|
23 | CC common |
---|
24 | CC all the common defined in opa |
---|
25 | CC |
---|
26 | CC |
---|
27 | CC OUTPUT : : no |
---|
28 | CC ------ |
---|
29 | CC |
---|
30 | CC MODIFICATIONS: |
---|
31 | CC -------------- |
---|
32 | CC original : O. Aumont (2004) |
---|
33 | CC---------------------------------------------------------------------- |
---|
34 | CC parameters and commons |
---|
35 | CC ====================== |
---|
36 | CDIR$ NOLIST |
---|
37 | USE oce_trc |
---|
38 | USE trp_trc |
---|
39 | USE sms |
---|
40 | IMPLICIT NONE |
---|
41 | #include "domzgr_substitute.h90" |
---|
42 | CDIR$ LIST |
---|
43 | CC---------------------------------------------------------------------- |
---|
44 | CC local declarations |
---|
45 | CC ================== |
---|
46 | INTEGER ji, jj, jk |
---|
47 | REAL xlim1,xlim2,xlim3,xlim4 |
---|
48 | REAL xconctemp,xconctemp2,xconctempn,xconctempn2 |
---|
49 | C |
---|
50 | C Tuning of the iron concentration to a minimum |
---|
51 | C level that is set to the detection limit |
---|
52 | C ------------------------------------- |
---|
53 | C |
---|
54 | |
---|
55 | trn(:,:,:,jpfer)=max(trn(:,:,:,jpfer),1.E-11) |
---|
56 | C |
---|
57 | C Computation of a variable Ks for iron on diatoms |
---|
58 | C taking into account that increasing biomass is |
---|
59 | C made of generally bigger cells |
---|
60 | C ------------------------------------------------ |
---|
61 | C |
---|
62 | DO jk=1,jpkm1 |
---|
63 | DO jj=1,jpj |
---|
64 | DO ji=1,jpi |
---|
65 | xconctemp=max(0.,trn(ji,jj,jk,jpdia)-5E-7) |
---|
66 | xconctemp2=min(5.E-7,trn(ji,jj,jk,jpdia)) |
---|
67 | xconctempn=max(0.,trn(ji,jj,jk,jpphy)-1E-6) |
---|
68 | xconctempn2=min(1.E-6,trn(ji,jj,jk,jpphy)) |
---|
69 | concdfe(ji,jj,jk)=(xconctemp2*conc3+0.4E-9* |
---|
70 | . xconctemp)/(xconctemp2+xconctemp+rtrn) |
---|
71 | concdfe(ji,jj,jk)=max(conc3,concdfe(ji,jj,jk)) |
---|
72 | concnfe(ji,jj,jk)=(xconctempn2*conc2+0.08E-9* |
---|
73 | . xconctempn)/(xconctempn2+xconctempn+rtrn) |
---|
74 | concnfe(ji,jj,jk)=max(conc2,concnfe(ji,jj,jk)) |
---|
75 | END DO |
---|
76 | END DO |
---|
77 | END DO |
---|
78 | C |
---|
79 | DO jk = 1,jpkm1 |
---|
80 | DO jj = 1,jpj |
---|
81 | DO ji = 1,jpi |
---|
82 | C |
---|
83 | C Michaelis-Menten Limitation term for nutrients |
---|
84 | C Small flagellates |
---|
85 | C ----------------------------------------------- |
---|
86 | C |
---|
87 | xnanono3(ji,jj,jk)=trn(ji,jj,jk,jpno3)*concnnh4 |
---|
88 | & /(conc0*concnnh4+concnnh4*trn(ji,jj,jk,jpno3)+ |
---|
89 | & conc0*trn(ji,jj,jk,jpnh4)) |
---|
90 | xnanonh4(ji,jj,jk)=trn(ji,jj,jk,jpnh4)*conc0 |
---|
91 | & /(conc0*concnnh4+concnnh4*trn(ji,jj,jk,jpno3)+ |
---|
92 | & conc0*trn(ji,jj,jk,jpnh4)) |
---|
93 | xlim1=xnanono3(ji,jj,jk)+xnanonh4(ji,jj,jk) |
---|
94 | xlim2=trn(ji,jj,jk,jppo4)/(trn(ji,jj,jk,jppo4)+concnnh4) |
---|
95 | xlim3=trn(ji,jj,jk,jpfer)/(trn(ji,jj,jk,jpfer) |
---|
96 | & +concnfe(ji,jj,jk)) |
---|
97 | xlimphy(ji,jj,jk)=min(xlim1,xlim2,xlim3) |
---|
98 | xlim1=trn(ji,jj,jk,jpnh4)/(concnnh4+trn(ji,jj,jk,jpnh4)) |
---|
99 | xlim3=trn(ji,jj,jk,jpfer)/(trn(ji,jj,jk,jpfer)+conc2) |
---|
100 | xlim4=trn(ji,jj,jk,jpdoc)/(trn(ji,jj,jk,jpdoc)+xkdoc2) |
---|
101 | xlimbac(ji,jj,jk)=min(xlim1,xlim2,xlim3)*xlim4 |
---|
102 | C |
---|
103 | END DO |
---|
104 | END DO |
---|
105 | END DO |
---|
106 | C |
---|
107 | DO jk = 1,jpkm1 |
---|
108 | DO jj = 1,jpj |
---|
109 | DO ji = 1,jpi |
---|
110 | C |
---|
111 | C Michaelis-Menten Limitation term for nutrients |
---|
112 | C Diatoms |
---|
113 | C ---------------------------------------------- |
---|
114 | C |
---|
115 | xdiatno3(ji,jj,jk)=trn(ji,jj,jk,jpno3)*concdnh4 |
---|
116 | & /(conc1*concdnh4+concdnh4*trn(ji,jj,jk,jpno3)+ |
---|
117 | & conc1*trn(ji,jj,jk,jpnh4)) |
---|
118 | xdiatnh4(ji,jj,jk)=trn(ji,jj,jk,jpnh4)*conc1 |
---|
119 | & /(conc1*concdnh4+concdnh4*trn(ji,jj,jk,jpno3)+ |
---|
120 | & conc1*trn(ji,jj,jk,jpnh4)) |
---|
121 | |
---|
122 | xlim1=xdiatno3(ji,jj,jk)+xdiatnh4(ji,jj,jk) |
---|
123 | xlim2=trn(ji,jj,jk,jppo4)/(trn(ji,jj,jk,jppo4)+concdnh4) |
---|
124 | xlim3=trn(ji,jj,jk,jpsil)/(trn(ji,jj,jk,jpsil)+xksi(ji,jj)) |
---|
125 | xlim4=trn(ji,jj,jk,jpfer)/(trn(ji,jj,jk,jpfer) |
---|
126 | & +concdfe(ji,jj,jk)) |
---|
127 | xlimdia(ji,jj,jk)=min(xlim1,xlim2,xlim3,xlim4) |
---|
128 | xlim4=trn(ji,jj,jk,jpfer)/(trn(ji,jj,jk,jpfer)+conc3) |
---|
129 | xlimdia2(ji,jj,jk)=min(xlim1,xlim2,xlim3,xlim4) |
---|
130 | C |
---|
131 | END DO |
---|
132 | END DO |
---|
133 | END DO |
---|
134 | C |
---|
135 | #endif |
---|
136 | RETURN |
---|
137 | END |
---|