Variable fevap of physiq is not used. Remove it from physiq and from the restart file. Remove the corresponding argument evap of pbl_surface. Use directly yqsurf instead of qairsol in pbl_surface.