correction for lsd model

hamidrezanorouzi · Oct 5, 2022 · 985e0a1 · 985e0a1
1 parent 82aecc6
commit 985e0a1
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/g_Prtcl_LSD_Model.f90 b/src/g_Prtcl_LSD_Model.f90
@@ -184,7 +184,7 @@ subroutine LSD_ContactForce_PP(this , ind , lnew )
 
  if( norm(ovlp_t) .gt. 0.0_RK ) then
 
- if( this%CF_Type == CFT_nLin_l ) then
+ if( this%CF_Type == CFT_LSD_l ) then
  ! limited overlap
  ftij = ftij*(ft_fric/ft)
  ovlp_t = (-1.0_RK) * (ftij/prop_ij%kt)