@@ -4057,14 +4057,6 @@ subroutine generic_WOMBATmid_update_from_source(tracer_list, Temp, Salt, &
40574057 + wombat% docremi(i,j,k) &
40584058 + wombat% zooresp(i,j,k) &
40594059 + wombat% mesresp(i,j,k) &
4060- + wombat% zooexcrphy(i,j,k) &
4061- + wombat% zooexcrdia(i,j,k) &
4062- + wombat% zooexcrdet(i,j,k) &
4063- + wombat% mesexcrphy(i,j,k) &
4064- + wombat% mesexcrdia(i,j,k) &
4065- + wombat% mesexcrdet(i,j,k) &
4066- + wombat% mesexcrbdet(i,j,k) &
4067- + wombat% mesexcrzoo(i,j,k) &
40684060 - wombat% phygrow(i,j,k) * wombat% phy_lnh4(i,j,k) / ( wombat% phy_lnit(i,j,k) + epsi ) &
40694061 - wombat% diagrow(i,j,k) * wombat% dia_lnh4(i,j,k) / ( wombat% dia_lnit(i,j,k) + epsi ) )
40704062
@@ -4245,6 +4237,14 @@ subroutine generic_WOMBATmid_update_from_source(tracer_list, Temp, Salt, &
42454237 + wombat% bdetremi(i,j,k) &
42464238 + wombat% phyresp(i,j,k) &
42474239 + wombat% diaresp(i,j,k) &
4240+ + wombat% zooexcrphy(i,j,k) &
4241+ + wombat% zooexcrdia(i,j,k) &
4242+ + wombat% zooexcrdet(i,j,k) &
4243+ + wombat% mesexcrphy(i,j,k) &
4244+ + wombat% mesexcrdia(i,j,k) &
4245+ + wombat% mesexcrdet(i,j,k) &
4246+ + wombat% mesexcrbdet(i,j,k) &
4247+ + wombat% mesexcrzoo(i,j,k) &
42484248 - wombat% docremi(i,j,k) )
42494249
42504250 ! Oxygen equation ! [molO2/kg]
@@ -4254,14 +4254,6 @@ subroutine generic_WOMBATmid_update_from_source(tracer_list, Temp, Salt, &
42544254 + wombat% docremi(i,j,k) * (1.0 - wombat% fdenitrif(i,j,k)) &
42554255 + wombat% zooresp(i,j,k) &
42564256 + wombat% mesresp(i,j,k) &
4257- + wombat% zooexcrphy(i,j,k) &
4258- + wombat% zooexcrdia(i,j,k) &
4259- + wombat% zooexcrdet(i,j,k) &
4260- + wombat% mesexcrphy(i,j,k) &
4261- + wombat% mesexcrdia(i,j,k) &
4262- + wombat% mesexcrdet(i,j,k) &
4263- + wombat% mesexcrbdet(i,j,k) &
4264- + wombat% mesexcrzoo(i,j,k) &
42654257 - wombat% phygrow(i,j,k) &
42664258 - wombat% diagrow(i,j,k) ) &
42674259 - 40 ./ 16 . * dtsb * ( wombat% ammox(i,j,k) )
@@ -4282,14 +4274,6 @@ subroutine generic_WOMBATmid_update_from_source(tracer_list, Temp, Salt, &
42824274 + wombat% docremi(i,j,k) &
42834275 + wombat% zooresp(i,j,k) &
42844276 + wombat% mesresp(i,j,k) &
4285- + wombat% zooexcrphy(i,j,k) &
4286- + wombat% zooexcrdia(i,j,k) &
4287- + wombat% zooexcrdet(i,j,k) &
4288- + wombat% mesexcrphy(i,j,k) &
4289- + wombat% mesexcrdia(i,j,k) &
4290- + wombat% mesexcrdet(i,j,k) &
4291- + wombat% mesexcrbdet(i,j,k) &
4292- + wombat% mesexcrzoo(i,j,k) &
42934277 - wombat% phygrow(i,j,k) &
42944278 - wombat% diagrow(i,j,k) &
42954279 - wombat% zooslopphy(i,j,k) * wombat% pic2poc(i,j,k) &
@@ -4304,14 +4288,6 @@ subroutine generic_WOMBATmid_update_from_source(tracer_list, Temp, Salt, &
43044288 + wombat% docremi(i,j,k) &
43054289 + wombat% zooresp(i,j,k) &
43064290 + wombat% mesresp(i,j,k) &
4307- + wombat% zooexcrphy(i,j,k) &
4308- + wombat% zooexcrdia(i,j,k) &
4309- + wombat% zooexcrdet(i,j,k) &
4310- + wombat% mesexcrphy(i,j,k) &
4311- + wombat% mesexcrdia(i,j,k) &
4312- + wombat% mesexcrdet(i,j,k) &
4313- + wombat% mesexcrbdet(i,j,k) &
4314- + wombat% mesexcrzoo(i,j,k) &
43154291 - wombat% phygrow(i,j,k) &
43164292 - wombat% diagrow(i,j,k) &
43174293 - wombat% zooslopphy(i,j,k) * wombat% pic2poc(i,j,k) &
@@ -4328,14 +4304,6 @@ subroutine generic_WOMBATmid_update_from_source(tracer_list, Temp, Salt, &
43284304 + wombat% docremi(i,j,k) &
43294305 + wombat% zooresp(i,j,k) &
43304306 + wombat% mesresp(i,j,k) &
4331- + wombat% zooexcrphy(i,j,k) &
4332- + wombat% zooexcrdia(i,j,k) &
4333- + wombat% zooexcrdet(i,j,k) &
4334- + wombat% mesexcrphy(i,j,k) &
4335- + wombat% mesexcrdia(i,j,k) &
4336- + wombat% mesexcrdet(i,j,k) &
4337- + wombat% mesexcrbdet(i,j,k) &
4338- + wombat% mesexcrzoo(i,j,k) &
43394307 - wombat% phygrow(i,j,k) * wombat% phy_lnh4(i,j,k) / ( wombat% phy_lnit(i,j,k) + epsi ) &
43404308 - wombat% diagrow(i,j,k) * wombat% dia_lnh4(i,j,k) / ( wombat% dia_lnit(i,j,k) + epsi ) ) &
43414309 + dtsb * ( 0.0 &
0 commit comments