Changeset 3978 for trunk

Timestamp:

Nov 27, 2025, 5:24:15 PM (6 days ago)

Author:

debatzbr

Message:

Pluto PCM: Initialize gas profiles for clouds in 1D.
BBT

File:

• trunk/LMDZ.PLUTO/libf/phypluto/dyn1d/rcm1d.F (modified) (5 diffs)

trunk/LMDZ.PLUTO/libf/phypluto/dyn1d/rcm1d.F

@@ -22 +22 @@
       use time_phylmdz_mod, only: daysec, dtphys, diagfi_output_rate,
      &                            nday
-      use callkeys_mod, only: tracer, specOLR,pceil,haze
+      use callkeys_mod, only: tracer, specOLR,pceil,haze,
+     &                        callmufi,callmuclouds
       USE comvert_mod, ONLY: ap,bp,aps,bps,pa,preff, sig,
      &                       presnivs,pseudoalt,scaleheight
@@ -79 +80 @@
       INTEGER lecttsoil     ! lecture of tsoil from proftsoil
       INTEGER lecthaze      ! lecture of haze from profhaze
+      INTEGER lectC2H2      ! lecture of gases from profC2H2
+      INTEGER lectC2H6      ! lecture of gases from profC2H6
+      INTEGER lectC4H2      ! lecture of gases from profC4H2
+      INTEGER lectC6H6      ! lecture of gases from profC6H6
+      INTEGER lectHCN       ! lecture of gases from profHCN
       REAL day              ! date during the run
       INTEGER day_step      ! number of time steps per day
@@ -96 +102 @@
       integer :: i_co_ice=0      ! tracer index of co ice
       integer :: i_co_gas=0      ! tracer index of co gas
+      integer :: i_C2H2_mugas=0  ! tracer index of C2H2 gas
+      integer :: i_C2H6_mugas=0  ! tracer index of C2H6 gas
+      integer :: i_C4H2_mugas=0  ! tracer index of C4H2 gas
+      integer :: i_C6H6_mugas=0  ! tracer index of C6H6 gas
+      integer :: i_HCN_mugas=0   ! tracer index of HCN gas
       integer :: i_prec_haze=0   ! tracer index of haze
       integer :: i_haze=0  ! tracer index of haze
@@ -343 +354 @@
            elseif (tname(iq)=="co_gas") then
              i_co_gas=iq
+           elseif (tname(iq)=="C2H2_mugas") then
+             i_C2H2_mugas=iq
+           elseif (tname(iq)=="C2H6_mugas") then
+             i_C2H6_mugas=iq
+           elseif (tname(iq)=="C4H2_mugas") then
+             i_C4H2_mugas=iq
+           elseif (tname(iq)=="C6H6_mugas") then
+             i_C6H6_mugas=iq
+           elseif (tname(iq)=="HCN_mugas") then
+             i_HCN_mugas=iq
            elseif (tname(iq)=="haze") then
              i_haze=iq
@@ -909 +930 @@
       endif
       endif
-! Initialize cloud fraction and oceanic ice !AF24: removed
-! Initialize slab ocean !AF24: removed
-! Initialize chemical species !AF24: removed photochem
+
+!     Initialize gas profiles for clouds
+!     ------------------------------------------
+      if (callmufi.and.callmuclouds) then
+        lectC2H2 = 0 ! default value for lectC2H2
+        call getin("lectC2H2",lectC2H2)
+        
+        if (lectC2H2 == 1) then
+          OPEN(15,file='profC2H2',status='old',form='formatted',err=501)
+          DO iq = 1, nq
+            if (iq.eq.i_C2H2_mugas) then
+              DO ilayer=1,nlayer
+                READ (15,*) q(ilayer,iq)
+              ENDDO
+            endif
+          ENDDO
+          GOTO 601
+501       STOP'Problem with profC2H2 file'
+601       CONTINUE
+          CLOSE(15)
+        endif
+
+        lectC2H6 = 0 ! default value for lectC2H6
+        call getin("lectC2H6",lectC2H6)
+        
+        if (lectC2H6 == 1) then
+          OPEN(15,file='profC2H6',status='old',form='formatted',err=502)
+          DO iq = 1, nq
+            if (iq.eq.i_C2H6_mugas) then
+              DO ilayer=1,nlayer
+                READ (15,*) q(ilayer,iq)
+              ENDDO
+            endif
+          ENDDO
+          GOTO 602
+502       STOP'Problem with profC2H6 file'
+602       CONTINUE
+          CLOSE(15)
+        endif
+
+        lectC4H2 = 0 ! default value for lectC4H2
+        call getin("lectC4H2",lectC4H2)
+        
+        if (lectC4H2 == 1) then
+          OPEN(15,file='profC4H2',status='old',form='formatted',err=503)
+          DO iq = 1, nq
+            if (iq.eq.i_C4H2_mugas) then
+              DO ilayer=1,nlayer
+                READ (15,*) q(ilayer,iq)
+              ENDDO
+            endif
+          ENDDO
+          GOTO 603
+503       STOP'Problem with profC4H2 file'
+603       CONTINUE
+          CLOSE(15)
+        endif
+
+        lectC6H6 = 0 ! default value for lectC6H6
+        call getin("lectC6H6",lectC6H6)
+        
+        if (lectC6H6 == 1) then
+          OPEN(15,file='profC6H6',status='old',form='formatted',err=504)
+          DO iq = 1, nq
+            if (iq.eq.i_C6H6_mugas) then
+              DO ilayer=1,nlayer
+                READ (15,*) q(ilayer,iq)
+              ENDDO
+            endif
+          ENDDO
+          GOTO 604
+504       STOP'Problem with profC6H6 file'
+604       CONTINUE
+          CLOSE(15)
+        endif
+
+        lectHCN = 0 ! default value for lectHCN
+        call getin("lectHCN",lectHCN)
+        
+        if (lectHCN == 1) then
+          OPEN(15,file='profHCN',status='old',form='formatted',err=505)
+          DO iq = 1, nq
+            if (iq.eq.i_HCN_mugas) then
+              DO ilayer=1,nlayer
+                READ (15,*) q(ilayer,iq)
+              ENDDO
+            endif
+          ENDDO
+          GOTO 605
+505       STOP'Problem with profHCN file'
+605       CONTINUE
+          CLOSE(15)
+        endif
+      endif ! end of callmufi.and.callmuclouds