BUGFIX: line 160 - divisiion by dy instead of dx --This line, and those below, will be ignored-- M concoutput.f90