Context Navigation

← Previous Revision
Next Revision →
Normal
Revision Log

readaerosol_optic_ecrad.F90

Last change on this file was 5202, checked in by Laurent Fairhead, 3 months ago
Updating cirrus branch to trunk revision 5171
File size: 16.2 KB

Rev	Line
[4853]	1	! $Id: readaerosol_optic_ecrad.F90
	2	!
	3	SUBROUTINE readaerosol_optic_ecrad(debut, aerosol_couple, ok_alw, ok_volcan, &
	4	flag_aerosol, flag_bc_internal_mixture, itap, rjourvrai, &
	5	pdtphys, pplay, paprs, t_seri, rhcl, presnivs, &
[5202]	6	tr_seri, mass_solu_aero, mass_solu_aero_pi, m_allaer)
[4853]	7	! tau_aero, piz_aero, cg_aero, &
	8	! tausum_aero, drytausum_aero, tau3d_aero )
	9
	10	! This routine will :
	11	! 1) recevie the aerosols(already read and interpolated) corresponding to flag_aerosol
	12	! 2) calculate the optical properties for the aerosols
	13	!
	14
	15	USE dimphy
	16	USE aero_mod
	17	USE phys_local_var_mod, only: sconcso4,sconcno3,sconcoa,sconcbc,sconcss,sconcdust, &
	18	concso4,concno3,concoa,concbc,concss,concdust,loadso4,loadoa,loadbc,loadss,loaddust, &
	19	loadno3,load_tmp1,load_tmp2,load_tmp3,load_tmp4,load_tmp5,load_tmp6,load_tmp7, &
[5202]	20	load_tmp8,load_tmp9,load_tmp10
[4853]	21
	22	USE infotrac_phy, ONLY: tracers, nqtot, nbtr
	23	USE YOMCST
	24
	25	IMPLICIT NONE
	26
	27	include "clesphys.h"
	28
	29	! Input arguments
	30	!****************************************************************************************
	31	LOGICAL, INTENT(IN) :: debut
	32	LOGICAL, INTENT(IN) :: aerosol_couple
	33	LOGICAL, INTENT(IN) :: ok_alw
	34	LOGICAL, INTENT(IN) :: ok_volcan
	35	INTEGER, INTENT(IN) :: flag_aerosol
	36	LOGICAL, INTENT(IN) :: flag_bc_internal_mixture
	37	INTEGER, INTENT(IN) :: itap
	38	REAL, INTENT(IN) :: rjourvrai
	39	REAL, INTENT(IN) :: pdtphys
	40	REAL, DIMENSION(klon,klev), INTENT(IN) :: pplay
	41	REAL, DIMENSION(klon,klev+1), INTENT(IN) :: paprs
	42	REAL, DIMENSION(klon,klev), INTENT(IN) :: t_seri
	43	REAL, DIMENSION(klon,klev), INTENT(IN) :: rhcl ! humidite relative ciel clair
	44	REAL, DIMENSION(klev), INTENT(IN) :: presnivs
	45	REAL, DIMENSION(klon,klev,nbtr), INTENT(IN) :: tr_seri ! concentration tracer
	46
	47	! Output arguments
	48	!****************************************************************************************
	49	REAL, DIMENSION(klon,klev), INTENT(OUT) :: mass_solu_aero ! Total mass for all soluble aerosols
	50	REAL, DIMENSION(klon,klev), INTENT(OUT) :: mass_solu_aero_pi ! -"- preindustrial values
[5202]	51	REAL, DIMENSION(klon,klev,naero_tot), INTENT(OUT) :: m_allaer
	52	! AI a passer par la suite en argument si besoin pour ecrad
	53	!REAL, DIMENSION(klon,klev,naero_tot), INTENT(OUT) :: m_allaer_pi !RAF
	54
[4853]	55	! REAL, DIMENSION(klon,klev,2,NSW), INTENT(OUT) :: tau_aero ! Aerosol optical thickness
	56	! REAL, DIMENSION(klon,klev,2,NSW), INTENT(OUT) :: piz_aero ! Single scattering albedo aerosol
	57	! REAL, DIMENSION(klon,klev,2,NSW), INTENT(OUT) :: cg_aero ! asymmetry parameter aerosol
	58	! REAL, DIMENSION(klon,nwave,naero_tot), INTENT(OUT) :: tausum_aero
	59	! REAL, DIMENSION(klon,naero_tot), INTENT(OUT) :: drytausum_aero
	60	! REAL, DIMENSION(klon,klev,nwave,naero_tot), INTENT(OUT) :: tau3d_aero
	61
	62	! Local variables
	63	!****************************************************************************************
	64	REAL, DIMENSION(klon) :: aerindex ! POLDER aerosol index
	65	REAL, DIMENSION(klon,klev) :: sulfacc ! SO4 accumulation concentration [ug/m3]
	66	REAL, DIMENSION(klon,klev) :: sulfcoarse ! SO4 coarse concentration [ug/m3]
	67	REAL, DIMENSION(klon,klev) :: bcsol ! BC soluble concentration [ug/m3]
	68	REAL, DIMENSION(klon,klev) :: bcins ! BC insoluble concentration [ug/m3]
	69	REAL, DIMENSION(klon,klev) :: pomsol ! POM soluble concentration [ug/m3]
	70	REAL, DIMENSION(klon,klev) :: pomins ! POM insoluble concentration [ug/m3]
	71	REAL, DIMENSION(klon,klev) :: cidust ! DUST aerosol concentration [ug/m3]
	72	REAL, DIMENSION(klon,klev) :: sscoarse ! SS Coarse concentration [ug/m3]
	73	REAL, DIMENSION(klon,klev) :: sssupco ! SS Super Coarse concentration [ug/m3]
	74	REAL, DIMENSION(klon,klev) :: ssacu ! SS Acumulation concentration [ug/m3]
	75	REAL, DIMENSION(klon,klev) :: nitracc ! nitrate accumulation concentration [ug/m3]
	76	REAL, DIMENSION(klon,klev) :: nitrcoarse ! nitrate coarse concentration [ug/m3]
	77	REAL, DIMENSION(klon,klev) :: nitrinscoarse ! nitrate insoluble coarse concentration [ug/m3]
	78	REAL, DIMENSION(klon,klev) :: sulfacc_pi
	79	REAL, DIMENSION(klon,klev) :: sulfcoarse_pi
	80	REAL, DIMENSION(klon,klev) :: bcsol_pi
	81	REAL, DIMENSION(klon,klev) :: bcins_pi
	82	REAL, DIMENSION(klon,klev) :: pomsol_pi
	83	REAL, DIMENSION(klon,klev) :: pomins_pi
	84	REAL, DIMENSION(klon,klev) :: cidust_pi
	85	REAL, DIMENSION(klon,klev) :: sscoarse_pi
	86	REAL, DIMENSION(klon,klev) :: sssupco_pi
	87	REAL, DIMENSION(klon,klev) :: ssacu_pi
	88	REAL, DIMENSION(klon,klev) :: nitracc_pi
	89	REAL, DIMENSION(klon,klev) :: nitrcoarse_pi
	90	REAL, DIMENSION(klon,klev) :: nitrinscoarse_pi
	91	REAL, DIMENSION(klon,klev) :: pdel, zrho
[5202]	92	REAL, DIMENSION(klon,klev,naero_tot) :: m_allaer_pi
[4853]	93
	94	integer :: id_ASBCM, id_ASPOMM, id_ASSO4M, id_ASMSAM, id_CSSO4M, id_CSMSAM, id_SSSSM
	95	integer :: id_CSSSM, id_ASSSM, id_CIDUSTM, id_AIBCM, id_AIPOMM, id_ASNO3M, id_CSNO3M, id_CINO3M
	96	INTEGER :: k, i, iq, itr
	97
	98	!--air density
	99	zrho(:,:)=pplay(:,:)/t_seri(:,:)/RD !--kg/m3
	100
	101	!****************************************************************************************
	102	! 1) Get aerosol mass
	103	!
	104	!****************************************************************************************
	105	!
	106	!
	107	IF (aerosol_couple) THEN !--we get aerosols from tr_seri array from INCA
	108	!
	109	!--copy fields from INCA tr_seri
	110	!--convert to ug m-3 unit for consistency with offline fields
	111	!
	112	itr = 0
	113	DO iq = 1,nqtot
	114	IF(.NOT. tracers(iq)%isInPhysics) CYCLE
	115	itr = itr+1
	116	SELECT CASE(trim(tracers(iq)%name))
	117	CASE ("ASBCM"); id_ASBCM = itr
	118	CASE ("ASPOMM"); id_ASPOMM = itr
	119	CASE ("ASSO4M"); id_ASSO4M = itr
	120	CASE ("ASMSAM"); id_ASMSAM = itr
	121	CASE ("CSSO4M"); id_CSSO4M = itr
	122	CASE ("CSMSAM"); id_CSMSAM = itr
	123	CASE ("SSSSM"); id_SSSSM = itr
	124	CASE ("CSSSM"); id_CSSSM = itr
	125	CASE ("ASSSM"); id_ASSSM = itr
	126	CASE ("CIDUSTM");id_CIDUSTM= itr
	127	CASE ("AIBCM"); id_AIBCM = itr
	128	CASE ("AIPOMM"); id_AIPOMM = itr
	129	CASE ("ASNO3M"); id_ASNO3M = itr
	130	CASE ("CSNO3M"); id_CSNO3M = itr
	131	CASE ("CINO3M"); id_CINO3M = itr
	132	END SELECT
	133	END DO
	134
	135	bcsol(:,:) = tr_seri(:,:,id_ASBCM) zrho(:,:)1.e9 ! ASBCM
	136	pomsol(:,:) = tr_seri(:,:,id_ASPOMM) zrho(:,:)1.e9 ! ASPOMM
	137	sulfacc(:,:) = (tr_seri(:,:,id_ASSO4M)+tr_seri(:,:,id_ASMSAM))zrho(:,:)1.e9 ! ASSO4M (=SO4) + ASMSAM (=MSA)
	138	sulfcoarse(:,:) = (tr_seri(:,:,id_CSSO4M)+tr_seri(:,:,id_CSMSAM))zrho(:,:)1.e9 ! CSSO4M (=SO4) + CSMSAM (=MSA)
	139	sssupco(:,:) = tr_seri(:,:,id_SSSSM) zrho(:,:)1.e9 ! SSSSM
	140	sscoarse(:,:) = tr_seri(:,:,id_CSSSM) zrho(:,:)1.e9 ! CSSSM
	141	ssacu(:,:) = tr_seri(:,:,id_ASSSM) zrho(:,:)1.e9 ! ASSSM
	142	cidust(:,:) = tr_seri(:,:,id_CIDUSTM) zrho(:,:)1.e9 ! CIDUSTM
	143	bcins(:,:) = tr_seri(:,:,id_AIBCM) zrho(:,:)1.e9 ! AIBCM
	144	pomins(:,:) = tr_seri(:,:,id_AIPOMM) zrho(:,:)1.e9 ! AIPOMM
	145	nitracc(:,:) = tr_seri(:,:,id_ASNO3M) zrho(:,:)1.e9 ! ASNO3M
	146	nitrcoarse(:,:) = tr_seri(:,:,id_CSNO3M) zrho(:,:)1.e9 ! CSNO3M
	147	nitrinscoarse(:,:)= tr_seri(:,:,id_CINO3M) zrho(:,:)1.e9 ! CINO3M
	148	!
	149	bcsol_pi(:,:) = 0.0 ! ASBCM pre-ind
	150	pomsol_pi(:,:) = 0.0 ! ASPOMM pre-ind
	151	sulfacc_pi(:,:) = 0.0 ! ASSO4M (=SO4) + ASMSAM (=MSA) pre-ind
	152	sulfcoarse_pi(:,:) = 0.0 ! CSSO4M (=SO4) + CSMSAM (=MSA) pre-ind
	153	sssupco_pi(:,:) = 0.0 ! SSSSM pre-ind
	154	sscoarse_pi(:,:) = 0.0 ! CSSSM pre-ind
	155	ssacu_pi(:,:) = 0.0 ! ASSSM pre-ind
	156	cidust_pi(:,:) = 0.0 ! CIDUSTM pre-ind
	157	bcins_pi(:,:) = 0.0 ! AIBCM pre-ind
	158	pomins_pi(:,:) = 0.0 ! AIPOMM pre-ind
	159	nitracc_pi(:,:) = 0.0 ! ASNO3M pre-ind
	160	nitrcoarse_pi(:,:) = 0.0 ! CSNO3M pre-ind
	161	nitrinscoarse_pi(:,:)= 0.0 ! CINO3M
	162	!
	163	ELSE !--not aerosol_couple
	164	!
	165	! Read and interpolate sulfate
	166	IF ( flag_aerosol .EQ. 1 .OR. flag_aerosol .EQ. 6 .OR. flag_aerosol .EQ. 7 ) THEN
	167
	168	CALL readaerosol_interp(id_ASSO4M_phy, itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, sulfacc, sulfacc_pi,loadso4)
	169	ELSE
	170	sulfacc(:,:) = 0. ; sulfacc_pi(:,:) = 0.
	171	loadso4=0.
	172	ENDIF
	173
	174	! Read and interpolate bcsol and bcins
	175	IF ( flag_aerosol .EQ. 2 .OR. flag_aerosol .EQ. 6 .OR. flag_aerosol .EQ. 7 ) THEN
	176
	177	! Get bc aerosol distribution
	178	CALL readaerosol_interp(id_ASBCM_phy, itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, bcsol, bcsol_pi, load_tmp1 )
	179	CALL readaerosol_interp(id_AIBCM_phy, itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, bcins, bcins_pi, load_tmp2 )
	180	loadbc(:)=load_tmp1(:)+load_tmp2(:)
	181	ELSE
	182	bcsol(:,:) = 0. ; bcsol_pi(:,:) = 0.
	183	bcins(:,:) = 0. ; bcins_pi(:,:) = 0.
	184	loadbc=0.
	185	ENDIF
	186
	187	! Read and interpolate pomsol and pomins
	188	IF ( flag_aerosol .EQ. 3 .OR. flag_aerosol .EQ. 6 .OR. flag_aerosol .EQ. 7 ) THEN
	189
	190	CALL readaerosol_interp(id_ASPOMM_phy, itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, pomsol, pomsol_pi, load_tmp3)
	191	CALL readaerosol_interp(id_AIPOMM_phy, itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, pomins, pomins_pi, load_tmp4)
	192	loadoa(:)=load_tmp3(:)+load_tmp4(:)
	193	ELSE
	194	pomsol(:,:) = 0. ; pomsol_pi(:,:) = 0.
	195	pomins(:,:) = 0. ; pomins_pi(:,:) = 0.
	196	loadoa=0.
	197	ENDIF
	198
	199	! Read and interpolate csssm, ssssm, assssm
	200	IF (flag_aerosol .EQ. 4 .OR. flag_aerosol .EQ. 6 .OR. flag_aerosol .EQ. 7 ) THEN
	201
	202	CALL readaerosol_interp(id_SSSSM_phy ,itap, pdtphys,rjourvrai, &
	203	debut, pplay, paprs, t_seri, sssupco, sssupco_pi, load_tmp5)
	204	CALL readaerosol_interp(id_CSSSM_phy ,itap, pdtphys,rjourvrai, &
	205	debut, pplay, paprs, t_seri, sscoarse,sscoarse_pi, load_tmp6)
	206	CALL readaerosol_interp(id_ASSSM_phy ,itap, pdtphys,rjourvrai, &
	207	debut, pplay, paprs, t_seri, ssacu, ssacu_pi, load_tmp7)
	208	loadss(:)=load_tmp5(:)+load_tmp6(:)+load_tmp7(:)
	209	ELSE
	210	sscoarse(:,:) = 0. ; sscoarse_pi(:,:) = 0.
	211	ssacu(:,:) = 0. ; ssacu_pi(:,:) = 0.
	212	sssupco(:,:) = 0. ; sssupco_pi = 0.
	213	loadss=0.
	214	ENDIF
	215
	216	! Read and interpolate cidustm
	217	IF (flag_aerosol .EQ. 5 .OR. flag_aerosol .EQ. 6 .OR. flag_aerosol .EQ. 7 ) THEN
	218
	219	CALL readaerosol_interp(id_CIDUSTM_phy, itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, cidust, cidust_pi, loaddust)
	220
	221	ELSE
	222	cidust(:,:) = 0. ; cidust_pi(:,:) = 0.
	223	loaddust=0.
	224	ENDIF
	225	!
	226	! Read and interpolate asno3m, csno3m, cino3m
	227	IF (flag_aerosol .EQ. 6 .OR. flag_aerosol .EQ. 7 ) THEN
	228
	229	CALL readaerosol_interp(id_ASNO3M_phy, itap, pdtphys, rjourvrai, &
	230	debut, pplay, paprs, t_seri, nitracc, nitracc_pi, load_tmp8)
	231	CALL readaerosol_interp(id_CSNO3M_phy, itap, pdtphys, rjourvrai, &
	232	debut, pplay, paprs, t_seri, nitrcoarse, nitrcoarse_pi, load_tmp9)
	233	CALL readaerosol_interp(id_CINO3M_phy, itap, pdtphys, rjourvrai, &
	234	debut, pplay, paprs, t_seri, nitrinscoarse, nitrinscoarse_pi, load_tmp10)
	235	loadno3(:)=load_tmp8(:)+load_tmp9(:)+load_tmp10(:)
	236
	237	ELSE
	238	nitracc(:,:) = 0.0 ; nitracc_pi(:,:) = 0.0
	239	nitrcoarse(:,:) = 0.0 ; nitrcoarse_pi(:,:) = 0.0
	240	nitrinscoarse(:,:) = 0.0 ; nitrinscoarse_pi(:,:)= 0.0
	241	loadno3(:)=0.0
	242	ENDIF
	243	!
	244	! CSSO4M is set to 0 as not reliable
	245	sulfcoarse(:,:) = 0.0 ! CSSO4M (=SO4) + CSMSAM (=MSA)
	246	sulfcoarse_pi(:,:) = 0.0 ! CSSO4M (=SO4) + CSMSAM (=MSA) pre-ind
	247
	248	ENDIF !--not aerosol_couple
	249
	250	!
	251	! Store all aerosols mixing ratios in one variable for radiation scheme (unit kg/kg for ECRAD)
	252	! present-day values
	253	m_allaer(:,:,id_ASBCM_phy) = bcsol(:,:) /1.e9/zrho(:,:) ! ASBCM
	254	m_allaer(:,:,id_ASPOMM_phy) = pomsol(:,:) /1.e9/zrho(:,:) ! ASPOMM
	255	m_allaer(:,:,id_ASSO4M_phy) = sulfacc(:,:) /1.e9/zrho(:,:) ! ASSO4M (= SO4)
	256	m_allaer(:,:,id_CSSO4M_phy) = sulfcoarse(:,:) /1.e9/zrho(:,:) ! CSSO4M
	257	m_allaer(:,:,id_SSSSM_phy) = sssupco(:,:) /1.e9/zrho(:,:) ! SSSSM
	258	m_allaer(:,:,id_CSSSM_phy) = sscoarse(:,:) /1.e9/zrho(:,:) ! CSSSM
	259	m_allaer(:,:,id_ASSSM_phy) = ssacu(:,:) /1.e9/zrho(:,:) ! ASSSM
	260	m_allaer(:,:,id_CIDUSTM_phy)= cidust(:,:) /1.e9/zrho(:,:) ! CIDUSTM
	261	m_allaer(:,:,id_AIBCM_phy) = bcins(:,:) /1.e9/zrho(:,:) ! AIBCM
	262	m_allaer(:,:,id_ASNO3M_phy) = nitracc(:,:) /1.e9/zrho(:,:) ! ASNO3M
	263	m_allaer(:,:,id_CSNO3M_phy) = nitrcoarse(:,:) /1.e9/zrho(:,:) ! CSNO3M
	264	m_allaer(:,:,id_CINO3M_phy) = nitrinscoarse(:,:)/1.e9/zrho(:,:) ! CINO3M
	265	m_allaer(:,:,id_AIPOMM_phy) = pomins(:,:) /1.e9/zrho(:,:) ! AIPOMM
	266	m_allaer(:,:,id_STRAT_phy) = 0.0
	267
	268	! pre-industrial (pi) values
	269	m_allaer_pi(:,:,id_ASBCM_phy) = bcsol_pi(:,:) /1.e9/zrho(:,:) ! ASBCM pre-ind
	270	m_allaer_pi(:,:,id_ASPOMM_phy) = pomsol_pi(:,:) /1.e9/zrho(:,:) ! ASPOMM pre-ind
	271	m_allaer_pi(:,:,id_ASSO4M_phy) = sulfacc_pi(:,:) /1.e9/zrho(:,:) ! ASSO4M (= SO4) pre-ind
	272	m_allaer_pi(:,:,id_CSSO4M_phy) = sulfcoarse_pi(:,:) /1.e9/zrho(:,:) ! CSSO4M pre-ind
	273	m_allaer_pi(:,:,id_SSSSM_phy) = sssupco_pi(:,:) /1.e9/zrho(:,:) ! SSSSM pre-ind
	274	m_allaer_pi(:,:,id_CSSSM_phy) = sscoarse_pi(:,:) /1.e9/zrho(:,:) ! CSSSM pre-ind
	275	m_allaer_pi(:,:,id_ASSSM_phy) = ssacu_pi(:,:) /1.e9/zrho(:,:) ! ASSSM pre-ind
	276	m_allaer_pi(:,:,id_CIDUSTM_phy)= cidust_pi(:,:) /1.e9/zrho(:,:) ! CIDUSTM pre-ind
	277	m_allaer_pi(:,:,id_AIBCM_phy) = bcins_pi(:,:) /1.e9/zrho(:,:) ! AIBCM pre-ind
	278	m_allaer_pi(:,:,id_ASNO3M_phy) = nitracc_pi(:,:) /1.e9/zrho(:,:) ! ASNO3M pre-ind
	279	m_allaer_pi(:,:,id_CSNO3M_phy) = nitrcoarse_pi(:,:) /1.e9/zrho(:,:) ! CSNO3M pre-ind
	280	m_allaer_pi(:,:,id_CINO3M_phy) = nitrinscoarse_pi(:,:)/1.e9/zrho(:,:) ! CINO3M pre-ind
	281	m_allaer_pi(:,:,id_AIPOMM_phy) = pomins_pi(:,:) /1.e9/zrho(:,:) ! AIPOMM pre-ind
	282	m_allaer_pi(:,:,id_STRAT_phy) = 0.0
	283
	284	!
	285	! Calculate the total mass of all soluble aersosols (in unit ug /m3)
	286	! to be revisited for AR6
	287	mass_solu_aero(:,:) = sulfacc(:,:) + bcsol(:,:) + pomsol(:,:) + nitracc(:,:) + ssacu(:,:)
	288	mass_solu_aero_pi(:,:) = sulfacc_pi(:,:) + bcsol_pi(:,:) + pomsol_pi(:,:) + nitracc_pi(:,:) + ssacu_pi(:,:)
	289
	290	!****************************************************************************************
	291	! 2) Calculate optical properties for the aerosols
	292	!
	293	!****************************************************************************************
	294	DO k = 1, klev
	295	DO i = 1, klon
	296	pdel(i,k) = paprs(i,k) - paprs (i,k+1)
	297	ENDDO
	298	ENDDO
	299
	300	! Diagnostics calculation for CMIP5 protocol unit (in unit kg/m3)
	301	sconcso4(:) =m_allaer(:,1,id_ASSO4M_phy)*1.e-9
	302	sconcno3(:) =(m_allaer(:,1,id_ASNO3M_phy)+m_allaer(:,1,id_CSNO3M_phy)+m_allaer(:,1,id_CINO3M_phy))*1.e-9
	303	sconcoa(:) =(m_allaer(:,1,id_ASPOMM_phy)+m_allaer(:,1,id_AIPOMM_phy))*1.e-9
	304	sconcbc(:) =(m_allaer(:,1,id_ASBCM_phy)+m_allaer(:,1,id_AIBCM_phy))*1.e-9
	305	sconcss(:) =(m_allaer(:,1,id_ASSSM_phy)+m_allaer(:,1,id_CSSSM_phy)+m_allaer(:,1,id_SSSSM_phy))*1.e-9
	306	sconcdust(:) =m_allaer(:,1,id_CIDUSTM_phy)*1.e-9
	307	concso4(:,:) =m_allaer(:,:,id_ASSO4M_phy)*1.e-9
	308	concno3(:,:) =(m_allaer(:,:,id_ASNO3M_phy)+m_allaer(:,:,id_CSNO3M_phy)+m_allaer(:,:,id_CINO3M_phy))*1.e-9
	309	concoa(:,:) =(m_allaer(:,:,id_ASPOMM_phy)+m_allaer(:,:,id_AIPOMM_phy))*1.e-9
	310	concbc(:,:) =(m_allaer(:,:,id_ASBCM_phy)+m_allaer(:,:,id_AIBCM_phy))*1.e-9
	311	concss(:,:) =(m_allaer(:,:,id_ASSSM_phy)+m_allaer(:,:,id_CSSSM_phy)+m_allaer(:,:,id_SSSSM_phy))*1.e-9
	312	concdust(:,:)=m_allaer(:,:,id_CIDUSTM_phy)*1.e-9
	313
	314	END SUBROUTINE readaerosol_optic_ecrad

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format