E:/turkuaz/krylov/codes/newton/1.2b/main.f90

Go to the documentation of this file.

program main

use physicalParameters
use postProcessor
use myFunctions
implicit none

integer i

real(8), allocatable :: X(:)


! initialize

call initializePhysics()

!print*,maxNewtonIter

!stop







allocate(X(N))




call initialGuess(X,n)





call newton(X,n)










!POST-PROCESSING




if (printSolution.EQV.ON) then
        
        ! hepsi tek bir module içinde contains ile yaazýlabilir
        ! local variabller ile çalýþýp mesh olayý ayarlanmalý
        ! bunlar yoksa boþ dönecek
        
        !CALL postProcesser()   !<---------- will be the main routine
        !call anExtraRoutineLikeCalculateFlux  !<--- will be already written in the post processor so no error will be generated
        CALL visualizer(X,N)

end if








if (displaySolution.EQV.ON) then
        print*,''
        print*,'solution:'
        print*,''
        do i=1,n
                print*,X(i)
        end do
end if








end program main

---