In procedure phytrac, separate dummy argument coefh(klon, klev) into coefh(klon, 2:klev) and cdragh(klon).