Changeset 4127 for trunk/LMDZ.GENERIC

Timestamp:

Mar 13, 2026, 10:58:37 AM (4 weeks ago)

Author:

aslmd

Message:

minor fortran coding error in the Rayleigh scattering routine and more major bug in missing albedo initialization in Mesoscale mode

Location:

trunk/LMDZ.GENERIC/libf/phygeneric

Files:

• physiq_mod.F90 (modified) (1 diff)
• rad_correlatedk_rayleigh_scattering_opacity.F90 (modified) (5 diffs)

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

@@ -620 +620 @@
 
          day_ini = pday
+         DO nw=1,L_NSPECTV
+            albedo(:,nw)=albedodat(:)
+         ENDDO
 #endif
 

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

@@ -48 +48 @@
       real, intent(out) :: tauray(L_LEVELS,L_NSPECTV)
       real*8 wl,wn
-      integer N,Nfine,ifine,igas,k
-      parameter(Nfine=500.0)
+      integer N,ifine,igas,k
+      integer, parameter :: Nfine = 500
       real*8 :: Fk ! King factor for the depolarization
       real*8 :: ng(L_LEVELS) ! real refractive index 
@@ -131 +131 @@
          endif
          tauvari(igas,:) = 0.
+         tauconsti(igas,:) = 0.
       enddo
       
@@ -149 +150 @@
          bwidth = (10000.0/BWNV(N)) - (10000.0/BWNV(N+1))
          do ifine=1,Nfine
-            wl=bstart+dble(ifine)*bwidth/Nfine
-            wn=BWNV(N)+dble(ifine)*(BWNV(N+1)-BWNV(N))/Nfine
-
-            tauvar(:)=0.0
+            wl=bstart+dble(ifine)*bwidth/dble(Nfine)
+            wn=BWNV(N)+dble(ifine)*(BWNV(N+1)-BWNV(N))/dble(Nfine)
+
+            tauvar(:)=0.0d0
             do igas=1,ngasmx
                if (maxval(mass_frac(igas,:)).ge.1e-2) then
@@ -319 +320 @@
                  ! pmid*scalep -> mbar to Pa
                  ! muvar/1000 -> g/mol to kg/mol
-               
                  tauvar(:)=tauvar(:)+tauconsti(igas,:)*tauvari(igas,:)
-                 
+          
                endif !greater than 0.01
 
@@ -327 +327 @@
 
             call rad_blackbody_planck_law_wavelength(dble(wl*1e-6),dble(tstellar),df)
-            df=df*bwidth/Nfine
+            df=df*bwidth/dble(Nfine)
             tauwei=tauwei+df
             tausum(:)=tausum(:)+tauvar(:)*df