AI1 = 0.D0
AR1H = DCP*AR1-DSP*AI1
AI1 = DCP*AI1+DSP*AR1
C2(IK,LC) = AI1*CH2(IK,IP)
DS2 = AI1
AI2 = AI1
CH2(IK,LC) = AI1*C2(IK,IP)