diff --git a/modulefiles/module_base.hera b/modulefiles/module_base.hera
index 32d9c72d5a..b77fbe6e86 100644
--- a/modulefiles/module_base.hera
+++ b/modulefiles/module_base.hera
@@ -36,4 +36,7 @@ module load z/1.2.11
 
 #Load from emc.nemspara
 module use -a /scratch1/NCEPDEV/nems/emc.nemspara/soft/modulefiles
-module load esmf/8.0.0
+module load netcdf_parallel/4.7.4.release
+module load esmf/8.1.0bs27
+module load hdf5_parallel/1.10.6.release
+
diff --git a/modulefiles/module_base.orion b/modulefiles/module_base.orion
index 4c790db468..2f78059706 100644
--- a/modulefiles/module_base.orion
+++ b/modulefiles/module_base.orion
@@ -50,7 +50,7 @@ module load prod_util/1.2.0
 ##
 module use /apps/contrib/NCEPLIBS/lib/modulefiles
 module load netcdfp/4.7.4.release
-module load esmflocal/8.1.0.21bs.release
+module load esmflocal/8.1.0.27bs.release
 module load post-intel-sandybridge/8.0.5
 
 module load wgrib/2.0.8
diff --git a/sorc/build_fv3_coupled.sh b/sorc/build_fv3_coupled.sh
index d33992c74c..5aee3e911e 100755
--- a/sorc/build_fv3_coupled.sh
+++ b/sorc/build_fv3_coupled.sh
@@ -22,6 +22,6 @@ if [ $target = orion ]; then target=orion.intel ; fi
 cd fv3_coupled.fd/
 FV3=$( pwd -P )/FV3
 cd tests/
-./compile.sh "$FV3" "$target" "CCPP=Y SUITES=FV3_GFS_2017_coupled,FV3_GFS_2017_satmedmf_coupled,FV3_GFS_v15p2_coupled MOM6=Y CICE=Y WW3=Y CMEPS=Y" 1
+./compile.sh "$FV3" "$target" "CCPP=Y SUITES=FV3_GFS_2017_coupled,FV3_GFS_2017_satmedmf_coupled,FV3_GFS_v15p2_coupled MOM6=Y CICE6=Y WW3=Y CMEPS=Y" 1
 mv -f fcst_1.exe ../NEMS/exe/nems_fv3_ccpp_mom6_cice5_ww3.x
 
diff --git a/sorc/checkout.sh b/sorc/checkout.sh
index 98fcd2c3af..431a3afba8 100755
--- a/sorc/checkout.sh
+++ b/sorc/checkout.sh
@@ -16,7 +16,7 @@ rm -f ${topdir}/checkout-fv3_coupled.log
 if [[ ! -d fv3_coupled.fd ]] ; then
     git clone https://github.com/ufs-community/ufs-s2s-model fv3_coupled.fd >> ${topdir}/checkout-fv3_coupled.log 2>&1
     cd fv3_coupled.fd
-    git checkout d7b7081f410eb84902a43c7479e14d44dfed0144 
+    git checkout s2s_prototype5.0
     git submodule update --init --recursive
     cd ${topdir}
 else
diff --git a/ush/forecast_postdet.sh b/ush/forecast_postdet.sh
index a991d28590..3b4acb0e26 100755
--- a/ush/forecast_postdet.sh
+++ b/ush/forecast_postdet.sh
@@ -746,6 +746,7 @@ CICE_postdet()
         echo "Link CICE fixed files"
         $NLN -sf $FIXcice/${ice_grid_file} $DATA/
         $NLN -sf $FIXcice/${ice_kmt_file} $DATA/
+        $NLN -sf $FIXcice/$MESHICE $DATA/
 }
 
 CICE_nml()
diff --git a/ush/nems.configure.med_atm_ocn_ice_wav.IN b/ush/nems.configure.med_atm_ocn_ice_wav.IN
index 239ab765fe..703ae4b89e 100644
--- a/ush/nems.configure.med_atm_ocn_ice_wav.IN
+++ b/ush/nems.configure.med_atm_ocn_ice_wav.IN
@@ -31,7 +31,6 @@ OCN_attributes::
   DumpFields = @[DumpFields]
   ProfileMemory = false
   OverwriteSlice = true
-  dbug_flag = 0
 ::
 
 # ICE #
@@ -42,7 +41,10 @@ ICE_attributes::
   DumpFields = @[DumpFields]
   ProfileMemory = false
   OverwriteSlice = true
-  dbug_flag = 0 
+  mesh_ice = @[MESHICE]
+  stop_n = @[FHMAX]
+  stop_option = nhours
+  stop_ymd = -999
 ::
 
 # WAV #
@@ -109,4 +111,6 @@ ALLCOMP_attributes::
       restart_n = @[restart_interval_hours]
       restart_option = nhours
       restart_ymd = -999
+      dbug_flag = 0
+      use_coldstart = .true.
 ::
diff --git a/ush/nems.configure.medcold_atm_ocn_ice_wav.IN b/ush/nems.configure.medcold_atm_ocn_ice_wav.IN
index 7e8629bf3c..cd011889e3 100644
--- a/ush/nems.configure.medcold_atm_ocn_ice_wav.IN
+++ b/ush/nems.configure.medcold_atm_ocn_ice_wav.IN
@@ -31,7 +31,6 @@ OCN_attributes::
   DumpFields = @[DumpFields]
   ProfileMemory = false
   OverwriteSlice = true
-  dbug_flag = 0
 ::
 
 # ICE #
@@ -42,7 +41,10 @@ ICE_attributes::
   DumpFields = @[DumpFields]
   ProfileMemory = false
   OverwriteSlice = true
-  dbug_flag = 0  
+  mesh_ice = @[MESHICE]
+  stop_n = @[FHMAX]
+  stop_option = nhours
+  stop_ymd = -999
 ::
 
 # WAV #
@@ -109,4 +111,5 @@ ALLCOMP_attributes::
       restart_n = 1
       restart_option = nhours
       restart_ymd = -999
+      dbug_flag = 0
 ::
diff --git a/ush/nems_configure.sh b/ush/nems_configure.sh
index 840cb685d7..d6ed6a5e5f 100755
--- a/ush/nems_configure.sh
+++ b/ush/nems_configure.sh
@@ -74,6 +74,8 @@ fi
 if [ $cplice = .true. ]; then
 	sed -i -e "s;@\[ice_model\];$ICE_model;g" tmp1
 	sed -i -e "s;@\[ice_petlist_bounds\];$ice_petlist_bounds;g" tmp1
+        sed -i -e "s;@\[MESHICE\];$MESHICE;g" tmp1
+        sed -i -e "s;@\[FHMAX\];$FHMAX_GFS;g" tmp1
 fi
 if [ $cplchem = .true. ]; then
 	sed -i -e "s;@\[chem_model\];gsd;g" tmp1
diff --git a/ush/parsing_namelists_CICE.sh b/ush/parsing_namelists_CICE.sh
index 26e7527cda..66edff5251 100644
--- a/ush/parsing_namelists_CICE.sh
+++ b/ush/parsing_namelists_CICE.sh
@@ -4,371 +4,594 @@ CICE_namelists(){
 
 cat > ice_in <<eof  
 &setup_nml
-    days_per_year  = 365
-  , use_leap_years = .true.
-  , year_init      = $year
-  , istep0         = $steps
-  , dt             = $ICETIM
-  , npt            = $npt
-  , ndtd           = 1
-  , runtype        = '$RUNTYPE' 
-  , runid          = 'cpcice' 
-  , ice_ic         = '$iceic'
-  , restart        = .true.
-  , restart_ext    = .false.
-  , use_restart_time = $USE_RESTART_TIME
-  , restart_format = 'nc'
-  , lcdf64         = .false.
-  , restart_dir    = './restart/'
-  , restart_file   = 'iced'
-  , pointer_file   = './restart/ice.restart_file'
-  , dumpfreq       = '$dumpfreq'
-  , dumpfreq_n     =  $dumpfreq_n
-  , dump_last      = .false.  
-  , diagfreq       = 6
-  , diag_type      = 'file'
-  , diag_file      = 'ice_diag.d'
-  , print_global   = .true.
-  , print_points   = .true.
-  , latpnt(1)      =  90.
-  , lonpnt(1)      =   0.
-  , latpnt(2)      = -65.
-  , lonpnt(2)      = -45.
-  , dbug           = .false.
-  , histfreq       = 'm','d','h','x','x'
-  , histfreq_n     =  0 , 0 , 6 , 1 , 1
-  , hist_avg       = $cice_hist_avg
-  , history_dir    = './history/'
-  , history_file   = 'iceh'
-  , write_ic       = .true.
-  , incond_dir     = './history/'
-  , incond_file    = 'iceh_ic'
+   days_per_year  = 365
+   use_leap_years = .true.
+   year_init      = $year
+   istep0         = $steps
+   dt             = $ICETIM
+   npt            = $npt
+   ndtd           = 1
+   runtype        = '$RUNTYPE' 
+   runid          = 'cpcice' 
+   ice_ic         = '$iceic'
+   restart        = .true.
+   restart_ext    = .false.
+   use_restart_time = $USE_RESTART_TIME
+   restart_format = 'nc'
+   lcdf64         = .false.
+   numin          = 21
+   numax          = 89
+   restart_dir    = './RESTART/'
+   restart_file   = 'iced'
+   pointer_file   = './ice.restart_file'
+   dumpfreq       = '$dumpfreq'
+   dumpfreq_n     =  $dumpfreq_n
+   dump_last      = .false.  
+   bfbflag        = 'off'
+   diagfreq       = 6
+   diag_type      = 'file'
+   diag_file      = 'ice_diag.d'
+   print_global   = .true.
+   print_points   = .true.
+   latpnt(1)      =  90.
+   lonpnt(1)      =   0.
+   latpnt(2)      = -65.
+   lonpnt(2)      = -45.
+   dbug           = .false.
+   histfreq       = 'm','d','h','x','x'
+   histfreq_n     =  0 , 0 , 6 , 1 , 1
+   hist_avg       = $cice_hist_avg
+   history_dir    = './history/'
+   history_file   = 'iceh'
+   write_ic       = .true.
+   incond_dir     = './history/'
+   incond_file    = 'iceh_ic'
+   version_name   = 'CICE_6.0.2'
 /
 
 &grid_nml
-    grid_format  = 'nc'
-  , grid_type    = 'displaced_pole'
-  , grid_file    = '$ice_grid_file'
-  , kmt_file     = '$ice_kmt_file'
-  , kcatbound    = 0
-/
-
-&domain_nml
-    nprocs = 48 
-  , processor_shape   = 'slenderX2'
-  , distribution_type = 'cartesian'
-  , distribution_wght = 'latitude'
-  , ew_boundary_type  = 'cyclic'
-  , ns_boundary_type  = 'tripole'
-  , maskhalo_dyn      = .false.
-  , maskhalo_remap    = .false.
-  , maskhalo_bound    = .false.
+   grid_format  = 'nc'
+   grid_type    = 'tripole'
+   grid_file    = '$ice_grid_file'
+   kmt_file     = '$ice_kmt_file'
+   use_bathymetry = .false.
+   kcatbound    = 0
+   ncat         = 5
+   nfsd         = 1
+   nilyr        = 7
+   nslyr        = 1
+   nblyr        = 1
+   nfsd         = 1
 /
 
 &tracer_nml
-    tr_iage      = .true.
-  , restart_age  = .false.
-  , tr_FY        = .true.
-  , restart_FY   = .false.
-  , tr_lvl       = .true.
-  , restart_lvl  = .false.
-  , tr_pond_cesm = .false.
-  , restart_pond_cesm = .false.
-  , tr_pond_topo = .false.
-  , restart_pond_topo = .false.
-  , tr_pond_lvl  = $tr_pond_lvl
-  , restart_pond_lvl  = $restart_pond_lvl
-  , tr_aero      = .false.
-  , restart_aero = .false.
+   n_aero       = 1
+   n_zaero      = 0
+   n_algae      = 0
+   n_doc        = 0
+   n_dic        = 0
+   n_don        = 0
+   n_fed        = 0
+   n_fep        = 0
+   tr_iage      = .true.
+   restart_age  = .false.
+   tr_FY        = .true.
+   restart_FY   = .false.
+   tr_lvl       = .true.
+   restart_lvl  = .false.
+   tr_pond_cesm = .false.
+   restart_pond_cesm = .false.
+   tr_pond_topo = .false.
+   restart_pond_topo = .false.
+   tr_pond_lvl  = $tr_pond_lvl
+   restart_pond_lvl  = $restart_pond_lvl
+   tr_aero      = .false.
+   restart_aero = .false.
+   tr_fsd       = .false.
+   restart_fsd  = .false.
 /
 
 &thermo_nml
-    kitd              = 1
-  , ktherm            = 1
-  , conduct           = 'MU71'
-  , a_rapid_mode      =  0.5e-3
-  , Rac_rapid_mode    =    10.0
-  , aspect_rapid_mode =     1.0
-  , dSdt_slow_mode    = -5.0e-8
-  , phi_c_slow_mode   =    0.05
-  , phi_i_mushy       =    0.85
+   kitd              = 1
+   ktherm            = 1
+   conduct           = 'MU71'
+   a_rapid_mode      =  0.5e-3
+   Rac_rapid_mode    =    10.0
+   aspect_rapid_mode =     1.0
+   dSdt_slow_mode    = -5.0e-8
+   phi_c_slow_mode   =    0.05
+   phi_i_mushy       =    0.85
+   sw_redist         = .true.
 /
 
 &dynamics_nml
-    kdyn            = 1
-  , ndte            = 120
-  , revised_evp     = .false.
-  , advection       = 'remap'
-  , kstrength       = 1
-  , krdg_partic     = 1
-  , krdg_redist     = 1
-  , mu_rdg          = 3
+   kdyn            = 1
+   ndte            = 120
+   revised_evp     = .false.
+   kevp_kernel     = 0
+   brlx            = 300.0
+   arlx            = 300.0
+   ssh_stress      = 'coupled'
+   advection       = 'remap'
+   kstrength       = 1
+   krdg_partic     = 1
+   krdg_redist     = 1
+   mu_rdg          = 3
+   Cf              = 17.
+   Ktens           = 0.
+   e_ratio         = 2.
+   basalstress     = .false.
+   k1              = 8.
+   coriolis        = 'latitude'
+   kridge          = 1
+   ktransport      = 1
 /
 
 &shortwave_nml
-    shortwave       = 'dEdd'
-  , albedo_type     = 'default'
-  , albicev         = 0.78
-  , albicei         = 0.36
-  , albsnowv        = 0.98
-  , albsnowi        = 0.70 
-  , ahmax           = 0.3
-  , R_ice           = 0.
-  , R_pnd           = 0.
-  , R_snw           = 1.5
-  , dT_mlt          = 1.5
-  , rsnw_mlt        = 1500.
+   shortwave       = 'dEdd'
+   albedo_type     = 'default'
+   albicev         = 0.78
+   albicei         = 0.36
+   albsnowv        = 0.98
+   albsnowi        = 0.70 
+   ahmax           = 0.3
+   R_ice           = 0.
+   R_pnd           = 0.
+   R_snw           = 1.5
+   dT_mlt          = 1.5
+   rsnw_mlt        = 1500.
+   kalg            = 0.0
 /
 
 &ponds_nml
-    hp1             = 0.01
-  , hs0             = 0.
-  , hs1             = 0.03
-  , dpscale         = 1.0e-3
-  , frzpnd          = 'hlid'
-  , snowinfil       = .true.
-  , rfracmin        = 0.15
-  , rfracmax        = 1.
-  , pndaspect       = 0.8
+   hp1             = 0.01
+   hs0             = 0.
+   hs1             = 0.03
+   dpscale         = 1.0e-3
+   frzpnd          = 'hlid'
+   rfracmin        = 0.15
+   rfracmax        = 1.
+   pndaspect       = 0.8
 /
 
-&zbgc_nml
-    tr_brine        = .false.
-  , restart_hbrine  = .false.
-  , skl_bgc         = .false.
-  , bgc_flux_type   = 'Jin2006'
-  , restart_bgc     = .false.
-  , restore_bgc     = .false.
-  , bgc_data_dir    = 'unknown_bgc_data_dir'
-  , sil_data_type   = 'default'
-  , nit_data_type   = 'default'
-  , tr_bgc_C_sk     = .false.
-  , tr_bgc_chl_sk   = .false.
-  , tr_bgc_Am_sk    = .false.
-  , tr_bgc_Sil_sk   = .false.
-  , tr_bgc_DMSPp_sk = .false.
-  , tr_bgc_DMSPd_sk = .false.
-  , tr_bgc_DMS_sk   = .false.
-  , phi_snow        = 0.5
+&forcing_nml
+   formdrag        = .false.
+   atmbndy         = 'default'
+   calc_strair     = .true.
+   calc_Tsfc       = .true.
+   highfreq        = .false.
+   natmiter        = 5
+   ustar_min       = 0.0005
+   emissivity      = 0.95
+   fbot_xfer_type  = 'constant'
+   update_ocn_f    = $FRAZIL_FWSALT
+   l_mpond_fresh   = .false.
+   tfrz_option     = 'linear_salt'
+   restart_coszen  = .true.
+   oceanmixed_ice  = .false.
+   wave_spec_type  = 'none'
+   wave_spec_file  = 'unknown_wave_spec_file'
+   nfreq           = 25
+   restore_ice     = .false.
+   restore_ocn     = .false.
+   trestore        =  90
+   precip_units    = 'mm_per_month'
+   default_season  = 'winter'
+   atm_data_type   = 'ncar'
+   ocn_data_type   = 'default'
+   bgc_data_type   = 'default'
+   fe_data_type    = 'default'
+   ice_data_type   = 'default'
+   fyear_init      = 1997
+   ycycle          = 1
+   atm_data_format = 'nc'
+   atm_data_dir    = './INPUT/gx3_forcing_fields.nc'
+   bgc_data_dir    = 'unknown_bgc_data_dir'
+   ocn_data_format = 'bin'
+   ocn_data_dir    = '/unknown_ocn_data_dir'
+   oceanmixed_file = 'unknown_oceanmixed_file'
 /
 
-&forcing_nml
-    formdrag        = .false.
-  , atmbndy         = 'default'
-  , fyear_init      = 1997
-  , ycycle          = 1
-  , atm_data_format = 'bin'
-  , atm_data_type   = 'none'
-  , atm_data_dir    = '/home/Fei.Liu/noscrub/lanl_cice_data/'
-  , calc_strair     = .true.
-  , calc_Tsfc       = .true.
-  , precip_units    = 'mm_per_month'
-  , ustar_min       = 0.0005
-  , update_ocn_f    = $FRAZIL_FWSALT
-  , oceanmixed_ice  = .false.
-  , ocn_data_format = 'bin'
-  , sss_data_type   = 'default'
-  , sst_data_type   = 'default'
-  , ocn_data_dir    = 'unknown_ocn_data_dir'
-  , oceanmixed_file = 'unknown_oceanmixed_file'
-  , restore_sst     = .false.
-  , trestore        =  90
-  , restore_ice     = .false.
+&domain_nml
+   nprocs = $ICEPETS
+   nx_global         = $NX_GLB
+   ny_global         = $NY_GLB
+   block_size_x      = $(( 2 * ( $NX_GLB / $ICEPETS ) ))
+   block_size_y      = $(( $NY_GLB / 2 ))
+   max_blocks        = -1
+   processor_shape   = 'slenderX2'
+   distribution_type = 'cartesian'
+   distribution_wght = 'latitude'
+   ew_boundary_type  = 'cyclic'
+   ns_boundary_type  = 'tripole'
+   maskhalo_dyn      = .false.
+   maskhalo_remap    = .false.
+   maskhalo_bound    = .false.
+/
+
+&zbgc_nml
+   tr_brine        = .false.
+   restart_hbrine  = .false.
+   tr_zaero        = .false.
+   modal_aero      = .false.
+   skl_bgc         = .false.
+   z_tracers       = .false.
+   dEdd_algae      = .false.
+   solve_zbgc      = .false.
+   bgc_flux_type   = 'Jin2006'
+   restore_bgc     = .false.
+   restart_bgc     = .false.
+   scale_bgc       = .false.
+   solve_zsal      = .false.
+   restart_zsal    = .false.
+   tr_bgc_Nit      = .true.
+   tr_bgc_C        = .true.
+   tr_bgc_chl      = .false.
+   tr_bgc_Am       = .true.
+   tr_bgc_Sil      = .true.
+   tr_bgc_DMS      = .false.
+   tr_bgc_PON      = .true.
+   tr_bgc_hum      = .true.
+   tr_bgc_DON      = .false.
+   tr_bgc_Fe       = .true. 
+   grid_o          = 0.006
+   grid_o_t        = 0.006
+   l_sk            = 0.024
+   grid_oS         = 0.0
+   l_skS           = 0.028
+   phi_snow        = -0.3
+   initbio_frac    = 0.8
+   frazil_scav     = 0.8  
+   ratio_Si2N_diatoms = 1.8                         
+   ratio_Si2N_sp      = 0.0
+   ratio_Si2N_phaeo   = 0.0
+   ratio_S2N_diatoms  = 0.03  
+   ratio_S2N_sp       = 0.03 
+   ratio_S2N_phaeo    = 0.03
+   ratio_Fe2C_diatoms = 0.0033
+   ratio_Fe2C_sp      = 0.0033
+   ratio_Fe2C_phaeo   = 0.1
+   ratio_Fe2N_diatoms = 0.023 
+   ratio_Fe2N_sp      = 0.023
+   ratio_Fe2N_phaeo   = 0.7
+   ratio_Fe2DON       = 0.023
+   ratio_Fe2DOC_s     = 0.1
+   ratio_Fe2DOC_l     = 0.033
+   fr_resp            = 0.05
+   tau_min            = 5200.0
+   tau_max            = 173000.0
+   algal_vel          = 0.0000000111
+   R_dFe2dust         = 0.035
+   dustFe_sol         = 0.005
+   chlabs_diatoms     = 0.03
+   chlabs_sp          = 0.01
+   chlabs_phaeo       = 0.05
+   alpha2max_low_diatoms = 0.8
+   alpha2max_low_sp      = 0.67
+   alpha2max_low_phaeo   = 0.67
+   beta2max_diatoms   = 0.018
+   beta2max_sp        = 0.0025
+   beta2max_phaeo     = 0.01
+   mu_max_diatoms     = 1.44
+   mu_max_sp          = 0.851
+   mu_max_phaeo       = 0.851
+   grow_Tdep_diatoms  = 0.06
+   grow_Tdep_sp       = 0.06
+   grow_Tdep_phaeo    = 0.06
+   fr_graze_diatoms   = 0.0
+   fr_graze_sp        = 0.1
+   fr_graze_phaeo     = 0.1
+   mort_pre_diatoms   = 0.007
+   mort_pre_sp        = 0.007
+   mort_pre_phaeo     = 0.007
+   mort_Tdep_diatoms  = 0.03
+   mort_Tdep_sp       = 0.03
+   mort_Tdep_phaeo    = 0.03
+   k_exude_diatoms    = 0.0
+   k_exude_sp         = 0.0
+   k_exude_phaeo      = 0.0
+   K_Nit_diatoms      = 1.0
+   K_Nit_sp           = 1.0
+   K_Nit_phaeo        = 1.0
+   K_Am_diatoms       = 0.3
+   K_Am_sp            = 0.3
+   K_Am_phaeo         = 0.3
+   K_Sil_diatoms      = 4.0
+   K_Sil_sp           = 0.0
+   K_Sil_phaeo        = 0.0
+   K_Fe_diatoms       = 1.0
+   K_Fe_sp            = 0.2
+   K_Fe_phaeo         = 0.1
+   f_don_protein      = 0.6
+   kn_bac_protein     = 0.03
+   f_don_Am_protein   = 0.25
+   f_doc_s            = 0.4
+   f_doc_l            = 0.4
+   f_exude_s          = 1.0
+   f_exude_l          = 1.0
+   k_bac_s            = 0.03
+   k_bac_l            = 0.03
+   T_max              = 0.0
+   fsal               = 1.0
+   op_dep_min         = 0.1
+   fr_graze_s         = 0.5
+   fr_graze_e         = 0.5
+   fr_mort2min        = 0.5
+   fr_dFe             = 0.3
+   k_nitrif           = 0.0
+   t_iron_conv        = 3065.0
+   max_loss           = 0.9
+   max_dfe_doc1       = 0.2
+   fr_resp_s          = 0.75
+   y_sk_DMS           = 0.5
+   t_sk_conv          = 3.0
+   t_sk_ox            = 10.0
+   algaltype_diatoms  = 0.0
+   algaltype_sp       = 0.5
+   algaltype_phaeo    = 0.5
+   nitratetype        = -1.0
+   ammoniumtype       = 1.0
+   silicatetype       = -1.0
+   dmspptype          = 0.5
+   dmspdtype          = -1.0
+   humtype            = 1.0
+   doctype_s          = 0.5
+   doctype_l          = 0.5
+   dontype_protein    = 0.5
+   fedtype_1          = 0.5
+   feptype_1          = 0.5
+   zaerotype_bc1      = 1.0
+   zaerotype_bc2      = 1.0
+   zaerotype_dust1    = 1.0
+   zaerotype_dust2    = 1.0
+   zaerotype_dust3    = 1.0
+   zaerotype_dust4    = 1.0
+   ratio_C2N_diatoms  = 7.0
+   ratio_C2N_sp       = 7.0
+   ratio_C2N_phaeo    = 7.0
+   ratio_chl2N_diatoms= 2.1
+   ratio_chl2N_sp     = 1.1
+   ratio_chl2N_phaeo  = 0.84
+   F_abs_chl_diatoms  = 2.0
+   F_abs_chl_sp       = 4.0
+   F_abs_chl_phaeo    = 5.0
+   ratio_C2N_proteins = 7.0
 /
 
 &icefields_nml
-    f_tmask         = .true.
-  , f_tarea         = .true.
-  , f_uarea         = .true.
-  , f_dxt           = .false.
-  , f_dyt           = .false.
-  , f_dxu           = .false.
-  , f_dyu           = .false.
-  , f_HTN           = .false.
-  , f_HTE           = .false.
-  , f_ANGLE         = .true.
-  , f_ANGLET        = .true.
-  , f_NCAT          = .true.
-  , f_VGRDi         = .false.
-  , f_VGRDs         = .false.
-  , f_VGRDb         = .false.
-  , f_bounds        = .false.
-  , f_aice          = 'mdhxx' 
-  , f_hi            = 'mdhxx'
-  , f_hs            = 'mdhxx' 
-  , f_Tsfc          = 'mdhxx' 
-  , f_sice          = 'x' 
-  , f_uvel          = 'mdhxx' 
-  , f_vvel          = 'mdhxx' 
-  , f_fswdn         = 'mdhxx' 
-  , f_flwdn         = 'mdhxx'
-  , f_snow          = 'mdhxx' 
-  , f_snow_ai       = 'xxxxx' 
-  , f_rain          = 'mdhxx' 
-  , f_rain_ai       = 'xxxxx' 
-  , f_sst           = 'mdhxx' 
-  , f_sss           = 'mdhxx' 
-  , f_uocn          = 'x' 
-  , f_vocn          = 'x' 
-  , f_frzmlt        = 'x'
-  , f_fswfac        = 'x'
-  , f_fswabs        = 'mdhxx' 
-  , f_fswabs_ai     = 'xxxxx' 
-  , f_albsni        = 'mdhxx' 
-  , f_alvdr         = 'x'
-  , f_alidr         = 'x'
-  , f_albice        = 'x'
-  , f_albsno        = 'x'
-  , f_albpnd        = 'x'
-  , f_coszen        = 'x'
-  , f_flat          = 'mdhxx' 
-  , f_flat_ai       = 'xxxxx' 
-  , f_fsens         = 'mdhxx' 
-  , f_fsens_ai      = 'xxxxx' 
-  , f_flwup         = 'mdhxx' 
-  , f_flwup_ai      = 'xxxxx' 
-  , f_evap          = 'mdhxx' 
-  , f_evap_ai       = 'xxxxx' 
-  , f_Tair          = 'mdhxx' 
-  , f_Tref          = 'mdhxx' 
-  , f_Qref          = 'x'
-  , f_congel        = 'x' 
-  , f_frazil        = 'x' 
-  , f_snoice        = 'x' 
-  , f_dsnow         = 'x' 
-  , f_melts         = 'x'
-  , f_meltt         = 'mdhxx'
-  , f_meltb         = 'mdhxx'
-  , f_meltl         = 'x'
-  , f_fresh         = 'x'
-  , f_fresh_ai      = 'xxxxx'
-  , f_fsalt         = 'x'
-  , f_fsalt_ai      = 'xxxxx'
-  , f_fhocn         = 'x' 
-  , f_fhocn_ai      = 'xxxxx' 
-  , f_fswthru       = 'x' 
-  , f_fswthru_ai    = 'xxxxx' 
-  , f_fsurf_ai      = 'xxxxx'
-  , f_fcondtop_ai   = 'xxxxx'
-  , f_fmeltt_ai     = 'xxxxx' 
-  , f_strairx       = 'mdhxx' 
-  , f_strairy       = 'mdhxx' 
-  , f_strtltx       = 'x' 
-  , f_strtlty       = 'x' 
-  , f_strcorx       = 'x' 
-  , f_strcory       = 'x' 
-  , f_strocnx       = 'x' 
-  , f_strocny       = 'x' 
-  , f_strintx       = 'x' 
-  , f_strinty       = 'x'
-  , f_strength      = 'x'
-  , f_divu          = 'mdhxx'
-  , f_shear         = 'mdhxx'
-  , f_sig1          = 'x' 
-  , f_sig2          = 'x' 
-  , f_dvidtt        = 'mdhxx' 
-  , f_dvidtd        = 'mdhxx' 
-  , f_daidtt        = 'mdhxx'
-  , f_daidtd        = 'mdhxx' 
-  , f_mlt_onset     = 'mdhxx'
-  , f_frz_onset     = 'mdhxx'
-  , f_hisnap        = 'x'
-  , f_aisnap        = 'x'
-  , f_trsig         = 'x'
-  , f_icepresent    = 'x'
-  , f_iage          = 'x'
-  , f_FY            = 'x'
-  , f_aicen         = 'xxxxx'
-  , f_vicen         = 'xxxxx'
-  , f_Tinz          = 'x'
-  , f_Sinz          = 'x'
-  , f_Tsnz          = 'x'
-  , f_fsurfn_ai     = 'xxxxx'
-  , f_fcondtopn_ai  = 'xxxxx'
-  , f_fmelttn_ai    = 'xxxxx'
-  , f_flatn_ai      = 'xxxxx'
-  , f_s11           = 'mdhxx'
-  , f_s12           = 'mdhxx'
-  , f_s22           = 'mdhxx'
-  , f_yieldstress11 = 'mdhxx'
-  , f_yieldstress12 = 'mdhxx'
-  , f_yieldstress22 = 'mdhxx'
+   f_tmask         = .true.
+   f_blkmask       = .true.
+   f_tarea         = .true.
+   f_uarea         = .true.
+   f_dxt           = .false.
+   f_dyt           = .false.
+   f_dxu           = .false.
+   f_dyu           = .false.
+   f_HTN           = .false.
+   f_HTE           = .false.
+   f_ANGLE         = .true.
+   f_ANGLET        = .true.
+   f_NCAT          = .true.
+   f_VGRDi         = .false.
+   f_VGRDs         = .false.
+   f_VGRDb         = .false.
+   f_VGRDa         = .true.
+   f_bounds        = .false.
+   f_aice          = 'mdhxx' 
+   f_hi            = 'mdhxx'
+   f_hs            = 'mdhxx' 
+   f_Tsfc          = 'mdhxx' 
+   f_sice          = 'mdhxx' 
+   f_uvel          = 'mdhxx' 
+   f_vvel          = 'mdhxx'
+   f_uatm          = 'mdhxx'
+   f_vatm          = 'mdhxx'
+   f_fswdn         = 'mdhxx' 
+   f_flwdn         = 'mdhxx'
+   f_snowfrac      = 'x'
+   f_snow          = 'mdhxx' 
+   f_snow_ai       = 'x' 
+   f_rain          = 'mdhxx' 
+   f_rain_ai       = 'x' 
+   f_sst           = 'mdhxx' 
+   f_sss           = 'mdhxx' 
+   f_uocn          = 'mdhxx' 
+   f_vocn          = 'mdhxx' 
+   f_frzmlt        = 'mdhxx'
+   f_fswfac        = 'mdhxx'
+   f_fswint_ai     = 'x'
+   f_fswabs        = 'mdhxx' 
+   f_fswabs_ai     = 'x' 
+   f_albsni        = 'mdhxx' 
+   f_alvdr         = 'mdhxx'
+   f_alidr         = 'mdhxx'
+   f_alvdf         = 'mdhxx'
+   f_alidf         = 'mdhxx'
+   f_alvdr_ai      = 'x'
+   f_alidr_ai      = 'x'
+   f_alvdf_ai      = 'x'
+   f_alidf_ai      = 'x'
+   f_albice        = 'mdhxx'
+   f_albsno        = 'mdhxx'
+   f_albpnd        = 'mdhxx'
+   f_coszen        = 'mdhxx'
+   f_flat          = 'mdhxx' 
+   f_flat_ai       = 'x' 
+   f_fsens         = 'mdhxx' 
+   f_fsens_ai      = 'x' 
+   f_fswup         = 'x'
+   f_flwup         = 'mdhxx' 
+   f_flwup_ai      = 'x' 
+   f_evap          = 'mdhxx' 
+   f_evap_ai       = 'x' 
+   f_Tair          = 'mdhxx' 
+   f_Tref          = 'mdhxx' 
+   f_Qref          = 'mdhxx'
+   f_congel        = 'mdhxx' 
+   f_frazil        = 'mdhxx' 
+   f_snoice        = 'mdhxx' 
+   f_dsnow         = 'mdhxx' 
+   f_melts         = 'mdhxx'
+   f_meltt         = 'mdhxx'
+   f_meltb         = 'mdhxx'
+   f_meltl         = 'mdhxx'
+   f_fresh         = 'mdhxx'
+   f_fresh_ai      = 'x'
+   f_fsalt         = 'mdhxx'
+   f_fsalt_ai      = 'x'
+   f_fbot          = 'mdhxx'
+   f_fhocn         = 'mdhxx' 
+   f_fhocn_ai      = 'x' 
+   f_fswthru       = 'mdhxx' 
+   f_fswthru_ai    = 'x' 
+   f_fsurf_ai      = 'x'
+   f_fcondtop_ai   = 'x'
+   f_fmeltt_ai     = 'x' 
+   f_strairx       = 'mdhxx' 
+   f_strairy       = 'mdhxx' 
+   f_strtltx       = 'x' 
+   f_strtlty       = 'x' 
+   f_strcorx       = 'x' 
+   f_strcory       = 'x' 
+   f_strocnx       = 'mdhxx' 
+   f_strocny       = 'mdhxx' 
+   f_strintx       = 'x' 
+   f_strinty       = 'x'
+   f_taubx         = 'x'
+   f_tauby         = 'x'
+   f_strength      = 'x'
+   f_divu          = 'x'
+   f_shear         = 'x'
+   f_sig1          = 'x' 
+   f_sig2          = 'x'
+   f_sigP          = 'x' 
+   f_dvidtt        = 'x' 
+   f_dvidtd        = 'x' 
+   f_daidtt        = 'x'
+   f_daidtd        = 'x' 
+   f_dagedtt       = 'x'
+   f_dagedtd       = 'x'
+   f_mlt_onset     = 'mdhxx'
+   f_frz_onset     = 'mdhxx'
+   f_hisnap        = 'x'
+   f_aisnap        = 'x'
+   f_trsig         = 'x'
+   f_icepresent    = 'x'
+   f_iage          = 'x'
+   f_FY            = 'x'
+   f_aicen         = 'x'
+   f_vicen         = 'x'
+   f_vsnon         = 'x'
+   f_snowfracn     = 'x'
+   f_keffn_top     = 'x'
+   f_Tinz          = 'x'
+   f_Sinz          = 'x'
+   f_Tsnz          = 'x'
+   f_fsurfn_ai     = 'x'
+   f_fcondtopn_ai  = 'x'
+   f_fmelttn_ai    = 'x'
+   f_flatn_ai      = 'x'
+   f_fsensn_ai     = 'x'
 /
 
 &icefields_mechred_nml
-    f_alvl         = 'mdhxx'
-  , f_vlvl         = 'mdhxx'
-  , f_ardg         = 'mdhxx'
-  , f_vrdg         = 'mdhxx'
-  , f_dardg1dt     = 'x'
-  , f_dardg2dt     = 'x'
-  , f_dvirdgdt     = 'x'
-  , f_opening      = 'mdhxx'
-  , f_ardgn        = 'xxxxx'
-  , f_vrdgn        = 'xxxxx'
-  , f_dardg1ndt    = 'x'
-  , f_dardg2ndt    = 'x'
-  , f_dvirdgndt    = 'x'
-  , f_krdgn        = 'x'
-  , f_aparticn     = 'x'
-  , f_aredistn     = 'x'
-  , f_vredistn     = 'x'
-  , f_araftn       = 'x'
-  , f_vraftn       = 'x'
+   f_alvl         = 'x'
+   f_vlvl         = 'x'
+   f_ardg         = 'x'
+   f_vrdg         = 'x'
+   f_dardg1dt     = 'x'
+   f_dardg2dt     = 'x'
+   f_dvirdgdt     = 'x'
+   f_opening      = 'x'
+   f_ardgn        = 'x'
+   f_vrdgn        = 'x'
+   f_dardg1ndt    = 'x'
+   f_dardg2ndt    = 'x'
+   f_dvirdgndt    = 'x'
+   f_krdgn        = 'x'
+   f_aparticn     = 'x'
+   f_aredistn     = 'x'
+   f_vredistn     = 'x'
+   f_araftn       = 'x'
+   f_vraftn       = 'x'
 /
 
 &icefields_pond_nml
-    f_apondn       = 'xxxxx'
-  , f_apeffn       = 'xxxxx'
-  , f_hpondn       = 'xxxxx'
-  , f_apond        = 'mdhxx'
-  , f_hpond        = 'mdhxx'
-  , f_ipond        = 'mdhxx'
-  , f_apeff        = 'mdhxx'
-  , f_apond_ai     = 'xxxxx'
-  , f_hpond_ai     = 'xxxxx'
-  , f_ipond_ai     = 'xxxxx'
-  , f_apeff_ai     = 'xxxxx'
+   f_apondn       = 'x'
+   f_apeffn       = 'x'
+   f_hpondn       = 'x'
+   f_apond        = 'mdhxx'
+   f_hpond        = 'mdhxx'
+   f_ipond        = 'mdhxx'
+   f_apeff        = 'mdhxx'
+   f_apond_ai     = 'x'
+   f_hpond_ai     = 'x'
+   f_ipond_ai     = 'x'
+   f_apeff_ai     = 'x'
 /
 
 &icefields_bgc_nml
-    f_faero_atm    = 'x'
-  , f_faero_ocn    = 'x'
-  , f_aero         = 'x'
-  , f_fNO          = 'x'
-  , f_fNO_ai       = 'x'
-  , f_fNH          = 'x'
-  , f_fNH_ai       = 'x'
-  , f_fN           = 'x'
-  , f_fN_ai        = 'x'
-  , f_fSil         = 'x'
-  , f_fSil_ai      = 'x'
-  , f_bgc_N_sk     = 'x'
-  , f_bgc_C_sk     = 'x'
-  , f_bgc_chl_sk   = 'x'
-  , f_bgc_Nit_sk   = 'x'
-  , f_bgc_Am_sk    = 'x'
-  , f_bgc_Sil_sk   = 'x'
-  , f_bgc_DMSPp_sk = 'x'
-  , f_bgc_DMSPd_sk = 'x'
-  , f_bgc_DMS_sk   = 'x'
-  , f_bgc_Nit_ml   = 'x'
-  , f_bgc_Am_ml    = 'x'
-  , f_bgc_Sil_ml   = 'x'  
-  , f_bgc_DMSP_ml  = 'x'
-  , f_bTin         = 'x'
-  , f_bphi         = 'x' 
-  , f_fbri         = 'x'    
-  , f_hbri         = 'x'
-  , f_grownet      = 'x'
-  , f_PPnet        = 'x'
+   f_faero_atm    = 'x'
+   f_faero_ocn    = 'x'
+   f_aero         = 'x'
+   f_fbio         = 'x'
+   f_fbio_ai      = 'x'
+   f_zaero        = 'x'
+   f_bgc_S        = 'x'
+   f_bgc_N        = 'x'
+   f_bgc_C        = 'x'
+   f_bgc_DOC      = 'x'
+   f_bgc_DIC      = 'x'
+   f_bgc_chl      = 'x'
+   f_bgc_Nit      = 'x'
+   f_bgc_Am       = 'x'
+   f_bgc_Sil      = 'x'
+   f_bgc_DMSPp    = 'x'
+   f_bgc_DMSPd    = 'x'
+   f_bgc_DMS      = 'x'
+   f_bgc_DON      = 'x'
+   f_bgc_Fe       = 'x'
+   f_bgc_hum      = 'x'
+   f_bgc_PON      = 'x'
+   f_bgc_ml       = 'x'
+   f_upNO         = 'x'
+   f_upNH         = 'x'
+   f_bTin         = 'x'
+   f_bphi         = 'x'
+   f_iDi          = 'x'
+   f_iki          = 'x'
+   f_fbri         = 'x'
+   f_hbri         = 'x'
+   f_zfswin       = 'x'
+   f_bionet       = 'x'
+   f_biosnow      = 'x'
+   f_grownet      = 'x'
+   f_PPnet        = 'x'
+   f_algalpeak    = 'x'
+   f_zbgc_frac    = 'x'
 /
 
 &icefields_drag_nml
-    f_drag         = 'mdhxx'
-  , f_Cdn_atm      = 'mdhxx'
-  , f_Cdn_ocn      = 'mdhxx'
+   f_drag         = 'x'
+   f_Cdn_atm      = 'x'
+   f_Cdn_ocn      = 'x'
+/
+
+&icefields_fsd_nml
+   f_fsdrad       = 'x'
+   f_fsdperim     = 'x'
+   f_afsd         = 'x'
+   f_afsdn        = 'x'
+   f_dafsd_newi   = 'x'
+   f_dafsd_latg   = 'x'
+   f_dafsd_latm   = 'x'
+   f_dafsd_wave   = 'x'
+   f_dafsd_weld   = 'x'
+   f_wave_sig_ht  = 'x'
+   f_aice_ww      = 'x'
+   f_diam_ww      = 'x'
+   f_hice_ww      = 'x'
 /
 eof
 
-}
+} 
diff --git a/workflow/config/fcst.yaml b/workflow/config/fcst.yaml
index 36a8c7363c..93975c6de6 100644
--- a/workflow/config/fcst.yaml
+++ b/workflow/config/fcst.yaml
@@ -64,6 +64,12 @@ config_fcst:
     export NTHREADS_FV3="{doc.fv3_gfs_settings.layout.nth}"
     export memory_fcst="1024M"
 
+    #######################################################################
+    # COUPLING COMPONENTS
+    export OCN_model="{doc.ocn_settings.MOD}"
+    export ICE_model="{doc.ice_settings.MOD}"
+    export WAV_model="{doc.wave_settings.MOD}"
+
     #######################################################################
     # PE Mapping
     export ATMPETS="{doc.fv3_gfs_settings.ATMPETS}"
@@ -73,6 +79,10 @@ config_fcst:
     export WAVPETS="{doc.wave_settings.WAVPETS}"
     #
     #######################################################################    
+    # CICE parameters
+    export NX_GLB="{doc.ice_settings.NX_GLB}"
+    export NY_GLB="{doc.ice_settings.NY_GLB}"
+    export MESHICE="mesh.mx{doc.ice_settings.ICERES}.nc"
 
     if [[ "$machine" == "WCOSS_C" ]]; then
         export memory_fcst="1024M"
diff --git a/workflow/defaults/ice.yaml b/workflow/defaults/ice.yaml
index 2feccddf20..6e66fd76af 100644
--- a/workflow/defaults/ice.yaml
+++ b/workflow/defaults/ice.yaml
@@ -1,7 +1,9 @@
 # default settings for fv3 core
 
 ice_defaults: &ice_defaults
-  MOD: cice5
+  MOD: cice6
   ICERES: "025"
   CPL_ICEIC: CPC
   ICEPETS: 48
+  NX_GLB: 1440
+  NY_GLB: 1080