Changeset 4081 for trunk/LMDZ.GENERIC/libf/phygeneric/radcommon_h.F90

Timestamp:

Feb 24, 2026, 9:59:11 AM (11 days ago)

Author:

emillour

Message:

Generic PCM:
OpenMP bug fix in "rad_correlatedk_stellar_spectrum", all intermediate computations
should be done by all OpenMp? threads.
While at it did some cleanup:

• added some OMP threadprivate statements in radcommon_h.F90 (not alwyas necessary, but best practice is that all saved variables be threadprivate)
• turned rad_blackbody.F into a module and modernized routines.
• use clearer stategy (wrt OpenMP) in "rad_correlatedk_init_thermal.F90" and "rad_correlatedk_init_stellar.F90": have the master read in file and data and then broadcast to all cores.

EM

File:

• trunk/LMDZ.GENERIC/libf/phygeneric/radcommon_h.F90 (modified) (1 diff)

trunk/LMDZ.GENERIC/libf/phygeneric/radcommon_h.F90

@@ -63 +63 @@
 !
 
-      REAL*8 BWNI(L_NSPECTI+1), WNOI(L_NSPECTI), DWNI(L_NSPECTI), WAVEI(L_NSPECTI) !BWNI read by master in rad_correlatedk_init_thermal
-      REAL*8 BWNV(L_NSPECTV+1), WNOV(L_NSPECTV), DWNV(L_NSPECTV), WAVEV(L_NSPECTV) !BWNV read by master in rad_correlatedk_init_stellar
+      REAL*8 BWNI(L_NSPECTI+1) !BWNI read by master in rad_correlatedk_init_thermal
+!$OMP THREADPRIVATE(BWNI)
+      REAL*8 WNOI(L_NSPECTI), DWNI(L_NSPECTI), WAVEI(L_NSPECTI) 
+!$OMP THREADPRIVATE(WNOI,DWNI,WAVEI)
+      REAL*8 BWNV(L_NSPECTV+1) !BWNV read by master in rad_correlatedk_init_stellar
+!$OMP THREADPRIVATE(BWNV)
+      REAL*8 WNOV(L_NSPECTV), DWNV(L_NSPECTV), WAVEV(L_NSPECTV)
       REAL*8 STELLARF(L_NSPECTV)
-!$OMP THREADPRIVATE(WNOI,DWNI,WAVEI,&
-        !$OMP WNOV,DWNV,WAVEV,&
-        !$OMP STELLARF)
+!$OMP THREADPRIVATE(WNOV,DWNV,WAVEV,STELLARF)
 
       REAL*8 blami(L_NSPECTI+1)