Removing more commented out code.

Author: mwaxmonsky

src/physics/cam/diffusion_solver.F90

@@ -745,23 +745,15 @@ end function vd_lu_qdecomp
 
           ! Boundary layer thickness of "0._r8" signifies that the boundary
           ! condition is defined directly on the top interface.
-         !  decomp = fin_vol_lu_decomp(ztodt, p, &
-         !       coef_q_diff=kvh(:ncol,:)*dpidz_sq, &
-         !       upper_bndry=interface_boundary)
 
           if (.not. use_spcam) then
              dse(:ncol,:) = fin_vol_solve(ztodt, p, dse(:ncol,:), &
                                   ncol, pver, &
                                   coef_q_diff=kvh(:ncol,:)*dpidz_sq, &
                                   upper_bndry=interface_boundary, &
                                   l_cond=BoundaryData(dse_top(:ncol)))
-             !dse(:ncol,:) = dse(:ncol,:) + ddse
-            ! call decomp%left_div(dse(:ncol,:), &
-            !      l_cond=BoundaryData(dse_top(:ncol)))
           endif
 
-         !  call decomp%finalize()
-
           ! Calculate flux at top interface
 
           ! Modification : Why molecular diffusion does not work for dry static energy in all layers ?