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 12537 for NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/P4Z/p4zint.F90 – NEMO

Ignore:
Timestamp:
2020-03-11T16:02:54+01:00 (4 years ago)
Author:
aumont
Message:

Comments in routines have been revised and significantly augmented

File:
1 edited

Legend:

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

    r12349 r12537  
    11MODULE p4zint 
    2    !!====================================================================== 
     2   !!========================================================================= 
    33   !!                         ***  MODULE p4zint  *** 
    44   !! TOP :   PISCES interpolation and computation of various accessory fields 
    5    !!====================================================================== 
     5   !!========================================================================= 
    66   !! History :   1.0  !  2004-03 (O. Aumont) Original code 
    77   !!             2.0  !  2007-12  (C. Ethe, G. Madec)  F90 
     
    4343      ! Computation of phyto and zoo metabolic rate 
    4444      ! ------------------------------------------- 
     45      ! Generic temperature dependence (Eppley, 1972) 
    4546      tgfunc (:,:,:) = EXP( 0.063913 * tsn(:,:,:,jp_tem) ) 
     47      ! Temperature dependence of mesozooplankton (Buitenhuis et al. (2005)) 
    4648      tgfunc2(:,:,:) = EXP( 0.07608  * tsn(:,:,:,jp_tem) ) 
     49      ! Temperature dependence of picophytoplankton (Stawiarsky et al., 2016) 
    4750      tgfunc3(:,:,:) = EXP( 0.0825   * tsn(:,:,:,jp_tem) ) 
    4851 
    4952      ! Computation of the silicon dependant half saturation  constant for silica uptake 
    50       ! --------------------------------------------------- 
     53      ! This is based on an old study by Pondaven et al. (1998) 
     54      ! -------------------------------------------------------------------------------- 
    5155      DO ji = 1, jpi 
    5256         DO jj = 1, jpj 
     
    5660      END DO 
    5761      ! 
     62      ! At the end of each year, the half saturation constant for silica is  
     63      ! updated as this is based on the highest concentration reached over  
     64      ! the year 
     65      ! ------------------------------------------------------------------- 
    5866      IF( nday_year == nyear_len(1) ) THEN 
    5967         xksi   (:,:) = xksimax(:,:) 
Note: See TracChangeset for help on using the changeset viewer.