Changeset 9863 for NEMO/branches/2018/dev_r9838_ENHANCE04_MLF/src/TOP/PISCES/P2Z/p2zexp.F90

Timestamp:

2018-06-30T12:51:02+02:00 (6 years ago)

Author:

gm

Message:

#1911 (ENHANCE-04): simplified implementation of the Euler stepping at nit000

File:

• NEMO/branches/2018/dev_r9838_ENHANCE04_MLF/src/TOP/PISCES/P2Z/p2zexp.F90 (modified) (6 diffs)

NEMO/branches/2018/dev_r9838_ENHANCE04_MLF/src/TOP/PISCES/P2Z/p2zexp.F90

@@ -4 +4 @@
    !! TOP :   LOBSTER Compute loss of organic matter in the sediments
    !!======================================================================
-   !! History :    -   !  1999    (O. Aumont, C. Le Quere)  original code
-   !!              -   !  2001-05 (O. Aumont, E. Kestenare) add sediment computations
-   !!             1.0  !  2005-06 (A.-S. Kremeur) new temporal integration for sedpoc
-   !!             2.0  !  2007-12  (C. Deltel, G. Madec)  F90
-   !!             3.5  !  2012-03  (C. Ethe)  Merge PISCES-LOBSTER
-   !!----------------------------------------------------------------------
-   !!   p2z_exp        :  Compute loss of organic matter in the sediments
-   !!----------------------------------------------------------------------
-   USE oce_trc         !
-   USE trc
-   USE sms_pisces
-   USE p2zsed
-   USE lbclnk
-   USE prtctl_trc      ! Print control for debbuging
-   USE trd_oce
-   USE trdtrc
-   USE iom
+   !! History :   -   !  1999    (O. Aumont, C. Le Quere)  original code
+   !!             -   !  2001-05 (O. Aumont, E. Kestenare) add sediment computations
+   !!            1.0  !  2005-06 (A.-S. Kremeur) new temporal integration for sedpoc
+   !!            2.0  !  2007-12  (C. Deltel, G. Madec)  F90
+   !!            3.5  !  2012-03  (C. Ethe)  Merge PISCES-LOBSTER
+   !!----------------------------------------------------------------------
+   !!   p2z_exp       :  Compute loss of organic matter in the sediments
+   !!----------------------------------------------------------------------
+   USE oce_trc        !
+   USE trc            !
+   USE sms_pisces     !
+   USE p2zsed         !
+   USE lbclnk         !
+   USE prtctl_trc     ! Print control for debbuging
+   USE trd_oce        !
+   USE trdtrc         !
+   USE iom            !
 
    IMPLICIT NONE
@@ -30 +30 @@
 
    !
-   REAL(wp), ALLOCATABLE, DIMENSION(:,:)   ::   dminl     !: fraction of sinking POC released in sediments
-   REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) ::   dmin3     !: fraction of sinking POC released at each level
-   REAL(wp), ALLOCATABLE, DIMENSION(:,:)   ::   sedpocb   !: mass of POC in sediments
-   REAL(wp), ALLOCATABLE, DIMENSION(:,:)   ::   sedpocn   !: mass of POC in sediments
-   REAL(wp), ALLOCATABLE, DIMENSION(:,:)   ::   cmask     !: Coastal mask area
-   REAL(wp)                                ::   areacot   !: surface coastal area
+   REAL(wp), ALLOCATABLE, DIMENSION(:,:)   ::   dminl     ! fraction of sinking POC released in sediments
+   REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) ::   dmin3     ! fraction of sinking POC released at each level
+   REAL(wp), ALLOCATABLE, DIMENSION(:,:)   ::   sedpocb   ! mass of POC in sediments
+   REAL(wp), ALLOCATABLE, DIMENSION(:,:)   ::   sedpocn   ! mass of POC in sediments
+   REAL(wp), ALLOCATABLE, DIMENSION(:,:)   ::   cmask     ! Coastal mask area
+   REAL(wp)                                ::   areacot   ! surface coastal area
 
    !! * Substitutions
@@ -59 +59 @@
       !!              COLUMN BELOW THE SURFACE LAYER.
       !!---------------------------------------------------------------------
-      !!
-      INTEGER, INTENT( in ) ::   kt      ! ocean time-step index      
+      INTEGER, INTENT( in ) ::   kt   ! ocean time-step index      
       !!
       INTEGER  ::   ji, jj, jk, jl, ikt
       REAL(wp) ::   zgeolpoc, zfact, zwork, ze3t, zsedpocd, zmaskt
       REAL(wp), DIMENSION(jpi,jpj)   ::  zsedpoca
-      CHARACTER (len=25) :: charout
+      CHARACTER (len=25) ::   charout
       !!---------------------------------------------------------------------
       !
@@ -72 +71 @@
       IF( kt == nittrc000 )   CALL p2z_exp_init
 
-      zsedpoca(:,:) = 0.
-
-
-      ! VERTICAL DISTRIBUTION OF NEWLY PRODUCED BIOGENIC
-      ! POC IN THE WATER COLUMN
+      zsedpoca(:,:) = 0._wp
+
+
+      ! VERTICAL DISTRIBUTION OF NEWLY PRODUCED BIOGENIC POC IN THE WATER COLUMN
       ! (PARTS OF NEWLY FORMED MATTER REMAINING IN THE DIFFERENT
       ! LAYERS IS DETERMINED BY DMIN3 DEFINED IN sms_p2z.F90
@@ -93 +91 @@
    
 
-      zgeolpoc = 0.e0         !     Initialization
-      ! Release of nutrients from the "simple" sediment
-      DO jj = 2, jpjm1
+      zgeolpoc = 0._wp        !     Initialization
+      DO jj = 2, jpjm1           ! Release of nutrients from the "simple" sediment
          DO ji = fs_2, fs_jpim1
             ikt = mbkt(ji,jj) 
@@ -121 +118 @@
       ! Time filter and swap of arrays
       ! ------------------------------
-      IF( neuler == 0 .AND. kt == nittrc000 ) THEN        ! Euler time-stepping at first time-step
-        !                                             ! (only swap)
+      IF( l_1st_euler ) THEN        ! Euler time-stepping at first time-step   (only swap)
         sedpocn(:,:) = zsedpoca(:,:)
         !                                              
-      ELSE
+      ELSE                          ! Leap-Frog + Asselin filter
         !
         DO jj = 1, jpj
            DO ji = 1, jpi
-              zsedpocd = zsedpoca(ji,jj) - 2. * sedpocn(ji,jj) + sedpocb(ji,jj)      ! time laplacian on tracers
+              zsedpocd = zsedpoca(ji,jj) - 2. * sedpocn(ji,jj) + sedpocb(ji,jj)     ! time laplacian on tracers
               sedpocb(ji,jj) = sedpocn(ji,jj) + atfp * zsedpocd                     ! sedpocb <-- filtered sedpocn
-              sedpocn(ji,jj) = zsedpoca(ji,jj)                                       ! sedpocn <-- sedpoca
+              sedpocn(ji,jj) = zsedpoca(ji,jj)                                      ! sedpocn <-- sedpoca
            END DO
         END DO