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.
Changeset 775 for branches/dev_001_GM/NEMO/TOP_SRC/PISCES_SMS/p4zbio_std.h90 – NEMO

Ignore:
Timestamp:
2007-12-19T14:45:15+01:00 (16 years ago)
Author:
gm
Message:

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

File:
1 moved

Legend:

Unmodified
Added
Removed
  • branches/dev_001_GM/NEMO/TOP_SRC/PISCES_SMS/p4zbio_std.h90

    r774 r775  
    1 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
    2 CCC p4zbio : PISCES MODEL  - Standard parameterization 
    3 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
     1      !CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
     2      !CC p4zbio : PISCES MODEL  - Standard parameterization 
     3      !CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
    44 
    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         ! 
    1949      END DO 
    20  
    21       DO jk = 1,jpkm1 
    22         DO jj = 1,jpj 
    23           DO ji = 1,jpi 
    24 C      
    25 C     Evolution of Detritus 
    26 C     --------------------- 
    27 C      
    28           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 C     
    35 C     Evolution of rapid Detritus 
    36 C     --------------------- 
    37 C     
    38           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 C 
    44           END DO 
    45         END DO 
    46       END DO 
    47  
    48  
    49  
    50       DO jk = 1,jpkm1 
    51         DO jj = 1,jpj 
    52           DO ji = 1,jpi 
    53 C 
    54 C     Evolution of small biogenic Iron 
    55 C     -------------------------- 
    56 C 
    57           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 C 
    66 C     Evolution of big biogenic Iron 
    67 C     -------------------------- 
    68 C 
    69           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 DO 
    78         END DO 
    79       END DO 
    80  
Note: See TracChangeset for help on using the changeset viewer.