Changeset 11965 for branches

Timestamp:

2019-11-26T12:39:45+01:00 (4 years ago)

Author:

frrh

Message:

Save updates for special viscosity adjustments in S polar region.

File:

• branches/UKMO/dev_r5518_GO6_eORCA1_visc/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn_c3d.h90 (modified) (6 diffs)

branches/UKMO/dev_r5518_GO6_eORCA1_visc/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn_c3d.h90

@@ -205 +205 @@
       INTEGER ::   ifreq, il1, il2, ij, ii
       REAL(wp) ::   zahmeq, zcoff, zcoft, zmsk   ! local scalars
-      REAL(wp) ::   zemax , zemin, zeref, zahmm
+      REAL(wp) ::   zemax , zemin, zetmax, zefmax, zeref, zahmm, zemax75
       CHARACTER (len=15) ::   clexp
       INTEGER , POINTER, DIMENSION(:,:)  :: icof
+      INTEGER , POINTER, DIMENSION(:,:)  :: imsk
       REAL(wp), POINTER, DIMENSION(:  )  :: zcoef   
       REAL(wp), POINTER, DIMENSION(:,:)  :: zahm0
@@ -216 +217 @@
       !
       CALL wrk_alloc( jpi   , jpj   , icof  )
+      CALL wrk_alloc( jpi   , jpj   , imsk  )
       CALL wrk_alloc( jpk   ,         zcoef )
       CALL wrk_alloc( jpi   , jpj   , zahm0 )
@@ -253 +255 @@
          zahmeq = 5.0 * aht0
          zahmm  = min( 160000.0, ahm0)
-         zemax = MAXVAL ( e1t(:,:) * e2t(:,:), tmask(:,:,1) .GE. 0.5 )
-         zemin = MINVAL ( e1t(:,:) * e2t(:,:), tmask(:,:,1) .GE. 0.5 )
+         zemax = MAXVAL ( e1t(:,:) * e2t(:,:), ssmask(:,:) .GE. 0.5 )
+         zemin = MINVAL ( e1t(:,:) * e2t(:,:), ssmask(:,:) .GE. 0.5 )
+
          tempmask(:,:) = .FALSE.
          ! Pre calculate mask for zeref since embedding the following
          ! term in the MAXVAL operation offends the Cray compiler for no 
          ! justifiable reason under certain conditions. 
-         tempmask(:,:) = (tmask(:,:,1) .GE. 0.5) .AND. (ABS(gphit(:,:)) .GT. 50.)
+         tempmask(:,:) = (ssmask(:,:) .GE. 0.5) .AND. (ABS(gphit(:,:)) .GT. 50.)
          zeref = MAXVAL ( e1t(:,:) * e2t(:,:), tempmask(:,:) )
  
@@ -321 +324 @@
       END DO
       ! f-point
-      icof(:,:) = icof(:,:) * tmask(:,:,1)
+      icof(:,:) = icof(:,:) * ssmask(:,:)
       DO jj = 1, jpjm1
          DO ji = 1, jpim1   ! NO vector opt.
-            zmsk = tmask(ji,jj+1,1) + tmask(ji+1,jj+1,1) + tmask(ji,jj,1) + tmask(ji,jj+1,1)
+            zmsk = ssmask(ji,jj+1) + ssmask(ji+1,jj+1) + ssmask(ji,jj) + ssmask(ji,jj+1)
             IF( zmsk == 0. ) THEN
                zcoff = 1.
@@ -387 +390 @@
       CALL lbc_lnk( ahm2, 'F', 1. )
 
+      IF( jp_cfg == 1 )   THEN         ! Limit AHM south of -75 (critical for Giant ice shelves with small e1/e2)
+         ! special orca1 treatment remove the grid size scaling. Need to restore it south of -75N
+         ! define max grid size south of -75
+         imsk(:,:) = 0
+         WHERE( gphit(:,:) < -75.0 ) imsk(:,:) = 1  
+         zemax75 = MAX( MAXVAL( e1t(:,:)*imsk(:,:) ), MAXVAL( e2t(:,:)*imsk(:,:) ) )
+         IF( lk_mpp )   CALL mpp_max(zemax75)
+         !
+         ! apply grid size scaling of ahm south of -75 (no change north of it)
+         DO jj = 1, jpj
+            DO ji = 1, jpi
+               IF ( gphit(ji,jj) < -75.0 ) THEN
+                  zetmax = MAX( e1t(ji,jj), e2t(ji,jj) )
+                  zefmax = MAX( e1f(ji,jj), e2f(ji,jj) )
+                  ahm1(ji,jj,:) = zetmax / zemax75 * ahm1(ji,jj,:)
+                  ahm2(ji,jj,:) = zefmax / zemax75 * ahm2(ji,jj,:)
+               END IF
+            END DO
+         END DO
+      END IF
+
 
       IF(lwp) THEN                    ! Control print
@@ -432 +456 @@
          WRITE(numout,*) '         ahm4 array level 1'
          CALL prihre(ahm4(:,:,1),jpi,jpj,1,jpi,1,1,jpj,1,1.e-3,numout)
-      ENDIF
-      !
+      ENDIF     
+      ! 
       IF(lwp .AND. lflush) CALL flush(numout)
       !
       CALL wrk_dealloc( jpi   , jpj   , icof  )
+      CALL wrk_dealloc( jpi   , jpj   , imsk  )
       CALL wrk_dealloc( jpk   ,         zcoef )
       CALL wrk_dealloc( jpi   , jpj   , zahm0 )