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

Timestamp:

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

Author:

mturbet

Message:

renaming and merging of radiative routines in phygeneric

File:

• trunk/LMDZ.GENERIC/libf/phygeneric/rad_correlatedk.F90 (moved) (moved from trunk/LMDZ.GENERIC/libf/phygeneric/callcorrk.F90) (26 diffs)

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

@@ -1 @@
-MODULE callcorrk_mod
+MODULE rad_correlatedk_mod
 
 IMPLICIT NONE
@@ -5 +5 @@
 CONTAINS
 
-      subroutine callcorrk(ngrid,nlayer,pq,nq,qsurf,zls,       &
+      subroutine rad_correlatedk(ngrid,nlayer,pq,nq,qsurf,zls,       &
           albedo,albedo_equivalent,emis,mu0,pplev,pplay,pt,    & 
           tsurf,fract,dist_star,aerosol,muvar,                 &
@@ -27 +27 @@
       use ioipsl_getin_p_mod, only: getin_p
       use gases_h, only: ngasmx
-      use radii_mod, only : su_aer_radii,co2_reffrad,h2o_reffrad,dust_reffrad,h2so4_reffrad,back2lay_reffrad
-      use aerosol_mod, only : iaero_co2,iaero_h2o,iaero_dust,iaero_h2so4, &
+      use aerosol_radius, only : aerosol_radius_init,aerosol_radius_co2,h2o_reffrad, &
+                    aerosol_radius_dust,aerosol_radius_h2so4,aerosol_radius_back2lay
+      use aerosol_global_variables , only : iaero_co2,iaero_h2o,iaero_dust,iaero_h2so4, &
                               iaero_back2lay, iaero_aurora,               &
                               iaero_venus1,iaero_venus2,iaero_venus2p,    &
                               iaero_venus3,iaero_venusUV
-      use aeropacity_mod, only: aeropacity
-      use aeroptproperties_mod, only: aeroptproperties
+      use aerosol_opacity_mod, only: aerosol_opacity
+      use aerosol_optical_properties_mod, only: aerosol_optical_properties
       use tracer_h, only: igcm_h2o_ice, igcm_h2o_vap, igcm_co2_ice
       use tracer_h, only: constants_epsi_generic
@@ -41 +42 @@
                               CLFvarying,tplanckmin,tplanckmax,global1d, &
                               generic_condensation, aerovenus, nvarlayer, varspec
-      use optcv_mod, only: optcv
-      use optci_mod, only: optci
-      use sfluxi_mod, only: sfluxi
-      use sfluxv_mod, only: sfluxv
-      use recombin_corrk_mod, only: corrk_recombin, call_recombin
+      use rad_correlatedk_opacities_stellar_mod, & 
+        only: rad_correlatedk_opacities_stellar
+      use rad_correlatedk_opacities_thermal_mod, &
+        only: rad_correlatedk_opacities_thermal
+      use rad_correlatedk_fluxes_thermal_mod, &
+        only: rad_correlatedk_fluxes_thermal
+      use rad_correlatedk_fluxes_stellar_mod, &
+        only: rad_correlatedk_fluxes_stellar
+      use rad_correlatedk_online_recombination_mod, &
+        only: corrk_recombin, &
+        rad_correlatedk_recombination_main
       use pindex_mod, only: pindex
       use generic_cloud_common_h, only: Psat_generic, epsi_generic
@@ -111 +118 @@
       REAL,INTENT(OUT) :: OSR_nu(ngrid,L_NSPECTV)        ! Outgoing SW radiation in each band (Normalized to the band width (W/m2/cm-1).
       REAL,INTENT(OUT) :: GSR_nu(ngrid,L_NSPECTV)        ! Surface SW radiation in each band (Normalized to the band width (W/m2/cm-1).
-      REAL,INTENT(OUT) :: tau_col(ngrid)                 ! Diagnostic from aeropacity.
+      REAL,INTENT(OUT) :: tau_col(ngrid)                 ! Diagnostic from aerosol_opacity.
       REAL,INTENT(OUT) :: albedo_equivalent(ngrid)       ! Spectrally Integrated Albedo. For Diagnostic. By MT2015
       REAL,INTENT(OUT) :: totcloudfrac(ngrid)            ! Column Fraction of clouds (%).
@@ -148 +155 @@
       REAL*8 tlevrad(L_LEVELS),plevrad(L_LEVELS)
 
-      ! Optical values for the optci/cv subroutines
+      ! Optical values for the rad_correlatedk_opacities_thermal/cv subroutines
       REAL*8 stel(L_NSPECTV),stel_fract(L_NSPECTV)
       ! NB: Arrays below are "save" to avoid reallocating them at every call
@@ -200 +207 @@
       character(len=10) :: tmp2
       character(len=100) :: message
-      character(len=10),parameter :: subname="callcorrk"
+      character(len=10),parameter :: subname="rad_correlatedk"
 
       ! For fixed water vapour profiles.
@@ -242 +249 @@
       if(firstcall) then
 
-        ! test on allocated necessary because of CLFvarying (two calls to callcorrk in physiq)
+        ! test on allocated necessary because of CLFvarying (two calls to rad_correlatedk in physiq)
         if(.not.allocated(QVISsQREF3d)) then
           allocate(QVISsQREF3d(ngrid,nlayer,L_NSPECTV,naerkind))
@@ -308 +315 @@
         endif
 
-         !!! ALLOCATED instances are necessary because of CLFvarying (strategy to call callcorrk twice in physiq...)
+         !!! ALLOCATED instances are necessary because of CLFvarying (strategy to call rad_correlatedk twice in physiq...)
          IF(.not.ALLOCATED(QREFvis3d))THEN
            ALLOCATE(QREFvis3d(ngrid,nlayer,naerkind), stat=ok)
@@ -343 +350 @@
 #endif
 
-         call su_aer_radii(ngrid,nlayer,reffrad,nueffrad)
+         call aerosol_radius_init(ngrid,nlayer,reffrad,nueffrad)
          
          
@@ -351 +358 @@
 
       !this block is now done at firstcall of physiq_mod
-         ! print*, "callcorrk: Correlated-k data base folder:",trim(datadir)
+         ! print*, "rad_correlatedk: Correlated-k data base folder:",trim(datadir)
          ! call getin_p("corrkdir",corrkdir)
          ! print*, "corrkdir = ",corrkdir
@@ -359 +366 @@
          ! banddir=trim(adjustl(corrkdir))//'/'//trim(adjustl(banddir))
 
-         ! call setspi            ! Basic infrared properties.
-         ! call setspv            ! Basic visible properties.
-         ! call sugas_corrk       ! Set up gaseous absorption properties.
-         ! call suaer_corrk       ! Set up aerosol optical properties.
+         ! call rad_correlatedk_init_thermal            ! Basic infrared properties.
+         ! call rad_correlatedk_init_stellar            ! Basic visible properties.
+         ! call rad_correlatedk_read_opacity_tables        ! Set up gaseous absorption properties.
+         ! call rad_correlatedk_ini_aerosol       ! Set up aerosol optical properties.
         
 
-         ! now that L_NGAUSS has been initialized (by sugas_corrk)
+         ! now that L_NGAUSS has been initialized (by rad_correlatedk_read_opacity_tables )
          ! allocate related arrays
          if(.not.allocated(dtaui)) then
@@ -446 +453 @@
 
          if((igcm_h2o_vap.eq.0) .and. varactive .and. water)then
-            message='varactive in callcorrk but no h2o_vap tracer.'
+            message='varactive in rad_correlatedk but no h2o_vap tracer.'
             call abort_physic(subname,message,1)
          endif
@@ -485 +492 @@
 
          if ((iaer.eq.iaero_co2).and.tracer.and.(igcm_co2_ice.gt.0)) then ! Treat condensed co2 particles.
-            call co2_reffrad(ngrid,nlayer,nq,pq,reffrad(1,1,iaero_co2))
+            call aerosol_radius_co2(ngrid,nlayer,nq,pq,reffrad(1,1,iaero_co2))
 
             call planetwide_maxval(reffrad(:,:,iaero_co2),maxvalue)
@@ -517 +524 @@
          
          if(iaer.eq.iaero_dust)then
-            call dust_reffrad(ngrid,nlayer,reffrad(1,1,iaero_dust))
+            call aerosol_radius_dust(ngrid,nlayer,reffrad(1,1,iaero_dust))
             if (is_master) then
                print*,'Dust particle size = ',reffrad(1,1,iaer)/1.e-6,' um'
@@ -524 +531 @@
          
          if(iaer.eq.iaero_h2so4)then
-            call h2so4_reffrad(ngrid,nlayer,reffrad(1,1,iaero_h2so4))
+            call aerosol_radius_h2so4(ngrid,nlayer,reffrad(1,1,iaero_h2so4))
             if (is_master) then
                print*,'H2SO4 particle size =',reffrad(1,1,iaer)/1.e-6,' um'
@@ -531 +538 @@
          
           if(iaer.eq.iaero_back2lay)then
-            call back2lay_reffrad(ngrid,reffrad(1,1,iaero_back2lay),nlayer,pplev)
-         endif
-
-         !  For n-layer aerosol size set once for all at firstcall in su_aer_radii
+            call aerosol_radius_back2lay(ngrid,reffrad(1,1,iaero_back2lay),nlayer,pplev)
+         endif
+
+         !  For n-layer aerosol size set once for all at firstcall in aerosol_radius_init
 
 !         if(iaer.eq.iaero_aurora)then
@@ -549 +556 @@
 
       ! Get 3D aerosol optical properties.
-      call aeroptproperties(ngrid,nlayer,reffrad,nueffrad,         &
+      call aerosol_optical_properties(ngrid,nlayer,reffrad,nueffrad,         &
            QVISsQREF3d,omegaVIS3d,gVIS3d,                          &
            QIRsQREF3d,omegaIR3d,gIR3d,                             &
@@ -555 +562 @@
 
       ! Get aerosol optical depths.
-      call aeropacity(ngrid,nlayer,nq,pplay,pplev,pt,pq,zls,aerosol,      &
+      call aerosol_opacity(ngrid,nlayer,nq,pplay,pplev,pt,pq,zls,aerosol,      &
            reffrad,nueffrad,QREFvis3d,QREFir3d,                             & 
            tau_col,cloudfrac,totcloudfrac,clearsky)                
@@ -692 +699 @@
                        (pplev(ig,L_NLAYRAD-k)-pplev(ig,L_NLAYRAD-k+1))
                ! As 'aerosol' is at reference (visible) wavelenght we scale it as
-               ! it will be multplied by qxi/v in optci/v
+               ! it will be multplied by qxi/v in rad_correlatedk_opacities_thermal/v
                temp=aerosol(ig,L_NLAYRAD-k,iaer)/QREFvis3d(ig,L_NLAYRAD-k,iaer)
                tauaero(2*k+2,iaer)=max(temp*pweight,0.d0)
@@ -918 +925 @@
             qvar(1)=qvar(2)
 
-            write(*,*)trim(subname),' :Warning: reducing qvar in callcorrk.'
+            write(*,*)trim(subname),' :Warning: reducing qvar in rad_correlatedk.'
             write(*,*)trim(subname),' :Temperature profile no longer consistent ', &
                    'with saturated H2O. qsat=',satval
@@ -1018 +1025 @@
       ! -- JVO 20 : Also add a sanity test checking that tlevrad is
       !             within Planck function temperature boundaries,
-      !             which would cause gfluxi/sfluxi to crash.
+      !             which would cause rad_correlatedk_fluxes_solver_thermal/rad_correlatedk_fluxes_thermal to crash.
       do k=1,L_LEVELS
 
@@ -1127 +1134 @@
 ! Recombine reference corrk tables if needed - Added by JVO, 2020.
          if (corrk_recombin) then
-           call call_recombin(ig,nlayer,pq(ig,:,:),pplay(ig,:),pt(ig,:),qvar(:),tmid(:),pmid(:))
+           call rad_correlatedk_recombination_main(ig,nlayer,pq(ig,:,:),pplay(ig,:),pt(ig,:),qvar(:),tmid(:),pmid(:))
          endif
 ! ----------------------------------------------------------------
@@ -1149 +1156 @@
             endif
 
-            call optcv(dtauv,tauv,taucumv,plevrad,                 &
-                 qxvaer,qsvaer,gvaer,wbarv,cosbv,tauaero,          &
+            call rad_correlatedk_opacities_stellar(dtauv,    &
+                 tauv,taucumv,plevrad,                       &
+                 qxvaer,qsvaer,gvaer,wbarv,cosbv,tauaero,    &
                  tmid,pmid,taugsurf,qvar,muvarrad,fracvari)
 
-            call sfluxv(dtauv,tauv,taucumv,albv,dwnv,wbarv,cosbv,  &
-                 acosz,stel_fract,                                 &
-                 nfluxtopv,fluxtopvdn,nfluxoutv_nu,nfluxgndv_nu,   &
+            call rad_correlatedk_fluxes_stellar(dtauv,tauv,  & 
+                 taucumv,albv,dwnv,wbarv,cosbv,              &
+                 acosz,stel_fract,                           &
+                 nfluxtopv,fluxtopvdn,                       &
+                 nfluxoutv_nu,nfluxgndv_nu,                  &
                  fmnetv,fluxupv,fluxdnv,fzerov,taugsurf)
 
@@ -1192 +1202 @@
 !-----------------------------------------------------------------------
 
-         call optci(plevrad,tlevrad,dtaui,taucumi,                  &
-              qxiaer,qsiaer,giaer,cosbi,wbari,tauaero,tmid,pmid,    &
+         call rad_correlatedk_opacities_thermal(plevrad,          &
+              tlevrad,dtaui,taucumi,                              &
+              qxiaer,qsiaer,giaer,cosbi,wbari,tauaero,tmid,pmid,  &
               taugsurfi,qvar,muvarrad,fracvari)
 
-         call sfluxi(plevrad,tlevrad,dtaui,taucumi,ubari,albi,      &
-              wnoi,dwni,cosbi,wbari,nfluxtopi,nfluxtopi_nu,         & 
+         call rad_correlatedk_fluxes_thermal(plevrad,             &
+              tlevrad,dtaui,taucumi,ubari,albi,                   &
+              wnoi,dwni,cosbi,wbari,nfluxtopi,nfluxtopi_nu,       & 
               fmneti,fluxupi,fluxdni,fluxupi_nu,fzeroi,taugsurfi)
 
@@ -1346 +1358 @@
 
 
-    end subroutine callcorrk
-
-END MODULE callcorrk_mod
+    end subroutine rad_correlatedk
+
+END MODULE rad_correlatedk_mod