Changeset 3802

Timestamp:

Jun 10, 2025, 11:20:19 AM (4 days ago)

Author:

debatzbr

Message:

Pluto PCM: add some sanity checks to ensure haze stability
BBT

Location:

trunk/LMDZ.PLUTO/libf/phypluto

Files:

• aeropacity.F90 (modified) (1 diff)
• aeroptproperties.F90 (modified) (7 diffs)

trunk/LMDZ.PLUTO/libf/phypluto/aeropacity.F90

@@ -122 +122 @@
             ! Spherical aerosols
             sig = 0.2
-            dtau_aer(:,:,1) = m0as(:,:) * QREFvis3d(:,:,1) * pi * mp2m_rc_sph(:,:)**2 * exp(2*sig**2)
+            where ((m0as(:,:) >= 1e-8).and.(mp2m_rc_sph(:,:) >= 5e-9))
+               dtau_aer(:,:,1) = m0as(:,:) * QREFvis3d(:,:,1) * pi * mp2m_rc_sph(:,:)**2 * exp(2*sig**2)
+            elsewhere
+               dtau_aer(:,:,1) = 0d0
+            endwhere
             ! Fractal aerosols
             sig = 0.35
-            dtau_aer(:,:,2) = m0af(:,:) * QREFvis3d(:,:,2) * pi * mp2m_rc_fra(:,:)**2 * exp(2*sig**2)
+            where ((m0af(:,:) >= 1e-8).and.(mp2m_rc_fra(:,:) >= 1e-8))
+               dtau_aer(:,:,2) = m0af(:,:) * QREFvis3d(:,:,2) * pi * mp2m_rc_fra(:,:)**2 * exp(2*sig**2)
+            elsewhere
+               dtau_aer(:,:,2) = 0d0
+            endwhere
             
             ! write(*,*) 'dtau_as  :', MINVAL(dtau_aer(:,:,1)), '-', MAXVAL(dtau_aer(:,:,1))

trunk/LMDZ.PLUTO/libf/phypluto/aeroptproperties.F90

@@ -38 +38 @@
 !     ==============================================================
 
-!     Local variables
-!     ---------------
-
-
+!     INPUTS
+!     ------
+      INTEGER :: ngrid,nlayer
+!     Aerosol effective radius used for radiative transfer (meter)
+      REAL :: reffrad(ngrid,nlayer,naerkind)
+!     Aerosol effective variance used for radiative transfer (n.u.)
+      REAL,INTENT(IN) :: nueffrad(ngrid,nlayer,naerkind)
+
+!     OUTPUTS
+!     -------
+      REAL,INTENT(OUT) :: QVISsQREF3d(ngrid,nlayer,L_NSPECTV,naerkind)
+      REAL,INTENT(OUT) :: omegaVIS3d(ngrid,nlayer,L_NSPECTV,naerkind)
+      REAL,INTENT(OUT) :: gVIS3d(ngrid,nlayer,L_NSPECTV,naerkind)
+
+      REAL,INTENT(OUT) :: QIRsQREF3d(ngrid,nlayer,L_NSPECTI,naerkind)
+      REAL,INTENT(OUT) :: omegaIR3d(ngrid,nlayer,L_NSPECTI,naerkind)
+      REAL,INTENT(OUT) :: gIR3d(ngrid,nlayer,L_NSPECTI,naerkind)
+
+      REAL,INTENT(OUT) :: QREFvis3d(ngrid,nlayer,naerkind)
+      REAL,INTENT(OUT) :: QREFir3d(ngrid,nlayer,naerkind)
+
+!      REAL :: omegaREFvis3d(ngrid,nlayer,naerkind)
+!      REAL :: omegaREFir3d(ngrid,nlayer,naerkind)
 
 !     =============================================================
@@ -48 +67 @@
 
 !     Min. and max radius of the interpolation grid (in METERS)
-      REAL, PARAMETER :: refftabmin = 2e-8 !2e-8
+      REAL, PARAMETER :: refftabmin = 4e-9 !2e-8
 !      REAL, PARAMETER :: refftabmax = 35e-6
       REAL, PARAMETER :: refftabmax = 1e-3
@@ -95 +114 @@
       INTEGER :: ig,lg,chg
 
+
 !     Local saved variables
 !     ---------------------
-
 !     Radius axis of the interpolation grid
       REAL,SAVE :: refftab(refftabsize)
@@ -170 +189 @@
 !$OMP THREADPRIVATE(omegIRa,omegrefIRa,omegIRb,omegrefIRb,gIRa,gIRb)
 
+      ! Local variables:
+      !-----------------
       REAL :: radiusm
       REAL :: radiusr
-
-!     Inputs
-!     ------
-
-      INTEGER :: ngrid,nlayer
-!     Aerosol effective radius used for radiative transfer (meter)
-      REAL :: reffrad(ngrid,nlayer,naerkind)
-!     Aerosol effective variance used for radiative transfer (n.u.)
-      REAL,INTENT(IN) :: nueffrad(ngrid,nlayer,naerkind)
-
-!     Outputs
-!     -------
-
-      REAL,INTENT(OUT) :: QVISsQREF3d(ngrid,nlayer,L_NSPECTV,naerkind)
-      REAL,INTENT(OUT) :: omegaVIS3d(ngrid,nlayer,L_NSPECTV,naerkind)
-      REAL,INTENT(OUT) :: gVIS3d(ngrid,nlayer,L_NSPECTV,naerkind)
-
-      REAL,INTENT(OUT) :: QIRsQREF3d(ngrid,nlayer,L_NSPECTI,naerkind)
-      REAL,INTENT(OUT) :: omegaIR3d(ngrid,nlayer,L_NSPECTI,naerkind)
-      REAL,INTENT(OUT) :: gIR3d(ngrid,nlayer,L_NSPECTI,naerkind)
-
-      REAL,INTENT(OUT) :: QREFvis3d(ngrid,nlayer,naerkind)
-      REAL,INTENT(OUT) :: QREFir3d(ngrid,nlayer,naerkind)
-
-!      REAL :: omegaREFvis3d(ngrid,nlayer,naerkind)
-!      REAL :: omegaREFir3d(ngrid,nlayer,naerkind)
 
       REAL :: minrad    ! minimal radius in table .dat radiustab (outside 0)
@@ -260 +255 @@
       ENDIF ! of IF (first_allocate)
 
+
       DO iaer = 1, naerkind ! Loop on aerosol kind
+        !==================================================================
+        ! If there is one single particle size, optical
+        ! properties of the considered aerosol are homogeneous
         IF ( (nsize(iaer,1).EQ.1).AND.(nsize(iaer,2).EQ.1) ) THEN
-!==================================================================
-!       If there is one single particle size, optical
-!         properties of the considered aerosol are homogeneous
+
           DO lg = 1, nlayer
             DO ig = 1, ngrid
@@ -289 +286 @@
           endif
 
-        ELSE ! Varying effective radius and variance
+        !==================================================================
+        ! Various particule size
+        ! Varying effective radius and variance
+        ELSE
+
       DO idomain = 1, 2 ! Loop on visible or infrared channel
-!==================================================================
 !     1. Creating the effective radius and variance grid
 !     --------------------------------------------------
@@ -487 +487 @@
 !         2.1 Effective radius index and kx calculation
           var_tmp=reffrad(ig,lg,iaer)/refftabmin
-          var_tmp=log(var_tmp)*3.
+          if (var_tmp.GT.0.) then ! for debug
+            var_tmp=log(var_tmp)*3.
+          else
+            var_tmp=0.
+          endif
           var_tmp=var_tmp/logvratgrid+1.
           grid_i=floor(var_tmp)