404 |
& +viscAh_DLth(i,j)+viscAh_DSmg(i,j) |
& +viscAh_DLth(i,j)+viscAh_DSmg(i,j) |
405 |
#ifdef ALLOW_3D_VISCAH |
#ifdef ALLOW_3D_VISCAH |
406 |
& +viscAhDfld(i,j,k,bi,bj) |
& +viscAhDfld(i,j,k,bi,bj) |
407 |
|
#ifdef ALLOW_AUTODIFF |
408 |
|
& *viscFacAdj |
409 |
|
#endif |
410 |
#endif |
#endif |
411 |
viscAh_DMin(i,j)=MAX(viscAhGridMin*L2rdt,Uscl) |
viscAh_DMin(i,j)=MAX(viscAhGridMin*L2rdt,Uscl) |
412 |
viscAh_D(i,j)=MAX(viscAh_DMin(i,j),Alin) |
viscAh_D(i,j)=MAX(viscAh_DMin(i,j),Alin) |
418 |
& +viscA4_DLth(i,j)+viscA4_DSmg(i,j) |
& +viscA4_DLth(i,j)+viscA4_DSmg(i,j) |
419 |
#ifdef ALLOW_3D_VISCA4 |
#ifdef ALLOW_3D_VISCA4 |
420 |
& +viscA4Dfld(i,j,k,bi,bj) |
& +viscA4Dfld(i,j,k,bi,bj) |
421 |
|
#ifdef ALLOW_AUTODIFF |
422 |
|
& *viscFacAdj |
423 |
|
#endif |
424 |
#endif |
#endif |
425 |
viscA4_DMin(i,j)=MAX(viscA4GridMin*L4rdt,U4scl) |
viscA4_DMin(i,j)=MAX(viscA4GridMin*L4rdt,U4scl) |
426 |
viscA4_D(i,j)=MAX(viscA4_DMin(i,j),Alin) |
viscA4_D(i,j)=MAX(viscA4_DMin(i,j),Alin) |