Changeset 7123 for branches/UKMO/r5518_bdy_sponge_temp/NEMOGCM/NEMO/OPA_SRC/BDY/bdyini.F90

Timestamp:

2016-10-26T19:32:29+02:00 (8 years ago)

Author:

jcastill

Message:

Changes as in UKMO/dev_r5518_bdy_sponge_temp@5878

File:

• branches/UKMO/r5518_bdy_sponge_temp/NEMOGCM/NEMO/OPA_SRC/BDY/bdyini.F90 (modified) (4 diffs)

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

@@ -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
       !!
       NAMELIST/nambdy_index/ ctypebdy, nbdyind, nbdybeg, nbdyend
@@ -365 +366 @@
         IF(lwp) WRITE(numout,*) '      Width of relaxation zone = ', nn_rimwidth(ib_bdy)
         IF(lwp) WRITE(numout,*)
+        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
@@ -383 +392 @@
         ENDIF
      ENDIF
+
+     sponge_factor(:,:) = 1.0
 
       ! -------------------------------------------------
@@ -1104 +1115 @@
             END DO
          END DO 
+         ! Compute multiplier for diffusion for sponge layer
+         ! -------------------------------------------------
+         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