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 9169 for branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmort.F90 – NEMO

Ignore:
Timestamp:
2017-12-26T17:32:56+01:00 (6 years ago)
Author:
gm
Message:

dev_merge_2017: all SRC: finalize the removal of useless warning when reading namelist_cfg + remove all nn_closea + nn_msh replaced by a logical

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmort.F90

    r9124 r9169  
    77   !!             2.0  !  2007-12  (C. Ethe, G. Madec)  F90 
    88   !!---------------------------------------------------------------------- 
    9    !!   p4z_mort       :   Compute the mortality terms for phytoplankton 
    10    !!   p4z_mort_init  :   Initialize the mortality params for phytoplankton 
    11    !!---------------------------------------------------------------------- 
    12    USE oce_trc         !  shared variables between ocean and passive tracers 
    13    USE trc             !  passive tracers common variables  
    14    USE sms_pisces      !  PISCES Source Minus Sink variables 
    15    USE p4zprod         !  Primary productivity  
    16    USE p4zlim          !  Phytoplankton limitation terms 
    17    USE prtctl_trc      !  print control for debugging 
     9   !!   p4z_mort       : Compute the mortality terms for phytoplankton 
     10   !!   p4z_mort_init  : Initialize the mortality params for phytoplankton 
     11   !!---------------------------------------------------------------------- 
     12   USE oce_trc         ! shared variables between ocean and passive tracers 
     13   USE trc             ! passive tracers common variables  
     14   USE sms_pisces      ! PISCES Source Minus Sink variables 
     15   USE p4zprod         ! Primary productivity  
     16   USE p4zlim          ! Phytoplankton limitation terms 
     17   USE prtctl_trc      ! print control for debugging 
    1818 
    1919   IMPLICIT NONE 
     
    2323   PUBLIC   p4z_mort_init     
    2424 
    25    !! * Shared module variables 
    26    REAL(wp), PUBLIC :: wchl    !: 
    27    REAL(wp), PUBLIC :: wchld   !: 
    28    REAL(wp), PUBLIC :: wchldm  !: 
    29    REAL(wp), PUBLIC :: mprat   !: 
    30    REAL(wp), PUBLIC :: mprat2  !: 
     25   REAL(wp), PUBLIC ::   wchl     !: 
     26   REAL(wp), PUBLIC ::   wchld    !: 
     27   REAL(wp), PUBLIC ::   wchldm   !: 
     28   REAL(wp), PUBLIC ::   mprat    !: 
     29   REAL(wp), PUBLIC ::   mprat2   !: 
    3130 
    3231   !!---------------------------------------------------------------------- 
     
    3534   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    3635   !!---------------------------------------------------------------------- 
    37  
    3836CONTAINS 
    3937 
     
    4947      INTEGER, INTENT(in) ::   kt ! ocean time step 
    5048      !!--------------------------------------------------------------------- 
    51  
     49      ! 
    5250      CALL p4z_nano            ! nanophytoplankton 
    53  
     51      ! 
    5452      CALL p4z_diat            ! diatoms 
    55  
     53      ! 
    5654   END SUBROUTINE p4z_mort 
    5755 
     
    6563      !! ** Method  : - ??? 
    6664      !!--------------------------------------------------------------------- 
    67       INTEGER  :: ji, jj, jk 
    68       REAL(wp) :: zsizerat, zcompaph 
    69       REAL(wp) :: zfactfe, zfactch, zprcaca, zfracal 
    70       REAL(wp) :: ztortp , zrespp , zmortp  
    71       CHARACTER (len=25) :: charout 
     65      INTEGER  ::   ji, jj, jk 
     66      REAL(wp) ::   zsizerat, zcompaph 
     67      REAL(wp) ::   zfactfe, zfactch, zprcaca, zfracal 
     68      REAL(wp) ::   ztortp , zrespp , zmortp  
     69      CHARACTER (len=25) ::   charout 
    7270      !!--------------------------------------------------------------------- 
    7371      ! 
    7472      IF( ln_timing )   CALL timing_start('p4z_nano') 
    7573      ! 
    76       prodcal(:,:,:) = 0.  !: calcite production variable set to zero 
     74      prodcal(:,:,:) = 0._wp   ! calcite production variable set to zero 
    7775      DO jk = 1, jpkm1 
    7876         DO jj = 1, jpj 
     
    139137      !! ** Method  : - ??? 
    140138      !!--------------------------------------------------------------------- 
    141       INTEGER  ::  ji, jj, jk 
    142       REAL(wp) ::  zfactfe,zfactsi,zfactch, zcompadi 
    143       REAL(wp) ::  zrespp2, ztortp2, zmortp2 
    144       REAL(wp) ::  zlim2, zlim1 
    145       CHARACTER (len=25) :: charout 
     139      INTEGER  ::   ji, jj, jk 
     140      REAL(wp) ::   zfactfe,zfactsi,zfactch, zcompadi 
     141      REAL(wp) ::   zrespp2, ztortp2, zmortp2 
     142      REAL(wp) ::   zlim2, zlim1 
     143      CHARACTER (len=25) ::   charout 
    146144      !!--------------------------------------------------------------------- 
    147145      ! 
    148146      IF( ln_timing )   CALL timing_start('p4z_diat') 
    149147      ! 
    150  
    151148      !    Aggregation term for diatoms is increased in case of nutrient 
    152149      !    stress as observed in reality. The stressed cells become more 
     
    196193      END DO 
    197194      ! 
    198       IF(ln_ctl)   THEN  ! print mean trends (used for debugging) 
     195      IF(ln_ctl) THEN      ! print mean trends (used for debugging) 
    199196         WRITE(charout, FMT="('diat')") 
    200197         CALL prt_ctl_trc_info(charout) 
     
    214211      !! 
    215212      !! ** Method  :   Read the nampismort namelist and check the parameters 
    216       !!      called at the first timestep 
     213      !!              called at the first timestep 
    217214      !! 
    218215      !! ** input   :   Namelist nampismort 
     
    224221      !!---------------------------------------------------------------------- 
    225222      ! 
     223      IF(lwp) THEN 
     224         WRITE(numout,*)  
     225         WRITE(numout,*) 'p4z_mort_init : Initialization of phytoplankton mortality parameters' 
     226         WRITE(numout,*) '~~~~~~~~~~~~~' 
     227      ENDIF 
     228      ! 
    226229      REWIND( numnatp_ref )              ! Namelist nampismort in reference namelist : Pisces phytoplankton 
    227230      READ  ( numnatp_ref, namp4zmort, IOSTAT = ios, ERR = 901) 
    228 901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namp4zmort in reference namelist', lwp ) 
    229       ! 
     231901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namp4zmort in reference namelist', lwp ) 
    230232      REWIND( numnatp_cfg )              ! Namelist nampismort in configuration namelist : Pisces phytoplankton 
    231233      READ  ( numnatp_cfg, namp4zmort, IOSTAT = ios, ERR = 902 ) 
    232 902   IF( ios >  0 ) CALL ctl_nam ( ios , 'namp4zmort in configuration namelist', lwp ) 
    233       IF(lwm) WRITE ( numonp, namp4zmort ) 
     234902   IF( ios >  0 )   CALL ctl_nam ( ios , 'namp4zmort in configuration namelist', lwp ) 
     235      IF(lwm) WRITE( numonp, namp4zmort ) 
    234236      ! 
    235237      IF(lwp) THEN                         ! control print 
    236          WRITE(numout,*) ' ' 
    237          WRITE(numout,*) ' Namelist parameters for phytoplankton mortality, namp4zmort' 
    238          WRITE(numout,*) ' ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~' 
    239          WRITE(numout,*) '    quadratic mortality of phytoplankton      wchl      =', wchl 
    240          WRITE(numout,*) '    maximum quadratic mortality of diatoms    wchld     =', wchld 
    241          WRITE(numout,*) '    maximum quadratic mortality of diatoms    wchldm    =', wchldm 
    242          WRITE(numout,*) '    phytoplankton mortality rate              mprat     =', mprat 
    243          WRITE(numout,*) '    Diatoms mortality rate                    mprat2    =', mprat2 
     238         WRITE(numout,*) '   Namelist : namp4zmort' 
     239         WRITE(numout,*) '      quadratic mortality of phytoplankton        wchl   =', wchl 
     240         WRITE(numout,*) '      maximum quadratic mortality of diatoms      wchld  =', wchld 
     241         WRITE(numout,*) '      maximum quadratic mortality of diatoms      wchldm =', wchldm 
     242         WRITE(numout,*) '      phytoplankton mortality rate                mprat  =', mprat 
     243         WRITE(numout,*) '      Diatoms mortality rate                      mprat2 =', mprat2 
    244244      ENDIF 
    245245      ! 
Note: See TracChangeset for help on using the changeset viewer.