_RL Uscl,U4scl
U4scl=SQRT(KE(i,j))*L3*viscA4Re_max
U4scl=0.
viscA4_DMin(i,j)=MAX(viscA4GridMin*L4rdt,U4scl)
U4scl= SQRT(keZpt)*L3*viscA4Re_max
viscA4_ZMin(i,j)=MAX(viscA4GridMin*L4rdt,U4scl)