Changeset 13233 for NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/P4Z/p4zsed.F90

Timestamp:

2020-07-02T20:34:16+02:00 (4 years ago)

Author:

aumont

Message:

update of the PISCES comments

File:

• NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/P4Z/p4zsed.F90 (modified) (17 diffs)

NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/P4Z/p4zsed.F90

@@ -35 +35 @@
    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:  ) :: sdenit     !: Nitrate reduction in the sediments
    REAL(wp) :: r1_rday                  !: inverse of rday
-   LOGICAL, SAVE :: lk_sed
+   LOGICAL, SAVE :: lk_sed              !: Explicit sediment module
 
    !!----------------------------------------------------------------------
@@ -48 +48 @@
       !!                     ***  ROUTINE p4z_sed  ***
       !!
-      !! ** Purpose :   Compute loss of biogenic matter in the sediments. This
+      !! ** Purpose : Compute the loss of biogenic matter in the sediments. This
       !!              is by no way a real sediment model. The loss is simply 
       !!              computed from metamodels.
@@ -55 +55 @@
       !!              N2 fixation is modeled using an implicit approach
       !!
-      !! ** Method  : - ???
+      !! ** Method  : - Fluxes with the sediments are mainly modeled using
+      !!                statiscal metamodels.
       !!---------------------------------------------------------------------
       !
@@ -81 +82 @@
       IF( kt == nittrc000 .AND. knt == 1 )   THEN
           r1_rday  = 1. / rday
+          ! Configuration with an active two-way sediment module 
           IF (ln_sediment .AND. ln_sed_2way) THEN
              lk_sed = .TRUE.
@@ -109 +111 @@
          !                                              
          ALLOCATE( zironice(jpi,jpj) )
-         !                                              
+
          ! Compute the iron flux between sea ice and sea water
+         ! Simple parameterization assuming a fixed constant concentration in
+         ! sea-ice (icefeinput)
+         ! ------------------------------------------------------------------
          DO jj = 1, jpj
             DO ji = 1, jpi
@@ -133 +138 @@
          !                                              
          ALLOCATE( zsidep(jpi,jpj), zpdep(jpi,jpj,jpk), zirondep(jpi,jpj,jpk) )
+
          ! Iron, P and Si deposition at the surface
          ! Iron flux at the surface due to dust deposition. Solubility can be 
          ! be variable if ln_solub is set to true. In that case, solubility 
          ! has to be provided in the specific input file (read in p4zsbc)
+         ! mfrac is the mean iron relative weight content of dust
          ! ------------------------------------------------------------------
          IF( ln_solub ) THEN
@@ -143 +150 @@
             zirondep(:,:,1) = dustsolub  * dust(:,:) * mfrac * rfact2 / e3t_n(:,:,1) / 55.85 + 3.e-10 * r1_ryyss 
          ENDIF
+
          ! Si and P flux at the surface due to dust deposition. The content 
          ! and the solubility are hard coded
@@ -148 +156 @@
          zsidep(:,:)   = 8.8 * 0.075 * dust(:,:) * mfrac * rfact2 / e3t_n(:,:,1) / 28.1 
          zpdep (:,:,1) = 0.1 * 0.021 * dust(:,:) * mfrac * rfact2 / e3t_n(:,:,1) / 31. / po4r 
+
          ! Iron solubilization of particles in the water column
          ! dust in kg/m2/s ---> 1/55.85 to put in mol/Fe ;  wdust in m/d
@@ -159 +168 @@
          DO jk = 2, jpkm1
             zirondep(:,:,jk) = dust(:,:) * mfrac * zwdust * rfact2 * EXP( -gdept_n(:,:,jk) / (250. * wdust) )
-!            zpdep   (:,:,jk) = zirondep(:,:,jk) * 0.023
             zpdep   (:,:,jk) = zirondep(:,:,jk) * 0.38 / po4r
          END DO
+
          ! Solubilization of particles in the water column (Si, P, Fe)
          tra(:,:,1,jpsil) = tra(:,:,1,jpsil) + zsidep  (:,:)
@@ -197 +206 @@
             ENDDO
          ENDDO
+
          ! When prognostic ligands are activated, ligands are supplied 
          ! to the ocean by rivers. We assume that the amount of ligands
@@ -262 +272 @@
          ! ------------------------------------------------------
          IF( ln_ironsed ) THEN
-                            tra(:,:,:,jpfer) = tra(:,:,:,jpfer) + ironsed(:,:,:) * rfact2
+            tra(:,:,:,jpfer) = tra(:,:,:,jpfer) + ironsed(:,:,:) * rfact2
             !
             IF( lk_iomput .AND. knt == nrdttrc .AND. iom_use( "Ironsed" ) )   &
@@ -434 +444 @@
                   xdiano3 = trb(ji,jj,jk,jpno3) / ( concnno3 + trb(ji,jj,jk,jpno3) ) * (1. - xdianh4)
                   zlim = ( 1.- xdiano3 - xdianh4 )
+                  ! Nitrogen fixation is almost fully halted when the N 
+                  ! limitation term (xdiano3+xdianh4) is > 0.9
                   IF( zlim <= 0.1 )   zlim = 0.01
                   zfact = zlim * rfact2
@@ -454 +466 @@
                   xdiano3 = trb(ji,jj,jk,jpno3) / ( concnno3 + trb(ji,jj,jk,jpno3) ) * (1. - xdianh4)
                   zlim = ( 1.- xdiano3 - xdianh4 )
+
+                  ! Nitrogen fixation is almost fully halted when the N 
+                  ! limitation term (xdiano3+xdianh4) is > 0.9
                   IF( zlim <= 0.1 )   zlim = 0.01
                   zfact = zlim * rfact2
@@ -474 +489 @@
                DO ji = 1, jpi
                   zfact = nitrpot(ji,jj,jk) * nitrfix
+                  ! 1/3 of the diazotrophs growth is supposed to be excreted
+                  ! as NH4. 1/3 as DOC and the rest is routed POC and GOC as 
+                  ! a result of mortality by predation. Completely adhoc param 
                   tra(ji,jj,jk,jpnh4) = tra(ji,jj,jk,jpnh4) + zfact / 3.0
                   tra(ji,jj,jk,jptal) = tra(ji,jj,jk,jptal) + rno3 * zfact / 3.0
@@ -481 +499 @@
                   tra(ji,jj,jk,jpgoc) = tra(ji,jj,jk,jpgoc) + zfact * 1.0 / 3.0 * 1.0 / 3.0
                   tra(ji,jj,jk,jpoxy) = tra(ji,jj,jk,jpoxy) + ( o2ut + o2nit ) * zfact * 2.0 / 3.0 + o2nit * zfact / 3.0
+                  ! Fe/c of diazotrophs is assumed to be 30umol Fe/mol C
                   tra(ji,jj,jk,jpfer) = tra(ji,jj,jk,jpfer) - 30E-6 * zfact * 1.0 / 3.0
                   tra(ji,jj,jk,jpsfe) = tra(ji,jj,jk,jpsfe) + 30E-6 * zfact * 1.0 / 3.0 * 2.0 / 3.0
@@ -494 +513 @@
                DO ji = 1, jpi
                   zfact = nitrpot(ji,jj,jk) * nitrfix
+                  ! 1/3 of the diazotrophs growth is supposed to be excreted
+                  ! as NH4. 1/3 as DOC and the rest is routed POC and GOC as 
+                  ! a result of mortality by predation. Completely adhoc param 
                   tra(ji,jj,jk,jpnh4) = tra(ji,jj,jk,jpnh4) + zfact / 3.0
                   tra(ji,jj,jk,jptal) = tra(ji,jj,jk,jptal) + rno3 * zfact / 3.0
+                  ! N/P ratio of diazotrophs is supposed to be 46
                   tra(ji,jj,jk,jppo4) = tra(ji,jj,jk,jppo4) - 16.0 / 46.0 * zfact * ( 1.0 - 1.0 / 3.0 ) &
                   &                     * ztrpo4(ji,jj,jk) / (ztrpo4(ji,jj,jk) + ztrdop(ji,jj,jk) + rtrn)
@@ -510 +533 @@
                   tra(ji,jj,jk,jpgop) = tra(ji,jj,jk,jpgop) + 16.0 / 46.0 * zfact * 1.0 / 3.0 * 1.0 /3.0
                   tra(ji,jj,jk,jpoxy) = tra(ji,jj,jk,jpoxy) + ( o2ut + o2nit ) * zfact * 2.0 / 3.0 + o2nit * zfact / 3.0
+                  ! Fe/c of diazotrophs is assumed to be 30umol Fe/mol C
                   tra(ji,jj,jk,jpfer) = tra(ji,jj,jk,jpfer) - 30E-6 * zfact * 1.0 / 3.0 
                   tra(ji,jj,jk,jpsfe) = tra(ji,jj,jk,jpsfe) + 30E-6 * zfact * 1.0 / 3.0 * 2.0 / 3.0