Debugging Fortran 90 code in order to solve 2D heat equation optimisation problem.

I have written the coding implementation for a 2D heat equation problem. I need someone to debug the code in order for me to achieve an optimised temperature variation. The optimisation tool being used is IPOPT, I have already linked the optimisation code to the library. I need a good programmer to produce an output of optimised temperature variation. I have attached a pdf file describing the project. As well as three .f90 files which represent the coding implementation of the problem. The GNU represents the current un-optimised temperature variation. The aim is to produce a graph of the same format but to establish an OPTIMISED temperature variation. The .f file is the key file which i need to programmer to work on in order to achieve results. Currently, the code needs debugging. The time frame for this project is limited, i need this project done within the next 3 days.