Commit f21461ba authored by Philippe Klotz's avatar Philippe Klotz
parents 63c62f83 73eb32bd
!$Id: bio_fasham.F90,v 1.11 2007-01-06 11:49:15 kbk Exp $
!$id: bio_fasham.F90,v 1.11 2007-01-06 11:49:15 kbk Exp $
#include"cppdefs.h"
!-----------------------------------------------------------------------
!BOP
......@@ -214,12 +214,8 @@
write(10,901) k2
write(*,900) ' w_p = ',w_p
write(10,901) w_p
write(*,900) ' theta = ',theta
write(10,901) theta
write(*,900) ' w_pmin = ',w_pmin
write(10,901) w_pmin
write(*,900) ' w_pmax = ',w_pmax
write(10,901) w_pmax
write(*,900) ' gmax = ',gmax
write(10,901) gmax
900 format (a,f8.5)
901 format (f8.5)
......
......@@ -234,7 +234,6 @@
endif
900 format (a,f8.5)
901 format (f8.5)
! Conversion from day to second
......@@ -625,12 +624,12 @@
fac4 = hmratio**0.25
! Autotroph metabolic rate
amr1 = max(0.0,ca1*0.25*exp(Ea/(kBeV*T0**2)*(T(ci)/(1+T(ci)/T0))) )
amr2 = max(0.0,ca2*0.25*exp(Ea/(kBeV*T0**2)*(T(ci)/(1+T(ci)/T0)))/fac3)
amr1 = max(0d0,ca1*0.25*exp(Ea/(kBeV*T0**2)*(T(ci)/(1+T(ci)/T0))) )
amr2 = max(0d0,ca2*0.25*exp(Ea/(kBeV*T0**2)*(T(ci)/(1+T(ci)/T0)))/fac3)
! Heterotroph metabolic rate
hmr1 = max(0.0,ch1*0.25*exp(Eh/(kBeV*T0**2)*(T(ci)/(1+T(ci)/T0))) )
hmr2 = max(0.0,ch2*0.25*exp(Eh/(kBeV*T0**2)*(T(ci)/(1+T(ci)/T0)))/fac4)
hmr1 = max(0d0,ch1*0.25*exp(Eh/(kBeV*T0**2)*(T(ci)/(1+T(ci)/T0))) )
hmr2 = max(0d0,ch2*0.25*exp(Eh/(kBeV*T0**2)*(T(ci)/(1+T(ci)/T0)))/fac4)
else
amr1 = 1.0
amr2 = 1.0
......
......@@ -228,7 +228,6 @@
endif
900 format (a,f8.5)
901 format (f8.5)
! Conversion from day to second
......@@ -603,12 +602,12 @@
fac4 = hmratio**0.25
! Autotroph metabolic rate
amr1 = max(0.0,ca1*0.25*exp(Ea/(kBeV*T0**2)*(T(ci)/(1+T(ci)/T0))) )
amr2 = max(0.0,ca2*0.25*exp(Ea/(kBeV*T0**2)*(T(ci)/(1+T(ci)/T0)))/fac3)
amr1 = max(0d0,ca1*0.25*exp(Ea/(kBeV*T0**2)*(T(ci)/(1+T(ci)/T0))) )
amr2 = max(0d0,ca2*0.25*exp(Ea/(kBeV*T0**2)*(T(ci)/(1+T(ci)/T0)))/fac3)
! Heterotroph metabolic rate
hmr1 = max(0.0,ch1*0.25*exp(Eh/(kBeV*T0**2)*(T(ci)/(1+T(ci)/T0))) )
hmr2 = max(0.0,ch2*0.25*exp(Eh/(kBeV*T0**2)*(T(ci)/(1+T(ci)/T0)))/fac4)
hmr1 = max(0d0,ch1*0.25*exp(Eh/(kBeV*T0**2)*(T(ci)/(1+T(ci)/T0))) )
hmr2 = max(0d0,ch2*0.25*exp(Eh/(kBeV*T0**2)*(T(ci)/(1+T(ci)/T0)))/fac4)
else
amr1 = 1.0
amr2 = 1.0
......
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