Index: trunk/LMDZ.PLUTO/libf/phypluto/callcorrk.F90
===================================================================
--- trunk/LMDZ.PLUTO/libf/phypluto/callcorrk.F90	(revision 3889)
+++ trunk/LMDZ.PLUTO/libf/phypluto/callcorrk.F90	(revision 3890)
@@ -210,4 +210,5 @@
       integer Nfine,ifine
       parameter(Nfine=701)
+      character(len=100) :: file_path
       real,save :: levdat(Nfine),vmrdat(Nfine)
       REAL dtlw_hcn_c2h2(ngrid, nlayer) ! cooling rate (K/s) due to C2H2/HCN (diagnostic)
@@ -489,4 +490,14 @@
          endif ! end callmufi
       endif ! end radiative haze
+
+      ! If fixed profile of CH4 gas
+      IF (vmrch4_proffix) then
+         file_path=trim(datadir)//'/gas_prop/vmr_ch4.txt'
+         open(115,file=file_path,form='formatted')
+         do ifine=1,Nfine
+           read(115,*) levdat(ifine), vmrdat(ifine)
+         enddo
+         close(115)
+      ENDIF
 
       ! How much light do we get ?