In procedure clmain, do not compute ycoefh(:, 1) since this is a duplicate of cdragh. Thus, in procedure physiq, use cdragh as actual argument to phytrac.
trunk/Sources/phylmd/clmain.f
trunk/Sources/phylmd/physiq.f