Changeset 15803 for NEMO/branches/UKMO/NEMO_4.0.4_EAP_rheology_fix/src/ICE/icedyn_rhg_evp.F90

Timestamp:

2022-04-26T14:32:50+02:00 (2 years ago)

Author:

annkeen

Message:

Change default ice_cvg.nc output to be mean rather than maximum

File:

• NEMO/branches/UKMO/NEMO_4.0.4_EAP_rheology_fix/src/ICE/icedyn_rhg_evp.F90 (modified) (2 diffs)

NEMO/branches/UKMO/NEMO_4.0.4_EAP_rheology_fix/src/ICE/icedyn_rhg_evp.F90

@@ -985 +985 @@
       REAL(wp)          ::   zresm           ! local real 
       CHARACTER(len=20) ::   clname
-      REAL(wp), DIMENSION(jpi,jpj) ::   zres           ! check convergence
+      LOGICAL           ::   ll_maxcvg
+      REAL(wp), DIMENSION(jpi,jpj) ::   zres1           ! check convergence
+      REAL(wp), DIMENSION(jpi,jpj) ::   zres2           ! check convergence
+      REAL(wp), DIMENSION(2)         ::   ztmp
       !!----------------------------------------------------------------------
+      ll_maxcvg = .FALSE.
 
       ! create file
@@ -1015 +1019 @@
          zresm = 0._wp
       ELSE
-         DO jj = 1, jpj
-            DO ji = 1, jpi
-               zres(ji,jj) = MAX( ABS( pu(ji,jj) - pub(ji,jj) ) * umask(ji,jj,1), &
+         zresm = 0._wp
+         IF( ll_maxcvg ) THEN   ! error max over the domain
+           DO jj = 1, jpj
+              DO ji = 1, jpi
+                 zres1(ji,jj) = MAX( ABS( pu(ji,jj) - pub(ji,jj) ) * umask(ji,jj,1), &
                   &               ABS( pv(ji,jj) - pvb(ji,jj) ) * vmask(ji,jj,1) ) * zmsk15(ji,jj)
-            END DO
-         END DO
-         zresm = MAXVAL( zres )
-         CALL mpp_max( 'icedyn_rhg_evp', zresm )   ! max over the global domain
+                 zres2(ji,jj) = 0._wp
+              END DO
+           END DO
+           zresm = MAXVAL( zres1 )
+           CALL mpp_max( 'icedyn_rhg_evp', zresm )   ! max over the global domain
+         ELSE                    ! error averaged over the domain
+           DO jj = 1, jpj
+              DO ji = 1, jpi
+                 zres1(ji,jj) = MAX( ABS( pu(ji,jj) - pub(ji,jj) ) * umask(ji,jj,1), &
+                  &                 ABS( pv(ji,jj) - pvb(ji,jj) ) * vmask(ji,jj,1) ) * zmsk15(ji,jj)
+                 zres2(ji,jj) = zmsk15(ji,jj)
+              END DO
+           END DO
+           ztmp(1) = glob_sum( 'icedyn_rhg_evp', zres1 )
+           ztmp(2) = glob_sum( 'icedyn_rhg_evp', zres2 )
+           IF( ztmp(2) /= 0._wp )   zresm = ztmp(1) / ztmp(2)
+         ENDIF
       ENDIF