Changeset 10222 for NEMO/trunk/src/TOP/PISCES/P4Z/p4zsms.F90

Timestamp:

2018-10-25T11:42:23+02:00 (6 years ago)

Author:

aumont

Message:

update of the sediment module + enhancement, bug correction in PISCES generating O2 negative values

File:

• NEMO/trunk/src/TOP/PISCES/P4Z/p4zsms.F90 (modified) (3 diffs)

NEMO/trunk/src/TOP/PISCES/P4Z/p4zsms.F90

@@ -101 +101 @@
       ENDIF
       !
+      IF( ll_sbc ) CALL p4z_sbc( kt )   ! external sources of nutrients 
+      !
+#if ! defined key_sed_off
       CALL p4z_che              ! computation of chemical constants
       CALL p4z_int( kt )        ! computation of various rates for biogeochemistry
       !
-      IF( ll_sbc ) CALL p4z_sbc( kt )   ! external sources of nutrients 
-
       DO jnt = 1, nrdttrc          ! Potential time splitting if requested
          !
@@ -149 +150 @@
          END DO
       END IF
-      !
-      IF( lk_sed ) THEN 
+#endif
+      !
+      IF( ln_sediment ) THEN 
          !
          CALL sed_model( kt )     !  Main program of Sediment model
+         !
+         IF( ln_top_euler ) THEN
+            DO jn = jp_pcs0, jp_pcs1
+               trn(:,:,:,jn) = trb(:,:,:,jn)
+            END DO
+         ENDIF
          !
       ENDIF
@@ -352 +360 @@
       IF(lwp)  WRITE(numout,*)
 
-      IF( cn_cfg == "ORCA" .OR. cn_cfg == "orca" ) THEN
+      IF( cn_cfg == "ORCA" .OR. cn_cfg == "orca") THEN
          IF( .NOT. lk_c1d ) THEN      ! ORCA configuration (not 1D) !
             !                                                ! --------------------------- !