Alin=viscA4D+vg4*L4+Alth4
Alin=viscA4Z+vg4*L4+Alth4
_RL vg2,vg2Min,vg2Max,vg4,vg4Min,vg4Max
vg4=viscA4Grid*recip_dt