RACMSORG.inp for CAM-Chem

[Gwang-Jin Lee]

Hello,

I am trying to apply CAM-Chem output to initial and boundary value of WRF-Chem. I'm newbie of mozbc.

And I'm using RACM + MADE/SORGAM mechanism. And I met error messages with camchem-20220126185317231035.nc file when I ran mozbc:

chk_moz_vars: could not find SOA                              in  
 /scratch/e1329a01/WRF/JOBS/Cheongju28/emiss/camchem-20220126185317231035.nc
 NetCDF: Variable not found    
 
 fail to process netCDF file...
 in module_mozart_lib ... 

And I don't know how to map orgalk1i, orgalk1j, orgaro2i, orgaro2j, orgaro1i, orgaro1j, so4ai, so4aj, no3ai, no3aj, nu0 etc. to CAM-Chem variables.

(I attached ncdump results of camchem-20220126185317231035 and original RACMSORG.inp)

Best regards,

Gwang-Jin

[Héctor Jorquera González]

Hi Gwan

I am using this mapping, derived from two PDF docs:

spc_map = 'co -> CO',
          'o3 -> O3',
          'so2 -> SO2',
          'no -> NO',
          'no2 -> NO2',
          'h2o2 -> H2O2',
          'no3 -> NO3',
          'n2o5 -> N2O5',
          'hno3 -> HNO3',
          'nh3 -> NH3',
          'hcho -> CH2O',
          'pan -> PAN',
          'tpan -> 0.5*PAN',
          'tol -> TOLUENE',
          'xyl -> XYLENES',
          'hno4 -> HO2NO2',
          'eth -> C2H6',
          'hc3 -> C3H8',
          'hc5 -> BIGALK',
          'iso -> ISOP',
          'ald -> CH3CHO+MACR',
          'ket -> CH3COCH3+HYAC',
          'gly -> GLYOXAL',
          'mgly -> CH3COCHO',
          'csl -> CRESOL',
          'ol2 -> C2H4',
          'olt -> C3H6+MVK+.5*BIGENE',
          'oli -> .5*BIGENE',
          'onit -> ONITR',
          'ora2 -> CH3COOH',
          'ora1 -> HCOOH',
          'op1 -> CH3OOH',
          'sulf -> so4_a1+so4_a2+so4_a3;0.30173e6',
          'orgpai->0.1*pom_a4+0.1*soa1_a2+0.1*soa2_a2+0.1*soa3_a2+0.1*soa4_a2+0.1*soa5_a2+0.1*pom_a1+0.1*soa1_a1+0.1*soa2_a1+0.1*soa3_a1+0.1*soa4_a1+0.1*soa5_a1;0.414e9',
          'orgpaj->0.9*pom_a4+0.9*soa1_a2+0.9*soa2_a2+0.9*soa3_a2+0.9*soa4_a2+0.9*soa5_a2+0.9*pom_a1+0.9*soa1_a1+0.9*soa2_a1+0.9*soa3_a1+0.9*soa4_a1+0.9*soa5_a1;0.414e9',
          'eci->0.1*bc_a4+0.1*bc_a1;0.414e9',
          'ecj->0.9*bc_a4+0.9*bc_a1;0.414e9',
          'claj->0.6*ncl_a1;2e9',
          'naaj->0.4*ncl_a1;2e9',
          'seas->ncl_a1+ncl_a2+ncl_a3;2.e9',
          'so4ai->0.1*so4_a1+0.1*so4_a2+0.1*so4_a3;3.31e9',
          'so4aj->0.9*so4_a1+0.9*so4_a2+0.9*so4_a3;3.31e9',
          'no3ai->0.1*NH4;2.76e9',
          'no3aj->0.9*NH4;2.76e9'
          'orgalk1i->0.01*soa1_a2+0.01*soa2_a2+0.01*soa3_a2+0.01*soa4_a2+0.01*soa5_a2+0.01*soa1_a1+0.01*soa2_a1+0.01*soa3_a1+0.01*soa4_a1+0.01*soa5_a1;5.e9',
          'orgaro1i->0.01*soa1_a2+0.01*soa2_a2+0.01*soa3_a2+0.01*soa4_a2+0.01*soa5_a2+0.01*soa1_a1+0.01*soa2_a1+0.01*soa3_a1+0.01*soa4_a1+0.01*soa5_a1;5.e9',
          'orgaro2i->0.01*soa1_a2+0.01*soa2_a2+0.01*soa3_a2+0.01*soa4_a2+0.01*soa5_a2+0.01*soa1_a1+0.01*soa2_a1+0.01*soa3_a1+0.01*soa4_a1+0.01*soa5_a1;5.e9',
          'orgole1i->0.01*soa1_a2+0.01*soa2_a2+0.01*soa3_a2+0.01*soa4_a2+0.01*soa5_a2+0.01*soa1_a1+0.01*soa2_a1+0.01*soa3_a1+0.01*soa4_a1+0.01*soa5_a1;5.e9',
          'orgba1i->0.06*soa1_a2+0.06*soa2_a2+0.06*soa3_a2+0.06*soa4_a2+0.06*soa5_a2+0.06*soa1_a1+0.06*soa2_a1+0.06*soa3_a1+0.06*soa4_a1+0.06*soa5_a1;5.e9',
          'orgalk1j->0.09*soa1_a2+0.09*soa2_a2+0.09*soa3_a2+0.09*soa4_a2+0.09*soa5_a2+0.09*soa1_a1+0.09*soa2_a1+0.09*soa3_a1+0.09*soa4_a1+0.09*soa5_a1;5.e9',
          'orgaro1j->0.09*soa1_a2+0.09*soa2_a2+0.09*soa3_a2+0.09*soa4_a2+0.09*soa5_a2+0.09*soa1_a1+0.09*soa2_a1+0.09*soa3_a1+0.09*soa4_a1+0.09*soa5_a1;5.e9',
          'orgaro2j->0.09*soa1_a2+0.09*soa2_a2+0.09*soa3_a2+0.09*soa4_a2+0.09*soa5_a2+0.09*soa1_a1+0.09*soa2_a1+0.09*soa3_a1+0.09*soa4_a1+0.09*soa5_a1;5.e9',
          'orgole1j->0.09*soa1_a2+0.09*soa2_a2+0.09*soa3_a2+0.09*soa4_a2+0.09*soa5_a2+0.09*soa1_a1+0.09*soa2_a1+0.09*soa3_a1+0.09*soa4_a1+0.09*soa5_a1;5.e9',
          'orgba1j->0.54*soa1_a2+0.54*soa2_a2+0.54*soa3_a2+0.54*soa4_a2+0.54*soa5_a2+0.54*soa1_a1+0.54*soa2_a1+0.54*soa3_a1+0.54*soa4_a1+0.54*soa5_a1;5.e9',
          'p25i->0.02*dst_a3;1.17e9',
          'p25j->0.93*dst_a3;1.17e9',
          'soila->0.05*dst_a3;1.17e9',
          'antha->0.041*pom_a4+2.041*soa1_a2+2.041*soa2_a2+2.041*soa3_a2+2.041*soa4_a2+2.041*soa5_a2+0.041*pom_a1+2.041*soa1_a1+2.041*soa2_a1+2.041*soa3_a1+2.041*soa4_a1+2.041*soa5_a1;1.e9',
          'nu0->3.8e19*bc_a4+3.8e19*bc_a1+3.3e18*pom_a4+3.3e18*soa1_a2+3.3e18*pom_a1+3.3e18*soa2_a1+9.0e17*so4_a1+9.0e17*so4_a2+9.0e17*so4_a3+7.5e17*NH4+4e19*soa1_a1+4e19*soa1_a2;1.e-2'

Regards

Héctor

[Gwang-Jin Lee]

Dear Héctor,

I met some error message:

--------------------------------------

checking wrf variable antha
 len pom_a4 =            6    
 len soa1_a2 =            7    
 len soa2_a2 =            7    
 len soa3_a2 =            7    
 len soa4_a2 =            7    
 len soa5_a2 =            7    
 len pom_a1 =            6    
 len soa1_a1 =            7    
 len soa2_a1 =            7    
 len soa3_a1 =            7    
 len soa4_a1 =            7    
 len 2.041 =            5    
 chk_moz_vars: could not find 2.041                            in  

 /scratch/e1329a01/WRF/JOBS/Cheongju28/emiss/camchem-20220126185317231035.nc
 NetCDF: Variable not found    
 
 fail to process netCDF file...
 in module_mozart_lib ...

--------------------------------------

How should I modify mo_mozart_lib.f90 to solve this problem?

Thank you.

Best regards,

Gwang-Jin

2022년 1월 29일 토요일 오전 9시 7분 2초 UTC+9에 jorq...@ing.puc.cl님이 작성:

[Gwang-Jin Lee]

Dear Héctor,

I've solved the problem by modifiying the length of spc_map in main_bc_wrfchem.f90 to 188 

and the length mstring in mo_utils.f90:to 188.

Thank you. It works.

Best regards,

Gwang-Jin 

2022년 2월 3일 목요일 오후 2시 20분 11초 UTC+9에 Gwang-Jin Lee님이 작성:

[Ronald Opio]

Hi  Gwang-Jin. 

I wish to do something similar to yours. I kindly ask that please share with me your final file that you used as mozbc input.

With regards,

Ronald.