Changeset 10064 for branches/UKMO/CCRS_NUS_MC_STABLE/NEMOGCM/NEMO/OPA_SRC/BDY/bdyini.F90

Timestamp:

2018-08-23T14:34:13+02:00 (6 years ago)

Author:

csanchez

Message:

Added changes from Bijoy and Enda's for sponge layer, no tidal harm, and rimwidth=1

File:

• branches/UKMO/CCRS_NUS_MC_STABLE/NEMOGCM/NEMO/OPA_SRC/BDY/bdyini.F90 (modified) (6 diffs)

branches/UKMO/CCRS_NUS_MC_STABLE/NEMOGCM/NEMO/OPA_SRC/BDY/bdyini.F90

@@ -95 +95 @@
       REAL(wp), POINTER, DIMENSION(:,:)       ::   zfmask  ! temporary fmask array excluding coastal boundary condition (shlat)
 
-      !!
+      !! CEOD sponge
       NAMELIST/nambdy/ nb_bdy, ln_coords_file, cn_coords_file,                 &
          &             ln_mask_file, cn_mask_file, cn_dyn2d, nn_dyn2d_dta,     &
@@ -102 +102 @@
          &             cn_ice_lim, nn_ice_lim_dta,                           &
          &             rn_ice_tem, rn_ice_sal, rn_ice_age,                 &
-         &             ln_vol, nn_volctl, nn_rimwidth
+         &             ln_vol, nn_volctl, nn_rimwidth,                      &
+         &             ln_sponge, rn_sponge !CEOD 
       !!
       NAMELIST/nambdy_index/ ctypebdy, nbdyind, nbdybeg, nbdyend
@@ -414 +415 @@
         IF(lwp) WRITE(numout,*) '      Width of relaxation zone = ', nn_rimwidth(ib_bdy)
         IF(lwp) WRITE(numout,*)
-
+!CEOD
+        IF( ln_sponge(ib_bdy) ) THEN                     ! check sponge layer choice
+             IF(lwp) WRITE(numout,*) 'Sponge layer applied at open boundaries'
+             IF(lwp) WRITE(numout,*) 'Multiplier for diffusion in sponge layer :', rn_sponge
+             IF(lwp) WRITE(numout,*)
+           ELSE
+             IF(lwp) WRITE(numout,*) 'No volume correction applied at open boundaries'
+             IF(lwp) WRITE(numout,*)
+           ENDIF
+!
       ENDDO
 
@@ -432 +442 @@
         ENDIF
      ENDIF
+     sponge_factor(:,:) = 1.0 !CEOD
 
       ! -------------------------------------------------
@@ -1136 +1147 @@
             DO ib = 1, idx_bdy(ib_bdy)%nblen(igrd)
                nbr => idx_bdy(ib_bdy)%nbr(ib,igrd)
-               idx_bdy(ib_bdy)%nbw(ib,igrd) = 1.- TANH( FLOAT( nbr - 1 ) * 0.5 &
-                                            &  *(10./ FLOAT(nn_rimwidth(ib_bdy))) ) ! JGraham:modified for rim=15
-!               idx_bdy(ib_bdy)%nbw(ib,igrd) = 1.- TANH( FLOAT( nbr - 1 ) *0.5 )      ! tanh formulation
+!               idx_bdy(ib_bdy)%nbw(ib,igrd) = 1.- TANH( FLOAT( nbr - 1 ) *  0.5  &  
+!                                            &  *(10./ FLOAT(nn_rimwidth(ib_bdy))) ) ! JGraham:modified for rim=15
+               idx_bdy(ib_bdy)%nbw(ib,igrd) = 1.- TANH( FLOAT( nbr - 1 ) *0.5 )      ! tanh formulation !CBIJ MC config rim=1
 !               idx_bdy(ib_bdy)%nbw(ib,igrd) = (FLOAT(nn_rimwidth(ib_bdy)+1-nbr)/FLOAT(nn_rimwidth(ib_bdy)))**2.  ! quadratic
 !               idx_bdy(ib_bdy)%nbw(ib,igrd) =  FLOAT(nn_rimwidth(ib_bdy)+1-nbr)/FLOAT(nn_rimwidth(ib_bdy))       ! linear
@@ -1155 +1166 @@
             END DO
          END DO 
+      ! Compute multiplier for diffusion for sponge layer !CEOD
+           ! -------------------------------------------------
+           IF( ln_sponge(ib_bdy) ) THEN
+              igrd=1
+              DO ib = 1, idx_bdy(ib_bdy)%nblen(igrd)
+                 nbi => idx_bdy(ib_bdy)%nbi(ib,igrd)
+                 nbj => idx_bdy(ib_bdy)%nbj(ib,igrd)
+                 nbr => idx_bdy(ib_bdy)%nbr(ib,igrd)
+                 sponge_factor(nbi,nbj) = 1.0 + (rn_sponge-1.0) * ( 1.- TANH( FLOAT( nbr - 1 ) *0.5 ) )
+                 ! 0.5 factor says how fast tanh goes to 1 - reduce this to have
+                 ! effect on more of the rimwidth
+              END DO
+           ENDIF
 
       ENDDO