Changeset 11257

Timestamp:

2019-07-11T17:57:52+02:00 (5 years ago)

Author:

mathiot

Message:

add scaling depending on grid size for eorca1 south of -75N

File:

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

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

@@ -203 +203 @@
       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
@@ -213 +214 @@
       !
       CALL wrk_alloc( jpi   , jpj   , icof  )
+      CALL wrk_alloc( jpi   , jpj   , imsk  )
       CALL wrk_alloc( jpk   ,         zcoef )
       CALL wrk_alloc( jpi   , jpj   , zahm0 )
@@ -250 +252 @@
          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 )
          zeref = MAXVAL ( e1t(:,:) * e2t(:,:),   &
-             &   tmask(:,:,1) .GE. 0.5 .AND. ABS(gphit(:,:)) .GT. 50. )
+             &   ssmask(:,:) .GE. 0.5 .AND. ABS(gphit(:,:)) .GT. 50. )
  
          DO jj = 1, jpj
@@ -314 +316 @@
       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.
@@ -380 +382 @@
       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
@@ -428 +451 @@
       !
       CALL wrk_dealloc( jpi   , jpj   , icof  )
+      CALL wrk_dealloc( jpi   , jpj   , imsk  )
       CALL wrk_dealloc( jpk   ,         zcoef )
       CALL wrk_dealloc( jpi   , jpj   , zahm0 )