Changeset 13233 for NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/P4Z/p4zsms.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/p4zsms.F90 (modified) (4 diffs)

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

@@ -74 +74 @@
         IF( .NOT. ln_rsttr ) THEN
             CALL p4z_che            ! initialize the chemical constants
-            CALL ahini_for_at(hi)   !  set PH at kt=nit000
+            CALL ahini_for_at(hi)   ! set PH at kt=nit000
             t_oce_co2_flx_cum = 0._wp
         ELSE
@@ -84 +84 @@
       IF( ln_pisdmp .AND. MOD( kt - nn_dttrc, nn_pisdmp ) == 0 )   CALL p4z_dmp( kt )      ! Relaxation of some tracers
       !
-      rfact = r2dttrc
+      rfact = r2dttrc  ! time step of PISCES
       !
       IF( ( ln_top_euler .AND. kt == nittrc000 )  .OR. ( .NOT.ln_top_euler .AND. kt <= nittrc000 + nn_dttrc ) ) THEN
-         rfactr  = 1. / rfact
-         rfact2  = rfact / REAL( nrdttrc, wp )
-         rfact2r = 1. / rfact2
-         xstep = rfact2 / rday         ! Time step duration for biology
+         rfactr  = 1. / rfact  ! inverse of the time step
+         rfact2  = rfact / REAL( nrdttrc, wp )  ! time step of the biological SMS
+         rfact2r = 1. / rfact2  ! Inverse of the biological time step
+         xstep = rfact2 / rday         ! Time step duration for biology relative to a day
          IF(lwp) WRITE(numout,*) 
          IF(lwp) WRITE(numout,*) '    Passive Tracer  time step    rfact  = ', rfact, ' rdt = ', rdt
@@ -147 +147 @@
             tra(:,:,:,jn) = 0._wp
          END DO
-         !
+         ! Euler-forward temporal scheme
          IF( ln_top_euler ) THEN
             DO jn = jp_pcs0, jp_pcs1
@@ -167 +167 @@
          !
          CALL sed_model( kt )     !  Main program of Sediment model
-         !
+         ! Eulor forward temporal scheme
          IF( ln_top_euler ) THEN
             DO jn = jp_pcs0, jp_pcs1