Hello Sébastien and JMH,
I am attaching a zip file here containing all the files that I am having for ADR solution. After the last post, I have tried reducing the time step and refining the mesh further, although the solution has improved, but it is still not capturing the numerical solution.
Once again I have switched of the hydrodynamics, made the bed friction zero and imposed a constant u velocity of 0.03 m/s in a channel of about 10 KM length. The analytical and numerical solution are compared at x = 2.0 KM from the upstream end in the middle of the channel and at 5 m from the bed. The channel is 10 m deep.
In the attached zipped file you will also find sumit.f having routines where I am making the required changes. If you could kindly take a quick look and if you have any further suggestions I would really appreciate.
With kind reagrds,
Sumit