Index: LMDZ6/branches/LMDZ_ECRad/libf/phylmd/ecrad/radiation_aerosol_optics.F90 =================================================================== --- LMDZ6/branches/LMDZ_ECRad/libf/phylmd/ecrad/radiation_aerosol_optics.F90 (revision 4187) +++ LMDZ6/branches/LMDZ_ECRad/libf/phylmd/ecrad/radiation_aerosol_optics.F90 (revision 4188) @@ -89,4 +89,5 @@ & IAerosolClassUndefined, IAerosolClassIgnored, & & IAerosolClassHydrophobic, IAerosolClassHydrophilic + USE phys_local_var_mod, ONLY: rhcl integer, intent(in) :: nlev ! number of model levels @@ -187,5 +188,6 @@ end if - call gas%get(IH2O, IMassMixingRatio, h2o_mmr, istartcol=istartcol) +!AI juin 2022 + !call gas%get(IH2O, IMassMixingRatio, h2o_mmr, istartcol=istartcol) ! Loop over position @@ -195,6 +197,9 @@ ! saturation and the index to the relative-humidity index of ! hydrophilic-aerosol data - rh = h2o_mmr(jcol,jlev) / thermodynamics%h2o_sat_liq(jcol,jlev) - irh = ao%calc_rh_index(rh) +! AI juin 2022 +! rh = h2o_mmr(jcol,jlev) / thermodynamics%h2o_sat_liq(jcol,jlev) +! irh = ao%calc_rh_index(rh) + irh = ao%calc_rh_index(rhcl(jcol,jlev)) +! print*,'irh=',irh factor = ( thermodynamics%pressure_hl(jcol,jlev+1) & Index: LMDZ6/branches/LMDZ_ECRad/libf/phylmd/ecrad/radiation_scheme.F90 =================================================================== --- LMDZ6/branches/LMDZ_ECRad/libf/phylmd/ecrad/radiation_scheme.F90 (revision 4187) +++ LMDZ6/branches/LMDZ_ECRad/libf/phylmd/ecrad/radiation_scheme.F90 (revision 4188) @@ -13,5 +13,5 @@ SUBROUTINE RADIATION_SCHEME & ! Inputs - & (KIDIA, KFDIA, KLON, KLEV, KAEROLMDZ, NSW, & + & (KIDIA, KFDIA, KLON, KLEV, KAEROSOL, NSW, & & IDAY, TIME, & & PSOLAR_IRRADIANCE, & @@ -117,9 +117,10 @@ INTEGER(KIND=JPIM),INTENT(IN) :: KLEV ! Number of levels !INTEGER, INTENT(IN) :: KLON, KLEV -INTEGER(KIND=JPIM),INTENT(IN) :: KAEROLMDZ ! Number of aerosol types +!INTEGER(KIND=JPIM),INTENT(IN) :: KAEROLMDZ ! Number of aerosol types +INTEGER(KIND=JPIM),INTENT(IN) :: KAEROSOL INTEGER(KIND=JPIM),INTENT(IN) :: NSW ! Numbe of bands ! AI ATTENTION -INTEGER, PARAMETER :: KAEROSOL = 12 +!INTEGER, PARAMETER :: KAEROSOL = 12 ! *** Single-level fields @@ -313,6 +314,6 @@ if (lprint_input) then print*,'********** Verification des entrees *************' - print*,'KIDIA, KFDIA, KLON, KLEV, KAEROLMDZ, NSW =', & - KIDIA, KFDIA, KLON, KLEV, KAEROLMDZ, NSW + print*,'KIDIA, KFDIA, KLON, KLEV, KAEROSOL, NSW =', & + KIDIA, KFDIA, KLON, KLEV, KAEROSOL, NSW print*,'IDAY, TIME =', IDAY, TIME print*,'PSOLAR_IRRADIANCE =', PSOLAR_IRRADIANCE @@ -510,9 +511,9 @@ print*,'******** AEROSOLS (allocate + input) **************************************' -IF (NAERMACC > 0) THEN +!IF (NAERMACC > 0) THEN CALL aerosol%allocate(KLON, 1, KLEV, KAEROSOL) ! MACC climatology -ELSE - CALL aerosol%allocate(KLON, 1, KLEV, 6) ! Tegen climatology -ENDIF +!ELSE +! CALL aerosol%allocate(KLON, 1, KLEV, 6) ! Tegen climatology +!ENDIF ! Compute the dry mass of each layer neglecting humidity effects, in ! kg m-2, needed to scale some of the aerosol inputs Index: LMDZ6/branches/LMDZ_ECRad/libf/phylmd/ecrad/radiation_setup.F90 =================================================================== --- LMDZ6/branches/LMDZ_ECRad/libf/phylmd/ecrad/radiation_setup.F90 (revision 4187) +++ LMDZ6/branches/LMDZ_ECRad/libf/phylmd/ecrad/radiation_setup.F90 (revision 4188) @@ -238,8 +238,8 @@ ! hydrophobic black carbon, 11 ammonium sulphate, 12 inactive ! SO2 - rad_config%n_aerosol_types = 12 - if (lprint_setp) then - print*,'rad_config%n_aerosol_types =', rad_config%n_aerosol_types - endif +! rad_config%n_aerosol_types = 12 +! if (lprint_setp) then +! print*,'rad_config%n_aerosol_types =', rad_config%n_aerosol_types +! endif ! Indices to the aerosol optical properties in @@ -247,22 +247,22 @@ ! index hydrophilic aerosol types and positive numbers index ! hydrophobic aerosol types - rad_config%i_aerosol_type_map = 0 ! There can be up to 256 types +! rad_config%i_aerosol_type_map = 0 ! There can be up to 256 types ! if (lprint_setp) then ! print*,'rad_config%i_aerosol_type_map =', rad_config%i_aerosol_type_map ! endif - rad_config%i_aerosol_type_map(1:12) = (/ & - & -1, & ! Sea salt, size bin 1 (OPAC) - & -2, & ! Sea salt, size bin 2 (OPAC) - & -3, & ! Sea salt, size bin 3 (OPAC) - & 7, & ! Desert dust, size bin 1 (Woodward 2001) - & 8, & ! Desert dust, size bin 2 (Woodward 2001) - & 9, & ! Desert dust, size bin 3 (Woodward 2001) - & -4, & ! Hydrophilic organic matter (OPAC) - & 10, & ! Hydrophobic organic matter (OPAC) - & 11, & ! Black carbon (Boucher) - & 11, & ! Black carbon (Boucher) - & -5, & ! Ammonium sulphate (OPAC) - & 14 /) ! Stratospheric sulphate (hand edited from OPAC) +! rad_config%i_aerosol_type_map(1:12) = (/ & +! & -1, & ! Sea salt, size bin 1 (OPAC) +! & -2, & ! Sea salt, size bin 2 (OPAC) +! & -3, & ! Sea salt, size bin 3 (OPAC) +! & 7, & ! Desert dust, size bin 1 (Woodward 2001) +! & 8, & ! Desert dust, size bin 2 (Woodward 2001) +! & 9, & ! Desert dust, size bin 3 (Woodward 2001) +! & -4, & ! Hydrophilic organic matter (OPAC) +! & 10, & ! Hydrophobic organic matter (OPAC) +! & 11, & ! Black carbon (Boucher) +! & 11, & ! Black carbon (Boucher) +! & -5, & ! Ammonium sulphate (OPAC) +! & 14 /) ! Stratospheric sulphate (hand edited from OPAC) ! if (lprint_setp) then ! print*,'rad_config%i_aerosol_type_map =', rad_config%i_aerosol_type_map @@ -433,4 +433,20 @@ & 'photosynthetically active radiation, PAR') + rad_config%i_aerosol_type_map(1:13) = (/ & + & -1, & ! Sea salt, size bin 1 (OPAC) + & -2, & ! Sea salt, size bin 2 (OPAC) + & -3, & ! Sea salt, size bin 3 (OPAC) + & -4, & ! Hydrophilic organic matter (OPAC) + & -5, & ! Ammonium sulphate (OPAC) + & -6, & + & -7, & + & 1, & + & 2, & + & 3, & + & -8, & + & -9, & + & 4 /) ! Stratospheric sulphate (hand edited from OPAC) + rad_config%aerosol_optics_override_file_name = 'aerosol_optics_lmdz.nc' + ! IF (YRERAD%NAERMACC > 0) THEN ! With the MACC aerosol climatology we need to add in the Index: LMDZ6/branches/LMDZ_ECRad/libf/phylmd/ecrad/readaerosol_optic_ecrad.F90 =================================================================== --- LMDZ6/branches/LMDZ_ECRad/libf/phylmd/ecrad/readaerosol_optic_ecrad.F90 (revision 4188) +++ LMDZ6/branches/LMDZ_ECRad/libf/phylmd/ecrad/readaerosol_optic_ecrad.F90 (revision 4188) @@ -0,0 +1,311 @@ +! $Id: readaerosol_optic_ecrad.F90 +! +SUBROUTINE readaerosol_optic_ecrad(debut, aerosol_couple, ok_alw, ok_volcan, & + flag_aerosol, flag_bc_internal_mixture, itap, rjourvrai, & + pdtphys, pplay, paprs, t_seri, rhcl, presnivs, & + tr_seri, mass_solu_aero, mass_solu_aero_pi, & + tau_aero, piz_aero, cg_aero, & + tausum_aero, drytausum_aero, tau3d_aero ) + + ! This routine will : + ! 1) recevie the aerosols(already read and interpolated) corresponding to flag_aerosol + ! 2) calculate the optical properties for the aerosols + ! + + USE dimphy + USE aero_mod + USE phys_local_var_mod, only: sconcso4,sconcno3,sconcoa,sconcbc,sconcss,sconcdust, & + concso4,concno3,concoa,concbc,concss,concdust,loadso4,loadoa,loadbc,loadss,loaddust, & + loadno3,load_tmp1,load_tmp2,load_tmp3,load_tmp4,load_tmp5,load_tmp6,load_tmp7, & + load_tmp8,load_tmp9,load_tmp10,m_allaer + + USE infotrac_phy, ONLY: tracers, nqtot, nbtr + USE YOMCST + + IMPLICIT NONE + + include "clesphys.h" + + ! Input arguments + !**************************************************************************************** + LOGICAL, INTENT(IN) :: debut + LOGICAL, INTENT(IN) :: aerosol_couple + LOGICAL, INTENT(IN) :: ok_alw + LOGICAL, INTENT(IN) :: ok_volcan + INTEGER, INTENT(IN) :: flag_aerosol + LOGICAL, INTENT(IN) :: flag_bc_internal_mixture + INTEGER, INTENT(IN) :: itap + REAL, INTENT(IN) :: rjourvrai + REAL, INTENT(IN) :: pdtphys + REAL, DIMENSION(klon,klev), INTENT(IN) :: pplay + REAL, DIMENSION(klon,klev+1), INTENT(IN) :: paprs + REAL, DIMENSION(klon,klev), INTENT(IN) :: t_seri + REAL, DIMENSION(klon,klev), INTENT(IN) :: rhcl ! humidite relative ciel clair + REAL, DIMENSION(klev), INTENT(IN) :: presnivs + REAL, DIMENSION(klon,klev,nbtr), INTENT(IN) :: tr_seri ! concentration tracer + + ! Output arguments + !**************************************************************************************** + REAL, DIMENSION(klon,klev), INTENT(OUT) :: mass_solu_aero ! Total mass for all soluble aerosols + REAL, DIMENSION(klon,klev), INTENT(OUT) :: mass_solu_aero_pi ! -"- preindustrial values + REAL, DIMENSION(klon,klev,2,NSW), INTENT(OUT) :: tau_aero ! Aerosol optical thickness + REAL, DIMENSION(klon,klev,2,NSW), INTENT(OUT) :: piz_aero ! Single scattering albedo aerosol + REAL, DIMENSION(klon,klev,2,NSW), INTENT(OUT) :: cg_aero ! asymmetry parameter aerosol + REAL, DIMENSION(klon,nwave,naero_tot), INTENT(OUT) :: tausum_aero + REAL, DIMENSION(klon,naero_tot), INTENT(OUT) :: drytausum_aero + REAL, DIMENSION(klon,klev,nwave,naero_tot), INTENT(OUT) :: tau3d_aero + + ! Local variables + !**************************************************************************************** + REAL, DIMENSION(klon) :: aerindex ! POLDER aerosol index + REAL, DIMENSION(klon,klev) :: sulfacc ! SO4 accumulation concentration [ug/m3] + REAL, DIMENSION(klon,klev) :: sulfcoarse ! SO4 coarse concentration [ug/m3] + REAL, DIMENSION(klon,klev) :: bcsol ! BC soluble concentration [ug/m3] + REAL, DIMENSION(klon,klev) :: bcins ! BC insoluble concentration [ug/m3] + REAL, DIMENSION(klon,klev) :: pomsol ! POM soluble concentration [ug/m3] + REAL, DIMENSION(klon,klev) :: pomins ! POM insoluble concentration [ug/m3] + REAL, DIMENSION(klon,klev) :: cidust ! DUST aerosol concentration [ug/m3] + REAL, DIMENSION(klon,klev) :: sscoarse ! SS Coarse concentration [ug/m3] + REAL, DIMENSION(klon,klev) :: sssupco ! SS Super Coarse concentration [ug/m3] + REAL, DIMENSION(klon,klev) :: ssacu ! SS Acumulation concentration [ug/m3] + REAL, DIMENSION(klon,klev) :: nitracc ! nitrate accumulation concentration [ug/m3] + REAL, DIMENSION(klon,klev) :: nitrcoarse ! nitrate coarse concentration [ug/m3] + REAL, DIMENSION(klon,klev) :: nitrinscoarse ! nitrate insoluble coarse concentration [ug/m3] + REAL, DIMENSION(klon,klev) :: sulfacc_pi + REAL, DIMENSION(klon,klev) :: sulfcoarse_pi + REAL, DIMENSION(klon,klev) :: bcsol_pi + REAL, DIMENSION(klon,klev) :: bcins_pi + REAL, DIMENSION(klon,klev) :: pomsol_pi + REAL, DIMENSION(klon,klev) :: pomins_pi + REAL, DIMENSION(klon,klev) :: cidust_pi + REAL, DIMENSION(klon,klev) :: sscoarse_pi + REAL, DIMENSION(klon,klev) :: sssupco_pi + REAL, DIMENSION(klon,klev) :: ssacu_pi + REAL, DIMENSION(klon,klev) :: nitracc_pi + REAL, DIMENSION(klon,klev) :: nitrcoarse_pi + REAL, DIMENSION(klon,klev) :: nitrinscoarse_pi + REAL, DIMENSION(klon,klev) :: pdel, zrho +! REAL, DIMENSION(klon,klev,naero_tot) :: m_allaer + REAL, DIMENSION(klon,klev,naero_tot) :: m_allaer_pi !RAF + + integer :: id_ASBCM, id_ASPOMM, id_ASSO4M, id_ASMSAM, id_CSSO4M, id_CSMSAM, id_SSSSM + integer :: id_CSSSM, id_ASSSM, id_CIDUSTM, id_AIBCM, id_AIPOMM, id_ASNO3M, id_CSNO3M, id_CINO3M + INTEGER :: k, i, iq, itr + + !--air density + zrho(:,:)=pplay(:,:)/t_seri(:,:)/RD !--kg/m3 + + !**************************************************************************************** + ! 1) Get aerosol mass + ! + !**************************************************************************************** + ! + ! + IF (aerosol_couple) THEN !--we get aerosols from tr_seri array from INCA + ! + !--copy fields from INCA tr_seri + !--convert to ug m-3 unit for consistency with offline fields + ! + itr = 0 + DO iq = 1,nqtot + IF(.NOT. tracers(iq)%isInPhysics) CYCLE + itr = itr+1 + SELECT CASE(trim(tracers(iq)%name)) + CASE ("ASBCM"); id_ASBCM = itr + CASE ("ASPOMM"); id_ASPOMM = itr + CASE ("ASSO4M"); id_ASSO4M = itr + CASE ("ASMSAM"); id_ASMSAM = itr + CASE ("CSSO4M"); id_CSSO4M = itr + CASE ("CSMSAM"); id_CSMSAM = itr + CASE ("SSSSM"); id_SSSSM = itr + CASE ("CSSSM"); id_CSSSM = itr + CASE ("ASSSM"); id_ASSSM = itr + CASE ("CIDUSTM");id_CIDUSTM= itr + CASE ("AIBCM"); id_AIBCM = itr + CASE ("AIPOMM"); id_AIPOMM = itr + CASE ("ASNO3M"); id_ASNO3M = itr + CASE ("CSNO3M"); id_CSNO3M = itr + CASE ("CINO3M"); id_CINO3M = itr + END SELECT + END DO + + bcsol(:,:) = tr_seri(:,:,id_ASBCM) *zrho(:,:)*1.e9 ! ASBCM + pomsol(:,:) = tr_seri(:,:,id_ASPOMM) *zrho(:,:)*1.e9 ! ASPOMM + sulfacc(:,:) = (tr_seri(:,:,id_ASSO4M)+tr_seri(:,:,id_ASMSAM))*zrho(:,:)*1.e9 ! ASSO4M (=SO4) + ASMSAM (=MSA) + sulfcoarse(:,:) = (tr_seri(:,:,id_CSSO4M)+tr_seri(:,:,id_CSMSAM))*zrho(:,:)*1.e9 ! CSSO4M (=SO4) + CSMSAM (=MSA) + sssupco(:,:) = tr_seri(:,:,id_SSSSM) *zrho(:,:)*1.e9 ! SSSSM + sscoarse(:,:) = tr_seri(:,:,id_CSSSM) *zrho(:,:)*1.e9 ! CSSSM + ssacu(:,:) = tr_seri(:,:,id_ASSSM) *zrho(:,:)*1.e9 ! ASSSM + cidust(:,:) = tr_seri(:,:,id_CIDUSTM) *zrho(:,:)*1.e9 ! CIDUSTM + bcins(:,:) = tr_seri(:,:,id_AIBCM) *zrho(:,:)*1.e9 ! AIBCM + pomins(:,:) = tr_seri(:,:,id_AIPOMM) *zrho(:,:)*1.e9 ! AIPOMM + nitracc(:,:) = tr_seri(:,:,id_ASNO3M) *zrho(:,:)*1.e9 ! ASNO3M + nitrcoarse(:,:) = tr_seri(:,:,id_CSNO3M) *zrho(:,:)*1.e9 ! CSNO3M + nitrinscoarse(:,:)= tr_seri(:,:,id_CINO3M) *zrho(:,:)*1.e9 ! CINO3M + ! + bcsol_pi(:,:) = 0.0 ! ASBCM pre-ind + pomsol_pi(:,:) = 0.0 ! ASPOMM pre-ind + sulfacc_pi(:,:) = 0.0 ! ASSO4M (=SO4) + ASMSAM (=MSA) pre-ind + sulfcoarse_pi(:,:) = 0.0 ! CSSO4M (=SO4) + CSMSAM (=MSA) pre-ind + sssupco_pi(:,:) = 0.0 ! SSSSM pre-ind + sscoarse_pi(:,:) = 0.0 ! CSSSM pre-ind + ssacu_pi(:,:) = 0.0 ! ASSSM pre-ind + cidust_pi(:,:) = 0.0 ! CIDUSTM pre-ind + bcins_pi(:,:) = 0.0 ! AIBCM pre-ind + pomins_pi(:,:) = 0.0 ! AIPOMM pre-ind + nitracc_pi(:,:) = 0.0 ! ASNO3M pre-ind + nitrcoarse_pi(:,:) = 0.0 ! CSNO3M pre-ind + nitrinscoarse_pi(:,:)= 0.0 ! CINO3M + ! + ELSE !--not aerosol_couple + ! + ! Read and interpolate sulfate + IF ( flag_aerosol .EQ. 1 .OR. flag_aerosol .EQ. 6 .OR. flag_aerosol .EQ. 7 ) THEN + + CALL readaerosol_interp(id_ASSO4M_phy, itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, sulfacc, sulfacc_pi,loadso4) + ELSE + sulfacc(:,:) = 0. ; sulfacc_pi(:,:) = 0. + loadso4=0. + ENDIF + + ! Read and interpolate bcsol and bcins + IF ( flag_aerosol .EQ. 2 .OR. flag_aerosol .EQ. 6 .OR. flag_aerosol .EQ. 7 ) THEN + + ! Get bc aerosol distribution + CALL readaerosol_interp(id_ASBCM_phy, itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, bcsol, bcsol_pi, load_tmp1 ) + CALL readaerosol_interp(id_AIBCM_phy, itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, bcins, bcins_pi, load_tmp2 ) + loadbc(:)=load_tmp1(:)+load_tmp2(:) + ELSE + bcsol(:,:) = 0. ; bcsol_pi(:,:) = 0. + bcins(:,:) = 0. ; bcins_pi(:,:) = 0. + loadbc=0. + ENDIF + + ! Read and interpolate pomsol and pomins + IF ( flag_aerosol .EQ. 3 .OR. flag_aerosol .EQ. 6 .OR. flag_aerosol .EQ. 7 ) THEN + + CALL readaerosol_interp(id_ASPOMM_phy, itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, pomsol, pomsol_pi, load_tmp3) + CALL readaerosol_interp(id_AIPOMM_phy, itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, pomins, pomins_pi, load_tmp4) + loadoa(:)=load_tmp3(:)+load_tmp4(:) + ELSE + pomsol(:,:) = 0. ; pomsol_pi(:,:) = 0. + pomins(:,:) = 0. ; pomins_pi(:,:) = 0. + loadoa=0. + ENDIF + + ! Read and interpolate csssm, ssssm, assssm + IF (flag_aerosol .EQ. 4 .OR. flag_aerosol .EQ. 6 .OR. flag_aerosol .EQ. 7 ) THEN + + CALL readaerosol_interp(id_SSSSM_phy ,itap, pdtphys,rjourvrai, & + debut, pplay, paprs, t_seri, sssupco, sssupco_pi, load_tmp5) + CALL readaerosol_interp(id_CSSSM_phy ,itap, pdtphys,rjourvrai, & + debut, pplay, paprs, t_seri, sscoarse,sscoarse_pi, load_tmp6) + CALL readaerosol_interp(id_ASSSM_phy ,itap, pdtphys,rjourvrai, & + debut, pplay, paprs, t_seri, ssacu, ssacu_pi, load_tmp7) + loadss(:)=load_tmp5(:)+load_tmp6(:)+load_tmp7(:) + ELSE + sscoarse(:,:) = 0. ; sscoarse_pi(:,:) = 0. + ssacu(:,:) = 0. ; ssacu_pi(:,:) = 0. + sssupco(:,:) = 0. ; sssupco_pi = 0. + loadss=0. + ENDIF + + ! Read and interpolate cidustm + IF (flag_aerosol .EQ. 5 .OR. flag_aerosol .EQ. 6 .OR. flag_aerosol .EQ. 7 ) THEN + + CALL readaerosol_interp(id_CIDUSTM_phy, itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, cidust, cidust_pi, loaddust) + + ELSE + cidust(:,:) = 0. ; cidust_pi(:,:) = 0. + loaddust=0. + ENDIF + ! + ! Read and interpolate asno3m, csno3m, cino3m + IF (flag_aerosol .EQ. 6 .OR. flag_aerosol .EQ. 7 ) THEN + + CALL readaerosol_interp(id_ASNO3M_phy, itap, pdtphys, rjourvrai, & + debut, pplay, paprs, t_seri, nitracc, nitracc_pi, load_tmp8) + CALL readaerosol_interp(id_CSNO3M_phy, itap, pdtphys, rjourvrai, & + debut, pplay, paprs, t_seri, nitrcoarse, nitrcoarse_pi, load_tmp9) + CALL readaerosol_interp(id_CINO3M_phy, itap, pdtphys, rjourvrai, & + debut, pplay, paprs, t_seri, nitrinscoarse, nitrinscoarse_pi, load_tmp10) + loadno3(:)=load_tmp8(:)+load_tmp9(:)+load_tmp10(:) + + ELSE + nitracc(:,:) = 0.0 ; nitracc_pi(:,:) = 0.0 + nitrcoarse(:,:) = 0.0 ; nitrcoarse_pi(:,:) = 0.0 + nitrinscoarse(:,:) = 0.0 ; nitrinscoarse_pi(:,:)= 0.0 + loadno3(:)=0.0 + ENDIF + ! + ! CSSO4M is set to 0 as not reliable + sulfcoarse(:,:) = 0.0 ! CSSO4M (=SO4) + CSMSAM (=MSA) + sulfcoarse_pi(:,:) = 0.0 ! CSSO4M (=SO4) + CSMSAM (=MSA) pre-ind + + ENDIF !--not aerosol_couple + + ! + ! Store all aerosols mixing ratios in one variable for radiation scheme (unit kg/kg for ECRAD) + ! present-day values + m_allaer(:,:,id_ASBCM_phy) = bcsol(:,:) /1.e9/zrho(:,:) ! ASBCM + m_allaer(:,:,id_ASPOMM_phy) = pomsol(:,:) /1.e9/zrho(:,:) ! ASPOMM + m_allaer(:,:,id_ASSO4M_phy) = sulfacc(:,:) /1.e9/zrho(:,:) ! ASSO4M (= SO4) + m_allaer(:,:,id_CSSO4M_phy) = sulfcoarse(:,:) /1.e9/zrho(:,:) ! CSSO4M + m_allaer(:,:,id_SSSSM_phy) = sssupco(:,:) /1.e9/zrho(:,:) ! SSSSM + m_allaer(:,:,id_CSSSM_phy) = sscoarse(:,:) /1.e9/zrho(:,:) ! CSSSM + m_allaer(:,:,id_ASSSM_phy) = ssacu(:,:) /1.e9/zrho(:,:) ! ASSSM + m_allaer(:,:,id_CIDUSTM_phy)= cidust(:,:) /1.e9/zrho(:,:) ! CIDUSTM + m_allaer(:,:,id_AIBCM_phy) = bcins(:,:) /1.e9/zrho(:,:) ! AIBCM + m_allaer(:,:,id_ASNO3M_phy) = nitracc(:,:) /1.e9/zrho(:,:) ! ASNO3M + m_allaer(:,:,id_CSNO3M_phy) = nitrcoarse(:,:) /1.e9/zrho(:,:) ! CSNO3M + m_allaer(:,:,id_CINO3M_phy) = nitrinscoarse(:,:)/1.e9/zrho(:,:) ! CINO3M + m_allaer(:,:,id_AIPOMM_phy) = pomins(:,:) /1.e9/zrho(:,:) ! AIPOMM + m_allaer(:,:,id_STRAT_phy) = 0.0 + + ! pre-industrial (pi) values + m_allaer_pi(:,:,id_ASBCM_phy) = bcsol_pi(:,:) /1.e9/zrho(:,:) ! ASBCM pre-ind + m_allaer_pi(:,:,id_ASPOMM_phy) = pomsol_pi(:,:) /1.e9/zrho(:,:) ! ASPOMM pre-ind + m_allaer_pi(:,:,id_ASSO4M_phy) = sulfacc_pi(:,:) /1.e9/zrho(:,:) ! ASSO4M (= SO4) pre-ind + m_allaer_pi(:,:,id_CSSO4M_phy) = sulfcoarse_pi(:,:) /1.e9/zrho(:,:) ! CSSO4M pre-ind + m_allaer_pi(:,:,id_SSSSM_phy) = sssupco_pi(:,:) /1.e9/zrho(:,:) ! SSSSM pre-ind + m_allaer_pi(:,:,id_CSSSM_phy) = sscoarse_pi(:,:) /1.e9/zrho(:,:) ! CSSSM pre-ind + m_allaer_pi(:,:,id_ASSSM_phy) = ssacu_pi(:,:) /1.e9/zrho(:,:) ! ASSSM pre-ind + m_allaer_pi(:,:,id_CIDUSTM_phy)= cidust_pi(:,:) /1.e9/zrho(:,:) ! CIDUSTM pre-ind + m_allaer_pi(:,:,id_AIBCM_phy) = bcins_pi(:,:) /1.e9/zrho(:,:) ! AIBCM pre-ind + m_allaer_pi(:,:,id_ASNO3M_phy) = nitracc_pi(:,:) /1.e9/zrho(:,:) ! ASNO3M pre-ind + m_allaer_pi(:,:,id_CSNO3M_phy) = nitrcoarse_pi(:,:) /1.e9/zrho(:,:) ! CSNO3M pre-ind + m_allaer_pi(:,:,id_CINO3M_phy) = nitrinscoarse_pi(:,:)/1.e9/zrho(:,:) ! CINO3M pre-ind + m_allaer_pi(:,:,id_AIPOMM_phy) = pomins_pi(:,:) /1.e9/zrho(:,:) ! AIPOMM pre-ind + m_allaer_pi(:,:,id_STRAT_phy) = 0.0 + + ! + ! Calculate the total mass of all soluble aersosols (in unit ug /m3) + ! to be revisited for AR6 + mass_solu_aero(:,:) = sulfacc(:,:) + bcsol(:,:) + pomsol(:,:) + nitracc(:,:) + ssacu(:,:) + mass_solu_aero_pi(:,:) = sulfacc_pi(:,:) + bcsol_pi(:,:) + pomsol_pi(:,:) + nitracc_pi(:,:) + ssacu_pi(:,:) + + !**************************************************************************************** + ! 2) Calculate optical properties for the aerosols + ! + !**************************************************************************************** + DO k = 1, klev + DO i = 1, klon + pdel(i,k) = paprs(i,k) - paprs (i,k+1) + ENDDO + ENDDO + + ! Diagnostics calculation for CMIP5 protocol unit (in unit kg/m3) + sconcso4(:) =m_allaer(:,1,id_ASSO4M_phy)*1.e-9 + sconcno3(:) =(m_allaer(:,1,id_ASNO3M_phy)+m_allaer(:,1,id_CSNO3M_phy)+m_allaer(:,1,id_CINO3M_phy))*1.e-9 + sconcoa(:) =(m_allaer(:,1,id_ASPOMM_phy)+m_allaer(:,1,id_AIPOMM_phy))*1.e-9 + sconcbc(:) =(m_allaer(:,1,id_ASBCM_phy)+m_allaer(:,1,id_AIBCM_phy))*1.e-9 + sconcss(:) =(m_allaer(:,1,id_ASSSM_phy)+m_allaer(:,1,id_CSSSM_phy)+m_allaer(:,1,id_SSSSM_phy))*1.e-9 + sconcdust(:) =m_allaer(:,1,id_CIDUSTM_phy)*1.e-9 + concso4(:,:) =m_allaer(:,:,id_ASSO4M_phy)*1.e-9 + concno3(:,:) =(m_allaer(:,:,id_ASNO3M_phy)+m_allaer(:,:,id_CSNO3M_phy)+m_allaer(:,:,id_CINO3M_phy))*1.e-9 + concoa(:,:) =(m_allaer(:,:,id_ASPOMM_phy)+m_allaer(:,:,id_AIPOMM_phy))*1.e-9 + concbc(:,:) =(m_allaer(:,:,id_ASBCM_phy)+m_allaer(:,:,id_AIBCM_phy))*1.e-9 + concss(:,:) =(m_allaer(:,:,id_ASSSM_phy)+m_allaer(:,:,id_CSSSM_phy)+m_allaer(:,:,id_SSSSM_phy))*1.e-9 + concdust(:,:)=m_allaer(:,:,id_CIDUSTM_phy)*1.e-9 + +END SUBROUTINE readaerosol_optic_ecrad Index: LMDZ6/branches/LMDZ_ECRad/libf/phylmd/ecrad/setup_aerosol_optics_lmdz_m.F90 =================================================================== --- LMDZ6/branches/LMDZ_ECRad/libf/phylmd/ecrad/setup_aerosol_optics_lmdz_m.F90 (revision 4188) +++ LMDZ6/branches/LMDZ_ECRad/libf/phylmd/ecrad/setup_aerosol_optics_lmdz_m.F90 (revision 4188) @@ -0,0 +1,126 @@ +module setup_aerosol_optics_lmdz_m + + implicit none + +contains + + subroutine setup_aerosol_optics_lmdz(ao, file_name) + + ! Read aerosol optical properties. Note differences with + ! "radiation_aerosol_optics_data::setup_aerosol_optics": + + ! -- The input NetCDF file is not flat, it contains NetCDF groups. + + ! -- We do not define ao%ssa_mono_phobic, ao%g_mono_phobic, + ! ao%lidar_ratio_mono_phobic, ao%ssa_mono_philic, + ! ao%g_mono_philic, ao%lidar_ratio_mono_philic. They are not in + ! the input NetCDF file and they are not used by ECRad. + + ! -- We do not define ao%description_phobic_str and + ! ao%description_philic_str. We just leave the initialization + ! value, which is a blank. + + ! -- We have to cshift the shortwave fields because the the + ! shortwave bands are in ascending order in the NetCDF file while + ! they are not in ECRad. + + use radiation_aerosol_optics_data, only: aerosol_optics_type, & + IAerosolClassUndefined + use netcdf95, only: nf95_open, nf95_inq_grp_full_ncid, nf95_close, & + nf95_inq_dimid, nf95_inq_varid, nf95_inquire_dimension, & + nf95_get_var, nf95_gw_var + use netcdf, only: nf90_nowrite + + type(aerosol_optics_type), intent(out):: ao + + character(len=*), intent(in):: file_name + ! NetCDF file containing the aerosol optics data + + ! Local: + integer ncid, grpid, dimid, varid + + !----------------------------------------------------------------------- + + ao%use_hydrophilic = .true. + ao%use_monochromatic = .true. + print*,'file_name= ',file_name + call nf95_open(file_name, nf90_nowrite, ncid) + call nf95_inq_grp_full_ncid(ncid, "Hydrophilic", grpid) + call nf95_inq_dimid(grpid, "hur", dimid) + call nf95_inquire_dimension(grpid, dimid, nclen = ao%nrh) + allocate(ao%rh_lower(ao%nrh)) + call nf95_inq_varid(grpid, "hur_bounds", varid) + call nf95_get_var(grpid, varid, ao%rh_lower, count_nc = [1, ao%nrh]) + + ! Hydrophilic/LW_bands: + call nf95_inq_grp_full_ncid(ncid, "Hydrophilic/LW_bands", grpid) + call nf95_inq_varid(grpid, "asymmetry", varid) + call nf95_gw_var(grpid, varid, ao%g_lw_philic) + call nf95_inq_varid(grpid, "single_scat_alb", varid) + call nf95_gw_var(grpid, varid, ao%ssa_lw_philic) + call nf95_inq_varid(grpid, "mass_ext", varid) + call nf95_gw_var(grpid, varid, ao%mass_ext_lw_philic) + + ! Hydrophilic/SW_bands: + call nf95_inq_grp_full_ncid(ncid, "Hydrophilic/SW_bands", grpid) + call nf95_inq_varid(grpid, "asymmetry", varid) + call nf95_gw_var(grpid, varid, ao%g_sw_philic) + ao%g_sw_philic = cshift(ao%g_sw_philic, 1) + call nf95_inq_varid(grpid, "single_scat_alb", varid) + call nf95_gw_var(grpid, varid, ao%ssa_sw_philic) + ao%g_sw_philic = cshift(ao%ssa_sw_philic, 1) + call nf95_inq_varid(grpid, "mass_ext", varid) + call nf95_gw_var(grpid, varid, ao%mass_ext_sw_philic) + ao%g_sw_philic = cshift(ao%mass_ext_sw_philic, 1) + + ! Hydrophilic/Monochromatic: + call nf95_inq_grp_full_ncid(ncid, "Hydrophilic/Monochromatic", grpid) + call nf95_inq_varid(grpid, "mass_ext", varid) + call nf95_gw_var(grpid, varid, ao%mass_ext_mono_philic) + + ! Hydrophobic/LW_bands: + call nf95_inq_grp_full_ncid(ncid, "Hydrophobic/LW_bands", grpid) + call nf95_inq_varid(grpid, "asymmetry", varid) + call nf95_gw_var(grpid, varid, ao%g_lw_phobic) + call nf95_inq_varid(grpid, "single_scat_alb", varid) + call nf95_gw_var(grpid, varid, ao%ssa_lw_phobic) + call nf95_inq_varid(grpid, "mass_ext", varid) + call nf95_gw_var(grpid, varid, ao%mass_ext_lw_phobic) + + ! Hydrophobic/SW_bands: + call nf95_inq_grp_full_ncid(ncid, "Hydrophobic/SW_bands", grpid) + call nf95_inq_varid(grpid, "asymmetry", varid) + call nf95_gw_var(grpid, varid, ao%g_sw_phobic) + ao%g_sw_phobic = cshift(ao%g_sw_phobic, 1) + call nf95_inq_varid(grpid, "single_scat_alb", varid) + call nf95_gw_var(grpid, varid, ao%ssa_sw_phobic) + ao%g_sw_phobic = cshift(ao%ssa_sw_phobic, 1) + call nf95_inq_varid(grpid, "mass_ext", varid) + call nf95_gw_var(grpid, varid, ao%mass_ext_sw_phobic) + ao%g_sw_phobic = cshift(ao%mass_ext_sw_phobic, 1) + + ! Hydrophobic/Monochromatic: + call nf95_inq_grp_full_ncid(ncid, "Hydrophobic/Monochromatic", grpid) + call nf95_inq_varid(grpid, "mass_ext", varid) + call nf95_gw_var(grpid, varid, ao%mass_ext_mono_phobic) + + call nf95_close(ncid) + + ! Get array sizes + ao%n_bands_lw = size(ao%mass_ext_lw_phobic, 1) + ao%n_bands_sw = size(ao%mass_ext_sw_phobic, 1) + ao%n_mono_wl = size(ao%mass_ext_mono_phobic, 1) + ao%n_type_phobic = size(ao%mass_ext_lw_phobic, 2) + ao%n_type_philic = size(ao%mass_ext_lw_philic, 3) + + ! Allocate memory for mapping arrays + ao%ntype = ao%n_type_phobic + ao%n_type_philic + allocate(ao%iclass(ao%ntype)) + allocate(ao%itype(ao%ntype)) + + ao%iclass = IAerosolClassUndefined + ao%itype = 0 + + end subroutine setup_aerosol_optics_lmdz + +end module setup_aerosol_optics_lmdz_m Index: LMDZ6/branches/LMDZ_ECRad/libf/phylmd/phys_local_var_mod.F90 =================================================================== --- LMDZ6/branches/LMDZ_ECRad/libf/phylmd/phys_local_var_mod.F90 (revision 4187) +++ LMDZ6/branches/LMDZ_ECRad/libf/phylmd/phys_local_var_mod.F90 (revision 4188) @@ -24,4 +24,6 @@ REAL, SAVE, ALLOCATABLE :: tr_seri(:,:,:) !$OMP THREADPRIVATE(tr_seri) + REAL, SAVE, ALLOCATABLE :: rhcl(:,:) + !$OMP THREADPRIVATE(rhcl) REAL, SAVE, ALLOCATABLE :: d_t_dyn(:,:), d_q_dyn(:,:) !$OMP THREADPRIVATE(d_t_dyn, d_q_dyn) @@ -116,4 +118,7 @@ !$OMP THREADPRIVATE(d_ts, d_tr) +! aerosols + REAL, SAVE, ALLOCATABLE :: m_allaer (:,:,:) + !$OMP THREADPRIVATE(m_allaer) ! diagnostique pour le rayonnement REAL, SAVE, ALLOCATABLE :: topswad_aero(:), solswad_aero(:) ! diag @@ -608,4 +613,5 @@ l_mix(:,:,:)=0.;l_mixmin(:,:,:)=0.;tke_dissip(:,:,:)=0.;wprime(:,:,:)=0. ! doit etre initialse car pas toujours remplis + ALLOCATE(rhcl(klon,klev)) ALLOCATE(tr_seri(klon,klev,nbtr)) ALLOCATE(d_t_dyn(klon,klev),d_q_dyn(klon,klev)) @@ -647,4 +653,7 @@ ALLOCATE(d_u_lif(klon,klev),d_v_lif(klon,klev)) ALLOCATE(d_ts(klon,nbsrf), d_tr(klon,klev,nbtr)) + +! aerosols + ALLOCATE(m_allaer(klon,klev,naero_tot)) ! Special RRTM ALLOCATE(ZLWFT0_i(klon,klev+1),ZSWFT0_i(klon,klev+1),ZFLDN0(klon,klev+1)) @@ -937,4 +946,5 @@ DEALLOCATE(l_mixmin,l_mix, tke_dissip,wprime) + DEALLOCATE(rhcl) DEALLOCATE(tr_seri) DEALLOCATE(d_t_dyn,d_q_dyn) @@ -1046,5 +1056,6 @@ DEALLOCATE(solsw_aerop, solsw0_aerop) DEALLOCATE(topswcf_aerop, solswcf_aerop) - +!AI Aerosols + DEALLOCATE(m_allaer) !CK LW diagnostics DEALLOCATE(toplwad_aerop, sollwad_aerop) Index: LMDZ6/branches/LMDZ_ECRad/libf/phylmd/phys_output_write_mod.F90 =================================================================== --- LMDZ6/branches/LMDZ_ECRad/libf/phylmd/phys_output_write_mod.F90 (revision 4187) +++ LMDZ6/branches/LMDZ_ECRad/libf/phylmd/phys_output_write_mod.F90 (revision 4188) @@ -1568,4 +1568,5 @@ !--OLIVIER !This is warranted by treating INCA aerosols as offline aerosols +#ifndef CPP_ECRAD IF (flag_aerosol.GT.0) THEN IF (ALL(types_trac/='inca') .OR. config_inca=='aeNP') THEN @@ -1773,4 +1774,5 @@ ENDIF ENDIF + CALL histwrite_phy(o_lwcon, flwc) CALL histwrite_phy(o_iwcon, fiwc) @@ -1807,4 +1809,5 @@ #endif !solbnd end +#endif #endif Index: LMDZ6/branches/LMDZ_ECRad/libf/phylmd/physiq_mod.F90 =================================================================== --- LMDZ6/branches/LMDZ_ECRad/libf/phylmd/physiq_mod.F90 (revision 4187) +++ LMDZ6/branches/LMDZ_ECRad/libf/phylmd/physiq_mod.F90 (revision 4188) @@ -132,4 +132,5 @@ ! Variables locales pour effectuer les appels en serie t_seri,q_seri,ql_seri,qs_seri,u_seri,v_seri,tr_seri,rneb_seri, & + rhcl, & ! Dynamic tendencies (diagnostics) d_t_dyn,d_q_dyn,d_ql_dyn,d_qs_dyn,d_u_dyn,d_v_dyn,d_tr_dyn,d_rneb_dyn, & @@ -813,5 +814,5 @@ ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ! - REAL rhcl(klon,klev) ! humiditi relative ciel clair +! REAL rhcl(klon,klev) ! humiditi relative ciel clair REAL dialiq(klon,klev) ! eau liquide nuageuse REAL diafra(klon,klev) ! fraction nuageuse @@ -3969,5 +3970,5 @@ tausum_aero, tau3d_aero) ENDIF - ELSE ! RRTM radiation + ELSE IF (iflag_rrtm .EQ.1) THEN ! RRTM radiation IF (aerosol_couple .AND. config_inca == 'aero' ) THEN abort_message='config_inca=aero et rrtm=1 impossible' @@ -4035,5 +4036,19 @@ ! ENDIF + ELSE IF (iflag_rrtm .EQ.2) THEN ! ecrad RADIATION +#ifdef CPP_ECRAD + !--climatologies or INCA aerosols + CALL readaerosol_optic_ecrad( debut, aerosol_couple, ok_alw, ok_volcan, & + flag_aerosol, flag_bc_internal_mixture, itap, jD_cur-jD_ref, & + pdtphys, pplay, paprs, t_seri, rhcl, presnivs, & + tr_seri, mass_solu_aero, mass_solu_aero_pi, & + tau_aero_sw_rrtm, piz_aero_sw_rrtm, cg_aero_sw_rrtm, & + tausum_aero, drytausum_aero, tau3d_aero) +#else + abort_message='You should compile with -rad ecrad if running with iflag_rrtm=2' + CALL abort_physic(modname,abort_message,1) +#endif ENDIF + ELSE !--flag_aerosol = 0 tausum_aero(:,:,:) = 0. Index: LMDZ6/branches/LMDZ_ECRad/libf/phylmd/radlwsw_m.F90 =================================================================== --- LMDZ6/branches/LMDZ_ECRad/libf/phylmd/radlwsw_m.F90 (revision 4187) +++ LMDZ6/branches/LMDZ_ECRad/libf/phylmd/radlwsw_m.F90 (revision 4188) @@ -78,4 +78,5 @@ ! Besoin pour ECRAD de pctsrf, zmasq, longitude, altitude #ifdef CPP_ECRAD + USE phys_local_var_mod, ONLY: rhcl, m_allaer USE geometry_mod, ONLY: latitude, longitude USE phys_state_var_mod, ONLY: pctsrf @@ -323,5 +324,5 @@ ZQ_SNOW(klon,klev) ! Snow cloud mass mixing ratio (kg/kg) ? REAL(KIND=8) ZAEROSOL_OLD(KLON,6,KLEV), & ! - ZAEROSOL(KLON,KLEV,naero_tot) ! + ZAEROSOL(KLON,KLEV,naero_grp) ! ! OUTPUTS REAL(KIND=8) ZFLUX_DIR(klon), & ! Direct compt of surf flux into horizontal plane @@ -1181,7 +1182,7 @@ ! ! AI ATTENTION Aerosols A REVOIR -! DO i = 1, kdlon -! DO k = 1, kflev -! DO kk= 1, naero_tot + DO i = 1, kdlon + DO k = 1, kflev + DO kk= 1, naero_grp ! DO kk=1, NSW ! @@ -1194,8 +1195,9 @@ ! PCGA_NAT(i,kflev+1-k,kk)=cg_aero_sw_rrtm(i,k,1,kk) ! ZAEROSOL(i,kflev+1-k,kk)=m_allaer(i,k,kk) -! -! ENDDO -! ENDDO -! ENDDO + ZAEROSOL(i,kflev+1-k,kk)=m_allaer(i,k,kk) +! + ENDDO + ENDDO + ENDDO !-end OB ! @@ -1345,5 +1347,5 @@ CALL RADIATION_SCHEME & - & (ist, iend, klon, klev, naero_tot, NSW, & + & (ist, iend, klon, klev, naero_grp, NSW, & ! ??? naero_tot & day_cur, current_time, & Index: LMDZ6/branches/LMDZ_ECRad/libf/phylmd/setup_aerosol_optics_lmdz_m.F90 =================================================================== --- LMDZ6/branches/LMDZ_ECRad/libf/phylmd/setup_aerosol_optics_lmdz_m.F90 (revision 4187) +++ (revision ) @@ -1,125 +1,0 @@ -module setup_aerosol_optics_lmdz_m - - implicit none - -contains - - subroutine setup_aerosol_optics_lmdz(ao, file_name) - - ! Read aerosol optical properties. Note differences with - ! "radiation_aerosol_optics_data::setup_aerosol_optics": - - ! -- The input NetCDF file is not flat, it contains NetCDF groups. - - ! -- We do not define ao%ssa_mono_phobic, ao%g_mono_phobic, - ! ao%lidar_ratio_mono_phobic, ao%ssa_mono_philic, - ! ao%g_mono_philic, ao%lidar_ratio_mono_philic. They are not in - ! the input NetCDF file and they are not used by ECRad. - - ! -- We do not define ao%description_phobic_str and - ! ao%description_philic_str. We just leave the initialization - ! value, which is a blank. - - ! -- We have to cshift the shortwave fields because the the - ! shortwave bands are in ascending order in the NetCDF file while - ! they are not in ECRad. - - use radiation_aerosol_optics_data, only: aerosol_optics_type, & - IAerosolClassUndefined - use netcdf95, only: nf95_open, nf95_inq_grp_full_ncid, nf95_close, & - nf95_inq_dimid, nf95_inq_varid, nf95_inquire_dimension, & - nf95_get_var, nf95_gw_var - use netcdf, only: nf90_nowrite - - type(aerosol_optics_type), intent(out):: ao - - character(len=*), intent(in):: file_name - ! NetCDF file containing the aerosol optics data - - ! Local: - integer ncid, grpid, dimid, varid - - !----------------------------------------------------------------------- - - ao%use_hydrophilic = .true. - ao%use_monochromatic = .true. - call nf95_open(file_name, nf90_nowrite, ncid) - call nf95_inq_grp_full_ncid(ncid, "Hydrophilic", grpid) - call nf95_inq_dimid(grpid, "hur", dimid) - call nf95_inquire_dimension(grpid, dimid, nclen = ao%nrh) - allocate(ao%rh_lower(ao%nrh)) - call nf95_inq_varid(grpid, "hur_bounds", varid) - call nf95_get_var(grpid, varid, ao%rh_lower, count_nc = [1, ao%nrh]) - - ! Hydrophilic/LW_bands: - call nf95_inq_grp_full_ncid(ncid, "Hydrophilic/LW_bands", grpid) - call nf95_inq_varid(grpid, "asymmetry", varid) - call nf95_gw_var(grpid, varid, ao%g_lw_philic) - call nf95_inq_varid(grpid, "single_scat_alb", varid) - call nf95_gw_var(grpid, varid, ao%ssa_lw_philic) - call nf95_inq_varid(grpid, "mass_ext", varid) - call nf95_gw_var(grpid, varid, ao%mass_ext_lw_philic) - - ! Hydrophilic/SW_bands: - call nf95_inq_grp_full_ncid(ncid, "Hydrophilic/SW_bands", grpid) - call nf95_inq_varid(grpid, "asymmetry", varid) - call nf95_gw_var(grpid, varid, ao%g_sw_philic) - ao%g_sw_philic = cshift(ao%g_sw_philic, 1) - call nf95_inq_varid(grpid, "single_scat_alb", varid) - call nf95_gw_var(grpid, varid, ao%ssa_sw_philic) - ao%g_sw_philic = cshift(ao%ssa_sw_philic, 1) - call nf95_inq_varid(grpid, "mass_ext", varid) - call nf95_gw_var(grpid, varid, ao%mass_ext_sw_philic) - ao%g_sw_philic = cshift(ao%mass_ext_sw_philic, 1) - - ! Hydrophilic/Monochromatic: - call nf95_inq_grp_full_ncid(ncid, "Hydrophilic/Monochromatic", grpid) - call nf95_inq_varid(grpid, "mass_ext", varid) - call nf95_gw_var(grpid, varid, ao%mass_ext_mono_philic) - - ! Hydrophobic/LW_bands: - call nf95_inq_grp_full_ncid(ncid, "Hydrophobic/LW_bands", grpid) - call nf95_inq_varid(grpid, "asymmetry", varid) - call nf95_gw_var(grpid, varid, ao%g_lw_phobic) - call nf95_inq_varid(grpid, "single_scat_alb", varid) - call nf95_gw_var(grpid, varid, ao%ssa_lw_phobic) - call nf95_inq_varid(grpid, "mass_ext", varid) - call nf95_gw_var(grpid, varid, ao%mass_ext_lw_phobic) - - ! Hydrophobic/SW_bands: - call nf95_inq_grp_full_ncid(ncid, "Hydrophobic/SW_bands", grpid) - call nf95_inq_varid(grpid, "asymmetry", varid) - call nf95_gw_var(grpid, varid, ao%g_sw_phobic) - ao%g_sw_phobic = cshift(ao%g_sw_phobic, 1) - call nf95_inq_varid(grpid, "single_scat_alb", varid) - call nf95_gw_var(grpid, varid, ao%ssa_sw_phobic) - ao%g_sw_phobic = cshift(ao%ssa_sw_phobic, 1) - call nf95_inq_varid(grpid, "mass_ext", varid) - call nf95_gw_var(grpid, varid, ao%mass_ext_sw_phobic) - ao%g_sw_phobic = cshift(ao%mass_ext_sw_phobic, 1) - - ! Hydrophobic/Monochromatic: - call nf95_inq_grp_full_ncid(ncid, "Hydrophobic/Monochromatic", grpid) - call nf95_inq_varid(grpid, "mass_ext", varid) - call nf95_gw_var(grpid, varid, ao%mass_ext_mono_phobic) - - call nf95_close(ncid) - - ! Get array sizes - ao%n_bands_lw = size(ao%mass_ext_lw_phobic, 1) - ao%n_bands_sw = size(ao%mass_ext_sw_phobic, 1) - ao%n_mono_wl = size(ao%mass_ext_mono_phobic, 1) - ao%n_type_phobic = size(ao%mass_ext_lw_phobic, 2) - ao%n_type_philic = size(ao%mass_ext_lw_philic, 3) - - ! Allocate memory for mapping arrays - ao%ntype = ao%n_type_phobic + ao%n_type_philic - allocate(ao%iclass(ao%ntype)) - allocate(ao%itype(ao%ntype)) - - ao%iclass = IAerosolClassUndefined - ao%itype = 0 - - end subroutine setup_aerosol_optics_lmdz - -end module setup_aerosol_optics_lmdz_m