-
Notifications
You must be signed in to change notification settings - Fork 0
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Workplan for hooking up SUBDDs #64
Conversation
1636eb0
to
0a67bb6
Compare
#ifdef TRACERS_ON | ||
select case (iorw) | ||
case (ioread) | ||
call tracerIO(fid, 'read_dist') | ||
case (iowrite) | ||
call tracerIO(fid, 'write_dist') | ||
end select | ||
#endif | ||
#ifdef TRACERS_GC | ||
select case (iorw) | ||
case (ioread) | ||
call IO_CHEM(fid, 'read_dist') | ||
case (iowrite) | ||
call IO_CHEM(fid, 'write_dist') | ||
end select | ||
#endif | ||
#ifdef TRACERS_GC | ||
#elif TRACERS_ON | ||
select case (iorw) | ||
case (ioread) | ||
call IO_CHEM(fid, 'read_dist') | ||
call tracerIO(fid, 'read_dist') | ||
case (iowrite) | ||
call IO_CHEM(fid, 'write_dist') | ||
call tracerIO(fid, 'write_dist') | ||
end select | ||
#endif | ||
|
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
This looked to me like a copy-paste error (the same code appears twice).
call inc_subdd(subdd,k,sddarr2d) | ||
! NOTE: OH index is obtained using Ind_ function (from State_Chm_Mod) with | ||
! key 'OH' on L1426 of GEOS-Chem/GeosCore/diagnostics_mod.F90. | ||
case ('OH') | ||
call inc_subdd(subdd,k,ptropo) ! TODO: Output correct field, not ptropo | ||
! NOTE: NO index is obtained using Ind_ function (from State_Chm_Mod) with | ||
! key 'NO' on L282 of GEOS-Chem/GeosCore/hco_interface_gc_mod.F90. | ||
! This index is used on L2653 to access an entry of | ||
! State_Chm%Species. | ||
case ('NO') | ||
call inc_subdd(subdd,k,ptropo) ! TODO: Output correct field, not ptropo | ||
! NOTE: O3 appears in State_Diag with index 1 | ||
! (see L18820 of GEOS-Chem/Headers/state_diag_mod.F90) | ||
! NOTE: O3 index is obtained using Ind_ function (from State_Chm_Mod) with | ||
! key 'O3' on L284 of GEOS-Chem/GeosCore/hco_interface_gc_mod.F90. | ||
! This index is used on L2601 to access an entry of | ||
! State_Chm%Species. | ||
case ('O3') | ||
call inc_subdd(subdd,k,ptropo) ! TODO: Output correct field, not ptropo | ||
! NOTE: NO2 index is obtained using Ind_ function (from State_Chm_Mod) | ||
! with key 'NO2' on L283 of GEOS-Chem/GeosCore/hco_interface_gc_mod.F90. | ||
! This index is used on L2627 to access an entry of | ||
! State_Chm%Species. | ||
case ('NO2') | ||
call inc_subdd(subdd,k,ptropo) ! TODO: Output correct field, not ptropo | ||
! NOTE: CO index is obtained using Ind_ function (from State_Chm_Mod) with | ||
! key 'CO' on L936 of GEOS-Chem/GeosCore/carbon_gases_mod.F90. | ||
! This index is used on L219 to access an entry of | ||
! State_Chm%Species. | ||
case ('CO') | ||
call inc_subdd(subdd,k,ptropo) ! TODO: Output correct field, not ptropo | ||
! NOTE: CH4 index is obtained using Ind_ function (from State_Chm_Mod) | ||
! with key 'CH4' on L540 of GEOS-Chem/GeosCore/global_ch4_mod.F90. | ||
! This index is used on L549 to access an entry of | ||
! State_Chm%SpcData. | ||
case ('CH4') | ||
call inc_subdd(subdd,k,ptropo) ! TODO: Output correct field, not ptropo |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I guess these should be in an #ifdef TRACERS_GC
block?
We discussed in today's meeting that this is done more neatly in Lee's implementation in |
Closes #33.
Includes changes from #63 so that will need to be merged first.