Changeset 2611 for trunk/LMDZ.MARS/libf/aeronomars

Timestamp:

Jan 18, 2022, 3:32:49 PM (3 years ago)

Author:

romain.vande

Message:

LMDZ_MARS RV : Open_MP; Reading files in parallel for moldiff parametrisation
(callmoldiff = .true.)

File:

• trunk/LMDZ.MARS/libf/aeronomars/moldiff_red.F90 (modified) (5 diffs)

trunk/LMDZ.MARS/libf/aeronomars/moldiff_red.F90

@@ -6 +6 @@
 use geometry_mod, only: cell_area
 use planetwide_mod, only: planetwide_sumval
+USE mod_phys_lmdz_para, only: is_master,bcast
 
 implicit none
@@ -83 +84 @@
       integer,save :: il0  
 
+!$OMP THREADPRIVATE(qq,qnew,qint,FacMass,rhoK,rhokinit,wi,Wad,Uthermal,Lambdaexo,Hspecie,Mtot1,Mtot2,Mraf1,Mraf2,i_h2,i_h,i_d,il0)
+
 ! Tracer indexes in the GCM:
 !      integer,save :: g_co2=0
@@ -114 +117 @@
       real,save :: step
 
+!$OMP THREADPRIVATE(ncompdiff,gcmind,firstcall,dij,step)
 
 ! Initializations at first call
@@ -176 +180 @@
 
 ! find vertical index above which diffusion is computed
-
+      if(is_master) then
         do l=1,nlayer
         if (pplay(1,l) .gt. Pdiff) then
@@ -183 +187 @@
         enddo
         il0=il0+1
+      endif !(is_master)
+      CALL bcast(il0)
         print*,'vertical index for diffusion',il0,pplay(1,il0)