Changeset 9130 for branches/NERC

Timestamp:

2017-12-19T14:15:21+01:00 (7 years ago)

Author:

jpalmier

Message:

JPALM -- microboil-carb fail Check and report

Location:

branches/NERC/dev_r5518_GO6_Carb_Fail_from_GO6_8356/NEMOGCM/NEMO

Files:

• OPA_SRC/stpctl.F90 (modified) (1 diff)
• TOP_SRC/MEDUSA/trcbio_medusa.F90 (modified) (6 diffs)

branches/NERC/dev_r5518_GO6_Carb_Fail_from_GO6_8356/NEMOGCM/NEMO/OPA_SRC/stpctl.F90

@@ -157 +157 @@
       ztmin = 100.e0  ! arbitrary high minimum value
       DO jj = 2, jpjm1
-         DO ji = 1, jpi
+         DO ji = 2, jpim1
             IF( tmask(ji,jj,1) == 1) THEN
                ztmax = MAX(ztmax,tsn(ji,jj,1,jp_tem))     ! find local maximum

branches/NERC/dev_r5518_GO6_Carb_Fail_from_GO6_8356/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90

@@ -821 +821 @@
       !!==========================
       !! DIC Check
-      zmax =  -5.e0  ! arbitrary  low maximum value
+      zmax =  -5.0  ! arbitrary  low maximum value
       zmin =  4.0E4  ! arbitrary high minimum value
       DO jj = 2, jpjm1
-         DO ji = 1, jpi
+         DO ji = 2,jpim1
             IF( tmask(ji,jj,1) == 1) THEN
                zmax = MAX(zmax,zdic(ji,jj))     ! find local maximum
@@ -834 +834 @@
       IF( lk_mpp )   CALL mpp_min( zmin )       ! min over the global domain
       !
-      IF( zmax > 4.0E3) THEN  ! we've got a problem
+      IF( zmax .GT. 4.0E3) THEN  ! we've got a problem
          IF (lk_mpp) THEN
             CALL mpp_maxloc ( zdic(:,:),tmask(:,:,1), zmax, ii,ij )
          ELSE
-            ilocs = MAXLOC( zdic(:,:), mask = tmask(:,:,1) == 1.e0 )
+            ilocs = MAXLOC( zdic(:,:), mask = tmask(:,:,1) == 1. )
             ii = ilocs(1) + nimpp - 1
             ij = ilocs(2) + njmpp - 1
@@ -851 +851 @@
       ENDIF
       !
-      IF( zmin <= 0.0) THEN  ! we've got a problem
+      IF( zmin .LE. 0.0) THEN  ! we've got a problem
          IF (lk_mpp) THEN
             CALL mpp_minloc ( zdic(:,:),tmask(:,:,1), zmin, ii,ij )
          ELSE
-            ilocs = MINLOC( zdic(:,:), mask = tmask(:,:,1) == 1.e0 )
+            ilocs = MINLOC( zdic(:,:), mask = tmask(:,:,1) == 1. )
             ii = ilocs(1) + nimpp - 1
             ij = ilocs(2) + njmpp - 1
@@ -870 +870 @@
       !!==========================
       !! ALKALINITY Check
-      zmax =  -5.e0  ! arbitrary  low maximum value
+      zmax =  -5.0  ! arbitrary  low maximum value
       zmin =  4.0E4  ! arbitrary high minimum value
       DO jj = 2, jpjm1
-         DO ji = 1, jpi
+         DO ji = 2,jpim1
             IF( tmask(ji,jj,1) == 1) THEN
                zmax = MAX(zmax,zalk(ji,jj))     ! find local maximum
@@ -883 +883 @@
       IF( lk_mpp )   CALL mpp_min( zmin )       ! min over the global domain
       !
-      IF( zmax > 4.0E3) THEN  ! we've got a problem
+      IF( zmax .GT. 4.0E3) THEN  ! we've got a problem
          IF (lk_mpp) THEN
             CALL mpp_maxloc ( zalk(:,:),tmask(:,:,1), zmax, ii,ij )
          ELSE
-            ilocs = MAXLOC( zalk(:,:), mask = tmask(:,:,1) == 1.e0 )
+            ilocs = MAXLOC( zalk(:,:), mask = tmask(:,:,1) == 1. )
             ii = ilocs(1) + nimpp - 1
             ij = ilocs(2) + njmpp - 1
@@ -900 +900 @@
       ENDIF
       !
-      IF( zmin <= 0.0) THEN  ! we've got a problem
+      IF( zmin .LE. 0.0) THEN  ! we've got a problem
          IF (lk_mpp) THEN
             CALL mpp_minloc ( zalk(:,:),tmask(:,:,1), zmin, ii,ij )
          ELSE
-            ilocs = MINLOC( zalk(:,:), mask = tmask(:,:,1) == 1.e0 )
+            ilocs = MINLOC( zalk(:,:), mask = tmask(:,:,1) == 1. )
             ii = ilocs(1) + nimpp - 1
             ij = ilocs(2) + njmpp - 1