Skip to content

part of cam6_4_214: Write lat_rad, lon_rad to snapshot for CCPPization (CAM-SIMA) #1731

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.

Sign up for GitHub

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

Merged
1 commit merged into from
Oct 21, 2026
Merged

part of cam6_4_214: Write lat_rad, lon_rad to snapshot for CCPPization (CAM-SIMA) #1731

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
64 changes: 64 additions & 0 deletions src/control/cam_snapshot_common.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -28,6 +28,7 @@ module cam_snapshot_common
public :: cam_snapshot_ptend_outfld
public :: snapshot_type
public :: cam_state_snapshot_init
public :: cam_grid_snapshot_init
public :: cam_cnst_snapshot_init
public :: cam_tend_snapshot_init
public :: cam_ptend_snapshot_init
Expand All @@ -38,6 +39,7 @@ module cam_snapshot_common

private :: snapshot_addfld_nd
private :: state_snapshot_all_outfld
private :: grid_snapshot_all_outfld
private :: cnst_snapshot_all_outfld
private :: tend_snapshot_all_outfld
private :: cam_in_snapshot_all_outfld
Expand Down Expand Up @@ -72,6 +74,7 @@ module cam_snapshot_common
end type pbuf_info_type

integer :: nstate_var
integer :: ngrid_var
integer :: ncnst_var
integer :: ntend_var
integer :: ncam_in_var
Expand All @@ -82,6 +85,7 @@ module cam_snapshot_common

! Note the maximum number of variables for each type
type (snapshot_type) :: state_snapshot(30)
type (snapshot_type) :: grid_snapshot(2)
type (snapshot_type) :: cnst_snapshot(pcnst)
type (snapshot_type) :: tend_snapshot(6)
type (snapshot_type) :: cam_in_snapshot(pcnst+35) ! needs to be bigger than pcnst because cam_in%cflx is split by constituent and cam_in%dstflx by bin.
Expand Down Expand Up @@ -120,6 +124,9 @@ subroutine cam_snapshot_all_outfld(file_num, state, tend, cam_in, cam_out, pbuf)
! Write out all the state fields
call state_snapshot_all_outfld(lchnk, file_num, state)

! Write out the physics grid coordinates in radians
call grid_snapshot_all_outfld(lchnk, file_num)

! Write out all the constituent fields
call cnst_snapshot_all_outfld(lchnk, file_num, state%q)

Expand Down Expand Up @@ -288,6 +295,28 @@ subroutine cam_state_snapshot_init(cam_snapshot_before_num_in, cam_snapshot_afte

end subroutine cam_state_snapshot_init

subroutine cam_grid_snapshot_init(cam_snapshot_before_num, cam_snapshot_after_num)

!--------------------------------------------------------
! This subroutine does the addfld calls for the physics grid coordinates in radians.
! The history coordinates lat and lon are written in degrees, and converting them back to
! radians in CAM-SIMA differs from the values CAM physics used by one ulp at some columns,
! which seeds answer differences in every latitude-dependent scheme.
! The CAM-SIMA null dycore prefers these fields when they are present on the snapshot.
!--------------------------------------------------------

integer,intent(in) :: cam_snapshot_before_num, cam_snapshot_after_num

ngrid_var = 0

call snapshot_addfld( ngrid_var, grid_snapshot, cam_snapshot_before_num, cam_snapshot_after_num, &
'grid%lat_rad', 'lat_rad', 'radians', horiz_only)

call snapshot_addfld( ngrid_var, grid_snapshot, cam_snapshot_before_num, cam_snapshot_after_num, &
'grid%lon_rad', 'lon_rad', 'radians', horiz_only)

end subroutine cam_grid_snapshot_init

subroutine cam_cnst_snapshot_init(cam_snapshot_before_num, cam_snapshot_after_num)

!--------------------------------------------------------
Expand Down Expand Up @@ -787,6 +816,41 @@ subroutine snapshot_addfld(nddt_var, ddt_snapshot, cam_snapshot_before_num, cam_

end subroutine snapshot_addfld

subroutine grid_snapshot_all_outfld(lchnk, file_num)

use phys_grid, only: get_rlat_all_p, get_rlon_all_p

integer, intent(in) :: lchnk
integer, intent(in) :: file_num

real(r8) :: rcoord(pcols)
integer :: i

do i=1, ngrid_var

! Turn on the writing for only the requested tape (file_num)
call cam_history_snapshot_activate(trim(grid_snapshot(i)%standard_name), file_num)

! Select the coordinate which is being written
select case(grid_snapshot(i)%ddt_string)

case ('grid%lat_rad')
call get_rlat_all_p(lchnk, pcols, rcoord)

case ('grid%lon_rad')
call get_rlon_all_p(lchnk, pcols, rcoord)

case default
call endrun('ERROR in grid_snapshot_all_outfld: no case for '//trim(grid_snapshot(i)%ddt_string))

end select

call outfld(grid_snapshot(i)%standard_name, rcoord, pcols, lchnk)

end do

end subroutine grid_snapshot_all_outfld

subroutine state_snapshot_all_outfld(lchnk, file_num, state)

use physics_types, only: phys_te_idx, dyn_te_idx
Expand Down
2 changes: 2 additions & 0 deletions src/physics/cam/cam_snapshot.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -20,6 +20,7 @@ module cam_snapshot
use cam_logfile, only: iulog
use cam_snapshot_common, only: snapshot_type, cam_snapshot_deactivate, cam_snapshot_all_outfld, cam_snapshot_ptend_outfld
use cam_snapshot_common, only: snapshot_type, cam_state_snapshot_init, cam_cnst_snapshot_init, cam_tend_snapshot_init
use cam_snapshot_common, only: cam_grid_snapshot_init
use cam_snapshot_common, only: cam_ptend_snapshot_init, cam_in_snapshot_init, cam_out_snapshot_init
use cam_snapshot_common, only: cam_pbuf_snapshot_init, snapshot_addfld

Expand Down Expand Up @@ -64,6 +65,7 @@ subroutine cam_snapshot_init(cam_in_arr, cam_out_arr, pbuf, index)
if (cam_snapshot_before_num <= 0 .and. cam_snapshot_after_num <= 0) return ! No snapshot files are being requested

call cam_state_snapshot_init(cam_snapshot_before_num, cam_snapshot_after_num)
call cam_grid_snapshot_init(cam_snapshot_before_num, cam_snapshot_after_num)
call cam_cnst_snapshot_init(cam_snapshot_before_num, cam_snapshot_after_num)
call cam_tend_snapshot_init(cam_snapshot_before_num, cam_snapshot_after_num)
call cam_ptend_snapshot_init(cam_snapshot_after_num)
Expand Down
2 changes: 2 additions & 0 deletions src/physics/cam7/cam_snapshot.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -20,6 +20,7 @@ module cam_snapshot
use cam_logfile, only: iulog
use cam_snapshot_common, only: snapshot_type, cam_snapshot_deactivate, cam_snapshot_all_outfld, cam_snapshot_ptend_outfld
use cam_snapshot_common, only: snapshot_type, cam_state_snapshot_init, cam_cnst_snapshot_init, cam_tend_snapshot_init
use cam_snapshot_common, only: cam_grid_snapshot_init
use cam_snapshot_common, only: cam_ptend_snapshot_init, cam_in_snapshot_init, cam_out_snapshot_init
use cam_snapshot_common, only: cam_pbuf_snapshot_init, snapshot_addfld

Expand Down Expand Up @@ -64,6 +65,7 @@ subroutine cam_snapshot_init(cam_in_arr, cam_out_arr, pbuf, index)
if (cam_snapshot_before_num <= 0 .and. cam_snapshot_after_num <= 0) return ! No snapshot files are being requested

call cam_state_snapshot_init(cam_snapshot_before_num, cam_snapshot_after_num)
call cam_grid_snapshot_init(cam_snapshot_before_num, cam_snapshot_after_num)
call cam_cnst_snapshot_init(cam_snapshot_before_num, cam_snapshot_after_num)
call cam_tend_snapshot_init(cam_snapshot_before_num, cam_snapshot_after_num)
call cam_ptend_snapshot_init(cam_snapshot_after_num)
Expand Down
Loading