[341] | 1 | |
---|
[719] | 2 | CCC $Header$ |
---|
[341] | 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 --------------------------------------------------------------------------- |
---|
[333] | 6 | CDIR$ LIST |
---|
| 7 | SUBROUTINE p4zlim |
---|
[772] | 8 | #if defined key_top && defined key_pisces |
---|
[333] | 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 |
---|
[655] | 47 | REAL xlim1,xlim2,xlim3,xlim4,zno3,zferlim |
---|
[333] | 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 |
---|
[655] | 54 | do jk=1,jpk |
---|
| 55 | do jj=1,jpj |
---|
| 56 | do ji=1,jpi |
---|
| 57 | zno3=trn(ji,jj,jk,jpno3)*1E6 |
---|
| 58 | zferlim=max(1.5E-11*(zno3/40)**2,3E-12) |
---|
| 59 | zferlim=min(zferlim,1.5E-11) |
---|
| 60 | trn(ji,jj,jk,jpfer)=max(trn(ji,jj,jk,jpfer),zferlim) |
---|
| 61 | end do |
---|
| 62 | end do |
---|
| 63 | end do |
---|
[333] | 64 | C |
---|
| 65 | C Computation of a variable Ks for iron on diatoms |
---|
| 66 | C taking into account that increasing biomass is |
---|
| 67 | C made of generally bigger cells |
---|
| 68 | C ------------------------------------------------ |
---|
| 69 | C |
---|
| 70 | DO jk=1,jpkm1 |
---|
| 71 | DO jj=1,jpj |
---|
| 72 | DO ji=1,jpi |
---|
| 73 | xconctemp=max(0.,trn(ji,jj,jk,jpdia)-5E-7) |
---|
| 74 | xconctemp2=min(5.E-7,trn(ji,jj,jk,jpdia)) |
---|
| 75 | xconctempn=max(0.,trn(ji,jj,jk,jpphy)-1E-6) |
---|
| 76 | xconctempn2=min(1.E-6,trn(ji,jj,jk,jpphy)) |
---|
| 77 | concdfe(ji,jj,jk)=(xconctemp2*conc3+0.4E-9* |
---|
| 78 | . xconctemp)/(xconctemp2+xconctemp+rtrn) |
---|
| 79 | concdfe(ji,jj,jk)=max(conc3,concdfe(ji,jj,jk)) |
---|
| 80 | concnfe(ji,jj,jk)=(xconctempn2*conc2+0.08E-9* |
---|
| 81 | . xconctempn)/(xconctempn2+xconctempn+rtrn) |
---|
| 82 | concnfe(ji,jj,jk)=max(conc2,concnfe(ji,jj,jk)) |
---|
| 83 | END DO |
---|
| 84 | END DO |
---|
| 85 | END DO |
---|
| 86 | C |
---|
| 87 | DO jk = 1,jpkm1 |
---|
| 88 | DO jj = 1,jpj |
---|
| 89 | DO ji = 1,jpi |
---|
| 90 | C |
---|
| 91 | C Michaelis-Menten Limitation term for nutrients |
---|
| 92 | C Small flagellates |
---|
| 93 | C ----------------------------------------------- |
---|
| 94 | C |
---|
| 95 | xnanono3(ji,jj,jk)=trn(ji,jj,jk,jpno3)*concnnh4 |
---|
| 96 | & /(conc0*concnnh4+concnnh4*trn(ji,jj,jk,jpno3)+ |
---|
| 97 | & conc0*trn(ji,jj,jk,jpnh4)) |
---|
| 98 | xnanonh4(ji,jj,jk)=trn(ji,jj,jk,jpnh4)*conc0 |
---|
| 99 | & /(conc0*concnnh4+concnnh4*trn(ji,jj,jk,jpno3)+ |
---|
| 100 | & conc0*trn(ji,jj,jk,jpnh4)) |
---|
| 101 | xlim1=xnanono3(ji,jj,jk)+xnanonh4(ji,jj,jk) |
---|
| 102 | xlim2=trn(ji,jj,jk,jppo4)/(trn(ji,jj,jk,jppo4)+concnnh4) |
---|
| 103 | xlim3=trn(ji,jj,jk,jpfer)/(trn(ji,jj,jk,jpfer) |
---|
| 104 | & +concnfe(ji,jj,jk)) |
---|
| 105 | xlimphy(ji,jj,jk)=min(xlim1,xlim2,xlim3) |
---|
| 106 | xlim1=trn(ji,jj,jk,jpnh4)/(concnnh4+trn(ji,jj,jk,jpnh4)) |
---|
| 107 | xlim3=trn(ji,jj,jk,jpfer)/(trn(ji,jj,jk,jpfer)+conc2) |
---|
| 108 | xlim4=trn(ji,jj,jk,jpdoc)/(trn(ji,jj,jk,jpdoc)+xkdoc2) |
---|
| 109 | xlimbac(ji,jj,jk)=min(xlim1,xlim2,xlim3)*xlim4 |
---|
| 110 | C |
---|
| 111 | END DO |
---|
| 112 | END DO |
---|
| 113 | END DO |
---|
| 114 | C |
---|
| 115 | DO jk = 1,jpkm1 |
---|
| 116 | DO jj = 1,jpj |
---|
| 117 | DO ji = 1,jpi |
---|
| 118 | C |
---|
| 119 | C Michaelis-Menten Limitation term for nutrients |
---|
| 120 | C Diatoms |
---|
| 121 | C ---------------------------------------------- |
---|
| 122 | C |
---|
| 123 | xdiatno3(ji,jj,jk)=trn(ji,jj,jk,jpno3)*concdnh4 |
---|
| 124 | & /(conc1*concdnh4+concdnh4*trn(ji,jj,jk,jpno3)+ |
---|
| 125 | & conc1*trn(ji,jj,jk,jpnh4)) |
---|
| 126 | xdiatnh4(ji,jj,jk)=trn(ji,jj,jk,jpnh4)*conc1 |
---|
| 127 | & /(conc1*concdnh4+concdnh4*trn(ji,jj,jk,jpno3)+ |
---|
| 128 | & conc1*trn(ji,jj,jk,jpnh4)) |
---|
| 129 | |
---|
| 130 | xlim1=xdiatno3(ji,jj,jk)+xdiatnh4(ji,jj,jk) |
---|
| 131 | xlim2=trn(ji,jj,jk,jppo4)/(trn(ji,jj,jk,jppo4)+concdnh4) |
---|
| 132 | xlim3=trn(ji,jj,jk,jpsil)/(trn(ji,jj,jk,jpsil)+xksi(ji,jj)) |
---|
| 133 | xlim4=trn(ji,jj,jk,jpfer)/(trn(ji,jj,jk,jpfer) |
---|
| 134 | & +concdfe(ji,jj,jk)) |
---|
| 135 | xlimdia(ji,jj,jk)=min(xlim1,xlim2,xlim3,xlim4) |
---|
| 136 | C |
---|
| 137 | END DO |
---|
| 138 | END DO |
---|
| 139 | END DO |
---|
| 140 | C |
---|
| 141 | #endif |
---|
| 142 | RETURN |
---|
| 143 | END |
---|