- Timestamp:
- 2007-12-19T14:45:15+01:00 (16 years ago)
- File:
-
- 1 moved
Legend:
- Unmodified
- Added
- Removed
-
branches/dev_001_GM/NEMO/TOP_SRC/PISCES_SMS/p4zbio_std.h90
r774 r775 1 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC2 CCC p4zbio : PISCES MODEL - Standard parameterization3 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC1 !CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 2 !CC p4zbio : PISCES MODEL - Standard parameterization 3 !CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 4 4 5 DO jk = 1,jpkm1 6 DO jj = 1,jpj 7 DO ji = 1,jpi 8 C 9 C Evolution of DOC 10 C ---------------- 11 C 12 trn(ji,jj,jk,jpdoc) = trn(ji,jj,jk,jpdoc) 13 & +orem(ji,jj,jk)+excret2*prorca2(ji,jj,jk) 14 & +excret*prorca(ji,jj,jk)-olimi(ji,jj,jk)-denitr(ji,jj,jk) 15 & +grarem(ji,jj,jk)*(1.-sigma1)+grarem2(ji,jj,jk) 16 & *(1.-sigma2)-xaggdoc(ji,jj,jk)-xaggdoc2(ji,jj,jk) 17 END DO 18 END DO 5 DO jk = 1, jpkm1 6 7 ! Evolution of DOC 8 ! ---------------- 9 trn(:,:,jk,jpdoc) = trn(:,:,jk,jpdoc) + orem(:,:,jk) + excret2 * prorca2(:,:,jk) & 10 & + excret * prorca(:,:,jk) - olimi(:,:,jk) - denitr(:,:,jk) & 11 & + grarem(:,:,jk) * (1.-sigma1) + grarem2(:,:,jk) * (1.-sigma2) & 12 & - xaggdoc(:,:,jk) - xaggdoc2(:,:,jk) 13 14 ! Evolution of Detritus 15 ! --------------------- 16 trn(:,:,jk,jppoc) = trn(:,:,jk,jppoc) - grazpoc(:,:,jk) + grapoc(:,:,jk) - grazm (:,:,jk) & 17 & + respz (:,:,jk) - xagg (:,:,jk) + xaggdoc(:,:,jk) & 18 & + ( 1.-0.5 * zfracal(:,:,jk) ) * ( tortp(:,:,jk) + respp(:,:,jk) ) & 19 & + 0.5 * tortp2 (:,:,jk) & 20 & + orem2(:,:,jk) + tortz(:,:,jk) - orem(:,:,jk) 21 22 ! Evolution of rapid Detritus 23 ! --------------------------- 24 trn(:,:,jk,jpgoc) = trn(:,:,jk,jpgoc) + grapoc2(:,:,jk) + respp2(:,:,jk) + xagg (:,:,jk) & 25 & + tortz2 (:,:,jk) + respz2(:,:,jk) - orem2(:,:,jk) & 26 & + 0.5*zfracal(:,:,jk)*(respp(:,:,jk)+tortp(:,:,jk)) & 27 & + 0.5*tortp2(:,:,jk)+xaggdoc2(:,:,jk)-grazffe(:,:,jk) 28 29 ! Evolution of small biogenic Iron 30 ! -------------------------------- 31 trn(:,:,jk,jpsfe) = trn(:,:,jk,jpsfe) + unass * ( grazpf(:,:,jk) + grazsf(:,:,jk) ) & 32 & - grazpof(:,:,jk) - ( 1.- unass ) * grazmf(:,:,jk) & 33 & + ( 1.- 0.5 * zfracal(:,:,jk) ) * ( tortnf(:,:,jk) & 34 & + respnf(:,:,jk) ) + 0.5 * tortdf(:,:,jk) + ferat3 * & 35 & ( tortz(:,:,jk) + respz(:,:,jk) ) - ofer(:,:,jk) & 36 & + ofer2(:,:,jk) - xaggfe(:,:,jk) & 37 & + xscave(:,:,jk) * zdenom1(:,:,jk) 38 39 ! Evolution of big biogenic Iron 40 ! ------------------------------ 41 trn(:,:,jk,jpbfe) = trn(:,:,jk,jpbfe) + unass2 * ( graznf (:,:,jk) + grazf(:,:,jk) + grazfff(:,:,jk) & 42 & + grazpof(:,:,jk) + ferat3 * grazz(:,:,jk) ) & 43 & + ferat3 * ( tortz2 (:,:,jk) + respz2(:,:,jk) ) & 44 & - ofer2(:,:,jk) & 45 & + 0.5 * zfracal(:,:,jk) * ( respnf(:,:,jk) + tortnf(:,:,jk) ) & 46 & + 0.5 * tortdf (:,:,jk) + respdf(:,:,jk) + xaggfe(:,:,jk) & 47 & + xbactfer(:,:,jk) - grazfff(:,:,jk) + xscave(:,:,jk) * zdenom2(:,:,jk) 48 ! 19 49 END DO 20 21 DO jk = 1,jpkm122 DO jj = 1,jpj23 DO ji = 1,jpi24 C25 C Evolution of Detritus26 C ---------------------27 C28 trn(ji,jj,jk,jppoc) = trn(ji,jj,jk,jppoc)29 & -grazpoc(ji,jj,jk)+grapoc(ji,jj,jk)-grazm(ji,jj,jk)30 & +respz(ji,jj,jk)-xagg(ji,jj,jk)+xaggdoc(ji,jj,jk)31 & +(1.-0.5*zfracal(ji,jj,jk))*(tortp(ji,jj,jk)32 & +respp(ji,jj,jk))+0.5*tortp2(ji,jj,jk)33 & +orem2(ji,jj,jk)+tortz(ji,jj,jk)-orem(ji,jj,jk)34 C35 C Evolution of rapid Detritus36 C ---------------------37 C38 trn(ji,jj,jk,jpgoc) = trn(ji,jj,jk,jpgoc)39 & +grapoc2(ji,jj,jk)+respp2(ji,jj,jk)+xagg(ji,jj,jk)40 & +tortz2(ji,jj,jk)+respz2(ji,jj,jk)-orem2(ji,jj,jk)41 & +0.5*zfracal(ji,jj,jk)*(respp(ji,jj,jk)+tortp(ji,jj,jk))42 & +0.5*tortp2(ji,jj,jk)+xaggdoc2(ji,jj,jk)-grazffe(ji,jj,jk)43 C44 END DO45 END DO46 END DO47 48 49 50 DO jk = 1,jpkm151 DO jj = 1,jpj52 DO ji = 1,jpi53 C54 C Evolution of small biogenic Iron55 C --------------------------56 C57 trn(ji,jj,jk,jpsfe) = trn(ji,jj,jk,jpsfe)58 & +unass*(grazpf(ji,jj,jk)+grazsf(ji,jj,jk))59 & -grazpof(ji,jj,jk)-(1.-unass)*grazmf(ji,jj,jk)60 & +(1.-0.5*zfracal(ji,jj,jk))*(tortnf(ji,jj,jk)61 & +respnf(ji,jj,jk))+0.5*tortdf(ji,jj,jk)+ferat3*62 & (tortz(ji,jj,jk)+respz(ji,jj,jk))-ofer(ji,jj,jk)63 & +ofer2(ji,jj,jk)-xaggfe(ji,jj,jk)64 & +xscave(ji,jj,jk)*zdenom1(ji,jj,jk)65 C66 C Evolution of big biogenic Iron67 C --------------------------68 C69 trn(ji,jj,jk,jpbfe) = trn(ji,jj,jk,jpbfe)70 & +unass2*(graznf(ji,jj,jk)+grazf(ji,jj,jk)+grazfff(ji,jj,jk)71 & +grazpof(ji,jj,jk)+ferat3*grazz(ji,jj,jk))+ferat3*72 & (tortz2(ji,jj,jk)+respz2(ji,jj,jk))-ofer2(ji,jj,jk)73 & +0.5*zfracal(ji,jj,jk)*(respnf(ji,jj,jk)+tortnf(ji,jj,jk))74 & +0.5*tortdf(ji,jj,jk)+respdf(ji,jj,jk)+xaggfe(ji,jj,jk)75 & +xbactfer(ji,jj,jk)-grazfff(ji,jj,jk)+xscave(ji,jj,jk)76 & *zdenom2(ji,jj,jk)77 END DO78 END DO79 END DO80
Note: See TracChangeset
for help on using the changeset viewer.