Changeset 3936 for trunk/LMDZ.PLUTO/libf/phypluto/tracer_h.F90

Timestamp:

Oct 25, 2025, 3:21:13 PM (2 months ago)

Author:

tbertrand

Message:

PLUTO PCM : correcting a bug in hazecloud (wrong lyman alpha fluxes due to mu0 being negative during nighttime) + cleaning routines
TB

File:

• trunk/LMDZ.PLUTO/libf/phypluto/tracer_h.F90 (modified) (3 diffs)

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

@@ -12 +12 @@
 
        integer, save :: nqtot  ! total number of tracers
-       integer, save :: nesp   ! number of species in the chemistry
-!$OMP THREADPRIVATE(nqtot,nesp)
+!$OMP THREADPRIVATE(nqtot)
 
        logical :: moderntracdef=.false. ! Standard or modern traceur.def
@@ -25 +24 @@
        real, save, allocatable :: rho_q(:)       ! tracer densities (kg.m-3)
        real, save, allocatable :: qext(:)        ! Single Scat. Extinction coeff at 0.67 um
-       real, save, allocatable :: alpha_lift(:)  ! saltation vertical flux/horiz flux ratio (m-1)
-       real, save, allocatable :: alpha_devil(:) ! lifting coeeficient by dust devil
        real, save, allocatable :: qextrhor(:)    ! Intermediate for computing opt. depth from q
 
@@ -41 +38 @@
        integer,save :: nmono
        real,save :: ref_r0        ! for computing reff=ref_r0*r0 (in log.n. distribution)
-!$OMP THREADPRIVATE(noms,mmol,aki,cpi,radius,rho_q,qext,alpha_lift,alpha_devil,qextrhor, &
+!$OMP THREADPRIVATE(noms,mmol,aki,cpi,radius,rho_q,qext,qextrhor, &
         !$OMP varian,r3n_q,rho_dust,rho_ice,rho_n2,lw_n2,ref_r0)
 
-       integer, save, allocatable :: is_chim(:) ! 1 if tracer used in chemistry, else 0
        integer, save, allocatable :: is_rad(:)  ! 1 if   ""    ""  in radiative transfer, else 0
-!$OMP THREADPRIVATE(is_chim,is_rad)
+!$OMP THREADPRIVATE(is_rad)
 
        integer, save, allocatable :: is_recomb(:)      ! 1 if tracer used in recombining scheme, else 0 (if 1, must have is_rad=1)