In clmain, local variable yfder was computed but not used. I think it was useful for coupling only. Variable fder_print of pbl_surface in LMDZ, which is output by LMDZ, corresponds to variable fder of physiq in LMDZ and LMDZE.