Index: trunk/LMDZ.MARS/libf/phymars/convadj.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/convadj.F	(revision 2584)
+++ trunk/LMDZ.MARS/libf/phymars/convadj.F	(revision 2586)
@@ -81,4 +81,7 @@
       LOGICAL vtest(ngrid),down
       logical,save :: firstcall=.true.
+
+
+!$OMP THREADPRIVATE(ico2,A,B,firstcall)
 
 !     for conservation test
Index: trunk/LMDZ.MARS/libf/phymars/nltecool.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/nltecool.F	(revision 2584)
+++ trunk/LMDZ.MARS/libf/phymars/nltecool.F	(revision 2586)
@@ -114,4 +114,7 @@
 
       logical    firstcall
+      
+!$OMP THREADPRIVATE(firstcall,ef1,ef2,co2vmr,n2covmr,o3pvmr,pnb)      
+      
       data       firstcall/.true./
       save       firstcall,ef1,ef2,co2vmr,n2covmr,o3pvmr,pnb
Index: trunk/LMDZ.MARS/libf/phymars/pbl_parameters.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/pbl_parameters.F	(revision 2584)
+++ trunk/LMDZ.MARS/libf/phymars/pbl_parameters.F	(revision 2586)
@@ -77,4 +77,7 @@
       REAL karman,nu
       DATA karman,nu/.41,0.001/
+
+!$OMP THREADPRIVATE(karman,nu)
+
       SAVE karman,nu
 
Index: trunk/LMDZ.MARS/libf/phymars/read_dust_scenario.F90
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/read_dust_scenario.F90	(revision 2584)
+++ trunk/LMDZ.MARS/libf/phymars/read_dust_scenario.F90	(revision 2586)
@@ -8,4 +8,5 @@
 use dust_param_mod, only: odpref, dustscaling_mode
 use planete_h, only: year_day
+USE mod_phys_lmdz_transfert_para, ONLY: bcast
 implicit none
 
@@ -80,4 +81,7 @@
      filename="dust_tes_MY26.nc"
    endif
+
+
+!$OMP MASTER
    
    ierr=nf90_open(trim(datadir)//"/"//trim(filename),NF90_NOWRITE,nid)
@@ -107,7 +111,14 @@
    ierr=nf90_inquire_dimension(nid,londim,len=lonlen)
 
+!$OMP END MASTER
+
+   call bcast(lonlen)
+   call bcast(latlen)
+   call bcast(timelen)
 
    if (.not.allocated(tautes)) allocate(tautes(lonlen,latlen,timelen))
    if (.not.allocated(lat)) allocate(lat(latlen), lon(lonlen), time(timelen))
+
+!$OMP MASTER
 
    ! "dustop" if loading visible extinction opacity
@@ -176,4 +187,11 @@
 
    ierr=nf90_close(nid)
+
+!$OMP END MASTER
+
+   call bcast(tautes)
+   call bcast(time)
+   call bcast(lat)
+   call bcast(lon)
 
 endif ! of if (firstcall)