Changeset 9073 for branches/NERC/dev_r5518_GO6_Carb_Fail_from_GO6_8356/NEMOGCM/NEMO/OPA_SRC/stpctl.F90

Timestamp:

2017-12-15T12:48:47+01:00 (6 years ago)

Author:

jpalmier

Message:

add all micro boil checks and securities

File:

• branches/NERC/dev_r5518_GO6_Carb_Fail_from_GO6_8356/NEMOGCM/NEMO/OPA_SRC/stpctl.F90 (modified) (2 diffs)

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

@@ -58 +58 @@
       INTEGER  ::   ii, ij, ik              ! temporary integers
       REAL(wp) ::   zumax, zsmin, zssh2     ! temporary scalars
+      REAL(wp) ::   ztmax, ztmin            ! temporary scalars
       INTEGER, DIMENSION(3) ::   ilocu      ! 
       INTEGER, DIMENSION(2) ::   ilocs      ! 
@@ -148 +149 @@
 9500  FORMAT (' kt=',i6,' min SSS: ',1pg11.4,', i j: ',2i5)
 
+! ====================================================================================================
+! ====================================================================================================
+      !                                              !AXY (25/10/17)
+      !                                              !* Test max/min limits of temperature
+      !                                              !  ----------------------------------
+      ztmax =  -5.e0  ! arbitrary  low maximum value
+      ztmin = 100.e0  ! arbitrary high minimum value
+      DO jj = 2, jpjm1
+         DO ji = 1, jpi
+            IF( tmask(ji,jj,1) == 1) THEN
+               ztmax = MAX(ztmax,tsn(ji,jj,1,jp_tem))     ! find local maximum
+               ztmin = MIN(ztmin,tsn(ji,jj,1,jp_tem))     ! find local minimum
+            ENDIF
+         END DO
+      END DO
+      IF( lk_mpp )   CALL mpp_max( ztmax )                ! max over the global domain
+      IF( lk_mpp )   CALL mpp_min( ztmin )                ! min over the global domain
+      !
+      IF( ztmax > 40.) THEN  ! we've got a problem
+         IF (lk_mpp) THEN
+            CALL mpp_maxloc ( tsn(:,:,1,jp_tem),tmask(:,:,1), ztmax, ii,ij )
+         ELSE
+            ilocs = MAXLOC( tsn(:,:,1,jp_tem), mask = tmask(:,:,1) == 1.e0 )
+            ii = ilocs(1) + nimpp - 1
+            ij = ilocs(2) + njmpp - 1
+         ENDIF
+         !
+         IF(lwp) THEN
+            WRITE(numout,*) 'stp_ctl:tracer anomaly: *****    WARNING     *****'
+            WRITE(numout,*) 'stp_ctl:tracer anomaly: sea surface temperature > 40C'
+            WRITE(numout,9600) kt, ztmax, ii, ij
+            WRITE(numout,*) 'stp_ctl:tracer anomaly: ***** END OF WARNING *****'
+         ENDIF
+      ENDIF
+      !
+      IF( ztmin < -3.) THEN  ! we've got a problem
+         IF (lk_mpp) THEN
+            CALL mpp_minloc ( tsn(:,:,1,jp_tem),tmask(:,:,1), ztmin, ii,ij )
+         ELSE
+            ilocs = MINLOC( tsn(:,:,1,jp_tem), mask = tmask(:,:,1) == 1.e0 )
+            ii = ilocs(1) + nimpp - 1
+            ij = ilocs(2) + njmpp - 1
+         ENDIF
+         !
+         IF(lwp) THEN
+            WRITE(numout,*) 'stp_ctl:tracer anomaly: *****    WARNING     *****'
+            WRITE(numout,*) 'stp_ctl:tracer anomaly: sea surface temperature < -3C'
+            WRITE(numout,9700) kt, ztmin, ii, ij
+            WRITE(numout,*) 'stp_ctl:tracer anomaly: ***** END OF WARNING *****'
+         ENDIF
+      ENDIF
+9600  FORMAT ('stp_ctl:tracer anomaly: kt=',i6,' max SST: ',f16.10,', i j: ',2i5)
+9700  FORMAT ('stp_ctl:tracer anomaly: kt=',i6,' min SST: ',f16.10,', i j: ',2i5)
+! ====================================================================================================
+! ====================================================================================================
       
       IF( lk_c1d )  RETURN          ! No log file in case of 1D vertical configuration