fcpertplus = fcppmem (i)
& fcref, fcpertplus, fcpertminus,
_RL fcpertplus, fcpertminus
fcpertplus = fc
print *, 'ph-check fcpertplus = ', fcpertplus
gfd = (fcpertplus-fcpertminus)
fcppmem ( ichknum ) = fcpertplus
& fcref, fcpertplus, fcpertminus