modify to enable BGC runs for Cycle 4. Next step is to set up restart…

… directory.

config.yaml

@@ -19,7 +19,7 @@ input: /g/data/ik11/inputs/access-om2/input_20201102/common_01deg_jra55
 submodels:
     - name: atmosphere
       model: yatm
-      exe: /g/data/ik11/inputs/access-om2/bin/yatm_44e8821.exe
+      exe: /g/data/ik11/inputs/access-om2/bin/yatm_0ab7295.exe
       input:
             - /g/data/ik11/inputs/access-om2/input_20201102/yatm_01deg
             - /g/data/qv56/replicas/input4MIPs/CMIP6/OMIP/MRI/MRI-JRA55-do-1-4-0/atmos/3hr/rsds/gr/v20190429
@@ -37,17 +37,20 @@ submodels:
 
     - name: ocean
       model: mom
-      exe: /g/data/ik11/inputs/access-om2/bin/fms_ACCESS-OM-BGC_af3a94d_libaccessom2_44e8821.x
+      exe: /g/data/ik11/inputs/access-om2/bin/fms_ACCESS-OM-BGC_1d9af9d_libaccessom2_0ab7295.x
       input:
+            - /g/data/ik11/inputs/access-om2/input_bgc_20211105/01deg
             - /g/data/ik11/inputs/access-om2/input_20201120/mom_01deg_150x100
             - /g/data/ik11/inputs/access-om2/input_20201102/mom_01deg
             - /g/data/ik11/inputs/access-om2/passive4-2020-11-13
       ncpus: 10461
 
     - name: ice
       model: cice5
-      exe: /g/data/ik11/inputs/access-om2/bin/cice_auscom_18x15.3600x2700_1682p_9864af2_libaccessom2_44e8821.exe
-      input: /g/data/ik11/inputs/access-om2/input_20201102/cice_01deg
+      exe: /g/data/ik11/inputs/access-om2/bin/cice_auscom_3600x2700_722p_3a5d05f_libaccessom2_0ab7295.exe
+      input: 
+            - /g/data/ik11/inputs/access-om2/input_bgc_20211105/01deg
+            - /g/data/ik11/inputs/access-om2/input_20201102/cice_01deg
       ncpus: 1682
 
 # Collation

ice/cice_in.nml

@@ -141,15 +141,15 @@
     restart_hbrine    = .false.
     restore_bgc       = .false.
     sil_data_type     = 'default'
-    skl_bgc           = .false.
+    skl_bgc           = .true.
     tr_bgc_am_sk      = .false.
     tr_bgc_c_sk       = .false.
     tr_bgc_chl_sk     = .false.
     tr_bgc_dms_sk     = .false.
     tr_bgc_dmspd_sk   = .false.
     tr_bgc_dmspp_sk   = .false.
     tr_bgc_sil_sk     = .false.
-    tr_brine          = .false.
+    tr_brine          = .true.
 /
 
 &forcing_nml
@@ -337,27 +337,27 @@
     f_bgc_dmsp_ml     = 'x'
     f_bgc_dmspd_sk    = 'x'
     f_bgc_dmspp_sk    = 'x'
-    f_bgc_n_sk        = 'x'
-    f_bgc_nit_ml      = 'x'
-    f_bgc_nit_sk      = 'x'
+    f_bgc_n_sk        = 'm'
+    f_bgc_nit_ml      = 'm'
+    f_bgc_nit_sk      = 'm'
     f_bgc_sil_ml      = 'x'
     f_bgc_sil_sk      = 'x'
     f_bphi            = 'x'
     f_btin            = 'x'
     f_faero_atm       = 'x'
     f_faero_ocn       = 'x'
-    f_fbri            = 'x'
+    f_fbri            = 'm'
     f_fn              = 'x'
-    f_fn_ai           = 'x'
+    f_fn_ai           = 'm'
     f_fnh             = 'x'
     f_fnh_ai          = 'x'
     f_fno             = 'x'
-    f_fno_ai          = 'x'
+    f_fno_ai          = 'm'
     f_fsil            = 'x'
     f_fsil_ai         = 'x'
     f_grownet         = 'x'
-    f_hbri            = 'x'
-    f_ppnet           = 'x'
+    f_hbri            = 'm'
+    f_ppnet           = 'm'
 /
 
 &icefields_drag_nml

ice/input_ice.nml

@@ -10,11 +10,11 @@
                         'runof_i', 'tair_i', 'qair_i', 'uwnd_i', 'vwnd_i',
                         'licalvf_i'
     fields_from_ocn   = 'sst_i', 'sss_i', 'ssu_i', 'ssv_i', 'sslx_i', 'ssly_i',
-                        'pfmice_i'
+                        'pfmice_i', 'ssn_i', 'ssalg_i'
     fields_to_ocn     = 'strsu_io', 'strsv_io', 'rain_io', 'snow_io', 'stflx_io',
                         'htflx_io', 'swflx_io', 'qflux_io', 'shflx_io',
                         'lwflx_io', 'runof_io', 'press_io', 'aice_io', 'melt_io',
-                        'form_io', 'licefw_io', 'licefh_io'
+                        'form_io', 'licefw_io', 'licefh_io', 'wnd10_io', 'nit_io', 'alg_io'
     gfdl_surface_flux = .true.
     ice_fwflux        = .true.
     ice_pressure_on   = .true.

namcouple

@@ -16,7 +16,7 @@
 ### 11 fields  atm -> ice
 ### 17 fields  ice -> ocn
 ###  7 fields  ocn -> ice
-   35
+   40
  $END
 #########################################################################
 # $JOBNAME
@@ -265,6 +265,21 @@ NA NA
 ##########
 licefh_io liceht -1 -1 0 i2o.nc IGNORED
 NA NA
+##########
+# Field 29 : 10 m winds
+########## 
+wnd10_io wnd_io -1 -1 0 i2o.nc IGNORED
+NA NA
+##########
+# Field 30 : ice-to-ocean nitrate flux [mmol m-2 s-1]
+##########
+nit_io iof_nit -1 -1 0 i2o.nc IGNORED
+NA NA
+##########
+# Field 31 : ice-to-ocean algae/phyto flux in nitrogen content [mmol m-2 s-1]
+##########
+alg_io iof_alg -1 -1 0 i2o.nc IGNORED
+NA NA
 ############################################################################
 #
 #                      OCEAN --->>> ICE
@@ -303,5 +318,15 @@ NA NA
 # Field 35 :  sea surface slope _y_ (m/m)
 ##########
 dssldy ssly_i -1 -1 0 o2i.nc IGNORED
+NA NA
+##########
+# Field 36 :  sea surface nitrate (mmol/m^3)
+##########
+n_surf ssn_i -1 -1 0 o2i.nc IGNORED
+NA NA
+##########
+# Field 37 :  sea surface algae/phytoplankton in nitrogen (mmol/m^3)
+##########
+alg_surf ssalg_i -1 -1 0 o2i.nc IGNORED
 NA NA
  $END