Changeset 12788


Ignore:
Timestamp:
2020-04-21T13:21:43+02:00 (6 months ago)
Author:
davestorkey
Message:

UKMO/NEMO_4.0.1_biharmonic_GM : New calculation for zmu with new rn_bgm_msc namelist parameter.

Location:
NEMO/branches/UKMO/NEMO_4.0.1_biharmonic_GM
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/UKMO/NEMO_4.0.1_biharmonic_GM/cfgs/SHARED/namelist_ref

    r12535 r12788  
    949949      !                       !  iso-neutral laplacian operator (ln_dynldf_iso=T) : 
    950950      rn_ahm_b      = 0.0         !  background eddy viscosity  [m2/s] 
    951    nn_bhm_ijk_t  = 0           !  space/time variation of bilaplacian GM coefficient : 
     951      nn_bhm_ijk_t  = 0           !  space/time variation of bilaplacian GM coefficient : 
    952952      !                             !  = 11  uniform coefficient in interior, zero top and bottom 
    953953      !                             !  = 12  linear tapering with dz to surface, zero at top and bottom 
     
    955955      rn_bhm_b      = 0.0         ! bilaplacian GM coefficient 
    956956      rn_bgmzcrit   = 0.0         ! critical depth for depth variation of bilaplacian GM coefficient 
     957      rn_bgm_msc    = 16.0        ! multiplicative parameter for Method of Stabilising Correction (BGM) 
    957958/ 
    958959!----------------------------------------------------------------------- 
  • NEMO/branches/UKMO/NEMO_4.0.1_biharmonic_GM/src/OCE/DYN/dynldf_lap_blp.F90

    r12777 r12788  
    146146      CALL dyn_ldf_lap( kt, tmask, fmask, pub, pvb, zulap, zvlap, 1 )    
    147147      ! 
    148       ! Calculate ratio of bhm and avm for stabilising correction terms 
     148      ! Calculate zmu for stabilising correction terms 
    149149      ! and add to avm to be included in the vertical diffusion calculation later. 
    150150      zmu(:,:,:) = 0.0_wp 
     
    152152         DO jj = 2, jpjm1 
    153153            DO ji = fs_2, jpim1   ! vector opt. 
    154                zmu(ji,jj,jk) = ( bhm(ji,jj,jk) / (avm(ji,jj,jk) + rsmall) ) * wmask(ji,jj,jk) 
     154               zmu(ji,jj,jk) = ( rn_bgm_msc * bhm(ji,jj,jk) / e1e2t(ji,jj) ) * wmask(ji,jj,jk) 
    155155               avm(ji,jj,jk) = avm(ji,jj,jk) + zmu(ji,jj,jk) 
    156156            ENDDO 
  • NEMO/branches/UKMO/NEMO_4.0.1_biharmonic_GM/src/OCE/LDF/ldfdyn.F90

    r12777 r12788  
    5555   REAL(wp), PUBLIC ::   rn_bhm_b              !: lateral bi-Laplacian GM background eddy viscosity  [m4/s] 
    5656   REAL(wp), PUBLIC ::   rn_bgmzcrit           !: critical depth (>=0) for linear tapering of bi-Lap. GM coef. to zero at surface [m]  
     57   REAL(wp), PUBLIC ::   rn_bgm_msc            !: multiplicative parameter for Method of Stabilising Correction: zero = no correction term. 
    5758   !                                    !!* Parameter to control the type of lateral viscous operator 
    5859   INTEGER, PARAMETER, PUBLIC ::   np_ERROR  =-10                       !: error in setting the operator 
     
    125126         &                 ln_dynldf_lev, ln_dynldf_hor, ln_dynldf_iso,   &   ! acting direction of the operator 
    126127         &                 nn_ahm_ijk_t , rn_Uv    , rn_Lv,   rn_ahm_b,   &   ! lateral eddy coefficient 
    127          &                 nn_bhm_ijk_t, rn_bhm_b, rn_bgmzcrit,           &  ! lateral bi-Lap. GM coefficient 
     128         &                 nn_bhm_ijk_t, rn_bhm_b, rn_bgmzcrit, rn_bgm_msc, & ! lateral bi-Lap. GM coefficient 
    128129         &                 rn_csmc      , rn_minfac    , rn_maxfac            ! Smagorinsky settings 
    129130      !!---------------------------------------------------------------------- 
     
    150151         WRITE(numout,*) '         bi-Laplacian GM operator             ln_dynldf_bgm = ', ln_dynldf_bgm 
    151152         WRITE(numout,*) '         critical depth for taper (if used)   rn_bgmzcrit   = ', rn_bgmzcrit 
     153         WRITE(numout,*) '         MSC parameter                        rn_bgm_msc    = ', rn_bgm_msc 
    152154         ! 
    153155         WRITE(numout,*) '      direction of action :' 
Note: See TracChangeset for help on using the changeset viewer.