Skip to content

Commit

Permalink
input to bc_surface from medusa added and Sed_input deleted
Browse files Browse the repository at this point in the history
  • Loading branch information
@a270105
a270105 committed Dec 13, 2024
1 parent 434f0e6 commit 0864237
Show file tree
Hide file tree
Showing 4 changed files with 83 additions and 24 deletions.
19 changes: 0 additions & 19 deletions src/gen_surface_forcing.F90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -1554,13 +1554,7 @@ SUBROUTINE sbc_do(partit, mesh)
icount = (/1,1/)
ncdata = 0.d0

! total_runoff = sum(runoff*area(1,:))*86400
! Does 'area' only contain values on one node? sum of area not equal total
! ocean surface area!
total_runoff = 8.76d5*86400
! if (mype==0) write(*,*) mype, 'total runoff (m3/day):', total_runoff
! if (mype==0) write(*,*) mype, 'runoff:', maxval(runoff),minval(runoff)
! if (mype==0) write(*,*) mype, 'surface area',maxval(area(1,:)),minval(area(1,:))

status=nf_open(sedfilename, nf_nowrite, ncid)
if(status.ne.nf_noerr) call handle_err(status)
Expand Down Expand Up @@ -1632,19 +1626,6 @@ SUBROUTINE sbc_do(partit, mesh)
lb_flux(:,n_lb) = -runoff*ncdata(n_lb)/total_runoff*lb_tscale
end do

! if (mype==0) write(*,*) mype, 'sum of surface area (m2)',
! sum(area(1,:))
! if (mype==0) write(*,*) mype, 'total ocean area (m2)', ocean_area
! if (mype==0) write(*,*) mype, 'DSi concentration in rivers',ncdata(4)/total_runoff
! if (mype==0) write(*,*) mype, 'DIC concentration in rivers',ncdata(2)/total_runoff
! if (mype==0) write(*,*) mype, 'Alk concentration in rivers',ncdata(3)/total_runoff
! if (mype==0) write(*,*) mype, 'DIN concentration in rivers',ncdata(1)/total_runoff
! if (mype==0) write(*,*) mype, 'DIN surface input:',minval(lb_flux(:,1)),maxval(lb_flux(:,1))
! if (mype==0) write(*,*) mype, 'DIC surface input:',minval(lb_flux(:,2)),maxval(lb_flux(:,2))
! if (mype==0) write(*,*) mype, 'Alk surface input:',minval(lb_flux(:,3)),maxval(lb_flux(:,3))
! if (mype==0) write(*,*) mype, 'DSi surface input:',minval(lb_flux(:,4)),maxval(lb_flux(:,4))
! if (mype==0) write(*,*) mype, 'DIC(calcite) surface input:',minval(lb_flux(:,5)),maxval(lb_flux(:,5))

end if ! add_loopback

else
Expand Down
1 change: 0 additions & 1 deletion src/int_recom/recom_extra.F90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -201,4 +201,3 @@ subroutine krill_resp(n, partit, mesh)
end if
endif
end subroutine krill_resp

8 changes: 8 additions & 0 deletions src/int_recom/recom_sinking.F90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -392,6 +392,14 @@ subroutine diff_ver_recom_expl(tr_num, tracers, partit, mesh)
bottom_flux = GlodecayBenthos(:,1) * Fe2N_benthos !*** DFe ***
CASE (1022)
bottom_flux = -GlodecayBenthos(:,2) * redO2C !*** O2 ***
CASE (1302)
if (ciso) then
bottom_flux = GlodecayBenthos(:,5) + GlodecayBenthos(:,6) !*** DIC_13 and Calc: DIC_13 ***
end if
CASE (1402)
if (ciso) then
bottom_flux = GlodecayBenthos(:,7) + GlodecayBenthos(:,8) !*** DIC_14 and Calc: DIC_14 ***
end if
CASE DEFAULT
if (partit%mype==0) then
write(*,*) 'check specified in boundary conditions'
Expand Down
79 changes: 75 additions & 4 deletions src/oce_ale_tracer.F90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -232,6 +232,11 @@ subroutine solve_tracers_ale(ice, dynamics, tracers, partit, mesh)
!!! !$ACC UPDATE DEVICE(tracers%work%fct_ttf_min, tracers%work%fct_ttf_max, tracers%work%fct_plus, tracers%work%fct_minus)
!$ACC UPDATE DEVICE (mesh%helem, mesh%hnode, mesh%hnode_new, mesh%zbar_3d_n, mesh%z_3d_n)
do tr_num=1, tracers%num_tracers

!YY: sinkflx needs to be reset at each time step
if(use_MEDUSA) then
SinkFlx = 0.0d0
endif

! do tracer AB (Adams-Bashfort) interpolation only for advectiv part
! needed
Expand Down Expand Up @@ -292,10 +297,20 @@ subroutine solve_tracers_ale(ice, dynamics, tracers, partit, mesh)
if (flag_debug .and. mype==0) print *, achar(27)//'[37m'//' --> call diff_tracers_ale'//achar(27)//'[0m'
call diff_tracers_ale(tr_num, dynamics, tracers, ice, partit, mesh)

! Radioactive decay of 14C and 39Ar
!___________________________________________________________________________
! Radioactive decay of 14C and 39Ar
if (tracers%data(tr_num)%ID == 14) tracers%data(tr_num)%values(:,:) = tracers%data(tr_num)%values(:,:) * exp(-decay14 * dt)
if (tracers%data(tr_num)%ID == 39) tracers%data(tr_num)%values(:,:) = tracers%data(tr_num)%values(:,:) * exp(-decay39 * dt)

!YY: C14 seems to be calculated both in fesom and recom
!YY: decay differently calculated???
#if defined(__recom)
! radioactive decay of 14C
if (ciso_14 .and. any(c14_tracer_id == tracers%data(tr_num)%ID)) then
tracers%data(tr_num)%values(:,:) = tracers%data(tr_num)%values(:,:) * (1 - lambda_14 * dt)
end if ! ciso & ciso_14
#endif

!___________________________________________________________________________
! relax to salt and temp climatology
if (flag_debug .and. mype==0) print *, achar(27)//'[37m'//' --> call relax_to_clim'//achar(27)//'[0m'
Expand Down Expand Up @@ -1614,16 +1629,35 @@ FUNCTION bc_surface(n, id, sval, nzmin, partit)
+ relax_salt(n) - real_salt_flux(n)*is_nonlinfs)
#if defined(__recom)
CASE (1001) ! DIN
if (use_MEDUSA .and. add_loopback) then ! OG: add is_MEDUSA_loopback flag is_MEDUSA_loopback flag * lb_flux(n,1)
bc_surface= dt*(AtmNInput(n) + RiverDIN2D(n) * is_riverinput &
+ ErosionTON2D(n) * is_erosioninput + lb_flux(n,1))
else
bc_surface= dt*(AtmNInput(n) + RiverDIN2D(n) * is_riverinput &
+ ErosionTON2D(n) * is_erosioninput)
end if

CASE (1002) ! DIC
if (use_MEDUSA .and. add_loopback) then
bc_surface= dt*(GloCO2flux_seaicemask(n) &
+ RiverDIC2D(n) * is_riverinput &
+ ErosionTOC2D(n) * is_erosioninput &
+ lb_flux(n,2) + lb_flux(n,5))
else
bc_surface= dt*(GloCO2flux_seaicemask(n) &
+ RiverDIC2D(n) * is_riverinput &
+ ErosionTOC2D(n) * is_erosioninput)
end if

CASE (1003) ! Alk
if (use_MEDUSA .and. add_loopback) then
bc_surface= dt*(virtual_alk(n) + relax_alk(n) &
+ RiverAlk2D(n) * is_riverinput &
+ lb_flux(n,3) + lb_flux(n,5)*2) !CaCO3:Alk burial=1:2
else
bc_surface= dt*(virtual_alk(n) + relax_alk(n) &
+ RiverAlk2D(n) * is_riverinput)
!bc_surface=0.0_WP
end if
CASE (1004:1010)
bc_surface=0.0_WP
CASE (1011) ! DON
Expand All @@ -1633,16 +1667,53 @@ FUNCTION bc_surface(n, id, sval, nzmin, partit)
CASE (1013:1017)
bc_surface=0.0_WP
CASE (1018) ! DSi
bc_surface=dt*(RiverDSi2D(n) * is_riverinput + ErosionTSi2D(n) * is_erosioninput)
if (use_MEDUSA .and. add_loopback) then
bc_surface=dt*(RiverDSi2D(n) * is_riverinput &
+ ErosionTSi2D(n) * is_erosioninput &
+ lb_flux(n,4))
else
bc_surface=dt*(RiverDSi2D(n) * is_riverinput + ErosionTSi2D(n) * is_erosioninput)
end if

CASE (1019) ! Fe
if (useRivFe) then
bc_surface= dt*(AtmFeInput(n) + RiverFe(n))
else
bc_surface= dt*AtmFeInput(n)
end if
CASE (1020:1021) ! Cal
bc_surface=0.0_WP
CASE (1022) ! OXY
bc_surface= dt*GloO2flux_seaicemask(n)
! bc_surface=0.0_WP
CASE (1023:1035)
CASE (1023:1033)
bc_surface=0.0_WP ! OG added bc for recom fields
CASE (1302) ! Before (1033) ! DIC_13
if (ciso) then
if (use_MEDUSA .and. add_loopback) then
bc_surface= dt*(GloCO2flux_seaicemask_13(n) &
+ lb_flux(n,6) + lb_flux(n,7))
else
bc_surface= dt*(GloCO2flux_seaicemask_13(n))
end if
else
bc_surface=0.0_WP
end if
CASE (1305:1321)
bc_surface=0.0_WP ! organic 13C
CASE (1402) ! Before (1034) ! DIC_14
if (ciso .and. ciso_14) then
if (use_MEDUSA .and. add_loopback .and. ciso_organic_14) then
bc_surface= dt*(GloCO2flux_seaicemask_14(n) &
+ lb_flux(n,8) + lb_flux(n,9))
else
bc_surface= dt*GloCO2flux_seaicemask_14(n)
end if
else
bc_surface=0.0_WP
end if
CASE (1405:1421)
bc_surface=0.0_WP ! organic 14C
#endif
CASE (101) ! apply boundary conditions to tracer ID=101
bc_surface= dt*(prec_rain(n))! - real_salt_flux(n)*is_nonlinfs)
Expand Down

0 comments on commit 0864237

Please sign in to comment.