Changeset 3360 for trunk/LMDZ.GENERIC/libf/phystd/newsedim.F

Timestamp:

Jun 5, 2024, 12:30:04 PM (6 months ago)

Author:

emillour

Message:

Generic PCM:
OpenMP fix in newsedim, add a missing THREADPRIVATE for "molrad" variable.
While at it make sedim and callsedim modules.
EM

File:

• trunk/LMDZ.GENERIC/libf/phystd/newsedim.F (modified) (4 diffs)

trunk/LMDZ.GENERIC/libf/phystd/newsedim.F

@@ +1 @@
+      MODULE newsedim_mod
+      
+      IMPLICIT NONE
+      
+      CONTAINS
+
       SUBROUTINE newsedim(ngrid,nlay,naersize,ptimestep,
      &  pplev,masse,epaisseur,pt,rd,rho,pqi,wq,iq)
@@ -4 +10 @@
       use ioipsl_getin_p_mod, only: getin_p
       use comcstfi_mod, only: r, g
-      use gases_h
+      use gases_h, only: gfrac, gnom, ngasmx
+      use gases_h, only: igas_CH4, igas_CO2, igas_H2, igas_H2O
+      use gases_h, only: igas_He, igas_N2
       use tracer_h, only : igcm_h2o_ice
       use watercommon_h, only: T_h2O_ice_liq,T_h2O_ice_clouds     
@@ -69 +77 @@
 c     Effective gas molecular radius (m)
       real,save :: molrad
+!$OMP THREADPRIVATE(molrad)
 
 c     local and saved variable
@@ -282 +291 @@
 c    &                wq(1,6),wq(1,7),pqi(1,6)
 
-      END
+      END SUBROUTINE newsedim
+      
+      END MODULE newsedim_mod