Alin=viscAhD+vg2*L2+Alth2
Alin=viscAhZ+vg2*L2+Alth2
_RL vg2,vg2Min,vg2Max,vg4,vg4Min,vg4Max
vg2=viscAhGrid*recip_dt