Changeset 9169 for branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zprod.F90

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:

• branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zprod.F90 (modified) (6 diffs)

branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zprod.F90

@@ -8 +8 @@
    !!             3.4  !  2011-05  (O. Aumont, C. Ethe) New parameterization of light limitation
    !!----------------------------------------------------------------------
-   !!   p4z_prod       :   Compute the growth Rate of the two phytoplanktons groups
-   !!   p4z_prod_init  :   Initialization of the parameters for growth
-   !!   p4z_prod_alloc :   Allocate variables for growth
+   !!   p4z_prod       : Compute the growth Rate of the two phytoplanktons groups
+   !!   p4z_prod_init  : Initialization of the parameters for growth
+   !!   p4z_prod_alloc : Allocate variables for growth
    !!----------------------------------------------------------------------
-   USE oce_trc         !  shared variables between ocean and passive tracers
-   USE trc             !  passive tracers common variables 
-   USE sms_pisces      !  PISCES Source Minus Sink variables
-   USE p4zlim          !  Co-limitations of differents nutrients
-   USE prtctl_trc      !  print control for debugging
-   USE iom             !  I/O manager
+   USE oce_trc         ! shared variables between ocean and passive tracers
+   USE trc             ! passive tracers common variables 
+   USE sms_pisces      ! PISCES Source Minus Sink variables
+   USE p4zlim          ! Co-limitations of differents nutrients
+   USE prtctl_trc      ! print control for debugging
+   USE iom             ! I/O manager
 
    IMPLICIT NONE
@@ -26 +26 @@
    PUBLIC   p4z_prod_alloc
 
-   !! * Shared module variables
-   LOGICAL , PUBLIC ::  ln_newprod      !:
-   REAL(wp), PUBLIC ::  pislopen         !:
-   REAL(wp), PUBLIC ::  pisloped        !:
-   REAL(wp), PUBLIC ::  xadap           !:
-   REAL(wp), PUBLIC ::  excretn          !:
-   REAL(wp), PUBLIC ::  excretd         !:
-   REAL(wp), PUBLIC ::  bresp           !:
-   REAL(wp), PUBLIC ::  chlcnm          !:
-   REAL(wp), PUBLIC ::  chlcdm          !:
-   REAL(wp), PUBLIC ::  chlcmin         !:
-   REAL(wp), PUBLIC ::  fecnm           !:
-   REAL(wp), PUBLIC ::  fecdm           !:
-   REAL(wp), PUBLIC ::  grosip          !:
+   LOGICAL , PUBLIC ::   ln_newprod   !:
+   REAL(wp), PUBLIC ::   pislopen     !:
+   REAL(wp), PUBLIC ::   pisloped     !:
+   REAL(wp), PUBLIC ::   xadap        !:
+   REAL(wp), PUBLIC ::   excretn      !:
+   REAL(wp), PUBLIC ::   excretd      !:
+   REAL(wp), PUBLIC ::   bresp        !:
+   REAL(wp), PUBLIC ::   chlcnm       !:
+   REAL(wp), PUBLIC ::   chlcdm       !:
+   REAL(wp), PUBLIC ::   chlcmin      !:
+   REAL(wp), PUBLIC ::   fecnm        !:
+   REAL(wp), PUBLIC ::   fecdm        !:
+   REAL(wp), PUBLIC ::   grosip       !:
 
    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   prmax    !: optimal production = f(temperature)
@@ -45 +44 @@
    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   quotad   !: proxy of N quota in diatomee
    
-   REAL(wp) :: r1_rday                !: 1 / rday
-   REAL(wp) :: texcretn               !: 1 - excretn 
-   REAL(wp) :: texcretd               !: 1 - excretd        
+   REAL(wp) ::   r1_rday    ! 1 / rday
+   REAL(wp) ::   texcretn   ! 1 - excretn 
+   REAL(wp) ::   texcretd   ! 1 - excretd        
 
    !!----------------------------------------------------------------------
@@ -65 +64 @@
       !! ** Method  : - ???
       !!---------------------------------------------------------------------
-      INTEGER, INTENT(in) :: kt, knt
+      INTEGER, INTENT(in) ::   kt, knt   !
       !
       INTEGER  ::   ji, jj, jk
@@ -475 +474 @@
          CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm)
      ENDIF
-     !
-     IF( ln_timing )  CALL timing_stop('p4z_prod')
-     !
+      !
+      IF( ln_timing )  CALL timing_stop('p4z_prod')
+      !
    END SUBROUTINE p4z_prod
 
@@ -492 +491 @@
       !! ** input   :   Namelist nampisprod
       !!----------------------------------------------------------------------
-      INTEGER :: ios                 ! Local integer output status for namelist read
+      INTEGER ::   ios   ! Local integer
       !
       NAMELIST/namp4zprod/ pislopen, pisloped, xadap, ln_newprod, bresp, excretn, excretd,  &
          &                 chlcnm, chlcdm, chlcmin, fecnm, fecdm, grosip
       !!----------------------------------------------------------------------
-
+      !
+      IF(lwp) THEN                         ! control print
+         WRITE(numout,*)
+         WRITE(numout,*) 'p4z_prod_init : phytoplankton growth'
+         WRITE(numout,*) '~~~~~~~~~~~~~'
+      ENDIF
+      !
       REWIND( numnatp_ref )              ! Namelist nampisprod in reference namelist : Pisces phytoplankton production
       READ  ( numnatp_ref, namp4zprod, IOSTAT = ios, ERR = 901)
-901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namp4zprod in reference namelist', lwp )
-
+901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namp4zprod in reference namelist', lwp )
       REWIND( numnatp_cfg )              ! Namelist nampisprod in configuration namelist : Pisces phytoplankton production
       READ  ( numnatp_cfg, namp4zprod, IOSTAT = ios, ERR = 902 )
-902   IF( ios >  0 ) CALL ctl_nam ( ios , 'namp4zprod in configuration namelist', lwp )
-      IF(lwm) WRITE ( numonp, namp4zprod )
+902   IF( ios >  0 )   CALL ctl_nam ( ios , 'namp4zprod in configuration namelist', lwp )
+      IF(lwm) WRITE( numonp, namp4zprod )
 
       IF(lwp) THEN                         ! control print
-         WRITE(numout,*) ' '
-         WRITE(numout,*) ' Namelist parameters for phytoplankton growth, namp4zprod'
-         WRITE(numout,*) ' ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~'
-         WRITE(numout,*) '    Enable new parame. of production (T/F)   ln_newprod    =', ln_newprod
-         WRITE(numout,*) '    mean Si/C ratio                           grosip       =', grosip
-         WRITE(numout,*) '    P-I slope                                 pislopen     =', pislopen
-         WRITE(numout,*) '    Acclimation factor to low light           xadap        =', xadap
-         WRITE(numout,*) '    excretion ratio of nanophytoplankton      excretn      =', excretn
-         WRITE(numout,*) '    excretion ratio of diatoms                excretd      =', excretd
+         WRITE(numout,*) '   Namelist : namp4zprod'
+         WRITE(numout,*) '      Enable new parame. of production (T/F)   ln_newprod    =', ln_newprod
+         WRITE(numout,*) '      mean Si/C ratio                           grosip       =', grosip
+         WRITE(numout,*) '      P-I slope                                 pislopen     =', pislopen
+         WRITE(numout,*) '      Acclimation factor to low light           xadap        =', xadap
+         WRITE(numout,*) '      excretion ratio of nanophytoplankton      excretn      =', excretn
+         WRITE(numout,*) '      excretion ratio of diatoms                excretd      =', excretd
          IF( ln_newprod )  THEN
-            WRITE(numout,*) '    basal respiration in phytoplankton        bresp        =', bresp
-            WRITE(numout,*) '    Maximum Chl/C in phytoplankton            chlcmin      =', chlcmin
+            WRITE(numout,*) '      basal respiration in phytoplankton        bresp        =', bresp
+            WRITE(numout,*) '      Maximum Chl/C in phytoplankton            chlcmin      =', chlcmin
          ENDIF
-         WRITE(numout,*) '    P-I slope  for diatoms                    pisloped     =', pisloped
-         WRITE(numout,*) '    Minimum Chl/C in nanophytoplankton        chlcnm       =', chlcnm
-         WRITE(numout,*) '    Minimum Chl/C in diatoms                  chlcdm       =', chlcdm
-         WRITE(numout,*) '    Maximum Fe/C in nanophytoplankton         fecnm        =', fecnm
-         WRITE(numout,*) '    Minimum Fe/C in diatoms                   fecdm        =', fecdm
+         WRITE(numout,*) '      P-I slope  for diatoms                    pisloped     =', pisloped
+         WRITE(numout,*) '      Minimum Chl/C in nanophytoplankton        chlcnm       =', chlcnm
+         WRITE(numout,*) '      Minimum Chl/C in diatoms                  chlcdm       =', chlcdm
+         WRITE(numout,*) '      Maximum Fe/C in nanophytoplankton         fecnm        =', fecnm
+         WRITE(numout,*) '      Minimum Fe/C in diatoms                   fecdm        =', fecdm
       ENDIF
       !