Hi everyone
Recently, I need to use MOSAIC scheme in CESM to finish my final project of one course. I have used CESM2.2.2 but I have never used chemsitry related models so I have many questions:
1. I do not know if I couple this scheme correctly.
At the beginning, I follow Use MOSAIC scheme in CESM/CAM-chem - CAM-Chem - wiki.ucar.edu to git clone MOSAIC (git clone GitHub - ESCOMP/CAM: Community Atmosphere Model -b zz_cam6_3_018_mosaic0.3 MOSAIC_CODE).
And then rename CESM2.2.2/components/cam to CESM2.2.2/components/cam_origin and copy mosaic_code to CESM2.2.2/components/ and rename it to cam.
Finally, I setup a FCHIST case and build
(I knew that mosaic in FHIST compset will not work. Because FCnudged needs to download so much data from MERRA2 so I use FCHIST to test.)
When running this case, I received an error in modal_aero_amicphys.F90:
5635 | if ( mwaer_from_advmass ) then
| 1
Error: IF clause at (1) requires a scalar LOGICAL expression
I review this code and I find in line 93: integer, public :: mwaer_from_advmass = .false. .
I modify integer to logical and this case could run but could not write cam output (cam.h0) successfully. It keeps running, the cam.h0 file keeps a fixed size and updates its status in real time (as indicated by modification time), yet no data has ever been written to it—only variable names.
Even though I use FCnudged to set up a new case, it is still so. In cesm.log, it says "imp_sol: time step 1800.000 failed to converge" and opened file cam.h0 and cam.r to write.
I want to know how I can debug this.
2. My final project is about decreasing NH3 to see how it impacts the SOA or ozone formation. I want to see the specific reactions from the source codes but after reading MZ272_TS1.1_simpleVBS_MOSAIC_20210528.in and chem_mech.doc, I do not find the information I want. I want to ask which files I should read to get the reactions related with NH3.
Looking forward to some suggestions and comments very much!
Recently, I need to use MOSAIC scheme in CESM to finish my final project of one course. I have used CESM2.2.2 but I have never used chemsitry related models so I have many questions:
1. I do not know if I couple this scheme correctly.
At the beginning, I follow Use MOSAIC scheme in CESM/CAM-chem - CAM-Chem - wiki.ucar.edu to git clone MOSAIC (git clone GitHub - ESCOMP/CAM: Community Atmosphere Model -b zz_cam6_3_018_mosaic0.3 MOSAIC_CODE).
And then rename CESM2.2.2/components/cam to CESM2.2.2/components/cam_origin and copy mosaic_code to CESM2.2.2/components/ and rename it to cam.
Finally, I setup a FCHIST case and build
(I knew that mosaic in FHIST compset will not work. Because FCnudged needs to download so much data from MERRA2 so I use FCHIST to test.)
When running this case, I received an error in modal_aero_amicphys.F90:
5635 | if ( mwaer_from_advmass ) then
| 1
Error: IF clause at (1) requires a scalar LOGICAL expression
I review this code and I find in line 93: integer, public :: mwaer_from_advmass = .false. .
I modify integer to logical and this case could run but could not write cam output (cam.h0) successfully. It keeps running, the cam.h0 file keeps a fixed size and updates its status in real time (as indicated by modification time), yet no data has ever been written to it—only variable names.
Even though I use FCnudged to set up a new case, it is still so. In cesm.log, it says "imp_sol: time step 1800.000 failed to converge" and opened file cam.h0 and cam.r to write.
I want to know how I can debug this.
2. My final project is about decreasing NH3 to see how it impacts the SOA or ozone formation. I want to see the specific reactions from the source codes but after reading MZ272_TS1.1_simpleVBS_MOSAIC_20210528.in and chem_mech.doc, I do not find the information I want. I want to ask which files I should read to get the reactions related with NH3.
Looking forward to some suggestions and comments very much!