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.

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

Add capability to run OAE experiments with an OmegaA threshold #291

Merged
merged 1 commit into from
Oct 20, 2023
Merged

Add capability to run OAE experiments with an OmegaA threshold #291

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
10 changes: 10 additions & 0 deletions cime_config/namelist_definition_blom.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -3836,6 +3836,16 @@
<desc>End year for ramp up in 'ramp' scenario</desc>
</entry>

<entry id="thrh_omegaa">
<type>real</type>
<category>bgcoafx</category>
<group>bgcoafx</group>
<values>
<value>-1.0</value>
</values>
<desc>OmegaA threshold for application of OA, negative = no threshold</desc>
</entry>

<!-- ========================= -->
<!-- namelist group: bgcparams -->
<!-- ========================= -->
Expand Down
4 changes: 3 additions & 1 deletion hamocc/mo_apply_oafx.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -70,8 +70,9 @@ subroutine apply_oafx(kpie,kpje,kpke,pddpo,omask,oafx)
!
!******************************************************************************
use mo_control_bgc, only: dtb,do_oalk
use mo_carbch, only: ocetra,oalkflx
use mo_param1_bgc, only: ialkali
use mo_carbch, only: ocetra,oalkflx,OmegaA
use mo_read_oafx, only: thrh_omegaa

implicit none

Expand All @@ -93,6 +94,7 @@ subroutine apply_oafx(kpie,kpje,kpke,pddpo,omask,oafx)
do j=1,kpje
do i=1,kpie
if (omask(i,j).gt.0.5) then
if (thrh_omegaa > 0.0 .and. OmegaA(i,j,1) > thrh_omegaa) cycle
oalkflx(i,j) = oafx(i,j)*dtb/365.
ocetra(i,j,1,ialkali)=ocetra(i,j,1,ialkali)+oalkflx(i,j)/pddpo(i,j,1)
endif
Expand Down
55 changes: 34 additions & 21 deletions hamocc/mo_read_oafx.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -74,14 +74,14 @@ module mo_read_oafx
implicit none

private
public :: ini_read_oafx,get_oafx,oalkscen,oalkfile
public :: ini_read_oafx,get_oafx,oalkscen,oalkfile,thrh_omegaa

character(len=128), save :: oalkscen =''
character(len=512), save :: oalkfile =''
real,allocatable, save :: oalkflx(:,:)
integer, save :: startyear,endyear
character(len=128), protected :: oalkscen =''
character(len=512), protected :: oalkfile =''
real,allocatable, protected :: oalkflx(:,:)
integer, protected :: startyear,endyear

real, parameter :: Pmol2kmol = 1.0e12
real, parameter :: Pmol2kmol = 1.0e12

! Parameter used in the definition of alkalinization scenarios not based on
! an input file. The following scenarios are defined in this module:
Expand All @@ -92,19 +92,26 @@ module mo_read_oafx
! Pmol ALK/yr-1 from year ramp_start to year ramp_end between
! latitude cdrmip_latmin and latitude cdrmip_latmax
!
real, protected :: addalk = 0.135 ! Pmol alkalinity/yr added in the
! scenarios. Read from namelist file
! to overwrite default value.
real, protected :: cdrmip_latmax = 70.0 ! Min and max latitude where
real, protected :: cdrmip_latmin = -60.0 ! alkalinity is added according
! to the CDRMIP protocol. Read from
! namelist file to overwrite default
! value.
integer, protected :: ramp_start = 2025 ! In 'ramp' scenario, start at
integer, protected :: ramp_end = 2035 ! 0 Pmol/yr in ramp_start, and max
! addalk Pmol/yr in ramp_end.
! Read from namelist file to
! overwrite default value.
! Values are read from namelist bgcoafx, which overwrites default values set
! here
real, protected :: addalk = 0.135 ! Pmol alkalinity/yr added in the
! scenarios.
real, protected :: cdrmip_latmax = 70.0 ! Min and max latitude where
real, protected :: cdrmip_latmin = -60.0 ! alkalinity is added according
! to the CDRMIP protocol.
integer, protected :: ramp_start = 2025 ! In 'ramp' scenario, start at
integer, protected :: ramp_end = 2035 ! 0 Pmol/yr at ramp_start, and
! arrive at addalk Pmol/yr in
! year ramp_end

! Parameter used for ALL alkalinization scenarios, read through namelist
! namelist bgcoafx, which overwrites default values set here
real, protected :: thrh_omegaa =-1.0 ! Limit the input of alkalinity by
! setting alkalinity-flux to zero
! for grid cells where Omegaa >
! thrh_omegaa (negative values mean
! no threshold considered)


logical, save :: lini = .false.

Expand Down Expand Up @@ -136,9 +143,9 @@ subroutine ini_read_oafx(kpie,kpje,pdlxp,pdlyp,pglat,omask)
!
!******************************************************************************
use mod_xc, only: xcsum,xchalt,mnproc,nbdy,ips
use mo_control_bgc, only: io_stdo_bgc,do_oalk,bgc_namelist,get_bgc_namelist
use mod_dia, only: iotype
use mod_nctools, only: ncfopn,ncgeti,ncfcls
use mo_control_bgc, only: io_stdo_bgc,do_oalk,bgc_namelist,get_bgc_namelist

implicit none

Expand All @@ -154,7 +161,7 @@ subroutine ini_read_oafx(kpie,kpje,pdlxp,pdlyp,pglat,omask)
real :: ztmp1(1-nbdy:kpie+nbdy,1-nbdy:kpje+nbdy)

namelist /bgcoafx/ oalkscen,oalkfile,addalk,cdrmip_latmax,cdrmip_latmin, &
& ramp_start,ramp_end
& ramp_start,ramp_end,thrh_omegaa

! Read parameters for alkalinization fluxes from namelist file
if(.not. allocated(bgc_namelist)) call get_bgc_namelist
Expand Down Expand Up @@ -248,6 +255,12 @@ subroutine ini_read_oafx(kpie,kpje,pdlxp,pdlyp,pglat,omask)
write(io_stdo_bgc,*)' ramping-up from ', ramp_start, ' to ', ramp_end
endif
endif

if(mnproc.eq.1 .and. thrh_omegaa > 0.0) then
write(io_stdo_bgc,*)' alkalinity flux will be limited by a threshold for Omega_a of ',thrh_omegaa
write(io_stdo_bgc,*)' '
endif


do j=1,kpje
do i=1,kpie
Expand Down
Loading