Commit 4e8d9489 authored by Markus Pahlow's avatar Markus Pahlow

remove trailing space

parent 93634994
......@@ -77,7 +77,7 @@
real diff_tx, diff_ty, diff_tz, zmax, cont, drho
# if defined O_benthic_denitrification
real bdeni, lno3, din15flag, sg_bdeni
# endif
# endif
# if defined O_isopycmix && defined O_gent_mcwilliams && !defined O_fct && !defined O_quicker
real adv_tziso, adv_txiso, adv_tyiso
# endif
......@@ -420,7 +420,7 @@
!-----------------------------------------------------------------------
! Initialize tracers.
! In the old code tracers were limited to positive values here
! this should not be done here! Only in npzd_src.F
! this should not be done here! Only in npzd_src.F
!-----------------------------------------------------------------------
tnpzd(1) = t(i,k,j,ipo4,taum1)
......@@ -503,14 +503,14 @@
fe_m = 4
! find fe data day index
call hunt (fe_x,fe_n,dayoyr,fe_jlo)
! initialize the fe data array at the right day
! initialize the fe data array at the right day
fe_k = min(max(fe_jlo-(fe_m-1)/2,1),fe_n+1-fe_m)
! interpolate the fe data, note this does not use the whole array
call polint (fe_x(fe_k),fe_y(fe_k)
& , fe_m,dayoyr,fe_conc,fe_dy)
# endif ! O_TMM (452)
fe_conc = max(fe_conc, 0.)
! calculate the fe limitation term
! calculate the fe limitation term
felimit = fe_conc/(kfe + fe_conc)
felimit_D = fe_conc/(kfe_D + fe_conc)
else
......@@ -607,7 +607,7 @@
# endif
# endif
& ,i,j,k,tlat(i,j))
! These are source/sink terms
! These are source/sink terms
snpzd = snpzd*rdtts(k)
expo = expo*rnbio(k)
# if defined O_npzd_iron
......@@ -735,8 +735,8 @@
# if defined O_benthic_denitrification
sg_bdeni = 0.
bdeni = 0.
# endif
# endif
!-----------------------------------------------------------------------
! calculate detritus at the bottom and remineralize
!-----------------------------------------------------------------------
......@@ -852,7 +852,7 @@
# if defined O_npzd_nitrogen || defined O_opem
# if defined O_benthic_denitrification
!-----------------------------------------------------------------------
! benthic denitrification model of Bohlen et al., 2012, GBC
! benthic denitrification model of Bohlen et al., 2012, GBC
! NO3 is removed out of bottom water nitrate.
! See Somes et al., 2012, BGS for additional details/results
!-----------------------------------------------------------------------
......@@ -860,24 +860,24 @@
no3flag = 0.5+sign(0.5,t(i,k,j,ino3,taum1)-trcmin)
din15flag = 1.0
lno3 = 0.5*tanh(t(i,k,j,ino3,taum1)*10 - 5.0)
bdeni = (0.06 + 0.19*0.99
& **(max(t(i,k,j,io2,taum1)*1000,trcmin)
& - max(t(i,k,j,ino3,taum1),trcmin)))
# if defined O_opem
& *max(expo_c,trcmin)
# else
& *max(expo,trcmin)*6.625
& *max(expo,trcmin)*6.625
# endif
bdeni = min(bdeni, expo)
bdeni = max(bdeni, 0.)
c if (k .ge. 1 .and. k .le. 10)
c if (k .ge. 1 .and. k .le. 10)
c & bdeni = sgbdfac*bdeni
bdeni = bdeni*(0.5 + lno3)*no3flag*din15flag
snpzd(5) = snpzd(5) + expo - bdeni
# else
# else
snpzd(5) = snpzd(5) + expo
# endif
# endif
......@@ -920,15 +920,15 @@ c & bdeni = sgbdfac*bdeni
rexpofe(i,k,j) = expofe
# endif
# endif
!-----------------------------------------------------------------------
! End Subgrid bathymetrie
! Additional output
! Additional output
!-----------------------------------------------------------------------
# if defined O_benthic_denitrification && O_save_npzd
rbdeni(i,k,jrow) = bdeni + sg_bdeni
# endif
# endif
!-----------------------------------------------------------------------
! set source/sink terms
!-----------------------------------------------------------------------
......@@ -1036,7 +1036,7 @@ c & bdeni = sgbdfac*bdeni
if (timavgperts .and. .not. euler2) then
ta_rnpp(i,k,jrow) = ta_rnpp(i,k,jrow) + rnpp(i,k,j)
ta_rgraz(i,k,jrow) = ta_rgraz(i,k,jrow) + rgraz(i,k,j)
ta_rgraz_Z(i,k,jrow) = ta_rgraz_Z(i,k,jrow)
ta_rgraz_Z(i,k,jrow) = ta_rgraz_Z(i,k,jrow)
& + rgraz_Z(i,k,j)
ta_rgraz_Det(i,k,jrow) = ta_rgraz_Det(i,k,jrow)
& + rgraz_Det(i,k,j)
......@@ -1118,7 +1118,7 @@ c & bdeni = sgbdfac*bdeni
& + rmorp_D(i,k,j)
ta_rnfix(i,k,jrow) = ta_rnfix(i,k,jrow) + rnfix(i,k,j)
# if defined O_benthic_denitrification
ta_rbdeni(i,k,jrow) = ta_rbdeni(i,k,jrow)
ta_rbdeni(i,k,jrow) = ta_rbdeni(i,k,jrow)
& + rbdeni(i,k,j)
# endif
# if defined O_npzzd
......@@ -1221,7 +1221,7 @@ c & bdeni = sgbdfac*bdeni
so2 = -src(i,k,j,ispo4)*redotp - rnfix(i,k,j)*1.25e-3
# endif
src(i,k,j,iso2) = so2*(0.5 + fo2)
# if defined O_npzd_nitrogen || defined O_opem
! add denitrification as source term for NO3
no3flag = 0.5+sign(0.5,t(i,k,j,ino3,taum1)-trcmin)
......@@ -1309,7 +1309,7 @@ c & bdeni = sgbdfac*bdeni
ta_rdeni(i,k,jrow) = ta_rdeni(i,k,jrow) + rdeni(i,k,j)
# endif
# if defined O_npzd_iron
ta_rremife(i,k,jrow) = ta_rremife(i,k,jrow)
ta_rremife(i,k,jrow) = ta_rremife(i,k,jrow)
& + rremife(i,k,j)
ta_rexpofe(i,k,jrow) = ta_rexpofe(i,k,jrow)
& + rexpofe(i,k,j)
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment