Changeset 4077 for trunk/LMDZ.GENERIC/libf/phygeneric/rad_correlatedk_opacities_stellar.F90

Timestamp:

Feb 18, 2026, 10:28:23 AM (7 weeks ago)

Author:

mturbet

Message:

renaming and merging of radiative routines in phygeneric

File:

• trunk/LMDZ.GENERIC/libf/phygeneric/rad_correlatedk_opacities_stellar.F90 (moved) (moved from trunk/LMDZ.GENERIC/libf/phygeneric/optcv.F90) (10 diffs)

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

@@ -1 @@
-MODULE optcv_mod
+MODULE rad_correlatedk_opacities_stellar_mod
 
 IMPLICIT NONE
@@ -5 +5 @@
 CONTAINS
 
-SUBROUTINE OPTCV(DTAUV,TAUV,TAUCUMV,PLEV,  &
+SUBROUTINE rad_correlatedk_opacities_stellar(DTAUV,TAUV,TAUCUMV,PLEV,  &
      QXVAER,QSVAER,GVAER,WBARV,COSBV,       &
      TAUAERO,TMID,PMID,TAUGSURF,QVAR,MUVAR,FRACVAR)
@@ -16 +16 @@
   use callkeys_mod, only: kastprof,continuum,graybody,callgasvis,varspec, &
                           rayleigh
-  use recombin_corrk_mod, only: corrk_recombin, gasv_recomb
+  use rad_correlatedk_online_recombination_mod, only: corrk_recombin, gasv_recomb
   use tpindex_mod, only: tpindex
-  use interpolate_continuum_mod, only: interpolate_continuum
-  use calc_rayleigh_mod, only: calc_rayleigh
+  use rad_correlatedk_continuum_interpolation_mod, only: rad_correlatedk_continuum_interpolation
+  use rad_correlatedk_rayleigh_scattering_opacity_mod, only: rad_correlatedk_rayleigh_scattering_opacity
 
   implicit none
@@ -138 +138 @@
      else
         dz(k) = dpr(k)*R*TMID(K)/(glat_ig*PMID(K))*mugaz/muvar(k)
-        U(k)  = Cmk*DPR(k)*mugaz/muvar(k)     ! only Cmk line in optci.F  
+        U(k)  = Cmk*DPR(k)*mugaz/muvar(k)     ! only Cmk line in rad_correlatedk_opacities_thermal.F  
             !JL13 the mugaz/muvar factor takes into account water meanmolecular weight if water is present
      endif
@@ -154 +154 @@
             !   but visible does not handle very well diffusion in first layer.
             !   The tauaero and tauray are thus set to 0 (a small value for rayleigh because the code crashes otherwise)
-            !   in the 4 first semilayers in optcv, but not optci.
+            !   in the 4 first semilayers in rad_correlatedk_opacities_stellar, but not rad_correlatedk_opacities_thermal.
             !   This solves random variations of the sw heating at the model top. 
   do iaer=1,naerkind
@@ -170 +170 @@
 
       if(rayleigh) then
-         call calc_rayleigh(QVAR,MUVAR,PMID,TMID,TAURAY)
+         call rad_correlatedk_rayleigh_scattering_opacity(QVAR,MUVAR,PMID,TMID,TAURAY)
       else
-         print*,'setspv: No Rayleigh scattering, check for NaN in output!'
+         print*,'rad_correlatedk_init_stellar: No Rayleigh scattering, check for NaN in output!'
          do NW=1,L_NSPECTV
             TAURAY(:,NW) = 1E-16
@@ -243 +243 @@
                      ((igas .eq. igas_CO2) .and. (jgas .eq. igas_CH4)) ) then
 
-                  call interpolate_continuum('',igas,jgas,'VI',nw,T_cont,p_cont,p_cross,dtemp,.false.)
+                  call rad_correlatedk_continuum_interpolation('',igas,jgas,'VI',nw,T_cont,p_cont,p_cross,dtemp,.false.)
 
                 endif
@@ -268 +268 @@
            
               ! JVO 2017 : added tmpk because the repeated calls to gasi/v increased dramatically
-              ! the execution time of optci/v -> ~ factor 2 on the whole radiative
+              ! the execution time of rad_correlatedk_opacities_thermal/v -> ~ factor 2 on the whole radiative
               ! transfer on the tested simulations !
 
@@ -334 +334 @@
             !JL18 It seems to be good to have aerosols in the first "radiative layer" of the gcm in the IR
             !   but not in the visible
-            !   The tauaero is thus set to 0 in the 4 first semilayers in optcv, but not optci.
+            !   The tauaero is thus set to 0 in the 4 first semilayers in rad_correlatedk_opacities_stellar, but not rad_correlatedk_opacities_thermal.
             !   This solves random variations of the sw heating at the model top. 
   do iaer=1,naerkind
@@ -407 +407 @@
 
 
-end subroutine optcv
-
-END MODULE optcv_mod
+end subroutine rad_correlatedk_opacities_stellar
+
+END MODULE rad_correlatedk_opacities_stellar_mod