source: branches/dev_001_GM/NEMO/TOP_SRC/PISCES_SMS/p4zbio_std.h90 @ 775

Last change on this file since 775 was 775, checked in by gm, 13 years ago

dev_001_GM - PISCES in F90 : encapsulation of all p4z…F files in module F90 + doctor norme for local variables - compilation OK

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 3.5 KB
Line 
1      !CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
2      !CC p4zbio : PISCES MODEL  - Standard parameterization
3      !CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
4
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         !
49      END DO
Note: See TracBrowser for help on using the repository browser.