c var.inc
	COMMON /ION_VAR/
     $	        psiin(njp1,nk2p1),psiis(njp1,nk2p1)