Changeset 11143 for NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser/src/OCE/stpctl.F90

Timestamp:

2019-06-19T18:15:12+02:00 (5 years ago)

Author:

agn

Message:

version w/o FK or local changes for NEMO repo

File:

• NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser/src/OCE/stpctl.F90 (modified) (5 diffs)

NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser/src/OCE/stpctl.F90

@@ -66 +66 @@
       INTEGER, DIMENSION(3)  ::   iu, is1, is2        ! min/max loc indices
       REAL(wp)               ::   zzz                 ! local real 
-      REAL(wp), DIMENSION(9) ::   zmax
-      LOGICAL                ::   ll_wrtstp, ll_colruns, ll_wrtruns
+      REAL(wp), DIMENSION(10) ::   zmax
+      LOGICAL                ::   ll_wrtstp, ll_colruns, ll_wrtruns, ll_isnan
       CHARACTER(len=20) :: clname
       !!----------------------------------------------------------------------
@@ -109 +109 @@
       ENDIF
       !
+      ll_isnan = ANY(ISNAN(tsn)) .OR. ANY(ISNAN(un))
+      IF (ll_isnan) nstop = nstop + 1
       !                                   !==  test of extrema  ==!
       IF( ll_wd ) THEN
@@ -124 +126 @@
          zmax(8) = MAXVAL(  ABS( wi(:,:,:) ) , mask = wmask(:,:,:) == 1._wp ) ! implicit vertical vel. max
          zmax(9) = MAXVAL(   Cu_adv(:,:,:)   , mask = tmask(:,:,:) == 1._wp ) !       cell Courant no. max
+      ENDIF
+      IF (ll_isnan) THEN
+         zmax(10) = 1._wp                                           ! stop indicator
+      ELSE
+         zmax(10) = 0._wp
       ENDIF
       !
@@ -147 +154 @@
       END IF
       !                                   !==  error handling  ==!
-      IF( ( ln_ctl .OR. lsomeoce ) .AND. (   &             ! have use mpp_max (because ln_ctl=.T.) or contains some ocean points
+      IF( ( (ln_ctl .OR. sn_cfctl%l_runstat) .OR. lsomeoce ) .AND. (   &             ! have use mpp_max (because ln_ctl=.T.) or contains some ocean points
          &  zmax(1) >   20._wp .OR.   &                    ! too large sea surface height ( > 20 m )
          &  zmax(2) >   10._wp .OR.   &                    ! too large velocity ( > 10 m/s)
@@ -153 +160 @@
          &  zmax(4) >= 100._wp .OR.   &                    ! too large sea surface salinity ( > 100 )
          &  zmax(4) <    0._wp .OR.   &                    ! too large sea surface salinity (keep this line for sea-ice)
-         &  ISNAN( zmax(1) + zmax(2) + zmax(3) ) ) ) THEN   ! NaN encounter in the tests
-         IF( lk_mpp .AND. ln_ctl ) THEN
+         &  zmax(10) >   0._wp  ) ) THEN   ! NaN encounter in the tests
+         ! &  ISNAN( zmax(1) + zmax(2) + zmax(3) ) ) ) THEN   ! NaN encounter in the tests
+         IF( lk_mpp .AND. (ln_ctl .OR. sn_cfctl%l_runstat)) THEN
             CALL mpp_maxloc( 'stpctl', ABS(sshn)        , ssmask(:,:)  , zzz, ih  )
             CALL mpp_maxloc( 'stpctl', ABS(un)          , umask (:,:,:), zzz, iu  )