Changeset 3585 for trunk/LMDZ.PLUTO


Ignore:
Timestamp:
Jan 20, 2025, 10:17:54 AM (9 hours ago)
Author:
debatzbr
Message:

Connecting microphysics to radiative transfer + miscellaneous cleans

Location:
trunk/LMDZ.PLUTO/libf/phypluto
Files:
21 edited

Legend:

Unmodified
Added
Removed

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

    r3572 r3585  
    55contains
    66
    7       Subroutine aeropacity(ngrid,nlayer,nq,pplay,pplev,pt, pq, &
    8          aerosol,reffrad,nueffrad, QREFvis3d,QREFir3d,tau_col)
     7      Subroutine aeropacity(ngrid,nlayer,nq,pplay,pplev,zzlev,pt,pq, &
     8         dtau_aer,reffrad,nueffrad, QREFvis3d,QREFir3d,tau_col)
    99
    1010       use radinc_h, only : L_TAUMAX,naerkind
    1111       use aerosol_mod, only: iaero_haze, i_haze
    12        USE tracer_h, only: noms,rho_n2,rho_ice,rho_q,mmol
     12       USE tracer_h, only: noms,rho_n2,rho_ice,rho_q,mmol,micro_indx
    1313       use comcstfi_mod, only: g, pi, mugaz, avocado
    1414       use geometry_mod, only: latitude
    15        use callkeys_mod, only: kastprof
     15       use callkeys_mod, only: kastprof, callmufi
     16       use mp2m_diagnostics
    1617       implicit none
    1718
    … …  
    2829!     update J.-B. Madeleine (2008)
    2930!     dust removal, simplification by Robin Wordsworth (2009)
    30 !     Generic n-layer aerosol - J. Vatant d'Ollone (2020)
    31 !     Radiative Generic Condensable Species - Lucas Teinturier (2022)
    3231!
    3332!     Input
    … …  
    4443!     Output
    4544!     ------
    46 !     aerosol            Aerosol optical depth in layer l, grid point ig
     45!     dtau_aer           Aerosol optical depth in layer l, grid point ig
    4746!     tau_col            Total column optical depth at grid point ig
    4847!
    … …  
    5453      REAL,INTENT(IN) :: pplay(ngrid,nlayer) ! mid-layer pressure (Pa)
    5554      REAL,INTENT(IN) :: pplev(ngrid,nlayer+1) ! inter-layer pressure (Pa)
     55      REAL,INTENT(IN) :: zzlev(ngrid,nlayer)   ! Altitude at the layer boundaries.
    5656      REAL,INTENT(IN) :: pq(ngrid,nlayer,nq) ! tracers (.../kg_of_air)
    5757      REAL,INTENT(IN) :: pt(ngrid,nlayer) ! mid-layer temperature (K)
    58       REAL,INTENT(OUT) :: aerosol(ngrid,nlayer,naerkind) ! aerosol optical depth
    59       REAL,INTENT(IN) :: reffrad(ngrid,nlayer,naerkind) ! aerosol effective radius
    60       REAL,INTENT(IN) :: nueffrad(ngrid,nlayer,naerkind) ! aerosol effective variance
     58      REAL,INTENT(OUT) :: dtau_aer(ngrid,nlayer,naerkind) ! aerosol optical depth
     59      REAL,INTENT(IN) :: reffrad(ngrid,nlayer,naerkind)   ! aerosol effective radius
     60      REAL,INTENT(IN) :: nueffrad(ngrid,nlayer,naerkind)  ! aerosol effective variance
    6161      REAL,INTENT(IN) :: QREFvis3d(ngrid,nlayer,naerkind) ! extinction coefficient in the visible
    6262      REAL,INTENT(IN) :: QREFir3d(ngrid,nlayer,naerkind)
    … …  
    6969      real dp_tropo(ngrid)
    7070      real dp_layer(ngrid)
     71
     72      ! Microphysical tracers
     73      real sig
     74      real m0as(ngrid,nlayer)
     75      real m0af(ngrid,nlayer)
    7176
    7277      INTEGER l,ig,iq,iaer,ia
    … …  
    110115
    111116      if (iaero_haze.ne.0) then
    112         iaer=iaero_haze
    113 !       1. Initialization
    114          aerosol(1:ngrid,1:nlayer,iaer)=0.0
    115 !       2. Opacity calculation
    116          DO ig=1, ngrid
    117                DO l=1,nlayer-1 ! to stop the rad tran bug
    118                   ! if fractal, radius doit etre equivalent sphere radius
    119                   aerosol0 =                         &
    120                        (  0.75 * QREFvis3d(ig,l,iaer) /        &
    121                        ( rho_q(i_haze) * reffrad(ig,l,iaer) )  ) *   &
    122                        ( pq(ig,l,i_haze) + 1.E-10 ) *         &
    123                        ( pplev(ig,l) - pplev(ig,l+1) ) / g
    124                   aerosol0           = max(aerosol0,1.e-10)
    125                   aerosol0           = min(aerosol0,L_TAUMAX)
    126                   aerosol(ig,l,iaer) = aerosol0
     117         if (callmufi) then
     118            ! Convert intensive microphysical tracers to extensive [m-2]
     119            m0as(:,:) = pq(:,:,micro_indx(1)) * (pplev(:,1:nlayer) - pplev(:,2:nlayer+1)) / g
     120            m0af(:,:) = pq(:,:,micro_indx(3)) * (pplev(:,1:nlayer) - pplev(:,2:nlayer+1)) / g
     121           
     122            ! Spherical aerosols
     123            sig = 0.2
     124            dtau_aer(:,:,1) = m0as(:,:) * QREFvis3d(:,:,1) * pi * mp2m_rc_sph(:,:)**2 * exp(2*sig**2)
     125            ! Fractal aerosols
     126            sig = 0.35
     127            dtau_aer(:,:,2) = m0af(:,:) * QREFvis3d(:,:,2) * pi * mp2m_rc_fra(:,:)**2 * exp(2*sig**2)
     128           
     129            ! write(*,*) 'dtau_as  :', MINVAL(dtau_aer(:,:,1)), '-', MAXVAL(dtau_aer(:,:,1))
     130            ! write(*,*) 'dtau_af  :', MINVAL(dtau_aer(:,:,2)), '-', MAXVAL(dtau_aer(:,:,2))
     131         
     132         else
     133            do iaer = 1, naerkind
     134               ! 1. Initialization
     135               dtau_aer(1:ngrid,1:nlayer,iaer)=0.0
     136               ! 2. Opacity calculation
     137               DO ig = 1, ngrid
     138                  DO l = 1, nlayer-1 ! to stop the rad tran bug
     139                     ! If fractal, radius doit etre equivalent sphere radius
     140                     ! Eq. 2.37 - Madeleine's PhD (2011).
     141                     aerosol0 =                         &
     142                        (  0.75 * QREFvis3d(ig,l,iaer) /        &
     143                        ( rho_q(i_haze) * reffrad(ig,l,iaer) )  ) *   &
     144                        ( pq(ig,l,i_haze) + 1.E-10 ) *         &
     145                        ( pplev(ig,l) - pplev(ig,l+1) ) / g
     146                     aerosol0            = max(aerosol0,1.e-10)
     147                     aerosol0            = min(aerosol0,L_TAUMAX)
     148                     dtau_aer(ig,l,iaer) = aerosol0
     149                  ENDDO
    127150               ENDDO
    128          ENDDO
    129          !QREF est le meme dans toute la colonne par def si size uniforme
    130          !print*, 'TB17: QREFvis3d=',QREFvis3d(1,:,1)
    131          !print*, 'TB17: rho_q=',rho_q(i_haze)
    132          !print*, 'TB17: reffrad=',reffrad(1,:,1)
    133          !print*, 'TB17: pq=',pq(1,:,i_haze)
    134          !print*, 'TB17: deltap=',pplev(1,1) - pplev(1,nlayer)
    135    end if ! if haze aerosols
     151               !QREF est le meme dans toute la colonne par def si size uniforme
     152               !print*, 'TB17: QREFvis3d=',QREFvis3d(1,:,1)
     153               !print*, 'TB17: rho_q=',rho_q(i_haze)
     154               !print*, 'TB17: reffrad=',reffrad(1,:,1)
     155               !print*, 'TB17: pq=',pq(1,:,i_haze)
     156               !print*, 'TB17: deltap=',pplev(1,1) - pplev(1,nlayer)
     157            enddo ! end iaer = 1, naerkind
     158         endif ! end callmufi
     159      endif ! if haze aerosols
    136160
    137161! --------------------------------------------------------------------------
    … …  
    142166      do l=1,nlayer
    143167         do ig=1,ngrid
    144             tau_col(ig) = tau_col(ig) + aerosol(ig,l,iaer)
     168            tau_col(ig) = tau_col(ig) + dtau_aer(ig,l,iaer)
    145169         end do
    146170      end do
    … …  
    150174      do l=1,nlayer
    151175         do iaer = 1, naerkind
    152             if(aerosol(ig,l,iaer).gt.1.e3)then
    153                print*,'WARNING: aerosol=',aerosol(ig,l,iaer)
     176            if(dtau_aer(ig,l,iaer).gt.1.e3)then
     177               print*,'WARNING: dtau_aer=',dtau_aer(ig,l,iaer)
    154178               print*,'at ig=',ig,',  l=',l,', iaer=',iaer
    155179               print*,'QREFvis3d=',QREFvis3d(ig,l,iaer)
    … …  
    164188         print*,'WARNING: tau_col=',tau_col(ig)
    165189         print*,'at ig=',ig
    166          print*,'aerosol=',aerosol(ig,:,:)
     190         print*,'dtau_aer=',dtau_aer(ig,:,:)
    167191         print*,'QREFvis3d=',QREFvis3d(ig,:,:)
    168192         print*,'reffrad=',reffrad(ig,:,:)

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

    r3353 r3585  
    279279              QREFvis3d(ig,lg,iaer)=QREFvis(iaer,1)
    280280              QREFir3d(ig,lg,iaer)=QREFir(iaer,1)
    281 !              omegaREFvis3d(ig,lg,iaer)=omegaREFvis(iaer,1)
    282 !              omegaREFir3d(ig,lg,iaer)=omegaREFir(iaer,1)
     281              !omegaREFvis3d(ig,lg,iaer)=omegaREFvis(iaer,1)
     282              !omegaREFir3d(ig,lg,iaer)=omegaREFir(iaer,1)
    283283            ENDDO
    284284          ENDDO
    285 
    286285
    287286          if (firstcall) then
    … …  
    297296      IF (firstcall) THEN
    298297
    299 !       1.1 Pi!
     298        ! 1.1 Pi
     299        !-------
    300300        pi = 2. * asin(1.e0)
    301301
    302 !       1.2 Effective radius
    303         refftab(1)    = refftabmin
     302        ! 1.2 Effective radius
     303        !---------------------
     304        refftab(1)           = refftabmin
    304305        refftab(refftabsize) = refftabmax
    305306
    … …  
    311312        enddo
    312313
    313 !       1.3 Effective variance
     314        ! 1.3 Effective variance
     315        !-----------------------
    314316        if(nuefftabsize.eq.1)then ! addded by RDW
    315317           print*,'Warning: no variance range in aeroptproperties'
    … …  
    344346      ENDIF ! of IF (firstcall)
    345347
    346 !       1.4 Radius middle point and range for Gauss integration
     348        ! 1.4 Radius middle point and range for Gauss integration
     349        !--------------------------------------------------------
    347350        radiusm=0.5*(radiustab(iaer,idomain,nsize(iaer,idomain)) + radiustab(iaer,idomain,1))
    348351        radiusr=0.5*(radiustab(iaer,idomain,nsize(iaer,idomain)) - radiustab(iaer,idomain,1))
    349352
    350 !       1.5 Interpolating data at the Gauss quadrature points:
     353        ! 1.5 Interpolating data at the Gauss quadrature points:
     354        !-------------------------------------------------------
    351355        DO gausind=1,ngau
    352356          drad=radiusr*radgaus(gausind)
    … …  
    829833!==================================================================
    830834
    831 
    832 
    833835      ENDDO ! idomain
    834836

  • trunk/LMDZ.PLUTO/libf/phypluto/callcorrk.F90

    r3572 r3585  
    77      subroutine callcorrk(ngrid,nlayer,pq,nq,qsurf,           &
    88          albedo,albedo_equivalent,emis,mu0,pplev,pplay,pt,    &
    9           zzlay,tsurf,fract,dist_star,aerosol,muvar,           &
     9          zzlay,zzlev,tsurf,fract,dist_star,dtau_aer,muvar,    &
    1010          dtlw,dtsw,fluxsurf_lw,                               &
    1111          fluxsurf_sw,fluxsurfabs_sw,fluxtop_lw,               &
    … …  
    3636                              optichaze,haze_radproffix,&
    3737                              methane,carbox,cooling,nlte,strobel,&
    38                               ch4fix,vmrch4_proffix,vmrch4fix
     38                              ch4fix,vmrch4_proffix,vmrch4fix,&
     39                              callmufi
    3940      use optcv_mod, only: optcv
    4041      use optci_mod, only: optci
    … …  
    4445      use planetwide_mod, only: planetwide_maxval, planetwide_minval
    4546      use radcommon_h, only: wavev,wavei
     47      use mp2m_diagnostics
    4648      implicit none
    4749
    … …  
    8183      REAL,INTENT(IN) :: pt(ngrid,nlayer)          ! Air temperature (K).
    8284      REAL,INTENT(IN) :: zzlay(ngrid,nlayer)       ! Mid-layer altitude
     85      REAL,INTENT(IN) :: zzlev(ngrid,nlayer)       ! Altitude at the layer boundaries.
    8386      REAL,INTENT(IN) :: tsurf(ngrid)              ! Surface temperature (K).
    8487      REAL,INTENT(IN) :: fract(ngrid)              ! Fraction of day.
    … …  
    8992
    9093      ! OUTPUT
    91       REAL,INTENT(OUT) :: aerosol(ngrid,nlayer,naerkind) ! Aerosol tau at reference wavelenght.
    92       REAL,INTENT(OUT) :: dtlw(ngrid,nlayer)             ! Heating rate (K/s) due to LW radiation.
    93       REAL,INTENT(OUT) :: dtsw(ngrid,nlayer)             ! Heating rate (K/s) due to SW radiation.
    94       REAL,INTENT(OUT) :: fluxsurf_lw(ngrid)             ! Incident LW flux to surf (W/m2).
    95       REAL,INTENT(OUT) :: fluxsurf_sw(ngrid)             ! Incident SW flux to surf (W/m2)
    96       REAL,INTENT(OUT) :: fluxsurfabs_sw(ngrid)          ! Absorbed SW flux by the surface (W/m2). By MT2015.
    97       REAL,INTENT(OUT) :: fluxtop_lw(ngrid)              ! Outgoing LW flux to space (W/m2).
    98       REAL,INTENT(OUT) :: fluxabs_sw(ngrid)              ! SW flux absorbed by the planet (W/m2).
    99       REAL,INTENT(OUT) :: fluxtop_dn(ngrid)              ! Incident top of atmosphere SW flux (W/m2).
    100       REAL,INTENT(OUT) :: OLR_nu(ngrid,L_NSPECTI)        ! Outgoing LW radiation in each band (Normalized to the band width (W/m2/cm-1).
    101       REAL,INTENT(OUT) :: OSR_nu(ngrid,L_NSPECTV)        ! Outgoing SW radiation in each band (Normalized to the band width (W/m2/cm-1).
    102       REAL,INTENT(OUT) :: GSR_nu(ngrid,L_NSPECTV)        ! Surface SW radiation in each band (Normalized to the band width (W/m2/cm-1).
    103       REAL,INTENT(OUT) :: tau_col(ngrid)                 ! Diagnostic from aeropacity.
    104       REAL,INTENT(OUT) :: albedo_equivalent(ngrid)       ! Spectrally Integrated Albedo. For Diagnostic. By MT2015
     94      REAL,INTENT(OUT) :: dtau_aer(ngrid,nlayer,naerkind) ! Aerosol tau at reference wavelenght.
     95      REAL,INTENT(OUT) :: dtlw(ngrid,nlayer)              ! Heating rate (K/s) due to LW radiation.
     96      REAL,INTENT(OUT) :: dtsw(ngrid,nlayer)              ! Heating rate (K/s) due to SW radiation.
     97      REAL,INTENT(OUT) :: fluxsurf_lw(ngrid)              ! Incident LW flux to surf (W/m2).
     98      REAL,INTENT(OUT) :: fluxsurf_sw(ngrid)              ! Incident SW flux to surf (W/m2)
     99      REAL,INTENT(OUT) :: fluxsurfabs_sw(ngrid)           ! Absorbed SW flux by the surface (W/m2). By MT2015.
     100      REAL,INTENT(OUT) :: fluxtop_lw(ngrid)               ! Outgoing LW flux to space (W/m2).
     101      REAL,INTENT(OUT) :: fluxabs_sw(ngrid)               ! SW flux absorbed by the planet (W/m2).
     102      REAL,INTENT(OUT) :: fluxtop_dn(ngrid)               ! Incident top of atmosphere SW flux (W/m2).
     103      REAL,INTENT(OUT) :: OLR_nu(ngrid,L_NSPECTI)         ! Outgoing LW radiation in each band (Normalized to the band width (W/m2/cm-1).
     104      REAL,INTENT(OUT) :: OSR_nu(ngrid,L_NSPECTV)         ! Outgoing SW radiation in each band (Normalized to the band width (W/m2/cm-1).
     105      REAL,INTENT(OUT) :: GSR_nu(ngrid,L_NSPECTV)         ! Surface SW radiation in each band (Normalized to the band width (W/m2/cm-1).
     106      REAL,INTENT(OUT) :: tau_col(ngrid)                  ! Diagnostic from aeropacity.
     107      REAL,INTENT(OUT) :: albedo_equivalent(ngrid)        ! Spectrally Integrated Albedo. For Diagnostic. By MT2015
    105108      REAL,INTENT(OUT) :: int_dtaui(ngrid,nlayer,L_NSPECTI) ! VI optical thickness of layers within narrowbands for diags ().
    106109      REAL,INTENT(OUT) :: int_dtauv(ngrid,nlayer,L_NSPECTV) ! IR optical thickness of layers within narrowbands for diags ().
    … …  
    186189!$OMP THREADPRIVATE(QREFvis3d,QREFir3d)
    187190
    188 
    189191      ! Miscellaneous :
    190       real*8  temp,temp1,temp2,pweight
     192      real*8  temp,temp1,temp2,pweight,sig
    191193      character(len=10) :: tmp1
    192194      character(len=10) :: tmp2
    … …  
    331333         if (is_master) call system('rm -f surf_vals_long.out')
    332334
    333          call su_aer_radii(ngrid,nlayer,reffrad,nueffrad)
    334 
    335 
    336335!--------------------------------------------------
    337336!             Set up correlated k
    … …  
    454453!     Effective radius and variance of the aerosols
    455454!--------------------------------------------------
     455      ! Radiative Hazes
    456456      if (optichaze) then
    457          do iaer=1,naerkind
    458             if ((iaer.eq.iaero_haze)) then
    459                call su_aer_radii(ngrid,nlayer,reffrad(1,1,iaer), &
    460                             nueffrad(1,1,iaer))
     457         if (callmufi) then
     458            ! Spherical aerosols
     459            sig = 0.2
     460            where (mp2m_rc_sph(:,:) > 1e-10)
     461               reffrad(:,:,1) = mp2m_rc_sph(:,:) * exp(5.*sig**2 / 2.)
     462            elsewhere
     463               reffrad(:,:,1) = 0d0     
     464            endwhere
     465            if (exp(sig**2) - 1 > 0.1) then
     466               nueffrad(:,:,1) = exp(sig**2) - 1
     467            else
     468               nueffrad(:,:,1) = 0.1   
    461469            endif
    462          end do !iaer=1,naerkind.
    463          if (haze_radproffix) then
     470            ! Fractal aerosols
     471            sig = 0.35
     472            where (mp2m_rc_fra(:,:) > 1e-10)
     473               reffrad(:,:,2) = mp2m_rc_fra(:,:) * exp(5.*sig**2 / 2.)
     474               elsewhere
     475               reffrad(:,:,2) = 0d0
     476               endwhere
     477            if (exp(sig**2) - 1 > 0.1) then
     478               nueffrad(:,:,2) = exp(sig**2) - 1
     479            else
     480               nueffrad(:,:,2) = 0.1   
     481            endif
     482           
     483         else
     484            do iaer=1,naerkind
     485               if ((iaer.eq.iaero_haze)) then
     486               call su_aer_radii(ngrid,nlayer,reffrad(1,1,iaer),nueffrad(1,1,iaer))
     487               endif
     488            end do ! iaer = 1, naerkind.
    464489            if (haze_radproffix) then
    465                 call haze_reffrad_fix(ngrid,nlayer,zzlay,&
    466                     reffrad,nueffrad)
    467             endif
    468 
    469             print*, 'haze_radproffix=T : fixed profile for haze rad'
    470          else
    471             print*,'reffrad haze:',reffrad(1,1,iaero_haze)
    472             print*,'nueff haze',nueffrad(1,1,iaero_haze)
    473          endif
    474       endif
    475 
     490               call haze_reffrad_fix(ngrid,nlayer,zzlay,reffrad,nueffrad)
     491               if (is_master) print*, 'haze_radproffix=T : fixed profile for haze rad'
     492            else
     493               if (is_master) print*,'reffrad haze:',reffrad(1,1,iaero_haze)
     494               if (is_master) print*,'nueff haze',nueffrad(1,1,iaero_haze)
     495            endif ! end haze_radproffix
     496         endif ! end callmufi
     497      endif ! end radiative haze
    476498
    477499      ! How much light do we get ?
    … …  
    489511
    490512         ! Get aerosol optical depths.
    491          call aeropacity(ngrid,nlayer,nq,pplay,pplev, pt,pq,aerosol,      &
     513         call aeropacity(ngrid,nlayer,nq,pplay,pplev,zzlev,pt,pq,dtau_aer,      &
    492514            reffrad,nueffrad,QREFvis3d,QREFir3d,                             &
    493515            tau_col)
    … …  
    665687               pweight=(pplay(ig,L_NLAYRAD-k)-pplev(ig,L_NLAYRAD-k+1))/   &
    666688                       (pplev(ig,L_NLAYRAD-k)-pplev(ig,L_NLAYRAD-k+1))
    667                ! As 'aerosol' is at reference (visible) wavelenght we scale it as
     689               ! As 'dtau_aer' is at reference (visible) wavelenght we scale it as
    668690               ! it will be multplied by qxi/v in optci/v
    669                temp=aerosol(ig,L_NLAYRAD-k,iaer)/QREFvis3d(ig,L_NLAYRAD-k,iaer)
     691               temp=dtau_aer(ig,L_NLAYRAD-k,iaer)/QREFvis3d(ig,L_NLAYRAD-k,iaer)
    670692               tauaero(2*k+2,iaer)=max(temp*pweight,0.d0)
    671693               tauaero(2*k+3,iaer)=max(temp-tauaero(2*k+2,iaer),0.d0)
    … …  
    9861008            print*,'fluxtop_dn=',fluxtop_dn(ig)
    9871009            print*,'acosz=',acosz
    988             print*,'aerosol=',aerosol(ig,:,:)
     1010            print*,'dtau_aer=',dtau_aer(ig,:,:)
    9891011            print*,'temp=   ',pt(ig,:)
    9901012            print*,'pplay=  ',pplay(ig,:)

  • trunk/LMDZ.PLUTO/libf/phypluto/callcorrk_pluto_mod.F90

    r3572 r3585  
    77    subroutine callcorrk_pluto(icount,ngrid,nlayer,pq,nq,qsurf,   &
    88          albedo,emis,mu0,pplev,pplay,pt,   &
    9           zzlay,tsurf,fract,dist_star,aerosol,    &
     9          zzlay,zzlev,tsurf,fract,dist_star,dtau_aer,    &
    1010          dtlw,dtsw,fluxsurf_lw,    &
    1111          fluxsurf_sw,fluxtop_lw,fluxtop_sw,fluxtop_dn, &
    … …  
    2626      use callkeys_mod, only: optichaze,ch4fix,cooling,methane,nlte,&
    2727                              strobel,vmrch4_proffix,specOLR,vmrch4fix,&
    28                               haze_radproffix
     28                              haze_radproffix,callmufi
    2929      use optcv_pluto_mod, only: optcv_pluto
    3030      use optci_pluto_mod, only: optci_pluto
    … …  
    3232      use sfluxv_pluto_mod, only: sfluxv_pluto
    3333      use mod_phys_lmdz_para, only : is_master
    34 
     34      use mp2m_diagnostics
    3535
    3636      implicit none
    … …  
    6161      INTEGER icount
    6262      INTEGER ngrid,nlayer
    63       REAL aerosol(ngrid,nlayer,naerkind) ! aerosol opacity tau
    64       REAL albedo(ngrid)                    ! SW albedo
    65       REAL emis(ngrid)                      ! LW emissivity
    66       REAL pplay(ngrid,nlayer)            ! pres. level in GCM mid of layer
    67       REAL zzlay(ngrid,nlayer)            ! altitude at the middle of the layers
    68       REAL pplev(ngrid,nlayer+1)          ! pres. level at GCM layer boundaries
    69 
     63      REAL dtau_aer(ngrid,nlayer,naerkind) ! aerosol opacity tau
     64      REAL albedo(ngrid)                   ! SW albedo
     65      REAL emis(ngrid)                     ! LW emissivity
     66      REAL pplay(ngrid,nlayer)             ! pres. level in GCM mid of layer
     67      REAL zzlay(ngrid,nlayer)             ! Altitude at the middle of the layers
     68      REAL zzlev(ngrid,nlayer)             ! Altitude at the layer boundaries.
     69      REAL pplev(ngrid,nlayer+1)           ! pres. level at GCM layer boundaries
     70 
    7071      REAL pt(ngrid,nlayer)               ! air temperature (K)
    7172      REAL tsurf(ngrid)                     ! surface temperature (K)
    … …  
    135136      INTEGER ig,l,k,nw,iaer,irad
    136137
     138      real*8 sig
     139
    137140      real fluxtoplanet
    138141      real*8 taugsurf(L_NSPECTV,L_NGAUSS-1)
    … …  
    234237         call setspv            ! basic visible properties
    235238
     239         !--------------------------------------------------
     240         !     Effective radius and variance of the aerosols
     241         !--------------------------------------------------
    236242         ! Radiative Hazes
    237243         if (optichaze) then
    238 
    239            ! AF24: TODO check duplicate suaer_corrk called from physiq_mod
    240           !  print*,'optichaze: starting suaer_corrk'
    241           !  call suaer_corrk       ! set up aerosol optical properties
    242           !  print*,'ending suaer_corrk'
    243 
    244            !--------------------------------------------------
    245            !     Effective radius and variance of the aerosols
    246            !--------------------------------------------------
    247            do iaer=1,naerkind
     244           if (callmufi) then
     245            ! Spherical aerosols
     246            sig = 0.2
     247            WHERE(mp2m_rc_sph(:,:) > 1e-10)
     248               reffrad(:,:,1) = mp2m_rc_sph(:,:) * exp(5.*sig**2 / 2.)
     249            ELSEWHERE
     250               reffrad(:,:,1) = 0d0         
     251            ENDWHERE
     252            nueffrad(:,:,1) = exp(sig**2) - 1
     253            ! Fractal aerosols
     254            sig = 0.35
     255            WHERE(mp2m_rc_fra(:,:) > 1e-10)
     256               reffrad(:,:,2) = mp2m_rc_fra(:,:) * exp(5.*sig**2 / 2.)
     257            ELSEWHERE
     258               reffrad(:,:,2) = 0d0
     259            ENDWHERE
     260            nueffrad(:,:,2) = exp(sig**2) - 1
     261           
     262           else
     263            do iaer=1,naerkind
    248264              if ((iaer.eq.iaero_haze)) then
    249                call su_aer_radii(ngrid,nlayer,reffrad(1,1,iaer),    &
    250                   nueffrad(1,1,iaer))
    251                 ! write(*,*) "Not supported yet"
    252                 ! STOP
     265              call su_aer_radii(ngrid,nlayer,reffrad(1,1,iaer),nueffrad(1,1,iaer))
    253266              endif
    254            end do !iaer=1,naerkind.
    255            if (haze_radproffix) then
    256               call haze_reffrad_fix(ngrid,nlayer,zzlay, &
    257                   reffrad,nueffrad)
     267            end do ! iaer = 1, naerkind.
     268            if (haze_radproffix) then
     269              call haze_reffrad_fix(ngrid,nlayer,zzlay,reffrad,nueffrad)
    258270              if (is_master) print*, 'haze_radproffix=T : fixed profile for haze rad'
    259            else
     271            else
    260272              if (is_master) print*,'reffrad haze:',reffrad(1,1,iaero_haze)
    261273              if (is_master) print*,'nueff haze',nueffrad(1,1,iaero_haze)
    262            endif
    263          endif ! radiative haze
     274            endif ! end haze_radproffix
     275           endif ! end callmufi
     276         endif ! end radiative haze
    264277
    265278         Cmk= 0.01 * 1.0 / (g * mugaz * 1.672621e-27) ! q_main=1.0 assumed
    … …  
    305318
    306319        ! Get aerosol optical depths.
    307         call aeropacity(ngrid,nlayer,nq,pplay,pplev, pt,pq,aerosol,      &
     320        call aeropacity(ngrid,nlayer,nq,pplay,pplev,zzlev,pt,pq,dtau_aer,      &
    308321            reffrad,nueffrad,QREFvis3d,QREFir3d,                             &
    309322            tau_col)
    … …  
    471484                   (pplev(ig,L_NLAYRAD-k)-pplev(ig,L_NLAYRAD-k+1))
    472485               if (QREFvis3d(ig,L_NLAYRAD-k,iaer).ne.0) then
    473                  temp=aerosol(ig,L_NLAYRAD-k,iaer)/ &
     486                 temp=dtau_aer(ig,L_NLAYRAD-k,iaer)/ &
    474487                   QREFvis3d(ig,L_NLAYRAD-k,iaer)
    475488               else

  • trunk/LMDZ.PLUTO/libf/phypluto/callsedim_pluto.F90

    r3353 r3585  
    33                     pq, pdqfi, pdqsed,pdqs_sed,nq,pphi)
    44
    5       use radinc_h, only : naerkind
    65      use tracer_h, only: igcm_ch4_ice,igcm_co_ice,radius,rho_q
    76      use comcstfi_mod, only:  g

  • trunk/LMDZ.PLUTO/libf/phypluto/ch4cloud.F90

    r3252 r3585  
    1        SUBROUTINE ch4cloud(ngrid,nlay,naersize, ptimestep,  &
     1       SUBROUTINE ch4cloud(ngrid,nlay,ptimestep,  &
    22                     pplev,pplay,pdpsrf,pzlev,pzlay,pt,pdt, &
    33                     pq,pdq,pdqcloud,pdqscloud,pdtcloud,  &
    … …  
    4747      real pq(ngrid,nlay,nq)     ! traceur (kg/kg)
    4848      real pdq(ngrid,nlay,nq)    ! tendance avant condensation  (kg/kg.s-1)
    49       integer naersize   ! nombre de traceurs radiativement actifs (=naerkind)
    5049      integer nq         ! nombre de traceurs
    5150

  • trunk/LMDZ.PLUTO/libf/phypluto/cocloud.F90

    r3252 r3585  
    1        SUBROUTINE cocloud(ngrid,nlay,naersize, ptimestep,   &
     1       SUBROUTINE cocloud(ngrid,nlay,ptimestep,   &
    22                     pplev,pplay,pdpsrf,pzlev,pzlay,pt,pdt,   &
    33                     pq,pdq,pdqcloud,pdqscloud,pdtcloud,  &
    … …  
    4444      real pq(ngrid,nlay,nq)     ! traceur (kg/kg)
    4545      real pdq(ngrid,nlay,nq)    ! tendance avant condensation  (kg/kg.s-1)
    46       integer naersize   ! nombre de traceurs radiativement actifs (=naerkind)
    4746      integer nq         ! nombre de traceurs
    4847

  • trunk/LMDZ.PLUTO/libf/phypluto/condense_n2.F90

    r3539 r3585  
    66      pdqc,pdicen2)
    77
    8   use radinc_h, only : naerkind
    98  use comgeomfi_h
    109  use comcstfi_mod, only: g, r, cpp, pi

  • trunk/LMDZ.PLUTO/libf/phypluto/datafile_mod.F90

    r3561 r3585  
    1919      character(LEN=100),save :: config_mufi ='datagcm/microphysics/config.cfg'
    2020!$OMP THREADPRIVATE(config_mufi)
     21      character(len=300),save :: aersprop_file
     22      character(len=300),save :: aerfprop_file
    2123
    2224      ! surfdir stores planetary topography, albedo, etc. (surface.nc files)

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

    r3572 r3585  
    7878
    7979  real dTstrat
    80   real,allocatable :: aerosol(:,:) ! aerosol tau (kg/kg)
     80  real,allocatable :: dtau_aer(:,:) ! aerosol tau (kg/kg)
    8181  real OLR_nu(1,L_NSPECTI)
    8282  real OSR_nu(1,L_NSPECTV)
    … …  
    320320  !write(*,*) 1,llm,nq,0,86400.0,1,1.0,latitude,longitude,cell_area,rad,g,r,cpp
    321321
    322   ! initialise naerkind (from callphys.def) and allocate aerosol(:,:)
     322  ! initialise naerkind (from callphys.def) and allocate dtau_aer(:,:)
    323323  naerkind=0 !default
    324324  call getin("naerkind",naerkind)
    325   allocate(aerosol(llm,naerkind))
    326   aerosol(:,:)=0
     325  allocate(dtau_aer(llm,naerkind))
     326  dtau_aer(:,:)=0
    327327
    328328  do iq=1,nq
    … …  
    374374          albedo_wv,albedo_equivalent,                    &
    375375          emis,mu0,plev,play,temp,                        &
    376           tsurf,fract,dist_star,aerosol,muvar,            &
     376          tsurf,fract,dist_star,dtau_aer,muvar,            &
    377377          dtlw,dtsw,fluxsurf_lw,fluxsurf_sw,              &
    378378          fluxsurfabs_sw,fluxtop_lw,                      &
    … …  
    385385     !     albedo_wv,'D',albedo_equivalent,'E',                    &
    386386     !     emis,'F',mu0,'G',plev,'H',play,'I',temp,'J',                        &
    387      !     tsurf,'K',fract,'L',dist_star,'M',aerosol,'N',muvar,'O',            &
     387     !     tsurf,'K',fract,'L',dist_star,'M',dtau_aer,'N',muvar,'O',            &
    388388     !     dtlw,'P',dtsw,'Q',fluxsurf_lw,'R',fluxsurf_sw,'S',              &
    389389     !     fluxsurfabs_sw,'T',fluxtop_lw,'U',                      &
    … …  
    426426  call callcorrk(1,nlayer,q,nq,qsurf,                          &
    427427       albedo_wv,albedo_equivalent,emis,mu0,plev,play,temp,    &
    428        tsurf,fract,dist_star,aerosol,muvar,                    &
     428       tsurf,fract,dist_star,dtau_aer,muvar,                    &
    429429       dtlw,dtsw,fluxsurf_lw,fluxsurf_sw,fluxsurfabs_sw,       &
    430430       fluxtop_lw, fluxabs_sw,fluxtop_dn,OLR_nu,OSR_nu,GSR_nu, &
    … …  
    437437     !     albedo_wv,'D',albedo_equivalent,'E',                    &
    438438     !     emis,'F',mu0,'G',plev,'H',play,'I',temp,'J',                        &
    439      !     tsurf,'K',fract,'L',dist_star,'M',aerosol,'N',muvar,'O',            &
     439     !     tsurf,'K',fract,'L',dist_star,'M',dtau_aer,'N',muvar,'O',            &
    440440     !     dtlw,'P',dtsw,'Q',fluxsurf_lw,'R',fluxsurf_sw,'S',              &
    441441     !     fluxsurfabs_sw,'T',fluxtop_lw,'U',                      &

  • trunk/LMDZ.PLUTO/libf/phypluto/dyn1d/rcm1d.F

    r3542 r3585  
    66      use mod_grid_phy_lmdz, only : regular_lonlat
    77      use infotrac, only: nqtot, tname
    8       use tracer_h, only: noms, is_condensable
     8      use tracer_h, only: noms
    99      use surfdat_h, only: albedodat, phisfi, dryness,
    1010     &                     zmea, zstd, zsig, zgam, zthe,

  • trunk/LMDZ.PLUTO/libf/phypluto/hazecloud.F90

    r3390 r3585  
    55!          zdqphot_ch4,zdqconv_prec,declin,zdqhaze_col)
    66
    7       use radinc_h, only : naerkind
    87      use comgeomfi_h
    98      use comcstfi_mod, only: pi, g

  • trunk/LMDZ.PLUTO/libf/phypluto/inifis_mod.F90

    r3572 r3585  
    1212  use radinc_h, only: ini_radinc_h, naerkind
    1313  use radcommon_h, only: ini_radcommon_h
    14   use radii_mod, only: radfixed, Nmix_n2
    15   use datafile_mod, only: datadir,config_mufi,hazeprop_file,hazerad_file,hazemmr_file,hazedens_file
     14  use radii_mod, only: radfixed
     15  use datafile_mod, only: datadir,hazeprop_file,hazerad_file,hazemmr_file,hazedens_file, &
     16                          config_mufi, aersprop_file, aerfprop_file
    1617  use comdiurn_h, only: sinlat, coslat, sinlon, coslon
    1718  use comgeomfi_h, only: totarea, totarea_planet
    … …  
    687688     call getin_p("config_mufi",config_mufi)
    688689     if (is_master) write(*,*)" config_mufi = ",config_mufi
     690
     691     if (is_master) write(*,*) "Spherical aerosol optical properties datafile"
     692     aersprop_file="optprop_rannou_r2-200nm_nu003.dat"  ! default file
     693     call getin_p("aersprop_file",aersprop_file)
     694     if (is_master) write(*,*) trim(rname)//" aersprop_file = ",trim(aersprop_file)
     695
     696     if (is_master) write(*,*) "Fractal aerosol optical properties datafile"
     697     aerfprop_file="optprop_rannou_fractal_r010nm_N1_1e4_d2.dat"  ! default file
     698     call getin_p("aerfprop_file",aerfprop_file)
     699     if (is_master) write(*,*) trim(rname)//" aerfprop_file = ",trim(aerfprop_file)
    689700 
    690701     if (is_master) write(*,*) "Use haze production from CH4 photolysis or production rate?"
    … …  
    12851296       call abort_physic(rname, 'if haze are on, tracers must be on!', 1)
    12861297     endif
     1298     if ((callmufi).and.(haze)) then
     1299       call abort_physic(rname, 'if haze are on, microphysics should be desactivated!', 1)
     1300     endif
     1301     if ((haze).and.(naerkind.gt.1)) then
     1302      call abort_physic(rname, 'if haze are on, naerkind must be equal to 1!', 1)
     1303     endif
     1304     if ((callmufi).and..not.(naerkind.gt.1)) then
     1305      call abort_physic(rname, 'if microphysics is on, naerkind must be > 1!', 1)
     1306     endif
     1307     if (.not.(callmufi.or.haze).and.(optichaze)) then
     1308      call abort_physic(rname, 'if microphysics and haze are off, optichaze must be desactivated!', 1)
     1309     endif
     1310     if ((callmufi.and.call_haze_prod_pCH4).and..not.(methane)) then
     1311      call abort_physic(rname, 'if haze production from CH4 photolysis is on, methane must be activated!', 1)
     1312     endif
    12871313     if (haze_proffix.and.sedimentation) then
    12881314         call abort_physic(rname, 'if haze profile is set, sedimentation must be deactivated', 1)

  • trunk/LMDZ.PLUTO/libf/phypluto/initracer.F90

    r3572 r3585  
    123123         ALLOCATE(is_recomb_qotf(nqtot))
    124124       ENDIF
    125        IF (.NOT. allocated(is_condensable)) allocate(is_condensable(nq)) !LT
    126        IF (.NOT. allocated(is_rgcs)) allocate(is_rgcs(nq)) !LT
    127        IF (.NOT. allocated(constants_mass)) allocate(constants_mass(nq))
    128        IF (.NOT. allocated(constants_delta_gasH)) allocate(constants_delta_gasH(nq))
    129        IF (.NOT. allocated(constants_Tref)) allocate(constants_Tref(nq))
    130        IF (.NOT. allocated(constants_Pref)) allocate(constants_Pref(nq))
    131        IF (.NOT. allocated(constants_RLVTT_generic)) allocate(constants_RLVTT_generic(nq))
    132        IF (.NOT. allocated(constants_RCPV_generic)) allocate(constants_RCPV_generic(nq))
    133125
    134126       !! initialization
    … …  
    148140       radius(:)=0.
    149141       qext(:)=0.
    150 
    151        ! For condensable tracers, by Lucas Teinturier and Noé Clément (2022)
    152 
    153        is_condensable(:)= 0
    154        is_rgcs(:) = 0
    155        constants_mass(:)=0
    156        constants_delta_gasH(:)=0
    157        constants_Tref(:)=0
    158        constants_Pref(:)=0
    159        constants_RLVTT_generic(:)=0
    160        constants_RCPV_generic(:)=0
    161142
    162143       rho_q(:) = 0. !need to be init here if we want to read it from modern traceur with get_tracdat
    … …  
    403384      if (is_master) close(407)
    404385
    405       ! Get specific data of condensable tracers
    406       do iq=1,nq
    407         if((is_condensable(iq)==1)) then
    408                 write(*,*) "There is a specie which is condensable, for generic condensation : ", noms(iq)
    409                 write(*,*) 'looking specie parameters for : ',noms(iq)(1:len(trim(noms(iq)))-4)
    410                 ! call specie_parameters_table(noms(iq)(1:len(trim(noms(iq)))-4))
    411                 ! constants_mass(iq)=m
    412                 ! constants_delta_gasH(iq)=delta_gasH
    413                 ! constants_Tref(iq)=Tref
    414                 ! constants_Pref(iq)=Pref
    415                 ! constants_RLVTT_generic(iq)=RLVTT_generic
    416                 ! constants_RCPV_generic(iq)=RCPV_generic
    417         else
    418                 write(*,*) "This tracer is not condensable, for generic condensation :  : ", noms(iq)
    419                 write(*,*) "We keep condensable constants at zero"
    420         endif !(is_condensable(iq)==1) .and. (index(noms(iq),"vap") .ne. 0))
    421       enddo ! iq=1,nq
    422 
    423386      ! Calculate number of species in the chemistry
    424387      nesp = sum(is_chim)
    425388      write(*,*) 'Number of species in the chemistry nesp = ',nesp
    426 
    427       ! Calculate number of generic tracers
    428       ngt = sum(is_condensable)
    429       write(*,*) 'Number of generic tracer is  ngt = ',ngt
    430 
    431       ! Calculate number of radiative generic condensable species
    432       n_rgcs = sum(is_rgcs)
    433       write(*,*)'Number of Radiative Generic Condensable Species is n_rgcs = ',n_rgcs
    434       if (n_rgcs> ngt/2) then
    435         write(*,*) 'You have more Radiative Generic Condensable Species than Generic Condensable Species'
    436         write(*,*)'This is not possible: check your Modern traceur.def'
    437         call abort_physic("initracer, issue with # of RGCS and GCS")
    438       endif
    439389
    440390      ! Calculate number of microphysical tracer
    … …  
    461411      write(*,*) "lw_n2 = ", lw_n2
    462412
    463       if (haze) then
     413      if (callmufi) then
     414        if (optichaze) then
     415          iaero_haze = 2
     416          write(*,*) 'Microphysical moment model'
     417          write(*,*) '--- number of haze aerosol = ', iaero_haze
     418        endif ! end optichaze
     419     
     420      elseif (haze) then
    464421        ! the sedimentation radius remains radius(igcm_haze)
    465422        if (fractal) then
    466            nmono=nb_monomer
     423          nmono=nb_monomer
    467424        else
    468            nmono=1
     425          nmono=1
    469426        endif ! end fractal
    470 
    471         ia=0
     427       
     428        ia = 0
    472429        if (optichaze) then
    473           ia=ia+1
    474           iaero_haze=ia
     430          ia = ia + 1
     431          iaero_haze = ia
    475432          write(*,*) '--- number of haze aerosol = ', iaero_haze
    476 
    477           block=0  ! Only one type of haze is active : the first one set in traceur.def
     433          block = 0 ! Only one type of haze is active : the first one set in traceur.def
    478434          do iq=1,nq
    479435            tracername=noms(iq)
    … …  
    484440              write(*,*) "i_haze=",i_haze
    485441              write(*,*) "Careful: if you set many haze traceurs in &
    486                 traceur.def,only ",tracername," will be radiatively active &
    487                 (first one in traceur.def)"
     442              traceur.def,only ",tracername," will be radiatively active &
     443              (first one in traceur.def)"
    488444            endif
    489           enddo
     445            enddo
    490446        endif ! end optichaze
    491447      endif ! end callmufi or haze
    … …  
    602558              is_recomb_qotf(iq)
    603559          end if
    604           !option is_condensable (LT)
    605           if (index(tracline,'is_condensable=') /=0) then
    606             read(tracline(index(tracline,'is_condensable=') &
    607               +len('is_condensable='):),*) is_condensable(iq)
    608             write(*,*) ' Parameter value (traceur.def) :'// &
    609               ' is_condensable=', is_condensable(iq)
    610           else
    611               write(*,*) ' Parameter value (default)     :'// &
    612               ' is_condensable=', is_condensable(iq)
    613           endif
    614560          !option radius
    615561          if (index(tracline,'radius=') .ne. 0) then
    … …  
    632578              ' rho=', rho_q(iq)
    633579          endif
    634           !option is_rgcs
    635           if (index(tracline,'is_rgcs') .ne. 0) then
    636             read(tracline(index(tracline,'is_rgcs=') &
    637               +len('is_rgcs='):),*) is_rgcs(iq)
    638             write(*,*)'Parameter value (traceur.def) :'// &
    639               'is_rgcs=',is_rgcs(iq)
    640           else
    641             write(*,*)'Parameter value (default) : '// &
    642               'is_rgcs = ',is_rgcs(iq)
    643           endif
    644580      end subroutine get_tracdat
    645581

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

    r3572 r3585  
    2121      use radcommon_h, only: sigma, glat, grav, BWNV, WNOI, DWNI, DWNV, WNOV
    2222      use suaer_corrk_mod, only: suaer_corrk
    23       use radii_mod, only: su_aer_radii,haze_reffrad_fix
    2423      use aerosol_mod, only: i_haze, haze_prof
    2524      use surfdat_h, only: phisfi, zmea, zstd, zsig, zgam, zthe, &
    … …  
    3534                          igcm_co_gas,igcm_co_ice,igcm_prec_haze,lw_n2,lw_ch4,lw_co,&
    3635                          alpha_lift, alpha_devil, qextrhor, &
    37                           nesp, is_chim, is_condensable
     36                          nesp, is_chim
    3837      use time_phylmdz_mod, only: ecritphy, iphysiq, nday
    3938      use phyetat0_mod, only: phyetat0,tab_cntrl_mod
    … …  
    286285!     for the "naerkind" optically active aerosols:
    287286
    288       real,save,allocatable :: aerosol(:,:,:) ! Aerosols
    289 !$OMP THREADPRIVATE(aerosol)
     287      real,save,allocatable :: dtau_aer(:,:,:) ! Aerosols
     288!$OMP THREADPRIVATE(dtau_aer)
    290289      real zh(ngrid,nlayer)               ! Potential temperature (K).
    291290      real pw(ngrid,nlayer)               ! Vertical velocity (m/s). (NOTE : >0 WHEN DOWNWARDS !!)
    … …  
    464463      real muvar(ngrid,nlayer+1) ! For Runaway Greenhouse 1D study. By RW
    465464
    466       real,save,allocatable :: reffcol(:,:)
    467 !$OMP THREADPRIVATE(reffcol)
    468 
    469465!  Non-oro GW tendencies
    470466      REAL d_u_hin(ngrid,nlayer), d_v_hin(ngrid,nlayer)
    … …  
    535531         ! allocate related local arrays
    536532         ! (need be allocated instead of automatic because of "naerkind")
    537          allocate(aerosol(ngrid,nlayer,naerkind))
    538          allocate(reffcol(ngrid,naerkind))
     533         allocate(dtau_aer(ngrid,nlayer,naerkind))
    539534
    540535#ifdef CPP_XIOS
    … …  
    10221017               if (oldplutocorrk) then
    10231018                  call callcorrk_pluto(icount,ngrid,nlayer,pq,nq,qsurf,          &
    1024                                albedo(:,1),emis,mu0,pplev,pplay,pt,                   &
    1025                                zzlay,tsurf,fract,dist_star,aerosol,              &
     1019                               albedo(:,1),emis,mu0,pplev,pplay,pt,              &
     1020                               zzlay,zzlev,tsurf,fract,dist_star,dtau_aer,       &
    10261021                               zdtlw,zdtsw,fluxsurf_lw,fluxsurf_sw,fluxtop_lw,   &
    10271022                               fluxabs_sw,fluxtop_dn,reffrad,tau_col,ptime,pday, &
    … …  
    10341029                call callcorrk(ngrid,nlayer,pq,nq,qsurf,  &
    10351030                              albedo,albedo_equivalent,emis,mu0,pplev,pplay,pt,   &
    1036                               zzlay,tsurf,fract,dist_star,aerosol,muvar,                &
     1031                              zzlay,zzlev,tsurf,fract,dist_star,dtau_aer,muvar,   &
    10371032                              zdtlw,zdtsw,fluxsurf_lw,fluxsurf_sw,                &
    10381033                              fluxsurfabs_sw,fluxtop_lw,                          &
    … …  
    14361431         zdqch4cloud(:,:,:)=0.
    14371432         if ((methane).and.(metcloud).and.(.not.fast)) THEN
    1438             call ch4cloud(ngrid,nlayer,naerkind,ptimestep,  &
     1433            call ch4cloud(ngrid,nlayer,ptimestep,  &
    14391434                      pplev,pplay,pdpsrf,zzlev,zzlay, pt,pdt,  &
    14401435                      pq,pdq,zdqch4cloud,zdqsch4cloud,zdtch4cloud,   &
    … …  
    14691464         zdqcocloud(:,:,:)=0.
    14701465         IF ((carbox).and.(monoxcloud).and.(.not.fast)) THEN
    1471             call cocloud(ngrid,nlayer,naerkind,ptimestep,   &
     1466            call cocloud(ngrid,nlayer,ptimestep,   &
    14721467                      pplev,pplay,pdpsrf,zzlev,zzlay, pt,pdt,  &
    14731468                      pq,pdq,zdqcocloud,zdqscocloud,zdtcocloud,   &
    … …  
    18451840               print*,'fluxtop_dn=',fluxtop_dn(ig)
    18461841               print*,'tau_col=',tau_col(ig)
    1847                print*,'aerosol=',aerosol(ig,:,:)
     1842               print*,'dtau_aer=',dtau_aer(ig,:,:)
    18481843               print*,'temp=   ',pt(ig,:)
    18491844               print*,'pplay=  ',pplay(ig,:)

  • trunk/LMDZ.PLUTO/libf/phypluto/prodhaze.F90

    r3247 r3585  
    33          flym_bot,flym_sol_bot,flym_ipm_bot,flym_sol,flym_ipm)
    44
    5       use radinc_h, only : naerkind
    65      use comgeomfi_h
    76      use comcstfi_mod, only: pi, g

  • trunk/LMDZ.PLUTO/libf/phypluto/radii_mod.F90

    r3572 r3585  
    44!  module to centralize the radii calculations for aerosols
    55!==================================================================
    6 
    7 !     N2 cloud properties (initialized in inifis)
    8       real,save :: Nmix_n2 ! Number mixing ratio of N2 ice particles
    9 !$OMP THREADPRIVATE(Nmix_n2)
    106
    117      ! flag to specify if we assume a constant fixed radius for particles
    … …  
    3026      use datafile_mod, only: hazerad_file
    3127      use aerosol_mod, only: iaero_haze, i_haze
    32       use tracer_h, only: radius, nqtot, is_rgcs, nmono
     28      use tracer_h, only: radius, nqtot, nmono
    3329      Implicit none
    3430
    … …  
    6763      USE tracer_h, only:rho_n2,nmono
    6864      use comcstfi_mod, only: pi
    69       use aerosol_mod, only: iaero_haze, i_haze
     65      use aerosol_mod, only: iaero_haze
    7066      use datafile_mod
    7167      Implicit none

  • trunk/LMDZ.PLUTO/libf/phypluto/spreadglacier_paleo.F90

    r3412 r3585  
    66      use comgeomfi_h
    77      use geometry_mod, only: latitude, longitude, cell_area
    8       use radinc_h, only : naerkind
    98      use tracer_h, only: igcm_n2,igcm_ch4_ice,igcm_co_ice
    109

  • trunk/LMDZ.PLUTO/libf/phypluto/spreadglacier_simple.F90

    r3473 r3585  
    22
    33      use geometry_mod, only: latitude, longitude, cell_area
    4       use radinc_h, only : naerkind
    54      use surfdat_h, only: phisfi
    65      use tracer_h, only: igcm_n2

  • trunk/LMDZ.PLUTO/libf/phypluto/suaer_corrk.F90

    r3572 r3585  
    1010      use radinc_h,    only: L_NSPECTI,L_NSPECTV,nsizemax,naerkind
    1111      use radcommon_h, only: blamv,blami,lamrefir,lamrefvis
    12       use datafile_mod, only: datadir, aerdir, hazeprop_file
     12      use datafile_mod, only: datadir, aerdir, &
     13                              hazeprop_file, aersprop_file, aerfprop_file
    1314
    1415      ! outputs
    … …  
    1718      use radcommon_h, only: qrefvis,qrefir,omegarefir !,omegarefvis
    1819      use aerosol_mod, only: iaero_haze
    19       use callkeys_mod, only: tplanet
     20      use callkeys_mod, only: tplanet, callmufi
    2021      use tracer_h, only: noms
    2122
    … …  
    129130      allocate(file_id(naerkind,2))
    130131
    131       do iaer=1,naerkind
    132          !     naerkind=1: haze
    133          if (iaer.eq.1) then
    134 
    135          ! Only one table of optical properties :
    136          if (is_master) write(*,*)'Suaer haze optical properties, using: ', &
    137                            TRIM(hazeprop_file)
    138 
    139          ! visible
    140          file_id(iaer,1)=TRIM(hazeprop_file)
    141          ! infrared
    142          file_id(iaer,2)=file_id(iaer,1)
    143 
    144          lamrefvis(iaer)=0.607E-6   ! reference wavelength for opacity vis pivot wavelength Cheng et al 2017
    145          lamrefir(iaer)=2.E-6   !  reference wavelength for opacity IR (in the LEISA range)
    146 
     132
     133      if (callmufi) then
     134         if (is_master) then
     135            write(*,*)'Suaer spherical aerosols optical properties, using: ', &
     136                       TRIM(aersprop_file)
     137            write(*,*)'Suaer fractal aerosols optical properties, using: ', &
     138                       TRIM(aerfprop_file)
    147139         endif
    148       enddo
    149 
    150       ! IF (naerkind .gt. 1) THEN ! AF24: ?
    151       !    print*,'naerkind = ',naerkind
    152       !    print*,'but we only have data for 1 type, exiting.'
    153       !    call abort
    154       ! ENDIF
     140         ! Visible
     141         file_id(1,1)=TRIM(aersprop_file)
     142         file_id(2,1)=TRIM(aerfprop_file)
     143         ! Infrared
     144         file_id(1,2)=file_id(1,1)
     145         file_id(2,2)=file_id(2,1)
     146
     147         do iaer=1,naerkind
     148            lamrefvis(iaer) = 0.607E-6 ! Reference wavelength for opacity vis pivot wavelength Cheng et al 2017
     149            lamrefir(iaer)  = 2.E-6    ! Reference wavelength for opacity IR (in the LEISA range)
     150         enddo ! end iaer=1,naerkind
     151         
     152      else
     153         do iaer=1,naerkind
     154            if (iaer.eq.1) then
     155               ! Only one table of optical properties :
     156               if (is_master) write(*,*)'Suaer haze optical properties, using: ', &
     157                                       TRIM(hazeprop_file)
     158               ! Visible
     159               file_id(iaer,1)=TRIM(hazeprop_file)
     160               ! Infrared
     161               file_id(iaer,2)=file_id(iaer,1)
     162
     163               lamrefvis(iaer) = 0.607E-6 ! Reference wavelength for opacity vis pivot wavelength Cheng et al 2017
     164               lamrefir(iaer)  = 2.E-6    ! Reference wavelength for opacity IR (in the LEISA range)
     165            endif ! end iaer.eq.1
     166
     167            if (naerkind.gt.1) then
     168               ! Visible
     169               file_id(iaer,1)=file_id(1,1)
     170               ! Infrared
     171               file_id(iaer,2)=file_id(1,1)
     172
     173               lamrefvis(iaer) = lamrefvis(1) ! Reference wavelength for opacity vis pivot wavelength Cheng et al 2017
     174               lamrefir(iaer)  = lamrefir(1)  ! Reference wavelength for opacity IR (in the LEISA range)
     175            endif ! end naerkind.gt.1
     176         enddo ! end iaer=1,naerkind
     177      endif ! end callmufi
    155178
    156179!------------------------------------------------------------------

  • trunk/LMDZ.PLUTO/libf/phypluto/tracer_h.F90

    r3572 r3585  
    1313       integer, save :: nqtot  ! total number of tracers
    1414       integer, save :: nesp   ! number of species in the chemistry
    15        integer, save :: ngt    ! number of generic tracers
    16        integer, save :: n_rgcs ! number of Radiative Generic Condensable Species
    17 !$OMP THREADPRIVATE(nqtot,nesp,ngt,n_rgcs)
     15!$OMP THREADPRIVATE(nqtot,nesp)
    1816
    1917       logical :: moderntracdef=.false. ! Standard or modern traceur.def
    … …  
    5452       integer, save, allocatable :: is_recomb_qotf(:) ! 1 if tracer recombination is done on-the-fly, else 0 (if 1, must have is_recomb_qset=0)
    5553!$OMP THREADPRIVATE(is_recomb,is_recomb_qset,is_recomb_qotf)
    56        integer, save, allocatable :: is_condensable(:)      ! 1 if tracer is generic, else 0 (added LT)
    57        integer,save,allocatable :: is_rgcs(:)               ! 1 if tracer is a radiative generic condensable specie, else 0 (added LT 2022)
    58 !$OMP THREADPRIVATE(is_condensable,is_rgcs)   !also added by LT
    5954       ! Lists of constants for condensable tracers
    60        real, save, allocatable :: constants_mass(:)                 ! molecular mass of the specie (g/mol)
    61        real, save, allocatable :: constants_delta_gasH(:)           ! Enthalpy of vaporization (J/mol)
    62        real, save, allocatable :: constants_Tref(:)                 ! Ref temperature for Clausis-Clapeyron (K)
    63        real, save, allocatable :: constants_Pref(:)                 ! Reference pressure for Clausius Clapeyron (Pa)
    64        real, save, allocatable :: constants_RLVTT_generic(:)        ! Latent heat of vaporization (J/kg)
    65        real, save, allocatable :: constants_RCPV_generic(:)         ! specific heat capacity of the tracer vapor at Tref
    66 !$OMP THREADPRIVATE(constants_mass,constants_delta_gasH,constants_Tref)
    67 !$OMP THREADPRIVATE(constants_Pref)
    68 !$OMP THREADPRIVATE(constants_RLVTT_generic,constants_RCPV_generic)
    6955
    7056! tracer indexes: these are initialized in initracer and should be 0 if the
Note: See TracChangeset for help on using the changeset viewer.