New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
p4zbio_kriest.h90 in branches/dev_001_GM/NEMO/TOP_SRC/PISCES_SMS – NEMO

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

Last change on this file since 775 was 775, checked in by gm, 16 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: 4.6 KB
Line 
1
2      !CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
3      !CC p4zbio : PISCES MODEL  - Kriest parameterization
4      !CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
5
6      DO jk = 1, jpkm1
7         DO jj = 1, jpj
8            DO ji = 1, jpi
9   
10               !     Evolution of DOC
11               !     ----------------     
12               trn(ji,jj,jk,jpdoc) = trn(ji,jj,jk,jpdoc) + orem(ji,jj,jk)                                  &
13                  &                                      + excret2 * prorca2(ji,jj,jk)                     &
14                  &                                      + excret  * prorca (ji,jj,jk) - olimi(ji,jj,jk)   &
15                  &                                      - denitr(ji,jj,jk)                                &
16                  &                                      + grarem (ji,jj,jk) * (1.-sigma1)                 &
17                  &                                      + grarem2(ji,jj,jk) * (1.-sigma2) - xaggdoc(ji,jj,jk)
18     
19               !     Evolution of Detritus
20               !     ---------------------
21               znumpoc = trn(ji,jj,jk,jpnum) / ( trn(ji,jj,jk,jppoc) + rtrn )
22               znumdoc = znumpoc
23
24               trn(ji,jj,jk,jppoc) = trn(ji,jj,jk,jppoc) - grazpoc(ji,jj,jk)                       &
25                  &                                      + grapoc (ji,jj,jk) + grapoc2(ji,jj,jk)   &
26                  &                                      - grazm  (ji,jj,jk) + tortz2 (ji,jj,jk)   &
27                  &                                      + respz  (ji,jj,jk) + respz2 (ji,jj,jk)   &
28                  &                                      + respp  (ji,jj,jk) + respp2 (ji,jj,jk)   &
29                  &                                      + tortp2 (ji,jj,jk) + tortz  (ji,jj,jk)   &
30                  &                                      + tortp  (ji,jj,jk) - orem   (ji,jj,jk)   &
31                  &                                      + xaggdoc(ji,jj,jk) - grazffe(ji,jj,jk)
32
33               !     Evolution of number of aggregates
34               !     ---------------------------------
35               trn(ji,jj,jk,jpnum) = trn(ji,jj,jk,jpnum) - xagg(ji,jj,jk)                                           &
36                  &                                      - ( orem(ji,jj,jk) + grazpoc(ji,jj,jk) ) * znumpoc         &
37                  &                                      + ( tortp(ji,jj,jk) ) * xkr_nnano                          &
38                  &                                      + ( tortp2(ji,jj,jk) + respp(ji,jj,jk) + tortz(ji,jj,jk)   &
39                  &                                        + grapoc(ji,jj,jk) - grazm(ji,jj,jk)                     &
40                  &                                        + respz(ji,jj,jk) ) * xkr_ndiat                          &
41                  &                                      + ( grapoc2(ji,jj,jk) + tortz2(ji,jj,jk)                   &
42                  &                                        + respz2 (ji,jj,jk) ) * xkr_nmeso                        &
43                  &                                      + respp2(ji,jj,jk) * xkr_naggr                             &
44                  &                                      + xaggdoc(ji,jj,jk) * znumdoc                              &
45                  &                                      - grazffe(ji,jj,jk) * znumpoc * wsbio4(ji,jj,jk)           &
46                  &                                      / ( wsbio3(ji,jj,jk) + rtrn )
47
48               trn(ji,jj,jk,jpnum) = MAX( trn(ji,jj,jk,jpnum), trn(ji,jj,jk,jppoc) / xkr_massp / xnumm(jk) )
49
50               trn(ji,jj,jk,jpnum) = MIN( trn(ji,jj,jk,jpnum), trn(ji,jj,jk,jppoc) / xkr_massp / 1.1       )
51
52               !     Evolution of biogenic Iron
53               !     --------------------------
54               trn(ji,jj,jk,jpsfe) = trn(ji,jj,jk,jpsfe) + unass * ( grazpf(ji,jj,jk) + grazsf(ji,jj,jk) )   &
55                  &                - ( 1.- unass2 ) * grazpof(ji,jj,jk) - ( 1.- unass ) * grazmf(ji,jj,jk)   &
56                  &                - ( 1.- unass2 ) * grazfff(ji,jj,jk) +       unass2  * ( graznf(ji,jj,jk)   &
57                  &                + grazf(ji,jj,jk) + ferat3 * grazz(ji,jj,jk) ) + ferat3                   &
58                  &                * (tortz2(ji,jj,jk)+respz2(ji,jj,jk)+tortz(ji,jj,jk)                      &
59                  &                + respz(ji,jj,jk) ) - ofer(ji,jj,jk) + ( respnf(ji,jj,jk)                 &
60                  &                + tortnf(ji,jj,jk) ) + tortdf(ji,jj,jk) + respdf(ji,jj,jk)                &
61                  &                + xbactfer(ji,jj,jk) + xscave(ji,jj,jk) * zdenom1(ji,jj,jk)
62            END DO
63         END DO
64      END DO
Note: See TracBrowser for help on using the repository browser.