work on FD-ACM
further work in progress (not a working code!!)
work in progress
first time working code in fortran