Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/module_initialize_les.F
===================================================================
--- trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/module_initialize_les.F	(revision 2018)
+++ trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/module_initialize_les.F	(revision 2020)
@@ -118,4 +118,5 @@
  REAL, DIMENSION(nl_max) :: prescribed_sw,prescribed_lw,prescribed_dyn
  REAL, DIMENSION(nl_max) :: hrsw,hrlw,hrdyn
+ REAL, DIMENSION(nl_max) :: lsf_dt,lsf_dq,lsfdt,lsfdq
  REAL, DIMENSION(nl_max) :: venus_hrdyn
  REAL, DIMENSION(nl_max) :: altitude
@@ -955,4 +956,17 @@
       ENDDO
       close(unit=20)
+    ENDIF
+
+    IF (planet.eq."generic") THEN
+      call read_lsf(lsfdt,lsfdq,nl_in)
+      open(unit=17,file="lsf.txt",action="write")
+      DO k=1,kte!-1
+        p_level = grid%znu(k)*(pd_surf - grid%p_top) + grid%p_top
+        lsf_dt = interp_0_log( lsfdt, pd_in, p_level, nl_in )
+        lsf_dq = interp_0_log( lsfdq, pd_in, p_level, nl_in )
+        write (17,*) lsf_dt(k),lsf_dq(k)
+      ENDDO
+    ENDIF
+
 
     open(unit=21,file="altitude.txt",action="write")
@@ -963,5 +977,4 @@
     ENDDO
     close(unit=21)
-    ENDIF
 
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
@@ -1317,3 +1330,31 @@
       end subroutine read_hr
 
+      subroutine read_lsf(dt,dq,n)
+      implicit none
+      integer n
+      real dt(n+1),dq(n+1)
+      logical end_of_file
+
+      integer k
+
+! first element is the surface
+
+      open(unit=12,file='input_lsf',form='formatted',status='old')
+      rewind(12)
+      end_of_file = .false.
+      k = 0
+      do while (.not. end_of_file)
+
+        read(12,*,end=103) dt(k+1),dq(k+1)
+        write(*,*) k,dt(k+1),dq(k+1)
+        k = k+1
+        go to 114
+ 103    end_of_file = .true.
+ 114    continue
+      enddo
+
+      close(unit=12,status = 'keep')
+
+      end subroutine read_lsf
+
 END MODULE module_initialize_ideal