Changeset 4077 for trunk/LMDZ.GENERIC/libf/phygeneric/aerosol_opacity.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/aerosol_opacity.F90 (moved) (moved from trunk/LMDZ.GENERIC/libf/phygeneric/aeropacity.F90) (9 diffs)

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

@@ -1 @@
-module aeropacity_mod
+module aerosol_opacity_mod
 
 implicit none
@@ -5 +5 @@
 contains
 
-      Subroutine aeropacity(ngrid,nlayer,nq,pplay,pplev,pt,pq,zls, &
+      Subroutine aerosol_opacity(ngrid,nlayer,nq,pplay,pplev,pt,pq,zls, &
          aerosol,reffrad,nueffrad, QREFvis3d,QREFir3d,tau_col, &
          cloudfrac,totcloudfrac,clearsky)
 
        use radinc_h, only : L_TAUMAX,naerkind
-       use aerosol_mod, only: iaero_nlay, iaero_generic, &
+       use aerosol_global_variables , only: iaero_nlay, iaero_generic, &
                               iaero_aurora, iaero_back2lay, iaero_co2, &
                               iaero_dust, iaero_h2o, iaero_h2so4, &
@@ -120 +120 @@
       IF (firstcall) THEN
         ia =0
-        write(*,*) "Tracers found in aeropacity:"
+        write(*,*) "Tracers found in aerosol_opacity:"
         do iq=1,nq
           tracername=noms(iq)
@@ -135 +135 @@
 
         if (noaero) then
-          print*, "No active aerosols found in aeropacity"
+          print*, "No active aerosols found in aerosol_opacity"
         else
           print*, "If you would like to use aerosols, make sure any old"
           print*, "start files are updated in newstart using the option"
           print*, "q=0"
-          write(*,*) "Active aerosols found in aeropacity:"
+          write(*,*) "Active aerosols found in aerosol_opacity:"
         endif
 
@@ -269 +269 @@
                   end do
 
-                  call abort_physic("aeropacity", "Something wrong happened on water ice liquid opacity calculation",1)
+                  call abort_physic("aerosol_opacity", "Something wrong happened on water ice liquid opacity calculation",1)
                endif
 
@@ -276 +276 @@
                do ig=1, ngrid
                   !do l=1,nlayer-1 ! to stop the rad tran bug
-                  do l=1,nlayer !JL18 if aerosols are present in the last layer we must account for them. Provides better upper boundary condition in the IR. They must however be put to zero in the sw (see optcv)
+                  do l=1,nlayer !JL18 if aerosols are present in the last layer we must account for them. Provides better upper boundary condition in the IR. They must however be put to zero in the sw (see rad_correlatedk_opacities_stellar)
                                 ! same correction should b-probably be done for other aerosol types.
                      aerosol(ig,l,iaer) =                                    & !modification by BC
@@ -294 +294 @@
                   do ig=1, ngrid
                      !do l=1,nlayer-1 ! to stop the rad tran bug
-                     do l=1,nlayer !JL18 if aerosols are present in the last layer we must account for them. Provides better upper boundary condition in the IR. They must however be put to zero in the sw (see optcv)
+                     do l=1,nlayer !JL18 if aerosols are present in the last layer we must account for them. Provides better upper boundary condition in the IR. They must however be put to zero in the sw (see rad_correlatedk_opacities_stellar)
                         CLFtot  = max(totcloudfrac(ig),0.01)
                         aerosol(ig,l,iaer)=aerosol(ig,l,iaer)/CLFtot
@@ -303 +303 @@
                   do ig=1, ngrid
                      !do l=1,nlayer-1 ! to stop the rad tran bug
-                     do l=1,nlayer !JL18 if aerosols are present in the last layer we must account for them. Provides better upper boundary condition in the IR. They must however be put to zero in the sw (see optcv)
+                     do l=1,nlayer !JL18 if aerosols are present in the last layer we must account for them. Provides better upper boundary condition in the IR. They must however be put to zero in the sw (see rad_correlatedk_opacities_stellar)
                         CLFtot  = CLFfixval
                         aerosol(ig,l,iaer)=aerosol(ig,l,iaer)/CLFtot
@@ -1092 +1092 @@
       ! end do
 
-    end subroutine aeropacity
+    end subroutine aerosol_opacity
       
-end module aeropacity_mod
+end module aerosol_opacity_mod