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.