In procedure clmain, separate coefh(klon, klev) into coefh(klon, 2:klev) and ycdragh(klon), coefm(klon, klev) into coefm(klon, 2:klev) and ycdragm(klon).