Changeset 12537 for NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/trcini_pisces.F90

Timestamp:

2020-03-11T16:02:54+01:00 (4 years ago)

Author:

aumont

Message:

Comments in routines have been revised and significantly augmented

File:

• NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/trcini_pisces.F90 (modified) (7 diffs)

NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/trcini_pisces.F90

@@ -3 +3 @@
    !!                         ***  MODULE trcini_pisces  ***
    !! TOP :   initialisation of the PISCES biochemical model
+   !!         This module is for LOBSTER, PISCES and PISCES-QUOTA
    !!======================================================================
    !! History :    -   !  1988-07  (E. Maier-Reiner) Original code
@@ -37 +38 @@
       !!
       !! ** Purpose :   Initialisation of the PISCES biochemical model
-      !!----------------------------------------------------------------------
-      !
+      !!                Allocation of the dynamic arrays
+      !!----------------------------------------------------------------------
+      !
+      ! Read the PISCES namelist
       CALL trc_nam_pisces
       !
@@ -54 +57 @@
       !! ** Purpose :   Initialisation of the PISCES biochemical model
       !!----------------------------------------------------------------------
-      USE p4zsms          ! Main P4Z routine
+      USE p4zsms          !  Main P4Z routine
       USE p4zche          !  Chemical model
       USE p4zsink         !  vertical flux of particulate matter due to sinking
@@ -111 +114 @@
       ierr = ierr +  p4z_lim_alloc()
       IF( ln_p4z ) THEN
+         ! PISCES part
          ierr = ierr +  p4z_prod_alloc()
          ierr = ierr +  p4z_meso_alloc()
       ELSE
+         ! PISCES-QUOTA part
          ierr = ierr +  p5z_lim_alloc()
          ierr = ierr +  p5z_prod_alloc()
@@ -128 +133 @@
 
       ! assign an index in trc arrays for each prognostic variables
+      ! This is based on the information read in the namelist_top
       DO jn = 1, jptra
         cltra = ctrcnm(jn) 
@@ -177 +183 @@
       ! Set biological ratios
       ! ---------------------
-      rno3    =  16._wp / 122._wp
-      po4r    =   1._wp / 122._wp
-      o2nit   =  32._wp / 122._wp
-      o2ut    = 133._wp / 122._wp
-      rdenit  =  ( ( o2ut + o2nit ) * 0.80 - rno3 - rno3 * 0.60 ) / rno3
-      rdenita =   3._wp /  5._wp
+      rno3    =  16._wp / 122._wp   ! C/N
+      po4r    =   1._wp / 122._wp   ! C/P
+      o2nit   =  32._wp / 122._wp   ! O2/C for nitrification
+      o2ut    = 133._wp / 122._wp   ! O2/C for ammonification
+      rdenit  =  ( ( o2ut + o2nit ) * 0.80 - rno3 - rno3 * 0.60 ) / rno3  ! Denitrification
+      rdenita =   3._wp /  5._wp    ! Denitrification
       IF( ln_p5z ) THEN
-         no3rat3 = no3rat3 / rno3
-         po4rat3 = po4rat3 / po4r
+         no3rat3 = no3rat3 / rno3   ! C/N ratio in zooplankton
+         po4rat3 = po4rat3 / po4r   ! C/P ratio in zooplankton
       ENDIF
 
@@ -250 +256 @@
       CALL p4z_opt_init          !  Optic: PAR in the water column
       IF( ln_p4z ) THEN
+         ! PISCES part
          CALL p4z_lim_init       !  co-limitations by the various nutrients
          CALL p4z_prod_init      !  phytoplankton growth rate over the global ocean.
       ELSE
+         ! PISCES-QUOTA part
          CALL p5z_lim_init       !  co-limitations by the various nutrients
          CALL p5z_prod_init      !  phytoplankton growth rate over the global ocean.