Modify fullchem #2374
Comments
yantosca
added
the
topic: Chemical Mechanisms
|
Thanks for writing @yangning-code. Maybe the best thing to do would be to add a line to the geos-chem/KPP/fullchem/commonIncludeVars.H Lines 112 to 213 in dd859ee and then to copy the values of BCPO into that in the SUBROUTINE fullChem_SetStateHet( I, J, L, &
id_SALA, id_SALAAL, id_SALC, &
id_SALCAL, Input_Opt, State_Chm, &
State_Met, H, RC )
!
! !USES:
!
USE ErrCode_Mod
USE GcKpp_Global
USE GcKpp_Parameters
USE PhysConstants, ONLY : AVO, PI
USE Input_Opt_Mod, ONLY : OptInput
USE rateLawUtilFuncs
USE State_Chm_Mod, ONLY : ChmState, Ind_
USE State_Met_Mod, ONLY : MetState
USE Aerosol_Mod, ONLY : BCPO ! <=== add this lineand then later on in that function add a line to copy values from !========================================================================
! Populate fields of the HetState object in gckpp_Global
!========================================================================
! Identify a box for debug printout within rate-law functions
H%debugBox = .FALSE.
! Constants (so that we can use these within KPP)
H%AVO = AVO
H%PI = PI
! Meteorology-related quantities
H%CldFr = MIN(MAX(State_Met%CLDF(I,J,L), 0.0_dp), 1.0_dp)
H%ClearFr = 1.0_dp - H%CldFr
H%QICE = State_Met%QI(I,J,L)
H%QLIQ = State_Met%QL(I,J,L)
H%vAir = State_Met%AIRVOL(I,J,L) * 1.0e6_dp
! Aerosol fields
H%nAeroType = State_Chm%nAeroType
H%aClArea = State_Chm%aClArea(I,J,L)
H%aClRadi = State_Chm%aClRadi(I,J,L)
H%aClVol = H%aClArea * H%aClRadi / 3.0_dp
H%AWATER(:) = State_Chm%IsorropAeroH2O(I,J,L,:)
H%xArea(1:NA) = State_Chm%AeroArea(I,J,L,1:NA)
H%xRadi(1:NA) = State_Chm%AeroRadi(I,J,L,1:NA)
H%xVol(1:NA) = H%xArea(1:NA) * H%xRadi(1:NA) / 3.0_dp
H%wetArea(1:NA) = State_Chm%WetAeroArea(I,J,L,1:NA)
H%xH2O(1:NA) = State_Chm%AeroH2O(I,J,L,1:NA) * 1.0e-6_dp
H%OMOC_POA = State_Chm%OMOC_POA(I,J)
H%OMOC_OPOA = State_Chm%OMOC_OPOA(I,J)
H%BCPO = BCPO(I,J,L) ! <=== add this line(note, 'State_Het Then make sure to pass State-Het to your rate law function and then you can use the BCPO. NOTE: There have been a LOT of changes since 13.3. Again, I would recommend using a newer version if possible. |
|
Because of the model version issue, I did not find the commonIncludeVars.H you mentioned. However, the fullchem.kpp in the version 13.3.2 I used has the same content as commonIncludeVars.H, so I added REAL(dp) :: BCPO Finally, I recompiled KPP, and then an error occurred when compiling the code again |
|
Thanks @yangning-code. Yes, the It looks like the code can't resolve the aerosol_mod.F90. The other thing you can do is to pass BCPO as an argument into the SetStateHet routine instead of trying to get it from aerosol_mod. |
|
It's a bit difficult for me, can you give me an example? |
Your name
Ning Yang
Your affiliation
JNU
Please provide a clear and concise description of your question or discussion topic.
I am using GEOSChemV13.3.2. I want to call the concentration of BCPO in GeosCore/aerosol_mod.F90 in KPP/fullchem/fullchem_RateLawFuncs.F90 so that I can do a mechanism update, but I don’t know how to add it. Can you give me some help?
The text was updated successfully, but these errors were encountered: