call reada(controlcard,'SCALSCP',scalscp,1.0d0)
call reada(controlcard,'CUTOFF',cutoff_corr,7.0d0)
call reada(controlcard,'DELT_CORR',delt_corr,0.5d0)
+ call reada(controlcard,'LIPSCALE',lipscale,1.0D0)
+c write (iout,*) "lipscale",lipscale
r0_corr=cutoff_corr-delt_corr
write (iout,*) "iparm",iparm," myparm",myparm