Changeset 4081 for trunk/LMDZ.GENERIC/libf/phygeneric/rad_blackbody.F

Timestamp:

Feb 24, 2026, 9:59:11 AM (5 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/rad_blackbody.F (modified) (1 diff)

trunk/LMDZ.GENERIC/libf/phygeneric/rad_blackbody.F

@@ +1 @@
+      module rad_blackbody_mod
+      
+      implicit none
+      
+      contains
+      
       subroutine rad_blackbody_planck_law_wavelength(blalong,blat,blae)
 
-      implicit double precision (a-h,o-z)
+      implicit none
+      
+      double precision, intent(in) :: blalong
+      double precision, intent(in) :: blat
+      double precision, intent(out) :: blae
+      
+!      double precision :: sigma,pi,c0,h,cbol,rind,c,c1,c2
 
       ! physical constants
-      sigma=5.670374D-8
-      pi=datan(1.d0)*4.d0
-      c0=2.997925d+08
-      h=6.62607d-34
-      cbol=1.380649d-23
-      rind=1.d0
-      c=c0/rind
-      c1=h*(c**2)
-      c2=h*c/cbol
+      double precision, parameter :: sigma=5.670374D-8
+      double precision, parameter :: pi=atan(1.d0)*4.d0
+      double precision, parameter :: c0=2.997925d+08
+      double precision, parameter :: h=6.62607d-34
+      double precision, parameter :: cbol=1.380649d-23
+      double precision, parameter :: rind=1.d0
+      double precision, parameter :: c=c0/rind
+      double precision, parameter :: c1=h*(c**2)
+      double precision, parameter :: c2=h*c/cbol
 
 
-      blae=2.d0*pi*c1/blalong**5/(dexp(c2/blalong/blat)-1.d0)
+      blae=2.d0*pi*c1/blalong**5/(exp(c2/blalong/blat)-1.d0)
 
-
-      return
-      end
+      end subroutine rad_blackbody_planck_law_wavelength
 
       subroutine rad_blackbody_planck_law_wavenumber(blalong,blat,blae)
 
-      implicit double precision (a-h,o-z)
+      implicit none
+
+      double precision, intent(in) :: blalong
+      double precision, intent(in) :: blat
+      double precision, intent(out) :: blae
+      
+!      double precision :: sigma,pi,c0,h,cbol,rind,c,c1,c2
 
       ! physical constants
-      sigma=5.670374D-8
-      pi=datan(1.d0)*4.d0
-      c0=2.997925d+08
-      h=6.62607d-34
-      cbol=1.380649d-23
-      rind=1.d0
-      c=c0/rind
-      c1=h*(c**2)
-      c2=h*c/cbol
+      double precision, parameter :: sigma=5.670374D-8
+      double precision, parameter :: pi=atan(1.d0)*4.d0
+      double precision, parameter :: c0=2.997925d+08
+      double precision, parameter :: h=6.62607d-34
+      double precision, parameter :: cbol=1.380649d-23
+      double precision, parameter :: rind=1.d0
+      double precision, parameter :: c=c0/rind
+      double precision, parameter :: c1=h*(c**2)
+      double precision, parameter :: c2=h*c/cbol
 
 
-      blae=2.d0*pi*c1*blalong**3/(dexp(c2*blalong/blat)-1.d0)
+      blae=2.d0*pi*c1*blalong**3/(exp(c2*blalong/blat)-1.d0)
 
-
-      return
-      end
+      end subroutine rad_blackbody_planck_law_wavenumber
+      
+      end module rad_blackbody_mod