Skip to content

Commit

Permalink
intmld and int100 bug fix: COSIMA/01deg_jra55_iaf#11 (comment)
Browse files Browse the repository at this point in the history
  • Loading branch information
hakase hayashida committed Feb 24, 2022
1 parent fa34a8d commit 398a0ae
Showing 1 changed file with 14 additions and 14 deletions.
28 changes: 14 additions & 14 deletions src/mom5/ocean_csiro_bgc/bio_v3.inc
View file
Original file line number Diff line number Diff line change
Expand Up @@ -423,26 +423,26 @@ do n = 1, instances !{
npp3d(i,j,k)=rdtts*npp3d(i,j,k)*grid%tmask(i,j,k)
if (Grid%zw(k) .le. mld(i,j)) then
adic_intmld(i,j) = adic_intmld(i,j) + biotr(i,k,id_adic) * thickness%dzt(i,j,k)
dic_intmld(i,j) = dic_intmld(i,j) + biotr(i,k,id_dic) * thickness%dzt(i,j,k)
o2_intmld(i,j) = o2_intmld(i,j) + biotr(i,k,id_o2) * thickness%dzt(i,j,k)
no3_intmld(i,j) = no3_intmld(i,j) + biotr(i,k,id_no3) * thickness%dzt(i,j,k)
fe_intmld(i,j) = fe_intmld(i,j) + biotr(i,k,id_fe) * thickness%dzt(i,j,k)
phy_intmld(i,j) = phy_intmld(i,j) + biotr(i,k,id_phy) * thickness%dzt(i,j,k)
det_intmld(i,j) = det_intmld(i,j) + biotr(i,k,id_det) * thickness%dzt(i,j,k)
adic_intmld(i,j)= adic_intmld(i,j)+ t_prog(ind_adic)%field(i,j,k,time%taum1)* thickness%dzt(i,j,k)
dic_intmld(i,j) = dic_intmld(i,j) + t_prog(ind_dic)%field(i,j,k,time%taum1) * thickness%dzt(i,j,k)
o2_intmld(i,j) = o2_intmld(i,j) + t_prog(ind_o2)%field(i,j,k,time%taum1) * thickness%dzt(i,j,k)
no3_intmld(i,j) = no3_intmld(i,j) + t_prog(ind_no3)%field(i,j,k,time%taum1) * thickness%dzt(i,j,k)
fe_intmld(i,j) = fe_intmld(i,j) + t_prog(ind_fe)%field(i,j,k,time%taum1) * thickness%dzt(i,j,k)
phy_intmld(i,j) = phy_intmld(i,j) + t_prog(ind_phy)%field(i,j,k,time%taum1) * thickness%dzt(i,j,k)
det_intmld(i,j) = det_intmld(i,j) + t_prog(ind_det)%field(i,j,k,time%taum1) * thickness%dzt(i,j,k)
pprod_gross_intmld(i,j) = pprod_gross_intmld(i,j) + pprod_gross(i,j,k) * thickness%dzt(i,j,k)
npp_intmld(i,j) = npp_intmld(i,j) + npp3d(i,j,k) * thickness%dzt(i,j,k)
radbio_intmld(i,j) = radbio_intmld(i,j) + radbio3d(i,j,k) * thickness%dzt(i,j,k)
endif
if (Grid%zw(k) .le. 100) then
adic_int100(i,j) = adic_int100(i,j) + biotr(i,k,id_adic) * thickness%dzt(i,j,k)
dic_int100(i,j) = dic_int100(i,j) + biotr(i,k,id_dic) * thickness%dzt(i,j,k)
o2_int100(i,j) = o2_int100(i,j) + biotr(i,k,id_o2) * thickness%dzt(i,j,k)
no3_int100(i,j) = no3_int100(i,j) + biotr(i,k,id_no3) * thickness%dzt(i,j,k)
fe_int100(i,j) = fe_int100(i,j) + biotr(i,k,id_fe) * thickness%dzt(i,j,k)
phy_int100(i,j) = phy_int100(i,j) + biotr(i,k,id_phy) * thickness%dzt(i,j,k)
det_int100(i,j) = det_int100(i,j) + biotr(i,k,id_det) * thickness%dzt(i,j,k)
adic_int100(i,j)= adic_int100(i,j)+ t_prog(ind_adic)%field(i,j,k,time%taum1)* thickness%dzt(i,j,k)
dic_int100(i,j) = dic_int100(i,j) + t_prog(ind_dic)%field(i,j,k,time%taum1) * thickness%dzt(i,j,k)
o2_int100(i,j) = o2_int100(i,j) + t_prog(ind_o2)%field(i,j,k,time%taum1) * thickness%dzt(i,j,k)
no3_int100(i,j) = no3_int100(i,j) + t_prog(ind_no3)%field(i,j,k,time%taum1) * thickness%dzt(i,j,k)
fe_int100(i,j) = fe_int100(i,j) + t_prog(ind_fe)%field(i,j,k,time%taum1) * thickness%dzt(i,j,k)
phy_int100(i,j) = phy_int100(i,j) + t_prog(ind_phy)%field(i,j,k,time%taum1) * thickness%dzt(i,j,k)
det_int100(i,j) = det_int100(i,j) + t_prog(ind_det)%field(i,j,k,time%taum1) * thickness%dzt(i,j,k)
pprod_gross_int100(i,j) = pprod_gross_int100(i,j) + pprod_gross(i,j,k) * thickness%dzt(i,j,k)
npp_int100(i,j) = npp_int100(i,j) + npp3d(i,j,k) * thickness%dzt(i,j,k)
radbio_int100(i,j) = radbio_int100(i,j) + radbio3d(i,j,k) * thickness%dzt(i,j,k)
Expand Down

0 comments on commit 398a0ae

Please sign in to comment.