[341] | 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 --------------------------------------------------------------------------- |
---|
[333] | 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 |
---|