Changeset 2616 for trunk/LMDZ.MARS/libf/phymars/callsedim_mod.F

Timestamp:

Jan 18, 2022, 4:28:31 PM (3 years ago)

Author:

romain.vande

Message:

LMDZ_MARS RV : Open_MP;
Put all the "save" variables as "!$OMP THREADPRIVATE" in phymars.
The code can now be tested, see README for more info

File:

• trunk/LMDZ.MARS/libf/phymars/callsedim_mod.F (modified) (4 diffs)

trunk/LMDZ.MARS/libf/phymars/callsedim_mod.F

@@ -102 +102 @@
 !                                    !   CCNs (m)
       real,save :: beta ! correction for the shape of the ice particles (cf. newsedim)
+      
+!$OMP THREADPRIVATE(beta)
+      
 c     for ice radius computation
       REAL Mo,No
@@ -128 +131 @@
       real,parameter :: rdimin=1.e-8 ! 1.e-7
       real,parameter :: rdimax=1.e-4
+      
+!$OMP THREADPRIVATE(rd,rdi)
 
 c       2) Second size distribution for the log-normal integration
@@ -138 +143 @@
       real sigma0
 
+!$OMP THREADPRIVATE(rr)
+      
 c       3) Other local variables used in doubleq
 
@@ -164 +171 @@
       INTEGER,SAVE :: iccnco2_meteor_mass ! index of tracer containing CCN number
       INTEGER,SAVE :: ico2_ice ! index of tracer containing CCN number
+      
+!$OMP THREADPRIVATE(idust_mass,idust_number,iccn_mass,iccn_number)
+!$OMP THREADPRIVATE(istormdust_mass,istormdust_number,itopdust_mass)
+!$OMP THREADPRIVATE(itopdust_number)
+!$OMP THREADPRIVATE(iccnco2_number,iccnco2_mass,iccnco2_h2o_number)
+!$OMP THREADPRIVATE(iccnco2_h2o_mass_ice)
+!$OMP THREADPRIVATE(iccnco2_h2o_mass_ccn,ico2_ice)
 
 
       LOGICAL,SAVE :: firstcall=.true.
+      
+!$OMP THREADPRIVATE(firstcall)