Changeset 3196 for trunk/LMDZ.PLUTO/libf

Timestamp:

Jan 31, 2024, 5:30:43 PM (11 months ago)

Author:

afalco

Message:

Pluto PCM:

Included N2 condensation from PLUTO.old & read haze aerosols.

AF

Location:

trunk/LMDZ.PLUTO/libf/phypluto

Files:

• aeroptproperties.F90 (modified) (12 diffs)
• physiq_mod.F90 (modified) (2 diffs)

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

@@ -38 +38 @@
 !     ==============================================================
 
-!     Local variables 
+!     Local variables
 !     ---------------
 
@@ -198 +198 @@
 !      REAL :: omegaREFvis3d(ngrid,nlayer,naerkind)
 !      REAL :: omegaREFir3d(ngrid,nlayer,naerkind)
+
+      REAL :: minrad    ! minimal radius in table .dat radiustab (outside 0)
+      REAL :: maxrad
 
 !     0. Allocate local saved arrays at firstcall
@@ -246 +249 @@
         ALLOCATE(qsqrefIRb(L_NSPECTI,ngau,naerkind))
         ALLOCATE(qrefIRb(ngau,naerkind))
-        
+
         ALLOCATE(omegIRa(L_NSPECTI,ngau,naerkind))
         ALLOCATE(omegrefIRa(ngau,naerkind))
@@ -253 +256 @@
         ALLOCATE(gIRa(L_NSPECTI,ngau,naerkind))
         ALLOCATE(gIRb(L_NSPECTI,ngau,naerkind))
-        
+
         first_allocate=.false.
       ENDIF ! of IF (first_allocate)
@@ -317 +320 @@
            enddo
         endif
+
+        ! check that radiustab (table.dat) cover the reffrad used for haze : iaer=1
+        minrad=min(MINVAL(radiustab(1,1,1:nsize(1,1))),MINVAL(radiustab(1,2,1:nsize(1,2))))
+        maxrad=min(MAXVAL(radiustab(1,1,1:nsize(1,1))),MAXVAL(radiustab(1,2,1:nsize(1,2))))
+        IF ((MINVAL(reffrad).LE.minrad).OR.(MAXVAL(reffrad).GE.maxrad)) then
+           WRITE(*,*) 'Warning: particle size in grid box #'
+           WRITE(*,*) ig,' is too large to be used by the '
+           WRITE(*,*) 'radiative transfer; please extend the '
+           WRITE(*,*) 'interpolation grid to larger grain sizes.'
+           WRITE(*,*) 'radiustab=',minrad,'-',maxrad
+           WRITE(*,*) 'reffrad=',MINVAL(reffrad),'-',MAXVAL(reffrad)
+           stop
+        ENDIF
 
         firstcall = .false.
@@ -412 +428 @@
             qsqrefVISb(m,gausind,iaer)=                         &
                     (1-kint)*QVISsQREF(m,iaer,radius_id) +      &
-                    kint*QVISsQREF(m,iaer,radius_id+1)    
+                    kint*QVISsQREF(m,iaer,radius_id+1)
             omegVISb(m,gausind,iaer)=                           &
                     (1-kint)*omegaVIS(m,iaer,radius_id) +       &
@@ -500 +516 @@
                     EXP(-dista(j,1,iaer,idomain,gausind) *            &
                     dista(j,1,iaer,idomain,gausind) *                 &
-                    0.5e0/sizedistk2)/(radiusm-drad)                  
+                    0.5e0/sizedistk2)/(radiusm-drad)
                   dista(j,1,iaer,idomain,gausind) =                   &
                     dista(j,1,iaer,idomain,gausind) /                 &
@@ -540 +556 @@
 !                     Convolution:
                       qextVISgrid(j,1,m,iaer) =              &
-                        qextVISgrid(j,1,m,iaer) +            & 
+                        qextVISgrid(j,1,m,iaer) +            &
                         weightgaus(gausind) *                &
                         (                                    &
@@ -609 +625 @@
                       (                                      &
                       omegrefVISb(gausind,iaer) *            &
-                      qrefVISb(gausind,iaer) *               & 
+                      qrefVISb(gausind,iaer) *               &
                       pi*radGAUSb(gausind,iaer,idomain) *    &
                       radGAUSb(gausind,iaer,idomain) *       &
@@ -622 +638 @@
 
                   qrefVISgrid(j,1,iaer)=qrefVISgrid(j,1,iaer) /          &
-                                normd(j,1,iaer,idomain)       
+                                normd(j,1,iaer,idomain)
                   qscatrefVISgrid(j,1,iaer)=qscatrefVISgrid(j,1,iaer) /  &
                                 normd(j,1,iaer,idomain)
@@ -736 +752 @@
                       )
                   ENDDO
- 
+
                   qrefIRgrid(j,1,iaer)=qrefIRgrid(j,1,iaer) /          &
                                 normd(j,1,iaer,idomain)
@@ -789 +805 @@
             omegaIR3d(ig,lg,m,iaer) =                              &
                         k1*omegIRgrid(grid_i,1,m,iaer) +           &
-                        k2*omegIRgrid(grid_i+1,1,m,iaer) 
-            gIR3d(ig,lg,m,iaer) =                                  & 
+                        k2*omegIRgrid(grid_i+1,1,m,iaer)
+            gIR3d(ig,lg,m,iaer) =                                  &
                         k1*gIRgrid(grid_i,1,m,iaer) +              &
                         k2*gIRgrid(grid_i+1,1,m,iaer)

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

@@ -1474 +1474 @@
          enddo
 
-         ! Generalised for arbitrary aerosols now. By LK
+         if (aerohaze) then
+            ! Generalised for arbitrary aerosols now. By LK
          reffcol(1:ngrid,1:naerkind)=0.0
             ! call n2_reffrad(ngrid,nlayer,nq,zq,reffrad(1,1,iaero_haze))
@@ -1496 +1497 @@
             ENDDO
          ENDDO
+         endif ! end of aerohaze
       endif ! end of 'tracer'