totalsmod.f

* © 2023. Triad National Security, LLC. All rights reserved.
* This program was produced under U.S. Government contract 89233218CNA000001 for Los Alamos National
* Laboratory (LANL), which is operated by Triad National Security, LLC for the U.S. Department of
* Energy/National Nuclear Security Administration. All rights in the program are reserved by Triad
* National Security, LLC, and the U.S. Department of Energy/National Nuclear Security Administration.
* The Government is granted for itself and others acting on its behalf a nonexclusive, paid-up,
* irrevocable worldwide license in this material to reproduce, prepare. derivative works, distribute
* copies to the public, perform publicly and display publicly, and to permit others to do so.
*This file is part of SuperNu.  SuperNu is released under the terms of the GNU GPLv3, see COPYING.
*Copyright (c) 2013-2022 Ryan T. Wollaeger and Daniel R. van Rossum.  All rights reserved.
      module totalsmod
c     ----------------
      implicit none
************************************************************************
* Miscellaneous domain-integrated quantities.
************************************************************************
c-- energy conservation check quantities
      real*8 :: tot_eext = 0d0 !time-integrated input energy from external source
      real*8 :: tot_emat = 0d0 !material energy
      real*8 :: tot_erad = 0d0 !census radiation energy
      real*8 :: tot_eout = 0d0 !energy escaped
      real*8 :: tot_evelo = 0d0 !total energy change to rad field from fluid
      real*8 :: tot_eerror = 0d0 !error in integral problem energy
      real*8 :: tot_esurf = 0d0
c-- initial external energy
      real*8 :: tot_eext0 = 0d0
c
c-- sources (domain integrated, not time integrated)
c-- gas_ sources (need to be collected)
      real*8 :: tot_sthermal = 0d0
      real*8 :: tot_smanufac = 0d0
      real*8 :: tot_sdecaygamma = 0d0
      real*8 :: tot_sdecaybeta = 0d0
c-- grd_ sources
      real*8 :: tot_sdeposgamma = 0d0
      real*8 :: tot_samp = 0d0
      real*8 :: tot_sanalvol = 0d0
      real*8 :: tot_sanalsurf = 0d0
c-- flx_ sources (negative)
      real*8 :: tot_sfluxgamma = 0d0
      real*8 :: tot_sflux = 0d0
c
      save
c
      contains
c
      subroutine totals_startup
c     -------------------------!{{{
      implicit none
************************************************************************
* The initial radiation field is estimated from the assumption
* that radiative work losses and energy deposition in the radiation
* field each recieve half of the total energy deposition.
* Once the radiation field is determined from this assumption, these
* values are adopted in the totals terms.
************************************************************************
c-- no radiation is assumed to have escaped
      tot_eout = 0d0
c-- energy decayed prior to first timestep
      tot_eext = tot_eext0
c-- exact balance
      tot_evelo = tot_eext-tot_erad-tot_emat !tot_emat is added to tot_eext in the first time step in gas_update
c!}}}
      end subroutine totals_startup
c
c
      subroutine totals_error
c     -----------------------!{{{
      implicit none
************************************************************************
* Check that all particle energy (weight) is accounted for from
* conservation in comoving quantities.
************************************************************************
      tot_eerror = (tot_eext-tot_evelo-tot_eout-tot_erad-tot_emat)/
     &  tot_eext
c!}}}
      end subroutine totals_error
c
      end module totalsmod
c vim: fdm=marker