Re: MOZBC aerosol mapping (CHEM=202)

646 views

Skip to first unread message

Rajesh Kumar

unread,

Nov 19, 2021, 12:45:09 PM11/19/21

to ummed singh, wrf-chem-fire_emiss

Dear Ummed,

Yes, this mapping is for the MAM4 aerosols and we need to modify the aerosols mapping for the MOSAIC aerosol scheme. It is Ok to provide total emissions of OC, BC, PM2.5, and PM10 through wrffirechemi* files to WRF-Chem even if you are using the MOSAIC scheme. In module_mosaic_addemiss.F, these total emissions are assigned to different MOSAIC bins. So, you can use the following mapping for aerosols in your fire_emis namelist:

'oc->1.4*OC;aerosol',

'bc->BC;aerosol',

'pm25->PM25 + -1.*OC + 1.*BC;aerosol'

'pm10->PM10 + -1.*PM25;aerosol'

I suggest using MODIS+VIIRS dataset if that is available for your target run year. You can stick to the NMVOC mapping from the example file.

Thanks,
Rajesh

On Fri, Nov 19, 2021 at 9:54 AM ummed singh <umme...@gmail.com> wrote:

Dear Sir,

I am title bit confused about mapping here for FINN2.4.

" &control
fire_directory = '/data16b/emmons/finn2/orig_txt_files/',
fire_filename(1) = 'MODVIRS_2019_03312020_take2_MOZ.txt',
mdlFilenm = 'ne30np4_091226_pentagons.nc',
mdlDir = '../grid_files/',
start_date = '2019-01-01',
end_date = '2019-01-05',
output_timing = 'daily',
Model = 'CAMSE',
resol = 'ne30',
EmisType = 'bb_surface',
defaultUnits = 'molecules/cm^2/s',
FinnVers = '2.2modvirs',

glb2fire_map = 'CO -> CO',
'CO2 -> CO2',
'CH4 -> CH4',
'NO -> NO',
'NO2 -> NO2',
'SO2 -> SO2',
'NH3 -> NH3',
'BIGALK -> BIGALK',
'BIGENE -> BIGENE',
'BZALD -> BZALD',
'CH3OH -> CH3OH',
'C2H5OH -> CH3CH2OH',
'CRESOL -> CRESOL',
'GLYALD -> GLYALD',
'HCOOH -> HCOOH',
'HONO -> HONO',
'HYAC -> HYAC',
'ISOP -> ISOP',
'MACR -> MACR',
'MEK -> MEK',
'MTERP -> APIN+BPIN+LIMON+MYRC',
'MVK -> MVK',
'PHENOL -> PHENOL',
'HCN->HCN',
'CH3CN -> CH3CN',
'TOLUENE -> TOLUENE',
'BENZENE -> BENZENE',
'XYLENES -> XYLENE',
'XYLOL -> XYLOL',
'pom_a4->1.4*OC;aerosol',
'bc_a4->BC;aerosol',
'num_bc_a4((particles/cm2/s)(molecules/mole)(g/kg))->5.60298303e18*BC;aerosol',
'num_pom_a4((particles/cm2/s)(molecules/mole)(g/kg))->1.33350996e19*OC;aerosol',
'SVOC->0.03251613*OC;aerosol',
'IVOC->0.04565217*C3H6+0.04782609*C3H8+0.03260870*C2H6+0.03043478*C2H4+0.06086957*BIGENE+0.0782609*BIGALK+0.06304348*CH3COCH3+0.0782609*MEK+0.04782609*CH3CHO+0.03260870*CH2O+0.0848*BENZENE+0.1*TOLUENE+0.1152*XYLENE'
/

My queries

1. Above mapping of species according to MAM4 so we should change it as wrfinput_do1 (chem=202) or is it work well?
like co=CO, co2=CO2..etc

2. Here mapping does not considered PM2.5,10, NO, CO, OC, BC and NMVOC ..etc. should we add these species in mapping?
https://www.acom.ucar.edu/Data/fire/data/finn2/README_FINNv22_November2020.pdf

3. Here I have download fir data for MODIS+VIIRS or only MODIS emission ..which is best for our run ? ...https://www.acom.ucar.edu/Data/fire/

With regards
Ummed Singh

On Fri, Nov 19, 2021 at 3:13 AM Rajesh Kumar <rku...@ucar.edu> wrote:
Dear Ummed,

You can look at o3_BXS and if you see positive definite values it means that species concentrations have been mapped correctly.

Thanks,
Rajesh

On Thu, Nov 18, 2021 at 7:08 AM ummed singh <umme...@gmail.com> wrote:
Dear Sir,

Thanks for your reply . Issue resolved.

Kindly tell me about the MOZBC update.

I have run real.exe with a nested domain . Therefore, MOZBC runs two times to update wrfinpute_d01 and d02.

Now, Which parameters will I check to find whether the wrfbdy_d01 is updated or not? or how ?

with regards
Ummed Singh

On Thu, Nov 18, 2021 at 7:12 PM Rajesh Kumar <rku...@ucar.edu> wrote:
Dear Ummed,

Could you please share the full output of mozbc run? The problem does not seem to be related to the installation.

Thanks,
Rajesh

On Wed, Nov 17, 2021 at 7:52 PM ummed singh <umme...@gmail.com> wrote:
Dear Sir,

I have download mozbc from this link

https://www.acom.ucar.edu/wrf-chem/download.shtml

which is compiled with serial netcdf3.6.3 and gfortran .....compilation shown ok ...

now it shows an error when run ..

" NetCDF: Unknown file format

fail to process netCDF file... "

If a new updated tool for mozbc is available for v4.3 then share the link. sir

With regards
Ummed Singh

On Wed, Nov 17, 2021 at 7:43 AM ummed singh <umme...@gmail.com> wrote:
Dear Sir,

I have checked rsl.error.0000 or rsl.out.0000 which showed real.exe successfully. then checked wrfinput_d02 ...does not show any error during ncdump -h.

what you suggest about these boundary conditions opt in the namelist, are it right? or need any modification.

&bdy_control
spec_bdy_width = 5, ! total number of rows for specified boundary value nudging
spec_zone = 1, ! number of points in specified zone
relax_zone = 4, ! number of points in relaxation zone
specified = .true., .false., .false., ! specified boundary condition
nested = .false., .true., .true., ! nested boundary conditions
/

with regards
Ummed SIngh

On Tue, Nov 16, 2021 at 11:36 PM Rajesh Kumar <rku...@ucar.edu> wrote:
Dear Ummed,

Can you try to ncview the file to make sure your wrfinput_d02 is created successfully by real.exe?

Rajesh

On Tue, Nov 16, 2021 at 10:57 AM ummed singh <umme...@gmail.com> wrote:
Dear Sir,

Thank you very much for the detailed explanation about speciation and its importance.

" I got this error during the mozbc run.

0.0000000E+00 5.6600001E-02 T T T T T T 1.000000
main: nspec = 80
Failed to open /home/IMDLRF/imdwrf/Atmos_AQ/WRF-4.3/test/em_real/wrfinput_d02
with open_option = 1
NetCDF: Unknown file format

fail to process netCDF file... "

But my real.exe run successfully.

With regards
Ummed Singh

On Tue, Nov 16, 2021 at 9:49 PM Rajesh Kumar <rku...@ucar.edu> wrote:
Dear Ummed,

Sulf in the gas-phase is actually the same as sulfate in aerosol phase but sulf represent sulfuric acid which is non-voaltile at atmospheric conditions. So, if you leave out sulf, you will lose all the sulfate mass.

Rajesh

On Tue, Nov 16, 2021 at 9:12 AM ummed singh <umme...@gmail.com> wrote:
Dear Sir,

If we do not consider or remove sulf from the gas phase which type of uncertainty may add in result.

With regards
Ummed Singh

On Tue, Nov 16, 2021 at 8:11 PM Rajesh Kumar <rku...@ucar.edu> wrote:
Dear Ummed,

Note that sulf is MOZCART is in gas-phase and the units are in ppmv. You can use the mapping you suggested but you need to convert the units from ug/kg to ppmv.

Rajesh

On Tue, Nov 16, 2021 at 7:31 AM ummed singh <umme...@gmail.com> wrote:
Dear Sir,

What about 'so4' mapping ?

can we map so4 for mozart as

'sulf -> 1.0* so4_a1+1.0*so4_a2+1.0*so4_a3'

to match the unit in both MOZART and wrf-chem ,...use. ..Alma method ... CAM-Chem aerosol species, including sulfates (so4_a1, so4_a2 and so4_a3), are in kg/kg, so then I have mapped sulfate using the following equation

ug/kg_dry_air = VMR * Mw(aero_mozart) / Mw(air) * 1e9 (https://www.acom.ucar.edu/wrf-chem/ConversionMozartWRFChem.pdf)

which resulted in

'sulf -> 3.31*so4_a1+3.31*so4_a2+3.31*so4_a3;1.e9'

Is this right ?

with regards
Ummed Singh

On Tue, Nov 16, 2021 at 6:36 PM Rajesh Kumar <rku...@ucar.edu> wrote:
Dear Ummed,

These mapping factors represent the mapping from the log-normal size distribution (MAM in CAM-Chem) to the sectional size distribution (MOSAIC in WRF-Chem). Alma Hodzic is a scientist at NCAR who has come up with these factors. I think if you look at the MAM and MOSAIC description papers and get the mean diameters for different sizes, you can derive these factors yourself.

Best,
Rajesh

On Mon, Nov 15, 2021 at 11:55 PM ummed singh <umme...@gmail.com> wrote:
Dear Sir,

My concern about the factors which used to multiply in the section bin. How they are selecting. and 2nd about so4 mapping in gas mapping (MOZART).

If we do not provide the "sulf" mapping. Does it not influence the secondary aerosol formation? and it also does not have NH4 gas species.

With regards
Ummed Singh

On Tue, Nov 16, 2021 at 12:10 PM Rajesh Kumar <rku...@ucar.edu> wrote:
Dear Ummed,

Thanks for sending the header. It looks like the
CAM-Chem output you are using is based on a newer version of the model that includes MAM representation of the aerosols. So, please use option 1 that you shared in your initial email.

Thanks,
Rajesh

Sent from my iPhone

On Nov 15, 2021, at 11:36 PM, ummed singh <umme...@gmail.com> wrote:

Dear Sir,

Find the whole record mention in camchem file

(base) [AEROIII@arclsv MOZBC]$ ncdump -h camchem-20211111231108812468.nc
netcdf camchem-20211111231108812468 {
dimensions:
time = UNLIMITED ; // (121 currently)
lev = 56 ;
lat = 48 ;
lon = 41 ;
chars = 8 ;
ilev = 57 ;
nbnd = 2 ;
variables:
float ALKNIT(time, lev, lat, lon) ;
ALKNIT:mdims = 1 ;
ALKNIT:units = "mol/mol" ;
ALKNIT:long_name = "ALKNIT concentration" ;
float BCARY(time, lev, lat, lon) ;
BCARY:mdims = 1 ;
BCARY:units = "mol/mol" ;
BCARY:long_name = "BCARY concentration" ;
float BENZENE(time, lev, lat, lon) ;
BENZENE:mdims = 1 ;
BENZENE:units = "mol/mol" ;
BENZENE:long_name = "BENZENE concentration" ;
float BIGALD(time, lev, lat, lon) ;
BIGALD:mdims = 1 ;
BIGALD:units = "mol/mol" ;
BIGALD:long_name = "BIGALD concentration" ;
float BIGALD1(time, lev, lat, lon) ;
BIGALD1:mdims = 1 ;
BIGALD1:units = "mol/mol" ;
BIGALD1:long_name = "BIGALD1 concentration" ;
float BIGALD2(time, lev, lat, lon) ;
BIGALD2:mdims = 1 ;
BIGALD2:units = "mol/mol" ;
BIGALD2:long_name = "BIGALD2 concentration" ;
float BIGALD3(time, lev, lat, lon) ;
BIGALD3:mdims = 1 ;
BIGALD3:units = "mol/mol" ;
BIGALD3:long_name = "BIGALD3 concentration" ;
float BIGALD4(time, lev, lat, lon) ;
BIGALD4:mdims = 1 ;
BIGALD4:units = "mol/mol" ;
BIGALD4:long_name = "BIGALD4 concentration" ;
float BIGALK(time, lev, lat, lon) ;
BIGALK:mdims = 1 ;
BIGALK:units = "mol/mol" ;
BIGALK:long_name = "BIGALK concentration" ;
float BIGENE(time, lev, lat, lon) ;
BIGENE:mdims = 1 ;
BIGENE:units = "mol/mol" ;
BIGENE:long_name = "BIGENE concentration" ;
float C2H2(time, lev, lat, lon) ;
C2H2:mdims = 1 ;
C2H2:units = "mol/mol" ;
C2H2:long_name = "C2H2 concentration" ;
float C2H4(time, lev, lat, lon) ;
C2H4:mdims = 1 ;
C2H4:units = "mol/mol" ;
C2H4:long_name = "C2H4 concentration" ;
float C2H5OH(time, lev, lat, lon) ;
C2H5OH:mdims = 1 ;
C2H5OH:units = "mol/mol" ;
C2H5OH:long_name = "C2H5OH concentration" ;
float C2H6(time, lev, lat, lon) ;
C2H6:mdims = 1 ;
C2H6:units = "mol/mol" ;
C2H6:long_name = "C2H6 concentration" ;
float C3H6(time, lev, lat, lon) ;
C3H6:mdims = 1 ;
C3H6:units = "mol/mol" ;
C3H6:long_name = "C3H6 concentration" ;
float C3H8(time, lev, lat, lon) ;
C3H8:mdims = 1 ;
C3H8:units = "mol/mol" ;
C3H8:long_name = "C3H8 concentration" ;
float CH2O(time, lev, lat, lon) ;
CH2O:mdims = 1 ;
CH2O:units = "mol/mol" ;
CH2O:long_name = "CH2O concentration" ;
float CH3CHO(time, lev, lat, lon) ;
CH3CHO:mdims = 1 ;
CH3CHO:units = "mol/mol" ;
CH3CHO:long_name = "CH3CHO concentration" ;
float CH3CN(time, lev, lat, lon) ;
CH3CN:mdims = 1 ;
CH3CN:units = "mol/mol" ;
CH3CN:long_name = "CH3CN concentration" ;
float CH3COCH3(time, lev, lat, lon) ;
CH3COCH3:mdims = 1 ;
CH3COCH3:units = "mol/mol" ;
CH3COCH3:long_name = "CH3COCH3 concentration" ;
float CH3COCHO(time, lev, lat, lon) ;
CH3COCHO:mdims = 1 ;
CH3COCHO:units = "mol/mol" ;
CH3COCHO:long_name = "CH3COCHO concentration" ;
float CH3COOH(time, lev, lat, lon) ;
CH3COOH:mdims = 1 ;
CH3COOH:units = "mol/mol" ;
CH3COOH:long_name = "CH3COOH concentration" ;
float CH3OH(time, lev, lat, lon) ;
CH3OH:mdims = 1 ;
CH3OH:units = "mol/mol" ;
CH3OH:long_name = "CH3OH concentration" ;
float CH3OOH(time, lev, lat, lon) ;
CH3OOH:mdims = 1 ;
CH3OOH:units = "mol/mol" ;
CH3OOH:long_name = "CH3OOH concentration" ;
float CH4(time, lev, lat, lon) ;
CH4:mdims = 1 ;
CH4:units = "mol/mol" ;
CH4:long_name = "CH4 concentration" ;
float CO(time, lev, lat, lon) ;
CO:mdims = 1 ;
CO:units = "mol/mol" ;
CO:long_name = "CO concentration" ;
float CRESOL(time, lev, lat, lon) ;
CRESOL:mdims = 1 ;
CRESOL:units = "mol/mol" ;
CRESOL:long_name = "CRESOL concentration" ;
float DMS(time, lev, lat, lon) ;
DMS:mdims = 1 ;
DMS:units = "mol/mol" ;
DMS:long_name = "DMS concentration" ;
float GLYOXAL(time, lev, lat, lon) ;
GLYOXAL:mdims = 1 ;
GLYOXAL:units = "mol/mol" ;
GLYOXAL:long_name = "GLYOXAL concentration" ;
float H2O(time, lev, lat, lon) ;
H2O:mdims = 1 ;
H2O:units = "mol/mol" ;
H2O:long_name = "water vapor concentration" ;
float H2O2(time, lev, lat, lon) ;
H2O2:mdims = 1 ;
H2O2:units = "mol/mol" ;
H2O2:long_name = "H2O2 concentration" ;
float HCN(time, lev, lat, lon) ;
HCN:mdims = 1 ;
HCN:units = "mol/mol" ;
HCN:long_name = "HCN concentration" ;
float HCOOH(time, lev, lat, lon) ;
HCOOH:mdims = 1 ;
HCOOH:units = "mol/mol" ;
HCOOH:long_name = "HCOOH concentration" ;
float HNO3(time, lev, lat, lon) ;
HNO3:mdims = 1 ;
HNO3:units = "mol/mol" ;
HNO3:long_name = "HNO3 concentration" ;
float HO2(time, lev, lat, lon) ;
HO2:mdims = 1 ;
HO2:units = "mol/mol" ;
HO2:long_name = "HO2 concentration" ;
float HO2NO2(time, lev, lat, lon) ;
HO2NO2:mdims = 1 ;
HO2NO2:units = "mol/mol" ;
HO2NO2:long_name = "HO2NO2 concentration" ;
float HONITR(time, lev, lat, lon) ;
HONITR:mdims = 1 ;
HONITR:units = "mol/mol" ;
HONITR:long_name = "HONITR concentration" ;
float HYAC(time, lev, lat, lon) ;
HYAC:mdims = 1 ;
HYAC:units = "mol/mol" ;
HYAC:long_name = "HYAC concentration" ;
float ISOP(time, lev, lat, lon) ;
ISOP:mdims = 1 ;
ISOP:units = "mol/mol" ;
ISOP:long_name = "ISOP concentration" ;
float ISOPNITA(time, lev, lat, lon) ;
ISOPNITA:mdims = 1 ;
ISOPNITA:units = "mol/mol" ;
ISOPNITA:long_name = "ISOPNITA concentration" ;
float ISOPNITB(time, lev, lat, lon) ;
ISOPNITB:mdims = 1 ;
ISOPNITB:units = "mol/mol" ;
ISOPNITB:long_name = "ISOPNITB concentration" ;
float IVOC(time, lev, lat, lon) ;
IVOC:mdims = 1 ;
IVOC:units = "mol/mol" ;
IVOC:long_name = "IVOC concentration" ;
float MACR(time, lev, lat, lon) ;
MACR:mdims = 1 ;
MACR:units = "mol/mol" ;
MACR:long_name = "MACR concentration" ;
float MEK(time, lev, lat, lon) ;
MEK:mdims = 1 ;
MEK:units = "mol/mol" ;
MEK:long_name = "MEK concentration" ;
float MPAN(time, lev, lat, lon) ;
MPAN:mdims = 1 ;
MPAN:units = "mol/mol" ;
MPAN:long_name = "MPAN concentration" ;
float MTERP(time, lev, lat, lon) ;
MTERP:mdims = 1 ;
MTERP:units = "mol/mol" ;
MTERP:long_name = "MTERP concentration" ;
float MVK(time, lev, lat, lon) ;
MVK:mdims = 1 ;
MVK:units = "mol/mol" ;
MVK:long_name = "MVK concentration" ;
float N2O(time, lev, lat, lon) ;
N2O:mdims = 1 ;
N2O:units = "mol/mol" ;
N2O:long_name = "N2O concentration" ;
float N2O5(time, lev, lat, lon) ;
N2O5:mdims = 1 ;
N2O5:units = "mol/mol" ;
N2O5:long_name = "N2O5 concentration" ;
float NH3(time, lev, lat, lon) ;
NH3:mdims = 1 ;
NH3:units = "mol/mol" ;
NH3:long_name = "NH3 concentration" ;
float NH4(time, lev, lat, lon) ;
NH4:mdims = 1 ;
NH4:units = "mol/mol" ;
NH4:long_name = "NH4 concentration" ;
float NO(time, lev, lat, lon) ;
NO:mdims = 1 ;
NO:units = "mol/mol" ;
NO:long_name = "NO concentration" ;
float NO2(time, lev, lat, lon) ;
NO2:mdims = 1 ;
NO2:units = "mol/mol" ;
NO2:long_name = "NO2 concentration" ;
float NO3(time, lev, lat, lon) ;
NO3:mdims = 1 ;
NO3:units = "mol/mol" ;
NO3:long_name = "NO3 concentration" ;
float NOA(time, lev, lat, lon) ;
NOA:mdims = 1 ;
NOA:units = "mol/mol" ;
NOA:long_name = "NOA concentration" ;
float O3(time, lev, lat, lon) ;
O3:mdims = 1 ;
O3:units = "mol/mol" ;
O3:long_name = "O3 concentration" ;
float O3S(time, lev, lat, lon) ;
O3S:mdims = 1 ;
O3S:units = "mol/mol" ;
O3S:long_name = "O3S concentration" ;
float OH(time, lev, lat, lon) ;
OH:mdims = 1 ;
OH:units = "mol/mol" ;
OH:long_name = "OH concentration" ;
float ONITR(time, lev, lat, lon) ;
ONITR:mdims = 1 ;
ONITR:units = "mol/mol" ;
ONITR:long_name = "ONITR concentration" ;
double P0 ;
P0:long_name = "reference pressure" ;
P0:units = "Pa" ;
float PAN(time, lev, lat, lon) ;
PAN:mdims = 1 ;
PAN:units = "mol/mol" ;
PAN:long_name = "PAN concentration" ;
float PBZNIT(time, lev, lat, lon) ;
PBZNIT:mdims = 1 ;
PBZNIT:units = "mol/mol" ;
PBZNIT:long_name = "PBZNIT concentration" ;
float PHENOL(time, lev, lat, lon) ;
PHENOL:mdims = 1 ;
PHENOL:units = "mol/mol" ;
PHENOL:long_name = "PHENOL concentration" ;
float PM25(time, lev, lat, lon) ;
PM25:mdims = 1 ;
PM25:units = "kg/m3" ;
PM25:long_name = "PM2.5 concentration" ;
float PS(time, lat, lon) ;
PS:units = "Pa" ;
PS:long_name = "Surface pressure" ;
float Q(time, lev, lat, lon) ;
Q:mdims = 1 ;
Q:units = "kg/kg" ;
Q:long_name = "Specific humidity" ;
float SO2(time, lev, lat, lon) ;
SO2:mdims = 1 ;
SO2:units = "mol/mol" ;
SO2:long_name = "SO2 concentration" ;
float SVOC(time, lev, lat, lon) ;
SVOC:mdims = 1 ;
SVOC:units = "mol/mol" ;
SVOC:long_name = "SVOC concentration" ;
float T(time, lev, lat, lon) ;
T:mdims = 1 ;
T:units = "K" ;
T:long_name = "Temperature" ;
float TERPNIT(time, lev, lat, lon) ;
TERPNIT:mdims = 1 ;
TERPNIT:units = "mol/mol" ;
TERPNIT:long_name = "TERPNIT concentration" ;
float TOLUENE(time, lev, lat, lon) ;
TOLUENE:mdims = 1 ;
TOLUENE:units = "mol/mol" ;
TOLUENE:long_name = "TOLUENE concentration" ;
float U(time, lev, lat, lon) ;
U:mdims = 1 ;
U:units = "m/s" ;
U:long_name = "Zonal wind" ;
float V(time, lev, lat, lon) ;
V:mdims = 1 ;
V:units = "m/s" ;
V:long_name = "Meridional wind" ;
float XYLENES(time, lev, lat, lon) ;
XYLENES:mdims = 1 ;
XYLENES:units = "mol/mol" ;
XYLENES:long_name = "XYLENES concentration" ;
float Z3(time, lev, lat, lon) ;
Z3:mdims = 1 ;
Z3:units = "m" ;
Z3:long_name = "Geopotential Height (above sea level)" ;
float bc_a1(time, lev, lat, lon) ;
bc_a1:mdims = 1 ;
bc_a1:units = "kg/kg" ;
bc_a1:long_name = "bc_a1 concentration" ;
float bc_a4(time, lev, lat, lon) ;
bc_a4:mdims = 1 ;
bc_a4:units = "kg/kg" ;
bc_a4:long_name = "bc_a4 concentration" ;
double ch4vmr(time) ;
ch4vmr:long_name = "ch4 volume mixing ratio" ;
double co2vmr(time) ;
co2vmr:long_name = "co2 volume mixing ratio" ;
int date(time) ;
date:long_name = "current date (YYYYMMDD)" ;
char date_written(time, chars) ;
int datesec(time) ;
datesec:long_name = "current seconds of current date" ;
float dst_a1(time, lev, lat, lon) ;
dst_a1:mdims = 1 ;
dst_a1:units = "kg/kg" ;
dst_a1:long_name = "dst_a1 concentration" ;
float dst_a2(time, lev, lat, lon) ;
dst_a2:mdims = 1 ;
dst_a2:units = "kg/kg" ;
dst_a2:long_name = "dst_a2 concentration" ;
float dst_a3(time, lev, lat, lon) ;
dst_a3:mdims = 1 ;
dst_a3:units = "kg/kg" ;
dst_a3:long_name = "dst_a3 concentration" ;
double f11vmr(time) ;
f11vmr:long_name = "f11 volume mixing ratio" ;
double f12vmr(time) ;
f12vmr:long_name = "f12 volume mixing ratio" ;
double gw(lat) ;
gw:_FillValue = -900. ;
gw:long_name = "latitude weights" ;
double hyai(ilev) ;
hyai:long_name = "hybrid A coefficient at layer interfaces" ;
double hyam(lev) ;
hyam:long_name = "hybrid A coefficient at layer midpoints" ;
double hybi(ilev) ;
hybi:long_name = "hybrid B coefficient at layer interfaces" ;
double hybm(lev) ;
hybm:long_name = "hybrid B coefficient at layer midpoints" ;
double ilev(ilev) ;
ilev:long_name = "hybrid level at interfaces (1000*(A+B))" ;
ilev:units = "hPa" ;
ilev:positive = "down" ;
ilev:standard_name = "atmosphere_hybrid_sigma_pressure_coordinate" ;
ilev:formula_terms = "a: hyai b: hybi p0: P0 ps: PS" ;
double lat(lat) ;
lat:_FillValue = -900. ;
lat:long_name = "latitude" ;
lat:units = "degrees_north" ;
double lev(lev) ;
lev:long_name = "hybrid level at midpoints (1000*(A+B))" ;
lev:units = "hPa" ;
lev:positive = "down" ;
lev:standard_name = "atmosphere_hybrid_sigma_pressure_coordinate" ;
lev:formula_terms = "a: hyam b: hybm p0: P0 ps: PS" ;
double lon(lon) ;
lon:_FillValue = -900. ;
lon:long_name = "longitude" ;
lon:units = "degrees_east" ;
int mdt ;
mdt:long_name = "timestep" ;
mdt:units = "s" ;
double n2ovmr(time) ;
n2ovmr:long_name = "n2o volume mixing ratio" ;
int nbdate ;
nbdate:long_name = "base date (YYYYMMDD)" ;
int nbsec ;
nbsec:long_name = "seconds of base date" ;
float ncl_a1(time, lev, lat, lon) ;
ncl_a1:mdims = 1 ;
ncl_a1:units = "kg/kg" ;
ncl_a1:long_name = "ncl_a1 concentration" ;
float ncl_a2(time, lev, lat, lon) ;
ncl_a2:mdims = 1 ;
ncl_a2:units = "kg/kg" ;
ncl_a2:long_name = "ncl_a2 concentration" ;
float ncl_a3(time, lev, lat, lon) ;
ncl_a3:mdims = 1 ;
ncl_a3:units = "kg/kg" ;
ncl_a3:long_name = "ncl_a3 concentration" ;
int ndbase ;
ndbase:long_name = "base day" ;
int ndcur(time) ;
ndcur:long_name = "current day (from base day)" ;
int nsbase ;
nsbase:long_name = "seconds of base day" ;
int nscur(time) ;
nscur:long_name = "current seconds of current day" ;
int nsteph(time) ;
nsteph:long_name = "current timestep" ;
float num_a1(time, lev, lat, lon) ;
num_a1:mdims = 1 ;
num_a1:units = " 1/kg" ;
num_a1:long_name = "num_a1 concentration" ;
float num_a2(time, lev, lat, lon) ;
num_a2:mdims = 1 ;
num_a2:units = " 1/kg" ;
num_a2:long_name = "num_a2 concentration" ;
float num_a3(time, lev, lat, lon) ;
num_a3:mdims = 1 ;
num_a3:units = " 1/kg" ;
num_a3:long_name = "num_a3 concentration" ;
float pom_a1(time, lev, lat, lon) ;
pom_a1:mdims = 1 ;
pom_a1:units = "kg/kg" ;
pom_a1:long_name = "pom_a1 concentration" ;
float pom_a4(time, lev, lat, lon) ;
pom_a4:mdims = 1 ;
pom_a4:units = "kg/kg" ;
pom_a4:long_name = "pom_a4 concentration" ;
float so4_a1(time, lev, lat, lon) ;
so4_a1:mdims = 1 ;
so4_a1:units = "kg/kg" ;
so4_a1:long_name = "so4_a1 concentration" ;
float so4_a2(time, lev, lat, lon) ;
so4_a2:mdims = 1 ;
so4_a2:units = "kg/kg" ;
so4_a2:long_name = "so4_a2 concentration" ;
float so4_a3(time, lev, lat, lon) ;
so4_a3:mdims = 1 ;
so4_a3:units = "kg/kg" ;
so4_a3:long_name = "so4_a3 concentration" ;
float soa1_a1(time, lev, lat, lon) ;
soa1_a1:mdims = 1 ;
soa1_a1:units = "kg/kg" ;
soa1_a1:long_name = "soa1_a1 concentration" ;
float soa1_a2(time, lev, lat, lon) ;
soa1_a2:mdims = 1 ;
soa1_a2:units = "kg/kg" ;
soa1_a2:long_name = "soa1_a2 concentration" ;
float soa2_a1(time, lev, lat, lon) ;
soa2_a1:mdims = 1 ;
soa2_a1:units = "kg/kg" ;
soa2_a1:long_name = "soa2_a1 concentration" ;
float soa2_a2(time, lev, lat, lon) ;
soa2_a2:mdims = 1 ;
soa2_a2:units = "kg/kg" ;
soa2_a2:long_name = "soa2_a2 concentration" ;
float soa3_a1(time, lev, lat, lon) ;
soa3_a1:mdims = 1 ;
soa3_a1:units = "kg/kg" ;
soa3_a1:long_name = "soa3_a1 concentration" ;
float soa3_a2(time, lev, lat, lon) ;
soa3_a2:mdims = 1 ;
soa3_a2:units = "kg/kg" ;
soa3_a2:long_name = "soa3_a2 concentration" ;
float soa4_a1(time, lev, lat, lon) ;
soa4_a1:mdims = 1 ;
soa4_a1:units = "kg/kg" ;
soa4_a1:long_name = "soa4_a1 concentration" ;
float soa4_a2(time, lev, lat, lon) ;
soa4_a2:mdims = 1 ;
soa4_a2:units = "kg/kg" ;
soa4_a2:long_name = "soa4_a2 concentration" ;
float soa5_a1(time, lev, lat, lon) ;
soa5_a1:mdims = 1 ;
soa5_a1:units = "kg/kg" ;
soa5_a1:long_name = "soa5_a1 concentration" ;
float soa5_a2(time, lev, lat, lon) ;
soa5_a2:mdims = 1 ;
soa5_a2:units = "kg/kg" ;
soa5_a2:long_name = "soa5_a2 concentration" ;
double sol_tsi(time) ;
sol_tsi:long_name = "total solar irradiance" ;
sol_tsi:units = "W/m2" ;
double time(time) ;
time:long_name = "time" ;
time:units = "days since 2000-03-01 00:00:00" ;
time:calendar = "gregorian" ;
time:bounds = "time_bnds" ;
double time_bnds(time, nbnd) ;
time_bnds:long_name = "time interval endpoints" ;
char time_written(time, chars) ;

// global attributes:
:Conventions = "CF-1.0" ;
:source = "CAM" ;
:case = "fmerra.2.1003.FCSD.f09.qfedcmip.56L.001.branch02" ;
:logname = "buchholz" ;
:host = "r6i6n33" ;
:initial_file = "/glade/p/cesmdata/cseg/inputdata/atm/cam/inic/fv/f.e20.FCSD.f09_f09_mg17.cesm2.1-exp002.001.cam.i.2005-01-01-00000_c180801.nc" ;
:topography_file = "/glade/p/cesmdata/cseg/inputdata/atm/cam/met/MERRA2/0.9x1.25/fv_0.9x1.25_nc3000_Nsw042_Nrs008_Co060_Fi001_ZR_sgh30_24km_GRNL_MERRA2_c171218.nc" ;
:model_doi_url = "https://doi.org/10.5065/D67H1H0V" ;
:time_period_freq = "hour_6" ;
:history = "Thu Nov 11 23:26:53 2021: /usr/local/nco-4.7.9/bin/ncrcat --ovr stage/525725.subset.fmerra.2.1003.FCSD.f09.qfedcmip.56L.001.branch02.cam.h1.2018-09-30-21600.nc stage/525725.subset.fmerra.2.1003.FCSD.f09.qfedcmip.56L.001.branch02.cam.h1.2018-10-01-21600.nc stage/525725.subset.fmerra.2.1003.FCSD.f09.qfedcmip.56L.001.branch02.cam.h1.2018-10-02-21600.nc stage/525725.subset.fmerra.2.1003.FCSD.f09.qfedcmip.56L.001.branch02.cam.h1.2018-10-03-21600.nc stage/525725.subset.fmerra.2.1003.FCSD.f09.qfedcmip.56L.001.branch02.cam.h1.2018-10-04-21600.nc stage/525725.subset.fmerra.2.1003.FCSD.f09.qfedcmip.56L.001.branch02.cam.h1.2018-10-05-21600.nc stage/525725.subset.fmerra.2.1003.FCSD.f09.qfedcmip.56L.001.branch02.cam.h1.2018-10-06-21600.nc stage/525725.subset.fmerra.2.1003.FCSD.f09.qfedcmip.56L.001.branch02.cam.h1.2018-10-07-21600.nc stage/525725.subset.fmerra.2.1003.FCSD.f09.qfedcmip.56L.001.branch02.cam.h1.2018-10-08-21600.nc stage/525725.subset.fmerra.2.1003.FCSD.f09.qfedcmip.56L.001.branch02.cam.h1.2018-10-09-21600.nc stage/525725.subset.fmerra.2.1003.FCSD.f09.qfedcmip.56L.001.branch02.cam.h1.2018-10-10-21600.nc stage/525725.subset.fmerra.2.1003.FCSD.f09.qfedcmip.56L.001.branch02.cam.h1.2018-10-11-21600.nc stage/525725.subset.fmerra.2.1003.FCSD.f09.qfedcmip.56L.001.branch02.cam.h1.2018-10-12-21600.nc stage/525725.subset.fmerra.2.1003.FCSD.f09.qfedcmip.56L.001.branch02.cam.h1.2018-10-13-21600.nc stage/525725.subset.fmerra.2.1003.FCSD.f09.qfedcmip.56L.001.branch02.cam.h1.2018-10-14-21600.nc stage/525725.subset.fmerra.2.1003.FCSD.f09.qfedcmip.56L.001.branch02.cam.h1.2018-10-15-21600.nc stage/525725.subset.fmerra.2.1003.FCSD.f09.qfedcmip.56L.001.branch02.cam.h1.2018-10-16-21600.nc stage/525725.subset.fmerra.2.1003.FCSD.f09.qfedcmip.56L.001.branch02.cam.h1.2018-10-17-21600.nc stage/525725.subset.fmerra.2.1003.FCSD.f09.qfedcmip.56L.001.branch02.cam.h1.2018-10-18-21600.nc stage/525725.subset.fmerra.2.1003.FCSD.f09.qfedcmip.56L.001.branch02.cam.h1.2018-10-19-21600.nc stage/525725.subset.fmerra.2.1003.FCSD.f09.qfedcmip.56L.001.branch02.cam.h1.2018-10-20-21600.nc stage/525725.subset.fmerra.2.1003.FCSD.f09.qfedcmip.56L.001.branch02.cam.h1.2018-10-21-21600.nc stage/525725.subset.fmerra.2.1003.FCSD.f09.qfedcmip.56L.001.branch02.cam.h1.2018-10-22-21600.nc stage/525725.subset.fmerra.2.1003.FCSD.f09.qfedcmip.56L.001.branch02.cam.h1.2018-10-23-21600.nc stage/525725.subset.fmerra.2.1003.FCSD.f09.qfedcmip.56L.001.branch02.cam.h1.2018-10-24-21600.nc stage/525725.subset.fmerra.2.1003.FCSD.f09.qfedcmip.56L.001.branch02.cam.h1.2018-10-25-21600.nc stage/525725.subset.fmerra.2.1003.FCSD.f09.qfedcmip.56L.001.branch02.cam.h1.2018-10-26-21600.nc stage/525725.subset.fmerra.2.1003.FCSD.f09.qfedcmip.56L.001.branch02.cam.h1.2018-10-27-21600.nc stage/525725.subset.fmerra.2.1003.FCSD.f09.qfedcmip.56L.001.branch02.cam.h1.2018-10-28-21600.nc stage/525725.subset.fmerra.2.1003.FCSD.f09.qfedcmip.56L.001.branch02.cam.h1.2018-10-29-21600.nc stage/525725.subset.fmerra.2.1003.FCSD.f09.qfedcmip.56L.001.branch02.cam.h1.2018-10-30-21600.nc stage/camchem-20211111231108812468.nc\n",
"Thu Nov 11 23:18:38 2021: ncks -d time,3,,1 -d lat,0.0,45.0 -d lon,50.0,100.0 /glade/collections/rda/data/ds313.7/2018/fmerra.2.1003.FCSD.f09.qfedcmip.56L.001.branch02.cam.h1.2018-09-30-21600.nc 525725.subset.fmerra.2.1003.FCSD.f09.qfedcmip.56L.001.branch02.cam.h1.2018-09-30-21600.nc" ;
:NCO = "netCDF Operators version 4.7.9 (Homepage = http://nco.sf.net, Code = http://github.com/nco/nco)" ;
}

with regards
Ummed Singh

On Tue, Nov 16, 2021 at 11:51 AM Rajesh Kumar <rku...@ucar.edu> wrote:
Dear Ummed,

Please send me the result of ncdump -h camchem-20211111231108812468.nc

Thanks,
Rajesh

Sent from my iPhone

On Nov 15, 2021, at 11:17 PM, ummed singh <umme...@gmail.com> wrote:

Dear Sir,

I am using MOZART/camchem output for MOZBC. but here I have a few dots about aerosol species mapping as per 2nd methods.

I have checked camchem-20211111231108812468.nc data then found that there is no data for these secondary aerosols in global output.

(base) [AEROIII@arclsv MOZBC]$ ncdump -v SOA camchem-20211111231108812468.nc ncdump: SOA: No such variable
(base) [AEROIII@arclsv MOZBC]$ ncdump -v OC1 camchem-20211111231108812468.nc
ncdump: OC1: No such variable
(base) [AEROIII@arclsv MOZBC]$ ncdump -v SO4 camchem-20211111231108812468.nc
ncdump: SO4: No such variable
(base) [AEROIII@arclsv MOZBC]$ ncdump -v NH4NO3 camchem-20211111231108812468.nc
ncdump: NH4NO3: No such variable
(base) [AEROIII@arclsv MOZBC]$ ncdump -v SA1 camchem-20211111231108812468.nc ncdump: SA1: No such variable
(base) [AEROIII@arclsv MOZBC]$ ncdump -v DUST1 camchem-20211111231108812468.nc
ncdump: DUST1: No such variable

........
as per MOZART speciation SO4 speciated as "slfo".... but camchem...nc data does not have so4.

What do you suggest in this condition?

With regards
Ummed Singh

On Tue, Nov 16, 2021 at 10:58 AM Rajesh Kumar <rku...@ucar.edu> wrote:
Dear Ummed,

These two methods are corresponding to two different models. The first is to use the WACCM output to provide boundary conditions to your domain and the second uses MOZART/CAM-Chem output. Which one did you download? You can send me a header of the global model output and I can tell you what model output you have.

Thanks,
Rajesh

PS: Are you done with the reviewers response? If not, I can work on it tomorrow morning.

Sent from my iPhone

On Nov 15, 2021, at 9:45 PM, ummed singh <umme...@gmail.com> wrote:

Dear sir,

I am a little bit confused about aerosol mapping In Mozart_Mosaic_aq 4bin. Here are two methods for mosaic mapping I have seen.

One is

'oc_a01->0.0093*pom_a1+0.7510*soa1_a2+0.0093*soa1_a1+0.7510*soa2_a2+0.0093*soa2_a1+0.7510*soa3_a2+0.0093*soa3_a1+0.7510*soa4_a2+0.0093*soa4_a1+0.7510*soa5_a2+0.0093*soa5_a1+0.1123*pom_a1+0.2376*soa1_a2+0.1123*soa1_a1+0.2376*soa2_a2+0.1123*soa2_a1+0.2376*soa3_a2+0.1123*soa3_a1+0.2376*soa4_a2+0.1123*soa4_a1+0.2376*soa5_a2+0.1123*soa5_a1;1.e9',
'oc_a02->0.3835*pom_a1+0.0113*soa1_a2+0.3835*soa1_a1+0.0133*soa2_a2+0.3835*soa2_a1+0.0113*soa3_a2+0.3838*soa3_a1+0.0113*soa4_a2+0.3838*soa4_a1+0.0113*soa5_a2+0.3838*soa5_a1+0.3783*pom_a1+0.0001*soa1_a2+0.3783*soa1_a1+0.0001*soa2_a2+0.3783*soa2_a1+0.0001*soa3_a2+0.3783*soa3_a1+0.0001*soa4_a2+0.3783*soa4_a1+0.0001*soa5_a2+0.3783*soa5_a1;1.e9',
'oc_a03->0.1077*pom_a1+0.0000*soa1_a2+0.1077*soa1_a1+0.0000*soa2_a2+0.1077*soa2_a1+0.0000*soa3_a2+0.1077*soa3_a1+0.0000*soa4_a2+0.1077*soa4_a1+0.0000*soa5_a2+0.1077*soa5_a1+0.0087*pom_a1+0.0000*soa1_a2+0.0087*soa1_a1+0.0000*soa2_a2+0.0087*soa2_a1+0.0000*soa3_a2+0.0087*soa3_a1+0.0000*soa4_a2+0.0087*soa4_a1+0.0000*soa5_a2+0.0087*soa5_a1;1.e9',
'oc_a04->0.0002*pom_a1+0.0000*soa1_a2+0.0002*soa1_a1+0.0000*soa2_a2+0.0002*soa2_a1+0.0000*soa3_a2+0.0002*soa3_a1+0.0000*soa4_a2+0.0002*soa4_a1+0.0000*soa5_a2+0.0002*soa5_a1+0.0000*pom_a1+0.0000*soa1_a2+0.0000*soa1_a1+0.0000*soa2_a2+0.0000*soa2_a1+0.0000*soa3_a2+0.0000*soa3_a1+0.0000*soa4_a2+0.0000*soa4_a1+0.0000*soa5_a2+0.0000*soa5_a1;1.e9',
'bc_a01->0.0093*bc_a1+0.0093*bc_a4+0.1123*bc_a1+0.1123*bc_a4;1.e9',
'bc_a02->0.3835*bc_a1+0.3835*bc_a4+0.3783*bc_a1+0.3783*bc_a4;1.e9',
'bc_a03->0.1077*bc_a1+0.1077*bc_a4+0.0087*bc_a1+0.0087*bc_a4;1.e9',
'bc_a04->0.0002*bc_a1+0.0002*bc_a4+0.0000*bc_a1+0.0000*bc_a4;1.e9',
'so4_a01->0.7510*so4_a2+0.0093*so4_a1+0.0000*so4_a3+0.2376*so4_a2+0.1123*so4_a1+0.0000*so4_a3;1.e9',
'so4_a02->0.0113*so4_a2+0.3835*so4_a1+0.0000*so4_a3+0.0001*so4_a2+0.3783*so4_a1+0.0002*so4_a3;1.e9',
'so4_a03->0.0000*so4_a2+0.1077*so4_a1+0.0061*so4_a3+0.0000*so4_a2+0.0087*so4_a1+0.0934*so4_a3;1.e9',
'so4_a04->0.0000*so4_a2+0.0002*so4_a1+0.4020*so4_a3+0.0000*so4_a2+0.0000*so4_a1+0.4983*so4_a3;1.e9',
'nh4_a01->0.1410*so4_a2+0.0033*so4_a1+0.0000*so4_a3+0.0446*so4_a2+0.0017*so4_a1+0.0000*so4_a3;1.e9',
'nh4_a02->0.0021*so4_a2+0.0210*so4_a1+0.0000*so4_a3+0.0000*so4_a2+0.0720*so4_a1+0.0000*so4_a3;1.e9',
'nh4_a03->0.0000*so4_a2+0.0202*so4_a1+0.0011*so4_a3+0.0000*so4_a2+0.0001*so4_a1+0.0175*so4_a3;1.e9',
'nh4_a04->0.0000*so4_a2+0.0000*so4_a1+0.0755*so4_a3+0.0000*so4_a2+0.0000*so4_a1+0.0755*so4_a3;1.e9',
'no3_a01->0.0000*so4_a2+0.0000*so4_a1+0.0000*so4_a3+0.0000*so4_a2+0.0000*so4_a1+0.0000*so4_a3;1.e9',
'no3_a02->0.0000*so4_a2+0.0000*so4_a1+0.0000*so4_a3+0.0000*so4_a2+0.0000*so4_a1+0.0000*so4_a3;1.e9',
'no3_a03->0.0000*so4_a2+0.0000*so4_a1+0.0000*so4_a3+0.0000*so4_a2+0.0000*so4_a1+0.0000*so4_a3;1.e9',
'no3_a04->0.0000*so4_a2+0.0000*so4_a1+0.0000*so4_a3+0.0000*so4_a2+0.0000*so4_a1+0.0000*so4_a3;1.e9',
'na_a01->0.2954*ncl_a2+0.0037*ncl_a1+0.0000*ncl_a3+0.0935*ncl_a2+0.0442*ncl_a1+0.0000*ncl_a3;1.e9',
'na_a02->0.0045*ncl_a2+0.1509*ncl_a1+0.0000*ncl_a3+0.0000*ncl_a2+0.1488*ncl_a1+0.0000*ncl_a3;1.e9',
'na_a03->0.0000*ncl_a2+0.0424*ncl_a1+0.0024*ncl_a3+0.0000*ncl_a2+0.0034*ncl_a1+0.0367*ncl_a3;1.e9',
'na_a04->0.0000*ncl_a2+0.0000*ncl_a1+0.1582*ncl_a3+0.0000*ncl_a2+0.0000*ncl_a1+0.1960*ncl_a3;1.e9',
'cl_a01->0.4555*ncl_a2+0.0056*ncl_a1+0.0000*ncl_a3+0.1441*ncl_a2+0.0681*ncl_a1+0.0000*ncl_a3;1.e9',
'cl_a02->0.0068*ncl_a2+0.2326*ncl_a1+0.0000*ncl_a3+0.0000*ncl_a2+0.2295*ncl_a1+0.0000*ncl_a3;1.e9',
'cl_a03->0.0000*ncl_a2+0.0654*ncl_a1+0.0037*ncl_a3+0.0000*ncl_a2+0.0055*ncl_a1+0.0567*ncl_a3;1.e9',
'cl_a04->0.0000*ncl_a2+0.0001*ncl_a1+0.2439*ncl_a3+0.0000*ncl_a2+0.0000*ncl_a1+0.3023*ncl_a3;1.e9',
'oin_a01->0.7510*dst_a2+0.0093*dst_a1+0.0000*dst_a3+0.2376*dst_a2+0.1123*dst_a1+0.0000*dst_a3;1.e9',
'oin_a02->0.0113*dst_a2+0.3835*dst_a1+0.0000*dst_a3+0.0001*dst_a2+0.3783*dst_a1+0.0002*dst_a3;1.e9',
'oin_a03->0.0000*dst_a2+0.1077*dst_a1+0.0061*dst_a3+0.0000*dst_a2+0.0087*dst_a1+0.0934*dst_a3;1.e9',
'oin_a04->0.0000*dst_a2+0.0002*dst_a1+0.4020*dst_a3+0.0000*dst_a2+0.0000*dst_a1+0.4983*dst_a3;1.e9',
'num_a01->0.9502*num_a2+0.2509*num_a1+0.0000*num_a3+0.0494*num_a2+0.4626*num_a1+0.0000*num_a3;1.0',
'num_a02->0.0004*num_a2+0.2470*num_a1+0.0007*num_a3+0.0000*num_a2+0.0377*num_a1+0.0232*num_a3;1.0',
'num_a03->0.0000*num_a2+0.0016*num_a1+0.1886*num_a3+0.0000*num_a2+0.0000*num_a1+0.4372*num_a3;1.0',
'num_a04->0.0000*num_a2+0.0000*num_a1+0.2935*num_a3+0.0000*num_a2+0.0000*num_a1+0.0566*num_a3;1.0',
/

2nd methods is

'oc_a01->9.60024e-2*OC1+9.60024e-2*OC2+1.152026*SOA;1.e9'
'oc_a02->.254356*OC1+.254356*OC2+3.05226*SOA;1.e9',
'oc_a03->6.00723e-2*OC1+6.00723e-2*OC2+.720865*SOA;1.e9',
'oc_a04->.00104*OC1+.00104*OC2+.012481*SOA;1.e9',
'bc_a01->6.5325E-01*CB1+6.5325E-01*CB2;0.413793103e9',
'bc_a02->2.5458E-01*CB1+2.5458E-01*CB2;0.413793103e9',
'bc_a03->4.0521E-03*CB1+4.0521E-03*CB2;0.413793103e9',
'bc_a04->1.6794E-06*CB1+1.6794E-06*CB2;0.413793103e9',
'so4_a01->2.4433e-2*SO4;3.31423e9',
'so4_a02->.47111*SO4;3.31423e9',
'so4_a03->.47857*SO4;3.31423e9',
'so4_a04->2.5727e-2*SO4;3.31423e9',
'no3_a01->2.4433e-2*NH4NO3;2.76186e9',
'no3_a02->.47111*NH4NO3;2.76186e9',
'no3_a03->.47857*NH4NO3;2.76186e9',
'no3_a04->2.5727e-2*NH4NO3;2.76186e9',
'nh4_a01->2.4433e-2*NH4;6.22876e8',
'nh4_a02->.47111*NH4;6.22876e8',
'nh4_a03->.47857*NH4;6.22876e8',
'nh4_a04->2.5727e-2*NH4;6.22876e8',
'na_a01->.2*SA1;1.e9',
'na_a02->.4*SA1;1.e9',
'na_a03->.2*SA1+.8*SA2;1.e9',
'na_a04->.8*SA3;1.e9',
'cl_a01->.3*SA1;1.e9',
'cl_a02->.6*SA1;1.e9',
'cl_a03->.3*SA1+1.2*SA2;1.e9',
'cl_a04->1.2*SA3;1.e9',
'oin_a01->0.5*[DUST1];1.e9',
'oin_a02->0.5*[DUST1]+0.24*[DUST2];1.e9',
'oin_a03->0.76*[DUST2]+1.0*[DUST3];1.e9',
'oin_a04->1.0*[DUST4];1.e9',
'num_a01-
>1.73e17*OC1+1.73e17*OC2+1.73e17*SOA+5.64e17*CB1+5.64e17*CB2+7.67e16*SO4+6.39
e16*NH4NO3+1.44e16*NH4+3.22e17*SA1+4.83e17*SA1+0.5*3.93e17*[DUST1];1.0',
'num_a02-
>1.71e16*OC1+1.71e16*OC2+1.71e16*SOA+9.91e15*CB1+9.91e15*CB2+6.06e16*SO4+5.05
e16*NH4NO3+1.14e16*NH4+2.90e16*SA1+4.36e16*SA1+0.5*1.17e16*[DUST1]+0.24*1.17e16
*[DUST2];1.0',
'num_a03-
>1.13e14*OC1+1.13e14*OC2+1.13e14*SOA+1.80e12*CB1+1.80e12*CB2+2.44e15*SO4+2.03
e15*NH4NO3+4.69e14*NH4+4.12e14*SA1+1.64e15*SA2+0.76*9.55e13*[DUST2]+9.55e13*[D
UST3];1.0',
'num_a04-
>4.09e10*OC1+4.09e10*OC2+4.09e10*SOA+3.42e07*CB1+3.42e07*CB2+3.25e12*SO4+2.71
e12*NH4NO3+6.12e11*NH4+7.60e13*SA3+1.49e12*[DUST4];1.0',

/

Please suggest to me,
Which one is best or why?

With regards
Ummed Singh

--
Rajesh Kumar, Project Scientist III

Research Applications Laboratory
National Center for Atmospheric Research
E-mail: rku...@ucar.edu
Phone: +1-303-497-1413
Web: https://staff.ucar.edu/users/rkumar

--
Rajesh Kumar, Project Scientist III

Research Applications Laboratory
National Center for Atmospheric Research
E-mail: rku...@ucar.edu
Phone: +1-303-497-1413
Web: https://staff.ucar.edu/users/rkumar

--
Rajesh Kumar, Project Scientist III

Research Applications Laboratory
National Center for Atmospheric Research
E-mail: rku...@ucar.edu
Phone: +1-303-497-1413
Web: https://staff.ucar.edu/users/rkumar

--
Rajesh Kumar, Project Scientist III

Research Applications Laboratory
National Center for Atmospheric Research
E-mail: rku...@ucar.edu
Phone: +1-303-497-1413
Web: https://staff.ucar.edu/users/rkumar

--
Rajesh Kumar, Project Scientist III

Research Applications Laboratory
National Center for Atmospheric Research
E-mail: rku...@ucar.edu
Phone: +1-303-497-1413
Web: https://staff.ucar.edu/users/rkumar

--
Rajesh Kumar, Project Scientist III

Research Applications Laboratory
National Center for Atmospheric Research
E-mail: rku...@ucar.edu
Phone: +1-303-497-1413
Web: https://staff.ucar.edu/users/rkumar

Rajesh Kumar, Project Scientist III

Research Applications Laboratory

National Center for Atmospheric Research

E-mail: rku...@ucar.edu

Phone: +1-303-497-1413

Web: https://staff.ucar.edu/users/rkumar

Sherin Hassan Bran

unread,

Feb 1, 2022, 3:57:53 AM2/1/22

to wrf-chem-fire_emiss, rkumar, umme...@gmail.com

Hi Rajesh,

Just a follow up question regarding the usage of FINNv2.4.

1) Species mapping which you shown above, BC should be subtracted from PM2.5 or it should be added ..??

2) As it in the website, FINNv2.4 were created for the MOZART T1 emissions, so does it work with MOZCART (chem_opt=112) or it only work with MOZCART T1 ..?? and if yes how about biomass_burn_opt (2 or 4) ..??