Ignore:
Timestamp:
2020-07-02T20:34:16+02:00 (8 months ago)
Author:
aumont
Message:

update of the PISCES comments

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/P4Z/p4zrem.F90

    r12682 r13233  
    9292      ! Computation of the mean bacterial concentration 
    9393      ! this parameterization has been deduced from a model version 
    94       ! that was modeling explicitely bacteria 
    95       ! ------------------------------------------------------- 
     94      ! that was modeling explicitely bacteria. This is a very old param  
     95      ! that will be very soon updated based on results from a much more 
     96      ! recent version of PISCES with bacteria. 
     97      ! ---------------------------------------------------------------- 
    9698      DO jk = 1, jpkm1 
    9799         DO jj = 1, jpj 
     
    111113      END DO 
    112114 
    113       IF( ln_p4z ) THEN 
     115      IF( ln_p4z ) THEN ! Standard PISCES code 
    114116         DO jk = 1, jpkm1 
    115117            DO jj = 1, jpj 
     
    152154            END DO 
    153155         END DO 
    154       ELSE 
     156      ELSE ! PISCES-QUOTA 
    155157         DO jk = 1, jpkm1 
    156158            DO jj = 1, jpj 
     
    231233 
    232234               ! Bacterial uptake of iron. No iron is available in DOC. So 
    233                ! Bacteries are obliged to take up iron from the water. Some 
     235               ! Bacteria are obliged to take up iron from the water. Some 
    234236               ! studies (especially at Papa) have shown this uptake to be significant 
    235237               ! --------------------------------------------------------------------- 
     
    262264         DO jj = 1, jpj 
    263265            DO ji = 1, jpi 
    264                ! Remineralization rate of BSi depedant on T and saturation 
     266 
     267               ! Remineralization rate of BSi dependent on T and saturation 
    265268               ! The parameterization is taken from Ridgwell et al. (2002)  
    266269               ! --------------------------------------------------------- 
     
    269272               zsatur2  = ( 1. + tsn(ji,jj,jk,jp_tem) / 400.)**37 
    270273               znusil   = 0.225  * ( 1. + tsn(ji,jj,jk,jp_tem) / 15.) * zsatur + 0.775 * zsatur2 * zsatur**9.25 
     274 
     275               ! Two fractions of bSi are considered : a labile one and a more 
     276               ! refractory one based on the commonly observed two step  
     277               ! dissolution of bSi (initial rapid dissolution followed by  
     278               ! more slowly dissolution). 
    271279               ! Computation of the vertical evolution of the labile fraction 
    272280               ! of bSi. This is computed assuming steady state. 
     281               ! -------------------------------------------------------------- 
    273282               IF ( gdept_n(ji,jj,jk) > zdep ) THEN 
    274283                  zfacsib(ji,jj,jk) = zfacsib(ji,jj,jk-1) * EXP( -0.5 * ( xsiremlab - xsirem )  & 
     
    280289               zsiremin = ( xsiremlab * zfacsi(ji,jj,jk) + xsirem * ( 1. - zfacsi(ji,jj,jk) ) ) * xstep * znusil 
    281290               zosil    = zsiremin * trb(ji,jj,jk,jpgsi) 
     291 
    282292               ! Update of the TRA arrays 
    283293               tra(ji,jj,jk,jpgsi) = tra(ji,jj,jk,jpgsi) - zosil 
Note: See TracChangeset for help on using the changeset viewer.