Index: /trunk/LMDZ.GENERIC/changelog.txt
===================================================================
--- /trunk/LMDZ.GENERIC/changelog.txt	(revision 4177)
+++ /trunk/LMDZ.GENERIC/changelog.txt	(revision 4178)
@@ -2265,2 +2265,5 @@
 unused module variables.
 While at it, put all thermcell_*.F90 routines in modules.
+
+== 10/04/2026 == EM
+Code cleanup: remove unused surface_nature.F and scopyi.F
Index: unk/LMDZ.GENERIC/libf/phygeneric/scopyi.F
===================================================================
--- /trunk/LMDZ.GENERIC/libf/phygeneric/scopyi.F	(revision 4177)
+++ 	(revision )
@@ -1,17 +1,0 @@
-      subroutine scopyi(n,sx,incx,sy,incy)
-c
-      IMPLICIT NONE
-c
-      integer n,incx,incy,ix,iy,i
-      integer sx((n-1)*incx+1),sy((n-1)*incy+1)
-c
-      iy=1
-      ix=1
-      do 10 i=1,n
-      sy(iy)=sx(ix)
-         ix=ix+incx
-         iy=iy+incy
-10    continue
-c
-      return
-      end
Index: unk/LMDZ.GENERIC/libf/phygeneric/surface_nature.F
===================================================================
--- /trunk/LMDZ.GENERIC/libf/phygeneric/surface_nature.F	(revision 4177)
+++ 	(revision )
@@ -1,64 +1,0 @@
-      SUBROUTINE surface_nature(ngrid,nq,obliquit,qsurf,qsurfliquid
-     &   ,qsurfsnow,rnat,oceanarea)
-
-      USE surfdat_h
-      USE comsoil_h
-      USE geometry_mod, ONLY: cell_area
-      USE tracer_h
-
-      IMPLICIT none
-
-!==================================================================
-!     
-!     Purpose
-!     -------
-!     Defines a few things
-!     
-!     Authors
-!     ------- 
-!     B. Charnay (2010)
-!     
-!     Called by
-!     ---------
-!     physiq.F
-!     
-!     Calls
-!     -----
-!     none
-!
-!     Notes
-!     -----
-!     rnat is terrain type: 0-ocean; 1-continent; 2-continental ice
-!     
-!==================================================================
-
-        integer ngrid,nq
-
-	REAL qsurf(ngrid,nq),ps(ngrid)
-	REAL qsurfliquid(ngrid)
-	REAL qsurfsnow(ngrid)
-	INTEGER iq, ig
-	INTEGER rnat(ngrid)
- 	REAL oceanarea
-	REAL obliquit
-
-	do ig=1,ngrid
-           rnat(ig)=1
-           dryness(ig)=1        !(coefficient for evaporation)
-           if (inertiedat(ig,1).gt.1E4) then
-              rnat(ig)=0
-           end if
-	end do
-
-! surface of all the oceans
-        
-        oceanarea=0.
-	do ig=1,ngrid
-           if (rnat(ig).eq.0)then
-              oceanarea=oceanarea+cell_area(ig)
-           end if
-	enddo
-
-        return
-        end
-