c**********************************************************************cc********************** sum two real arrays ***************************c        subroutine myreal_somme (c1,c2,c,nvecs)	        implicit  none        integer*4 nvecs,i        real*8    c1(nvecs),c2(nvecs),c(nvecs)        do 100 i = 1, nvecs           c(i) = c1(i) + c2(i) 100    continue        return        endc**********************************************************************cc********************** subtract two real arrays **********************c        subroutine myreal_diffr (c1,c2,c,nvecs)        implicit  none        integer*4 nvecs,i        real*8    c1(nvecs),c2(nvecs),c(nvecs)        do 100 i = 1, nvecs           c(i) = c1(i) - c2(i) 100    continue        return        endc**********************************************************************cc************* multiply an real array with a real constant ************c        subroutine myreal_mult (in,const,out,nvecs)        implicit  none        integer*4 nvecs,i        real*8    in(nvecs),const,out(nvecs)        do 100 i = 1, nvecs           out(i) = const*in(i) 100    continue        return        endc**********************************************************************cc**************** scalar sum product of two real arrays ***************c        subroutine myrescal (c1,c2,ps,nvecs)        implicit  none        integer*4 nvecs,i        real*8    c1(nvecs),c2(nvecs),ps        ps = 0.0        do 100 i = 1, nvecs           ps = ps + c1(i)*c2(i) 100    continue        return        endc**********************************************************************cc**********************************************************************c