Changeset 5878 for branches/UKMO/dev_r5518_bdy_sponge_temp/NEMOGCM/NEMO/OPA_SRC/BDY/bdy_oce.F90

Timestamp:

2015-11-11T19:11:01+01:00 (8 years ago)

Author:

deazer

Message:

Adds simple ramp of diffusion coefficient in BDY rim using TANH weighting

File:

• branches/UKMO/dev_r5518_bdy_sponge_temp/NEMOGCM/NEMO/OPA_SRC/BDY/bdy_oce.F90 (modified) (4 diffs)

branches/UKMO/dev_r5518_bdy_sponge_temp/NEMOGCM/NEMO/OPA_SRC/BDY/bdy_oce.F90

@@ -83 +83 @@
    LOGICAL                    ::   ln_mask_file             !: =T read bdymask from file
    LOGICAL                    ::   ln_vol                   !: =T volume correction             
+   LOGICAL, DIMENSION(jp_bdy) ::   ln_sponge                !: =T use sponge layer 
    !
    INTEGER                    ::   nb_bdy                   !: number of open boundary sets
@@ -88 +89 @@
    INTEGER                    ::   nn_volctl                !: = 0 the total volume will have the variability of the surface Flux E-P 
    !                                                        !  = 1 the volume will be constant during all the integration.
+   REAL(wp)                   ::   rn_sponge                !: multiplier of diffusion for sponge layer
    CHARACTER(len=20), DIMENSION(jp_bdy) ::   cn_dyn2d       ! Choice of boundary condition for barotropic variables (U,V,SSH)
    INTEGER, DIMENSION(jp_bdy)           ::   nn_dyn2d_dta   !: = 0 use the initial state as bdy dta ; 
@@ -118 +120 @@
    REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:), TARGET ::   bdyumask   !: Mask defining computational domain at U-points
    REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:), TARGET ::   bdyvmask   !: Mask defining computational domain at V-points
+   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:), TARGET ::   sponge_factor !: Multiplier for diffusion for sponge layer
 
    REAL(wp)                                    ::   bdysurftot !: Lateral surface of unstructured open boundary
-
    !!----------------------------------------------------------------------
    !! open boundary data variables
@@ -148 +150 @@
       !
       ALLOCATE( bdytmask(jpi,jpj) , bdyumask(jpi,jpj), bdyvmask(jpi,jpj),     &  
-         &      STAT=bdy_oce_alloc )
+         &      sponge_factor(jpi,jpj), STAT=bdy_oce_alloc )
       !
       ! Initialize masks