save_XY.f90

    subroutine save_XY
    
    include 'com.txt'
    
    open(unit = 1,file = 'Xc.txt')
    open(unit = 2,file = 'Yc.txt')
    
    do i = 1,Nx
        do i1 = 1,NumGLP
            write(1,*) Xc(i) + hx1*lambda(i1)
        end do
    end do
    
    do j = 1,Ny
        do j1 = 1,NumGLP
            write(2,*) Yc(j) + hy1*lambda(j1)
        end do
    end do
    
    close(1)
    close(2)
    
    end subroutine save_XY