Changeset 3917

Timestamp:

Sep 19, 2025, 12:17:56 PM (4 months ago)

Author:

tbertrand

Message:

Pluto PCM:
Adding Triton options (CH4 profile)
TB

Location:

trunk/LMDZ.PLUTO/libf/phypluto

Files:

• callcorrk.F90 (modified) (3 diffs)
• hazecloud.F90 (modified) (3 diffs)

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

@@ -37 +37 @@
                               methane,carbox,cooling,nlte,strobel,&
                               ch4fix,vmrch4_proffix,vmrch4fix,&
-                              callmufi
+                              callmufi,triton
       use optcv_mod, only: optcv
       use optci_mod, only: optci
@@ -209 +209 @@
       real eps_nlte_lw(ngrid,nlayer) ! CH4 NLTE efficiency factor for zdtsw
       integer Nfine,ifine
-      parameter(Nfine=701)
+      parameter(Nfine=600)
       character(len=100) :: file_path
       real,save :: levdat(Nfine),vmrdat(Nfine)
@@ -493 +493 @@
       ! If fixed profile of CH4 gas
       IF (vmrch4_proffix) then
+        IF (triton) then
+         file_path=trim(datadir)//'/gas_prop/vmr_ch4_triton.txt'
+        ELSE 
          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
+        open(115,file=file_path,form='formatted')
+        do ifine=1,Nfine
+          read(115,*) levdat(ifine), vmrdat(ifine)
+        enddo
+        close(115)
       ENDIF
 

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

@@ -8 +8 @@
       use tracer_h, only: igcm_haze, igcm_ch4_gas, igcm_n2, igcm_prec_haze, noms, mmol
       use geometry_mod, only: longitude, latitude ! in radians
-      use callkeys_mod, only: hazeconservch4, haze_ch4proffix, diurnal, tcon_ch4, k_ch4, ncratio_ch4
+      use callkeys_mod, only: hazeconservch4, haze_ch4proffix, diurnal, tcon_ch4, k_ch4, & 
+                              ncratio_ch4,triton
       use datafile_mod, only: datadir
 
@@ -90 +91 @@
 
       integer Nfine,ifine
-      parameter(Nfine=701)
+      parameter(Nfine=600)
       character(len=100) :: file_path
       real,save :: levdat(Nfine),vmrdat(Nfine)
@@ -167 +168 @@
       ! If fixed profile of CH4 gas
       if (haze_ch4proffix) then
+        if (triton) then
+         file_path=trim(datadir)//'/gas_prop/vmr_ch4_triton.txt'
+        else
          file_path=trim(datadir)//'/gas_prop/vmr_ch4.txt'
-         open(115,file=file_path,form='formatted')
-         do ifine=1,Nfine
+        endif
+        open(115,file=file_path,form='formatted')
+        do ifine=1,Nfine
            read(115,*) levdat(ifine), vmrdat(ifine)
-         enddo
-         close(115)
+        enddo
+        close(115)
       endif
       vmrch4(:,:)=0.