Changeset 3585 for trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90

Timestamp:

Jan 20, 2025, 10:17:54 AM (12 days ago)

Author:

debatzbr

Message:

Connecting microphysics to radiative transfer + miscellaneous cleans

File:

• trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90 (modified) (10 diffs)

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

@@ -21 +21 @@
       use radcommon_h, only: sigma, glat, grav, BWNV, WNOI, DWNI, DWNV, WNOV
       use suaer_corrk_mod, only: suaer_corrk
-      use radii_mod, only: su_aer_radii,haze_reffrad_fix
       use aerosol_mod, only: i_haze, haze_prof
       use surfdat_h, only: phisfi, zmea, zstd, zsig, zgam, zthe, &
@@ -35 +34 @@
                           igcm_co_gas,igcm_co_ice,igcm_prec_haze,lw_n2,lw_ch4,lw_co,&
                           alpha_lift, alpha_devil, qextrhor, &
-                          nesp, is_chim, is_condensable
+                          nesp, is_chim
       use time_phylmdz_mod, only: ecritphy, iphysiq, nday
       use phyetat0_mod, only: phyetat0,tab_cntrl_mod
@@ -286 +285 @@
 !     for the "naerkind" optically active aerosols:
 
-      real,save,allocatable :: aerosol(:,:,:) ! Aerosols
-!$OMP THREADPRIVATE(aerosol)
+      real,save,allocatable :: dtau_aer(:,:,:) ! Aerosols
+!$OMP THREADPRIVATE(dtau_aer)
       real zh(ngrid,nlayer)               ! Potential temperature (K).
       real pw(ngrid,nlayer)               ! Vertical velocity (m/s). (NOTE : >0 WHEN DOWNWARDS !!)
@@ -464 +463 @@
       real muvar(ngrid,nlayer+1) ! For Runaway Greenhouse 1D study. By RW
 
-      real,save,allocatable :: reffcol(:,:)
-!$OMP THREADPRIVATE(reffcol)
-
 !  Non-oro GW tendencies
       REAL d_u_hin(ngrid,nlayer), d_v_hin(ngrid,nlayer)
@@ -535 +531 @@
          ! allocate related local arrays
          ! (need be allocated instead of automatic because of "naerkind")
-         allocate(aerosol(ngrid,nlayer,naerkind))
-         allocate(reffcol(ngrid,naerkind))
+         allocate(dtau_aer(ngrid,nlayer,naerkind))
 
 #ifdef CPP_XIOS
@@ -1022 +1017 @@
                if (oldplutocorrk) then
                   call callcorrk_pluto(icount,ngrid,nlayer,pq,nq,qsurf,          &
-                               albedo(:,1),emis,mu0,pplev,pplay,pt,                   &
-                               zzlay,tsurf,fract,dist_star,aerosol,              &
+                               albedo(:,1),emis,mu0,pplev,pplay,pt,              &
+                               zzlay,zzlev,tsurf,fract,dist_star,dtau_aer,       &
                                zdtlw,zdtsw,fluxsurf_lw,fluxsurf_sw,fluxtop_lw,   &
                                fluxabs_sw,fluxtop_dn,reffrad,tau_col,ptime,pday, &
@@ -1034 +1029 @@
                 call callcorrk(ngrid,nlayer,pq,nq,qsurf,  &
                               albedo,albedo_equivalent,emis,mu0,pplev,pplay,pt,   &
-                              zzlay,tsurf,fract,dist_star,aerosol,muvar,                &
+                              zzlay,zzlev,tsurf,fract,dist_star,dtau_aer,muvar,   &
                               zdtlw,zdtsw,fluxsurf_lw,fluxsurf_sw,                &
                               fluxsurfabs_sw,fluxtop_lw,                          &
@@ -1436 +1431 @@
          zdqch4cloud(:,:,:)=0.
          if ((methane).and.(metcloud).and.(.not.fast)) THEN
-            call ch4cloud(ngrid,nlayer,naerkind,ptimestep,  &
+            call ch4cloud(ngrid,nlayer,ptimestep,  &
                       pplev,pplay,pdpsrf,zzlev,zzlay, pt,pdt,  &
                       pq,pdq,zdqch4cloud,zdqsch4cloud,zdtch4cloud,   &
@@ -1469 +1464 @@
          zdqcocloud(:,:,:)=0.
          IF ((carbox).and.(monoxcloud).and.(.not.fast)) THEN
-            call cocloud(ngrid,nlayer,naerkind,ptimestep,   &
+            call cocloud(ngrid,nlayer,ptimestep,   &
                       pplev,pplay,pdpsrf,zzlev,zzlay, pt,pdt,  &
                       pq,pdq,zdqcocloud,zdqscocloud,zdtcocloud,   &
@@ -1845 +1840 @@
                print*,'fluxtop_dn=',fluxtop_dn(ig)
                print*,'tau_col=',tau_col(ig)
-               print*,'aerosol=',aerosol(ig,:,:)
+               print*,'dtau_aer=',dtau_aer(ig,:,:)
                print*,'temp=   ',pt(ig,:)
                print*,'pplay=  ',pplay(ig,:)