Skip to content

Commit

Permalink
improves the readability of salt budget
Browse files Browse the repository at this point in the history
  • Loading branch information
jedwards4b committed Apr 28, 2023
1 parent d8ddff5 commit 962646a
Show file tree
Hide file tree
Showing 2 changed files with 6 additions and 4 deletions.
2 changes: 1 addition & 1 deletion cesm/driver/esm.F90
Original file line number Diff line number Diff line change
Expand Up @@ -796,7 +796,7 @@ subroutine esm_init_pelayout(driver, maxthreads, rc)
#ifndef NO_MPI2
use mpi , only : MPI_COMM_NULL, mpi_comm_size
#endif
use mct_mod , only : mct_world_init
use m_MCTWorld , only : mct_world_init => init

#ifdef MED_PRESENT
use med_internalstate_mod , only : med_id
Expand Down
8 changes: 5 additions & 3 deletions mediator/med_diag_mod.F90
Original file line number Diff line number Diff line change
Expand Up @@ -95,6 +95,8 @@ module med_diag_mod
character(*), parameter :: FA1 = "(' ',a12,6f15.8)"
character(*), parameter :: FA0r = "(' ',12x,8(6x,a8,1x))"
character(*), parameter :: FA1r = "(' ',a12,8f15.8)"
character(*), parameter :: FA0s = "(' ',12x,8(7x,a8,2x))"
character(*), parameter :: FA1s = "(' ',a12,8g18.8)"

! ---------------------------------
! C for component
Expand Down Expand Up @@ -2683,7 +2685,7 @@ subroutine med_diag_print_summary(data, ip, date, tod)
write(diagunit,*) ' '
write(diagunit,FAH) subname,'NET SALT BUDGET (kg/m2s): period = ',&
trim(budget_diags%periods(ip)%name), ': date = ',date,tod
write(diagunit,FA0r) ' atm',' lnd',' rof',' ocn',' ice nh',' ice sh',' glc',' *SUM* '
write(diagunit,FA0s) ' atm',' lnd',' rof',' ocn',' ice nh',' ice sh',' glc',' *SUM* '
do nf = f_salt_beg, f_salt_end
net_salt_atm = data(nf, c_atm_recv, ip) + data(nf, c_atm_send, ip)
net_salt_lnd = data(nf, c_lnd_recv, ip) + data(nf, c_lnd_send, ip)
Expand All @@ -2695,7 +2697,7 @@ subroutine med_diag_print_summary(data, ip, date, tod)
net_salt_tot = net_salt_atm + net_salt_lnd + net_salt_rof + net_salt_ocn + &
net_salt_ice_nh + net_salt_ice_sh + net_salt_glc

write(diagunit,FA1r) budget_diags%fields(nf)%name,&
write(diagunit,FA1s) budget_diags%fields(nf)%name,&
net_salt_atm, net_salt_lnd, net_salt_rof, net_salt_ocn, &
net_salt_ice_nh, net_salt_ice_sh, net_salt_glc, net_salt_tot
enddo
Expand All @@ -2718,7 +2720,7 @@ subroutine med_diag_print_summary(data, ip, date, tod)
sum_net_salt_tot = sum_net_salt_atm + sum_net_salt_lnd + sum_net_salt_rof + sum_net_salt_ocn + &
sum_net_salt_ice_nh + sum_net_salt_ice_sh + sum_net_salt_glc

write(diagunit,FA1r)' *SUM*',&
write(diagunit,FA1s)' *SUM*',&
sum_net_salt_atm, sum_net_salt_lnd, sum_net_salt_rof, sum_net_salt_ocn, &
sum_net_salt_ice_nh, sum_net_salt_ice_sh, sum_net_salt_glc, sum_net_salt_tot
end if
Expand Down

0 comments on commit 962646a

Please sign in to comment.