Changeset 9751 for NEMO/trunk/src/TOP/PISCES

Timestamp:

2018-06-06T15:00:58+02:00 (6 years ago)

Author:

nicolasmartin

Message:

Fix for OFF_PISCES ref cfg in DEBUG mode

Location:

NEMO/trunk/src/TOP/PISCES/P4Z

Files:

• p4zbio.F90 (modified) (1 diff)
• p4zlys.F90 (modified) (1 diff)
• p4zsms.F90 (modified) (1 diff)

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

@@ -99 +99 @@
       ENDIF
       !
-      CALL p4z_agg  ( kt, knt )     ! Aggregation of particles
+      CALL p4z_agg     ( kt, knt )     ! Aggregation of particles
       CALL p4z_rem     ( kt, knt )     ! remineralization terms of organic matter+scavenging of Fe
       CALL p4z_poc     ( kt, knt )     ! Remineralization of organic particles
-      IF( ln_ligand ) THEN
-        CALL p4z_ligand( kt, knt )
-      ENDIF
+      !
+      IF( ln_ligand )  &
+      & CALL p4z_ligand( kt, knt )
       !                                                             !
       IF(ln_ctl)   THEN  ! print mean trends (used for debugging)

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

@@ -123 +123 @@
 
       IF( lk_iomput .AND. knt == nrdttrc ) THEN
-         IF( iom_use( "PH"     ) ) CALL iom_put( "PH"    , -1. * LOG10( hi(:,:,:) )          * tmask(:,:,:) )
-         IF( iom_use( "CO3"    ) ) CALL iom_put( "CO3"   , zco3(:,:,:)    * 1.e+3            * tmask(:,:,:) )
-         IF( iom_use( "CO3sat" ) ) CALL iom_put( "CO3sat", zco3sat(:,:,:) * 1.e+3            * tmask(:,:,:) )
-         IF( iom_use( "DCAL"   ) ) CALL iom_put( "DCAL"  , zcaldiss(:,:,:) * 1.e+3 * rfact2r * tmask(:,:,:) )
+         IF( iom_use( "PH"     ) ) CALL iom_put( "PH"    , -1. * LOG10( MAX( hi(:,:,:), rtrn ) ) * tmask(:,:,:) )
+         IF( iom_use( "CO3"    ) ) CALL iom_put( "CO3"   , zco3(:,:,:)     * 1.e+3               * tmask(:,:,:) )
+         IF( iom_use( "CO3sat" ) ) CALL iom_put( "CO3sat", zco3sat(:,:,:)  * 1.e+3               * tmask(:,:,:) )
+         IF( iom_use( "DCAL"   ) ) CALL iom_put( "DCAL"  , zcaldiss(:,:,:) * 1.e+3 * rfact2r     * tmask(:,:,:) )
       ENDIF
       !

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

@@ -282 +282 @@
             CALL iom_get( numrtr, jpdom_autoglo, 'PH' , hi(:,:,:)  )
          ELSE
+            CALL p4z_che                              ! initialize the chemical constants
             CALL ahini_for_at(hi)
          ENDIF