Changeset 3299 for trunk/LMDZ.GENERIC/libf

Timestamp:

Apr 12, 2024, 4:48:27 PM (8 months ago)

Author:

emillour

Message:

Generic PCM:

• add some auxiliary functions in module tracer_h.F90: is_known_tracer(tname) returns ".true." if "tname" is a known tracer tracer_index(tname) returns the tracer index of tracer "tname"
• add the possibility of "volcanic outgasing" as point sources. controled by a "callvolcano" flag (default .false.) details about source location and injection details (frequency, length) need be provided in a "volcano.def" ASCII file read at runtime. (see routine "read_volcano" in volcano.F90)

AB+EM

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• callkeys_mod.F90 (modified) (1 diff)
• inifis_mod.F90 (modified) (1 diff)
• physiq_mod.F90 (modified) (4 diffs)
• tracer_h.F90 (modified) (1 diff)
• volcano.F90 (added)

trunk/LMDZ.GENERIC/libf/phystd/callkeys_mod.F90

@@ -6 +6 @@
       logical,save :: calladj,calltherm,co2cond,callsoil
 !$OMP THREADPRIVATE(calladj,calltherm,co2cond,callsoil)
+
+      ! do we have volcanoes (injecting tracers in atmosphere) ?
+      logical,save :: callvolcano
+!$OMP THREADPRIVATE(callvolcano)
+
       logical,save :: season,diurnal,tlocked,rings_shadow,lwrite
 !$OMP THREADPRIVATE(season,diurnal,tlocked,rings_shadow,lwrite)

trunk/LMDZ.GENERIC/libf/phystd/inifis_mod.F90

@@ -503 +503 @@
 !       ": ok_slab_heat_transp = ",ok_slab_heat_transp
 
+! Volcanoes (as point sources of tracers)
+      if (is_master) write(*,*) trim(rname)//": Erupt volcano ?"
+      callvolcano=.false. ! default value
+      call getin_p("callvolcano",callvolcano)
+      if (is_master) write(*,*) trim(rname)//": callvolcano = ",callvolcano
+      
+      ! Sanity check: volcano requires tracers
+      if (callvolcano.and.(.not.tracer)) then
+        call abort_physic(rname,&
+             " Error: volcano option requires using tracers",1)
+      endif
+
 ! Photochemistry and chemistry in the thermosphere
 

trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90

@@ -66 +66 @@
                               ok_slab_ocean, photochem, rings_shadow, rmean, &
                               season, sedimentation,generic_condensation, &
-                              specOLR, callthermos, &
+                              specOLR, callthermos, callvolcano, &
                               startphy_file, testradtimes, tlocked, &
                               tracer, UseTurbDiff, water, watercond, &
@@ -77 +77 @@
       use phys_state_var_mod
       use callcorrk_mod, only: callcorrk
-      use conduction_mod
+      use conduction_mod, only: conduction
+      use volcano_mod, only: volcano
       use pindex_mod, only: pindex
       use vdifc_mod, only: vdifc
@@ -358 +359 @@
       REAL,allocatable,save :: zdqschim(:,:)  ! Calchim_asis routine
 !$OMP THREADPRIVATE(zdqchim,zdqschim)
+      real zdqvolc(ngrid,nlayer,nq) ! injection by volcanoes (kg/kg_of_air/s)
 
       REAL array_zero1(ngrid)
@@ -1481 +1483 @@
       if (tracer) then
 
+  ! -----------------------------------------------------
+  !   VI.0. Volcanoes injecting tracers in the atmosphere
+  ! -----------------------------------------------------
+        if (callvolcano) then
+          call volcano(ngrid,nlayer,pplev,pu,pv,pt,zzlev,zdqvolc,nq,massarea,&
+                       zday,ptimestep)
+          pdq(1:ngrid,1:nlayer,1:nq) = pdq(1:ngrid,1:nlayer,1:nq) + &
+                                       zdqvolc(1:ngrid,1:nlayer,1:nq)
+        endif
+
   ! ---------------------
   !   VI.1. Water and ice

trunk/LMDZ.GENERIC/libf/phystd/tracer_h.F90

@@ -115 +115 @@
        !$OMP igcm_c2h3,igcm_c2h4,igcm_c2h6,igcm_ch2co,igcm_ch3co,igcm_hcaer)
 
+       contains
+       
+       ! some basic auxiliary functions for tracers
+       
+       logical function is_known_tracer(tracer_name)
+         ! checks if tracer_name is a know tracer
+         ! returns .true. if yes and .false. otherwise
+         implicit none
+         character(len=*),intent(in) :: tracer_name
+         integer :: iq
+         ! initialization
+         is_known_tracer=.false.
+         do iq=1,nqtot
+           if( trim(tracer_name) == trim(noms(iq)) ) then
+             is_known_tracer=.true.
+             exit
+           endif
+         enddo
+       end function is_known_tracer
+       
+       integer function tracer_index(tracer_name)
+         ! gets the index of the tracer tracer_name
+         ! returns index or 0 if unknown tracer
+         implicit none
+         character(len=*),intent(in) :: tracer_name
+         integer :: iq
+         ! initialization
+         tracer_index=0
+         do iq=1,nqtot
+           if( trim(tracer_name) == trim(noms(iq)) ) then
+             tracer_index=iq
+             exit
+           endif
+         enddo
+       end function tracer_index
+        
        end module tracer_h