Changeset 3585 for trunk/LMDZ.PLUTO/libf/phypluto/dyn1d/kcm1d.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/dyn1d/kcm1d.F90 (modified) (6 diffs)

trunk/LMDZ.PLUTO/libf/phypluto/dyn1d/kcm1d.F90

@@ -78 +78 @@
 
   real dTstrat
-  real,allocatable :: aerosol(:,:) ! aerosol tau (kg/kg)
+  real,allocatable :: dtau_aer(:,:) ! aerosol tau (kg/kg)
   real OLR_nu(1,L_NSPECTI)
   real OSR_nu(1,L_NSPECTV)
@@ -320 +320 @@
   !write(*,*) 1,llm,nq,0,86400.0,1,1.0,latitude,longitude,cell_area,rad,g,r,cpp
 
-  ! initialise naerkind (from callphys.def) and allocate aerosol(:,:)
+  ! initialise naerkind (from callphys.def) and allocate dtau_aer(:,:)
   naerkind=0 !default
   call getin("naerkind",naerkind)
-  allocate(aerosol(llm,naerkind))
-  aerosol(:,:)=0
+  allocate(dtau_aer(llm,naerkind))
+  dtau_aer(:,:)=0
 
   do iq=1,nq
@@ -374 +374 @@
           albedo_wv,albedo_equivalent,                    &
           emis,mu0,plev,play,temp,                        &
-          tsurf,fract,dist_star,aerosol,muvar,            &
+          tsurf,fract,dist_star,dtau_aer,muvar,            &
           dtlw,dtsw,fluxsurf_lw,fluxsurf_sw,              & 
           fluxsurfabs_sw,fluxtop_lw,                      &
@@ -385 +385 @@
      !     albedo_wv,'D',albedo_equivalent,'E',                    &
      !     emis,'F',mu0,'G',plev,'H',play,'I',temp,'J',                        &
-     !     tsurf,'K',fract,'L',dist_star,'M',aerosol,'N',muvar,'O',            &
+     !     tsurf,'K',fract,'L',dist_star,'M',dtau_aer,'N',muvar,'O',            &
      !     dtlw,'P',dtsw,'Q',fluxsurf_lw,'R',fluxsurf_sw,'S',              & 
      !     fluxsurfabs_sw,'T',fluxtop_lw,'U',                      &
@@ -426 +426 @@
   call callcorrk(1,nlayer,q,nq,qsurf,                          &
        albedo_wv,albedo_equivalent,emis,mu0,plev,play,temp,    &
-       tsurf,fract,dist_star,aerosol,muvar,                    &
+       tsurf,fract,dist_star,dtau_aer,muvar,                    &
        dtlw,dtsw,fluxsurf_lw,fluxsurf_sw,fluxsurfabs_sw,       &
        fluxtop_lw, fluxabs_sw,fluxtop_dn,OLR_nu,OSR_nu,GSR_nu, &
@@ -437 +437 @@
      !     albedo_wv,'D',albedo_equivalent,'E',                    &
      !     emis,'F',mu0,'G',plev,'H',play,'I',temp,'J',                        &
-     !     tsurf,'K',fract,'L',dist_star,'M',aerosol,'N',muvar,'O',            &
+     !     tsurf,'K',fract,'L',dist_star,'M',dtau_aer,'N',muvar,'O',            &
      !     dtlw,'P',dtsw,'Q',fluxsurf_lw,'R',fluxsurf_sw,'S',              & 
      !     fluxsurfabs_sw,'T',fluxtop_lw,'U',                      &