Changeset 348
- Timestamp:
- 09/15/21 09:32:46 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SOURCES/dragging_coulomb_friction_mod.f90
r334 r348 237 237 tauc_my(:,:)= cf*neffmy(:,:)*coef_sedim_my(:,:) 238 238 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???) 241 where (abs(uxbar(:,:)).gt.1) 242 betamx(:,:) = tauc_mx(:,:) * ( abs(uxbar(:,:))**(q_nolin-1.) ) / ( u0 ** q_nolin ) 241 243 elsewhere 242 244 betamx(:,:) = tauc_mx(:,:) / ( u0 ** q_nolin ) 243 245 endwhere 244 where (abs(uy (:,:,nz)).gt.1)245 betamy(:,:) = tauc_my(:,:) / ( abs(uy(:,:,nz))**(q_nolin-1.) ) / ( u0 ** q_nolin )246 where (abs(uybar(:,:)).gt.1) 247 betamy(:,:) = tauc_my(:,:) * ( abs(uybar(:,:))**(q_nolin-1.) ) / ( u0 ** q_nolin ) 246 248 elsewhere 247 249 betamy(:,:) = tauc_my(:,:) / ( u0 ** q_nolin )
Note: See TracChangeset
for help on using the changeset viewer.