Skip to content
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

Closed
wants to merge 5 commits into from
Closed

Conversation

jwallwork23
Copy link
Collaborator

@jwallwork23 jwallwork23 commented Dec 11, 2024

Closes #33.

Includes changes from #63 so that will need to be merged first.

  • Provide a workplan for where to hook up SUBDDs to GEOS-Chem diagnostics.
  • Turn SUBDDs back on in the rundeck.
  • Delete unused files that were added since E2.1.2.
  • Tidy up and fix some ifdefs.

@jwallwork23 jwallwork23 added the deliverable1 Associated with deliverable 1 label Dec 11, 2024
@jwallwork23 jwallwork23 self-assigned this Dec 11, 2024
Comment on lines -1019 to -1043
#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

Copy link
Collaborator Author

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).

@jwallwork23 jwallwork23 marked this pull request as ready for review December 20, 2024 15:32
@jwallwork23 jwallwork23 requested a review from ltmurray December 20, 2024 15:32
Comment on lines +1914 to +1950
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
Copy link
Collaborator Author

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?

@jwallwork23 jwallwork23 marked this pull request as draft December 20, 2024 16:54
@jwallwork23 jwallwork23 removed the request for review from ltmurray December 20, 2024 16:55
@jwallwork23
Copy link
Collaborator Author

We discussed in today's meeting that this is done more neatly in Lee's implementation in CHEM_DRV. (Search there for inc_subdd.) We should figure out why that isn't working rather than continuing with this approach.

@jwallwork23 jwallwork23 mentioned this pull request Dec 20, 2024
@jwallwork23 jwallwork23 added invalid This doesn't seem right bug Something isn't working and removed invalid This doesn't seem right labels Dec 20, 2024
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something isn't working deliverable1 Associated with deliverable 1
Projects
None yet
Development

Successfully merging this pull request may close these issues.

1 participant