Changeset 348


Ignore:
Timestamp:
09/15/21 09:32:46 (3 years ago)
Author:
aquiquet
Message:

Bug correction for Coulomb friction law

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/SOURCES/dragging_coulomb_friction_mod.f90

    r334 r348  
    237237tauc_my(:,:)= cf*neffmy(:,:)*coef_sedim_my(:,:) 
    238238 
    239 where (abs(ux(:,:,nz)).gt.1) 
    240    betamx(:,:) = tauc_mx(:,:) * ( abs(ux(:,:,nz))**(q_nolin-1.) ) / ( u0 ** q_nolin ) 
     239! ux/uy(:,:,nz) should be used but only uxbar/uybar are updated by diagno_L2 
     240! anyway: ux/uy(:,:,nz) are uxbar/uybar (as it should be???) 
     241where (abs(uxbar(:,:)).gt.1) 
     242   betamx(:,:) = tauc_mx(:,:) * ( abs(uxbar(:,:))**(q_nolin-1.) ) / ( u0 ** q_nolin ) 
    241243elsewhere 
    242244   betamx(:,:) = tauc_mx(:,:) / ( u0 ** q_nolin ) 
    243245endwhere 
    244 where (abs(uy(:,:,nz)).gt.1) 
    245    betamy(:,:) = tauc_my(:,:) / ( abs(uy(:,:,nz))**(q_nolin-1.) ) / ( u0 ** q_nolin ) 
     246where (abs(uybar(:,:)).gt.1) 
     247   betamy(:,:) = tauc_my(:,:) * ( abs(uybar(:,:))**(q_nolin-1.) ) / ( u0 ** q_nolin ) 
    246248elsewhere 
    247249   betamy(:,:) = tauc_my(:,:) / ( u0 ** q_nolin ) 
Note: See TracChangeset for help on using the changeset viewer.