Changeset 10064 for branches/UKMO/CCRS_NUS_MC_STABLE/NEMOGCM/NEMO/OPA_SRC/BDY/bdy_oce.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/bdy_oce.F90 (modified) (3 diffs)

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

@@ -90 +90 @@
    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  CEOD
    !
    INTEGER                    ::   nb_bdy                   !: number of open boundary sets
    INTEGER, DIMENSION(jp_bdy) ::   nn_rimwidth              !: boundary rim width for Flow Relaxation Scheme
    INTEGER                    ::   nn_volctl                !: = 0 the total volume will have the variability of the surface Flux E-P 
+   REAL(wp)                   ::   rn_sponge                !: multiplier of diffusion for sponge layer  CEOD
    !                                                        !  = 1 the volume will be constant during all the integration.
    CHARACTER(len=20), DIMENSION(jp_bdy) ::   cn_dyn2d       ! Choice of boundary condition for barotropic variables (U,V,SSH)
@@ -130 +132 @@
    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 CEOD
 
    REAL(wp)                                    ::   bdysurftot !: Lateral surface of unstructured open boundary
@@ -160 +163 @@
       !
       ALLOCATE( bdytmask(jpi,jpj) , bdyumask(jpi,jpj), bdyvmask(jpi,jpj),     &  
-         &      STAT=bdy_oce_alloc )
+         &      sponge_factor(jpi,jpj), STAT=bdy_oce_alloc )   !CEOD
       !
       ! Initialize masks