Error while runing CAM-chem model output for updating chemical boundary condition using mozbc

[Ganesan M]

I am trying run wrf chem model with RADM2/MADE/SORGAM chemitry and aerosol speciation (chem_opt = 11). Since the MOZART data has been replaced with CAM-chem output for recent years (2018), I tried updating the chemical boundary condition using CAM-chem model data in the mozbc utiliy, but I am getting errors and the mozbc is crashing without completing the process. Is there any changes i need to make to run the mozbc utility for RADM2/MADE/SORGAM chemistry option using CAM-chem model data?. Is there another way for me to update the chemical boundary condition as I specifically need to run the model for RADM2 chemistry option. I am attaching the error log and inp files for your reference, kindly help me with this issue.

[Stacy Walters]

Ganesan,

Would you please issue the following command and return the hdr.out

file to me as an email attachment?

ncdump -h hb0068 > hdr.out

Stacy

[Arthur Mizzi]

All/Stacy:

I download mozbc from the ACOM site yesterday.  I've built it on the NASA computer Pleiades and am using it to create IC/BCs from cam-chem output also downloaded from the CAM-Chem download site yesterday.  I'm experiencing the same suffix error as this user.  How was this problem solved?

Thanks,

Arthur

[Arthur Mizzi]

Never mind, I found the nameless setting that solves the problem.

[Ankan Sarkar]

Hello all,

I also want to run WRF-Chem using the RADM2-MADE/SORGAM chemistry option (chem_opt=11), just like Ganesh. I ran WRF-Chem successfully using MOZART4 data to include chemical IC/BC in the wrfbdy and wrfinput files. But, as he mentioned, the MOZART4 data is not available after January 21, 2018.  So, to simulate the recent periods, I have to use the CAM-Chem model output. I downloaded the CAM-Chem data, and I noticed that some variables are different as compared to MOZART4 dataset. In the RADM2SORG.inp file (which is provided with the Mozbc utility), the following mapping is given for aerosol species:

          'orgpai->2.03507e-2*OC1+2.03507e-2*OC2;1.e9',
          'orgpaj->0.37992*OC1+0.37992*OC2;1.e9',
          'antha->0.011199*OC1+0.011199*OC2+0.1344*SOA;1.e9',
          'eci->CB1+CB2;0.113571e9',
          'ecj->CB1+CB2;0.264125e9',
          'claj->1.2*SA1;1e9',
          'naaj->0.8*SA1;1e9',
          'seas->SA1+SA2+SA3;1.e9',
          'orgalk1i->.244208*SOA;1.e9',
          'orgalk1j->3.9595*SOA;1.e9',
          'orgaro2i->.244208*SOA;1.e9',
          'orgaro2j->3.9595*SOA;1.e9',
          'orgaro1i->0.248134*SOA;1.e9',
          'orgaro1j->4.669487*SOA;1.e9'
          'so4ai->1.5563e-3*SO4;3.31423e9',
          'so4aj->0.832*SO4;3.31423e9',
          'no3ai->1.5563e-3*NH4NO3;2.76186e9',
          'no3aj->0.832*NH4NO3;2.76186e9'
          'p25i->4.66064*[DUST1];1.e9',
          'p25j->4.66064*[DUST2];1.e9',
          'soila->4.66064*[DUST3];1.e9',
          'nu0->4.411755e14*OC1+4.411755e14*OC2+5.293920e15*SOA+4.3853645e15*CB1+4.3853645e15*CB2+4.3656995e15*SO4+1.034506e14*NH4NO3+3.954300e13*SA1+7.796950e13*[DUST1];2.387325e2'

But I checked the variables available in the CAM-Chem dataset, and there are no variables like OC1, OC2, SOA, CB1, CB2, SA1, SA2, SA3, SO4, NH4NO3, DUST1, DUST2, DUST3, etc. For your convenience, I attached the hdr.out file (generated after giving the command: 'ncdump -h').

Some similar variables, like dst_a1, dst_a2, and dst_a3 are available with CAM-Chem data. So, should I replace DUST1, DUST2, and 'DUST3 with these variables? For BC, there are two variables I can see: bc_a1 and bc_a4. So, are they the same as CB1 and CB2? Is any sample mapping for aerosol species available for the RADM2-MADE/SORGAM chemistry option while using the CAM-Chem model output for chemical IC/BC? Can anyone please provide me with the sample mapping?  That will be very helpful for me. Thank you for your time and consideration.

With regards,

Ankan