In procedure physiq, rename dsens and devap to dflux_t and dflux_q so they correspond to arguments with the same name in pbl_surface. (In LMDZ, dsens and devap are not used in physiq, the computation of fder is done inside pbl_surface.)