Changeset 3890

Timestamp:

Aug 12, 2025, 2:48:54 PM (28 hours ago)

Author:

debatzbr

Message:

Pluto PCM: Fix bug - read data file when vmrch4_proffix = .true.
BBT

File:

• trunk/LMDZ.PLUTO/libf/phypluto/callcorrk.F90 (modified) (2 diffs)

trunk/LMDZ.PLUTO/libf/phypluto/callcorrk.F90

@@ -210 +210 @@
       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 +490 @@
          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 ?