pass flux_bio to MOM for updating ocean bgc tracers within MOM instea…

…d of CICE. consequently, no need dz_ocean1.
hakaseh · Sep 16, 2021 · 05862df · 05862df
1 parent 08a1112
commit 05862df
Show file tree

Hide file tree

Showing 4 changed files with 7 additions and 14 deletions.
diff --git a/drivers/auscom/cpl_forcing_handler.F90 b/drivers/auscom/cpl_forcing_handler.F90
@@ -23,7 +23,7 @@ module cpl_forcing_handler
     use cpl_netcdf_setup
     use cpl_arrays_setup
     use ice_calendar, only: dt
-    use ice_zbgc_shared, only: ocean_bio,nlt_bgc_N,nlt_bgc_NO
+    use ice_zbgc_shared, only: ocean_bio,flux_bio,nlt_bgc_N,nlt_bgc_NO
 
 implicit none
 
@@ -812,9 +812,9 @@ subroutine get_i2o_fluxes
 !16 10m wind. To mask or not to mask?
   tiownd(:,:,:) = sqrt(uatm(:,:,:)**2 + vatm(:,:,:)**2)
 !17? sea surface nitrate updated after ice-water flux
-  tionit(:,:,:) = ocean_bio(:,:,nlt_bgc_NO,:)
+  tionit(:,:,:) = flux_bio(:,:,nlt_bgc_NO,:)*aice(:,:,:)
 !18 sea surface algae/phyto/nitrogen updated after ice-water flux
-  tioalg(:,:,:) = ocean_bio(:,:,nlt_bgc_N,:)
+  tioalg(:,:,:) = flux_bio(:,:,nlt_bgc_N,:)*aice(:,:,:)
 
 return
 end subroutine get_i2o_fluxes

diff --git a/source/ice_algae.F90 b/source/ice_algae.F90
@@ -198,7 +198,7 @@ subroutine skl_biogeochemistry (nx_block, ny_block,  &
                                       indxi,    indxj,     &
                                       nbtrcr,              &
                                       flux_bio, ocean_bio, &
-                                      dz_ocean1,     aicen,     &
+                                      aicen,     &
                                       meltb,    congel,    &
                                       fswthru,  first_ice, &
                                       trcrn,    grow_Cn)
@@ -215,7 +215,6 @@ subroutine skl_biogeochemistry (nx_block, ny_block,  &
          indxi, indxj    ! compressed indices for cells with aicen > puny
 
       real (kind=dbl_kind), intent(in) :: &
-         dz_ocean1, &  ! first ocean layer thickness (m)
          dt         ! time step 
 
       logical (kind=log_kind), dimension (nx_block,ny_block), &
@@ -466,9 +465,8 @@ subroutine skl_biogeochemistry (nx_block, ny_block,  &
 
          ! Uncomment to update ocean concentration
          ! Currently not coupled with ocean biogeochemistry
-         ! HH uncommented to enable ice-ocean BGC coupling for ACCESS-OM2
-         ocean_bio(i,j,nn) = ocean_bio(i,j,nn) &
-                           + flux_bio(i,j,nn)/dz_ocean1*aicen(i,j)
+         !ocean_bio(i,j,nn) = ocean_bio(i,j,nn) &
+         !                  + flux_bio(i,j,nn)/hmix(i,j)*aicen(i,j)
 
          if (cinit(ij,nn) < c0) then
               write(nu_diag,*) 'sk_bgc < 0 after algal fluxes, ij,nn,cinit,flux_bio',&

diff --git a/source/ice_zbgc.F90 b/source/ice_zbgc.F90
@@ -57,7 +57,7 @@ subroutine init_zbgc
          restore_bgc, skl_bgc, &
          tr_bgc_C_sk, tr_bgc_chl_sk, tr_bgc_Am_sk, tr_bgc_Sil_sk, &
          tr_bgc_DMSPp_sk, tr_bgc_DMSPd_sk, tr_bgc_DMS_sk, &
-         restart_bgc, restart_hbrine, phi_snow, bgc_flux_type, dz_ocean1
+         restart_bgc, restart_hbrine, phi_snow, bgc_flux_type
 
       !-----------------------------------------------------------------
       ! default values
@@ -80,7 +80,6 @@ subroutine init_zbgc
       restart_hbrine  = .false.  ! hbrine restart
       phi_snow        = p5       ! snow porosity
       bgc_flux_type   = 'Jin2006'! type of ocean-ice poston velocity ('constant')
-      dz_ocean1       = 2.3035 ! HH first layer thickness of ocean
 
       !-----------------------------------------------------------------
       ! read from input file
@@ -182,7 +181,6 @@ subroutine init_zbgc
       call broadcast_scalar(tr_bgc_DMSPp_sk,    master_task)
       call broadcast_scalar(tr_bgc_DMSPd_sk,    master_task)
       call broadcast_scalar(tr_bgc_DMS_sk,      master_task)
-      call broadcast_scalar(dz_ocean1,          master_task)
 
       if (skl_bgc) then
 
@@ -206,7 +204,6 @@ subroutine init_zbgc
          write(nu_diag,1010) ' tr_bgc_DMSPp_sk           = ', tr_bgc_DMSPp_sk
          write(nu_diag,1010) ' tr_bgc_DMSPd_sk           = ', tr_bgc_DMSPd_sk
          write(nu_diag,1010) ' tr_bgc_DMS_sk             = ', tr_bgc_DMS_sk
-         write(nu_diag,1010) ' dz_ocean1                 = ', dz_ocean1
 
       endif   ! master_task
 
@@ -726,7 +723,6 @@ subroutine biogeochemistry (dt, iblk)
                                          nbtrcr,                        &
                                          flux_bion(:,:,1:nbtrcr),       &
                                          ocean_bio(:,:,1:nbtrcr, iblk), &
-                                         dz_ocean1                    , &
                                          aicen    (:,:,        n,iblk), & 
                                          meltbn   (:,:,        n,iblk), &
                                          congeln  (:,:,        n,iblk), &

diff --git a/source/ice_zbgc_shared.F90 b/source/ice_zbgc_shared.F90
@@ -126,7 +126,6 @@ module ice_zbgc_shared
          thinS     = 0.05_dbl_kind     ! minimum ice thickness for brine
 
       real (kind=dbl_kind), public :: & 
-         dz_ocean1  = 2.3035_dbl_kind,  & ! first ocean layer thickness (m)
          phi_snow ,  &  ! porosity of snow
          flood_frac     ! fraction of ocean/meltwater that floods