DOI: 10.5281/zenodo.1308226
title: Cice6 Simulations: Quality Control For Community Based Sea Ice Model Development

description: CICE6 simulations: Quality Control for Community Based Sea Ice Model Development<br>
-----------------------------------------------------------------------------------

Data  from the CICE6 simulations by Elizabeth C. Hunke at Los Alamos National Laboratory described in the manuscript:

Roberts, Hunke, Allard, Bailey, Craig, Lemieux and Turner (2018): Quality Control for Community Based Sea Ice Model Development.

It was produced using the CICE sea ice model version 6.0.0.alpha with the baseline namelist:

&setup_nml<br>
    days_per_year  = 365<br>
    use_leap_years = .false.<br>
    year_init      = 1990<br>
    istep0         = 0<br>
    dt             = 3600.0<br>
    npt            = 43800<br>
    ndtd           = 1<br>
    runtype        = 'continue'<br>
    ice_ic         = '/usr/projects/climate/eclare/DATA/Consortium/CICE_data/ic/gx1/iced_gx1_v5.nc'<br>
    restart        = .true.<br>
    restart_ext    = .false.<br>
    use_restart_time = .true.<br>
    restart_format = 'nc'<br>
    lcdf64         = .false.<br>
    restart_dir    = './restart/'<br>
    restart_file   = 'iced'<br>
    pointer_file   = './ice.restart_file'<br>
    dumpfreq       = 'm'<br>
    dumpfreq_n     = 1<br>
    dump_last      = .false.<br>
    bfbflag        = .false.<br>
    diagfreq       = 24<br>
    diag_type      = 'stdout'<br>
    diag_file      = 'ice_diag.d'<br>
    print_global   = .true.<br>
    print_points   = .true.<br>
    latpnt(1)      =  90.<br>
    lonpnt(1)      =  0.<br>
    latpnt(2)      = -65.<br>
    lonpnt(2)      = -45.<br>
    dbug           = .false.<br>
    histfreq       = 'm','d','x','x','x'<br>
    histfreq_n     =  1 , 1 , 1 , 1 , 1<br>
    hist_avg       = .true.<br>
    history_dir    = './history/'<br>
    history_file   = 'iceh'<br>
    write_ic       = .true.<br>
    incond_dir     = './history/'<br>
    incond_file    = 'iceh_ic'<br>
/

&grid_nml<br>
    grid_format  = 'bin'<br>
    grid_type    = 'displaced_pole'<br>
    grid_file    = '/usr/projects/climate/eclare/DATA/Consortium/CICE_data/grid/gx1/grid_gx1.bin'<br>
    kmt_file     = '/usr/projects/climate/eclare/DATA/Consortium/CICE_data/grid/gx1/kmt_gx1.bin'<br>
    gridcpl_file = 'unknown_gridcpl_file'<br>
    kcatbound    = 0<br>
/

&domain_nml<br>
    nprocs = 32<br>
    processor_shape   = 'slenderX2'<br>
    distribution_type = 'cartesian'<br>
    distribution_wght = 'latitude'<br>
    ew_boundary_type  = 'cyclic'<br>
    ns_boundary_type  = 'open'<br>
    maskhalo_dyn      = .false.<br>
    maskhalo_remap    = .false.<br>
    maskhalo_bound    = .false.<br>
/

&tracer_nml<br>
    tr_iage      = .true.<br>
    restart_age  = .false.<br>
    tr_FY        = .true.<br>
    restart_FY   = .false.<br>
    tr_lvl       = .true.<br>
    restart_lvl  = .false.<br>
    tr_pond_cesm = .false.<br>
    restart_pond_cesm = .false.<br>
    tr_pond_topo = .false.<br>
    restart_pond_topo = .false.<br>
    tr_pond_lvl  = .true.<br>
    restart_pond_lvl  = .false.<br>
    tr_aero      = .false.<br>
    restart_aero = .false.<br>
/

&thermo_nml<br>
    kitd              = 1<br>
    ktherm            = 2<br>
    conduct           = 'bubbly'<br>
    a_rapid_mode      =  0.5e-3<br>
    Rac_rapid_mode    =    10.0<br>
    aspect_rapid_mode =     1.0<br>
    dSdt_slow_mode    = -5.0e-8<br>
    phi_c_slow_mode   =    0.05<br>
    phi_i_mushy       =    0.85<br>
/

&dynamics_nml<br>
    kdyn            = 1<br>
    ndte            = 120<br>
    revised_evp     = .false.<br>
    advection       = 'remap'<br>
    kstrength       = 1<br>
    krdg_partic     = 1<br>
    krdg_redist     = 1<br>
    mu_rdg          = 3<br>
    Cf              = 17.<br>
    Ktens           = 0.<br>
    e_ratio         = 2.<br>
    basalstress     = .false.<br>
/

&shortwave_nml<br>
    shortwave       = 'dEdd'<br>
    albedo_type     = 'default'<br>
    albicev         = 0.78<br>
    albicei         = 0.36<br>
    albsnowv        = 0.98<br>
    albsnowi        = 0.70 <br>
    ahmax           = 0.3<br>
    R_ice           = 0.<br>
    R_pnd           = 0.<br>
    R_snw           = 1.5<br>
    dT_mlt          = 1.5<br>
    rsnw_mlt        = 1500.<br>
    kalg            = 0.6<br>
/

&ponds_nml<br>
    hp1             = 0.01<br>
    hs0             = 0.<br>
    hs1             = 0.03<br>
    dpscale         = 1.e-3<br>
    frzpnd          = 'hlid'<br>
    rfracmin        = 0.15<br>
    rfracmax        = 1.<br>
    pndaspect       = 0.8<br>
/

&zbgc_nml<br>
    tr_brine        = .false.<br>
    restart_hbrine  = .false.<br>
    tr_zaero        = .false.<br>
    modal_aero      = .false.<br>
    skl_bgc         = .false.<br>
    z_tracers       = .false.<br>
    dEdd_algae      = .false.<br>
    solve_zbgc      = .false.<br>
    bgc_flux_type   = 'Jin2006'<br>
    restore_bgc     = .false.<br>
    restart_bgc     = .false.<br>
    scale_bgc       = .false.<br>
    solve_zsal      = .false.<br>
    restart_zsal    = .false.<br>
    bgc_data_dir    = '/uknown_bgc_data_dir'<br>
    sil_data_type   = 'default'<br>
    nit_data_type   = 'default'<br>
    fe_data_type    = 'default'<br>
    tr_bgc_Nit      = .true.<br>
    tr_bgc_C        = .true.<br>
    tr_bgc_chl      = .false.<br>
    tr_bgc_Am       = .true.<br>
    tr_bgc_Sil      = .true.<br>
    tr_bgc_DMS      = .false.<br>
    tr_bgc_PON      = .true.<br>
    tr_bgc_hum      = .true.<br>
    tr_bgc_DON      = .false.<br>
    tr_bgc_Fe       = .true. <br>
    grid_o          = 0.006<br>
    grid_o_t        = 0.006<br>
    l_sk            = 0.024<br>
    grid_oS         = 0.0<br>
    l_skS           = 0.028<br>
    phi_snow        = -0.3<br>
    initbio_frac    = 0.8<br>
    frazil_scav     = 0.8  <br>
    ratio_Si2N_diatoms = 1.8                         <br>
    ratio_Si2N_sp      = 0.0<br>
    ratio_Si2N_phaeo   = 0.0<br>
    ratio_S2N_diatoms  = 0.03  <br>
    ratio_S2N_sp       = 0.03 <br>
    ratio_S2N_phaeo    = 0.03<br>
    ratio_Fe2C_diatoms = 0.0033<br>
    ratio_Fe2C_sp      = 0.0033<br>
    ratio_Fe2C_phaeo   = 0.1<br>
    ratio_Fe2N_diatoms = 0.023 <br>
    ratio_Fe2N_sp      = 0.023<br>
    ratio_Fe2N_phaeo   = 0.7<br>
    ratio_Fe2DON       = 0.023<br>
    ratio_Fe2DOC_s     = 0.1<br>
    ratio_Fe2DOC_l     = 0.033<br>
    fr_resp            = 0.05<br>
    tau_min            = 5200.0<br>
    tau_max            = 173000.0<br>
    algal_vel          = 0.0000000111<br>
    R_dFe2dust         = 0.035<br>
    dustFe_sol         = 0.005<br>
    chlabs_diatoms     = 0.03<br>
    chlabs_sp          = 0.01<br>
    chlabs_phaeo       = 0.05<br>
    alpha2max_low_diatoms = 0.8<br>
    alpha2max_low_sp      = 0.67<br>
    alpha2max_low_phaeo   = 0.67<br>
    beta2max_diatoms   = 0.018<br>
    beta2max_sp        = 0.0025<br>
    beta2max_phaeo     = 0.01<br>
    mu_max_diatoms     = 1.44<br>
    mu_max_sp          = 0.851<br>
    mu_max_phaeo       = 0.851<br>
    grow_Tdep_diatoms  = 0.06<br>
    grow_Tdep_sp       = 0.06<br>
    grow_Tdep_phaeo    = 0.06<br>
    fr_graze_diatoms   = 0.0<br>
    fr_graze_sp        = 0.1<br>
    fr_graze_phaeo     = 0.1<br>
    mort_pre_diatoms   = 0.007<br>
    mort_pre_sp        = 0.007<br>
    mort_pre_phaeo     = 0.007<br>
    mort_Tdep_diatoms  = 0.03<br>
    mort_Tdep_sp       = 0.03<br>
    mort_Tdep_phaeo    = 0.03<br>
    k_exude_diatoms    = 0.0<br>
    k_exude_sp         = 0.0<br>
    k_exude_phaeo      = 0.0<br>
    K_Nit_diatoms      = 1.0<br>
    K_Nit_sp           = 1.0<br>
    K_Nit_phaeo        = 1.0<br>
    K_Am_diatoms       = 0.3<br>
    K_Am_sp            = 0.3<br>
    K_Am_phaeo         = 0.3<br>
    K_Sil_diatoms      = 4.0<br>
    K_Sil_sp           = 0.0<br>
    K_Sil_phaeo        = 0.0<br>
    K_Fe_diatoms       = 1.0<br>
    K_Fe_sp            = 0.2<br>
    K_Fe_phaeo         = 0.1<br>
    f_don_protein      = 0.6<br>
    kn_bac_protein     = 0.03<br>
    f_don_Am_protein   = 0.25<br>
    f_doc_s            = 0.4<br>
    f_doc_l            = 0.4<br>
    f_exude_s          = 1.0<br>
    f_exude_l          = 1.0<br>
    k_bac_s            = 0.03<br>
    k_bac_l            = 0.03<br>
    T_max              = 0.0<br>
    fsal               = 1.0<br>
    op_dep_min         = 0.1<br>
    fr_graze_s         = 0.5<br>
    fr_graze_e         = 0.5<br>
    fr_mort2min        = 0.5<br>
    fr_dFe             = 0.3<br>
    k_nitrif           = 0.0<br>
    t_iron_conv        = 3065.0<br>
    max_loss           = 0.9<br>
    max_dfe_doc1       = 0.2<br>
    fr_resp_s          = 0.75<br>
    y_sk_DMS           = 0.5<br>
    t_sk_conv          = 3.0<br>
    t_sk_ox            = 10.0<br>
    algaltype_diatoms  = 0.0<br>
    algaltype_sp       = 0.5<br>
    algaltype_phaeo    = 0.5<br>
    nitratetype        = -1.0<br>
    ammoniumtype       = 1.0<br>
    silicatetype       = -1.0<br>
    dmspptype          = 0.5<br>
    dmspdtype          = -1.0<br>
    humtype            = 1.0<br>
    doctype_s          = 0.5<br>
    doctype_l          = 0.5<br>
    dontype_protein    = 0.5<br>
    fedtype_1          = 0.5<br>
    feptype_1          = 0.5<br>
    zaerotype_bc1      = 1.0<br>
    zaerotype_bc2      = 1.0<br>
    zaerotype_dust1    = 1.0<br>
    zaerotype_dust2    = 1.0<br>
    zaerotype_dust3    = 1.0<br>
    zaerotype_dust4    = 1.0<br>
    ratio_C2N_diatoms  = 7.0<br>
    ratio_C2N_sp       = 7.0<br>
    ratio_C2N_phaeo    = 7.0<br>
    ratio_chl2N_diatoms= 2.1<br>
    ratio_chl2N_sp     = 1.1<br>
    ratio_chl2N_phaeo  = 0.84<br>
    F_abs_chl_diatoms  = 2.0<br>
    F_abs_chl_sp       = 4.0<br>
    F_abs_chl_phaeo    = 5.0<br>
    ratio_C2N_proteins = 7.0<br>
/

&forcing_nml<br>
    formdrag        = .false.<br>
    atmbndy         = 'default'<br>
    fyear_init      = 1990<br>
    ycycle          = 20<br>
    atm_data_format = 'bin'<br>
    atm_data_type   = 'LYq'<br>
    atm_data_dir    = '/usr/projects/climate/eclare/DATA/Consortium/CICE_data/forcing/gx1/COREII'<br>
    calc_strair     = .true.<br>
    highfreq        = .false.<br>
    natmiter        = 5<br>
    calc_Tsfc       = .true.<br>
    precip_units    = 'mm_per_sec'<br>
    ustar_min       = 0.0005<br>
    fbot_xfer_type  = 'constant'<br>
    update_ocn_f    = .false.<br>
    l_mpond_fresh   = .false.<br>
    tfrz_option     = 'mushy'<br>
    oceanmixed_ice  = .true.<br>
    ocn_data_format = 'nc'<br>
    sss_data_type   = 'ncar'<br>
    sst_data_type   = 'ncar'<br>
    ocn_data_dir    = '/usr/projects/climate/eclare/DATA/gx1v3/gx1v3/forcing/'<br>
    oceanmixed_file = 'oceanmixed_ice_depth.nc'<br>
    restore_sst     = .false.<br>
    trestore        =  90<br>
    restore_ice     = .false.<br>
/

&icefields_nml<br>
    f_tmask        = .true.<br>
    f_blkmask      = .true.<br>
    f_tarea        = .true.<br>
    f_uarea        = .true.<br>
    f_dxt          = .false.<br>
    f_dyt          = .false.<br>
    f_dxu          = .false.<br>
    f_dyu          = .false.<br>
    f_HTN          = .false.<br>
    f_HTE          = .false.<br>
    f_ANGLE        = .true.<br>
    f_ANGLET       = .true.<br>
    f_NCAT         = .true.<br>
    f_VGRDi        = .false.<br>
    f_VGRDs        = .false.<br>
    f_VGRDb        = .false.<br>
    f_VGRDa        = .true.<br>
    f_bounds       = .false.<br>
    f_aice         = 'm' <br>
    f_hi           = 'm'<br>
    f_hs           = 'm' <br>
    f_Tsfc         = 'm' <br>
    f_sice         = 'm' <br>
    f_uvel         = 'm' <br>
    f_vvel         = 'm' <br>
    f_uatm         = 'm' <br>
    f_vatm         = 'm' <br>
    f_fswdn        = 'm' <br>
    f_flwdn        = 'm'<br>
    f_snowfrac     = 'x'<br>
    f_snow         = 'x' <br>
    f_snow_ai      = 'm' <br>
    f_rain         = 'x' <br>
    f_rain_ai      = 'm' <br>
    f_sst          = 'm' <br>
    f_sss          = 'm' <br>
    f_uocn         = 'm' <br>
    f_vocn         = 'm' <br>
    f_frzmlt       = 'm'<br>
    f_fswfac       = 'm'<br>
    f_fswint_ai    = 'm'<br>
    f_fswabs       = 'x' <br>
    f_fswabs_ai    = 'm' <br>
    f_albsni       = 'm' <br>
    f_alvdr        = 'x'<br>
    f_alidr        = 'x'<br>
    f_alvdf        = 'x'<br>
    f_alidf        = 'x'<br>
    f_alvdr_ai     = 'x'<br>
    f_alidr_ai     = 'x'<br>
    f_alvdf_ai     = 'x'<br>
    f_alidf_ai     = 'x'<br>
    f_albice       = 'x'<br>
    f_albsno       = 'x'<br>
    f_albpnd       = 'x'<br>
    f_coszen       = 'x'<br>
    f_flat         = 'x' <br>
    f_flat_ai      = 'm' <br>
    f_fsens        = 'x' <br>
    f_fsens_ai     = 'm' <br>
    f_fswup        = 'x' <br>
    f_flwup        = 'x' <br>
    f_flwup_ai     = 'm' <br>
    f_evap         = 'x' <br>
    f_evap_ai      = 'm' <br>
    f_Tair         = 'm' <br>
    f_Tref         = 'x' <br>
    f_Qref         = 'x'<br>
    f_congel       = 'm' <br>
    f_frazil       = 'm' <br>
    f_snoice       = 'm' <br>
    f_dsnow        = 'x' <br>
    f_melts        = 'm'<br>
    f_meltt        = 'm'<br>
    f_meltb        = 'm'<br>
    f_meltl        = 'm'<br>
    f_fresh        = 'x'<br>
    f_fresh_ai     = 'm'<br>
    f_fsalt        = 'x'<br>
    f_fsalt_ai     = 'm'<br>
    f_fbot         = 'm'<br>
    f_fhocn        = 'x' <br>
    f_fhocn_ai     = 'm' <br>
    f_fswthru      = 'x' <br>
    f_fswthru_ai   = 'm' <br>
    f_fsurf_ai     = 'x'<br>
    f_fcondtop_ai  = 'x'<br>
    f_fmeltt_ai    = 'x' <br>
    f_strairx      = 'm' <br>
    f_strairy      = 'm' <br>
    f_strtltx      = 'x' <br>
    f_strtlty      = 'x' <br>
    f_strcorx      = 'x' <br>
    f_strcory      = 'x' <br>
    f_strocnx      = 'x' <br>
    f_strocny      = 'x' <br>
    f_strintx      = 'x' <br>
    f_strinty      = 'x'<br>
    f_taubx        = 'x'<br>
    f_tauby        = 'x'<br>
    f_strength     = 'm'<br>
    f_divu         = 'm'<br>
    f_shear        = 'm'<br>
    f_sig1         = 'm' <br>
    f_sig2         = 'm' <br>
    f_dvidtt       = 'm' <br>
    f_dvidtd       = 'm' <br>
    f_daidtt       = 'm'<br>
    f_daidtd       = 'm' <br>
    f_dagedtt      = 'm'<br>
    f_dagedtd      = 'm' <br>
    f_mlt_onset    = 'm'<br>
    f_frz_onset    = 'm'<br>
    f_hisnap       = 'd'<br>
    f_aisnap       = 'd'<br>
    f_trsig        = 'm'<br>
    f_icepresent   = 'm'<br>
    f_iage         = 'm'<br>
    f_FY           = 'x'<br>
    f_aicen        = 'x'<br>
    f_vicen        = 'x'<br>
    f_vsnon        = 'x'<br>
    f_snowfracn    = 'x'<br>
    f_keffn_top    = 'x'<br>
    f_Tinz         = 'x'<br>
    f_Sinz         = 'x'<br>
    f_Tsnz         = 'x'<br>
    f_fsurfn_ai    = 'x'<br>
    f_fcondtopn_ai = 'x'<br>
    f_fmelttn_ai   = 'x'<br>
    f_flatn_ai     = 'x'<br>
    f_fsensn_ai     = 'x'<br>
/

&icefields_mechred_nml<br>
    f_alvl         = 'm'<br>
    f_vlvl         = 'm'<br>
    f_ardg         = 'm'<br>
    f_vrdg         = 'm'<br>
    f_dardg1dt     = 'x'<br>
    f_dardg2dt     = 'x'<br>
    f_dvirdgdt     = 'x'<br>
    f_opening      = 'x'<br>
    f_ardgn        = 'x'<br>
    f_vrdgn        = 'x'<br>
    f_dardg1ndt    = 'x'<br>
    f_dardg2ndt    = 'x'<br>
    f_dvirdgndt    = 'x'<br>
    f_krdgn        = 'x'<br>
    f_aparticn     = 'x'<br>
    f_aredistn     = 'x'<br>
    f_vredistn     = 'x'<br>
    f_araftn       = 'x'<br>
    f_vraftn       = 'x'<br>
/

&icefields_pond_nml<br>
    f_apondn       = 'x'<br>
    f_apeffn       = 'x'<br>
    f_hpondn       = 'x'<br>
    f_apond        = 'm'<br>
    f_hpond        = 'm'<br>
    f_ipond        = 'm'<br>
    f_apeff        = 'm'<br>
    f_apond_ai     = 'm'<br>
    f_hpond_ai     = 'm'<br>
    f_ipond_ai     = 'm'<br>
    f_apeff_ai     = 'm'<br>
/

&icefields_bgc_nml<br>
    f_faero_atm    = 'x'<br>
    f_faero_ocn    = 'x'<br>
    f_aero         = 'x'  <br>
    f_fbio         = 'm'<br>
    f_fbio_ai      = 'm'<br>
    f_zaero        = 'x'<br>
    f_bgc_S        = 'm'<br>
    f_bgc_N        = 'm'<br>
    f_bgc_C        = 'x'<br>
    f_bgc_DOC      = 'm'<br>
    f_bgc_DIC      = 'x'<br>
    f_bgc_chl      = 'x'<br>
    f_bgc_Nit      = 'm'<br>
    f_bgc_Am       = 'm'<br>
    f_bgc_Sil      = 'm'<br>
    f_bgc_DMSPp    = 'x'<br>
    f_bgc_DMSPd    = 'x'<br>
    f_bgc_DMS      = 'x'<br>
    f_bgc_DON      = 'x'  <br>
    f_bgc_Fe       = 'm'  <br>
    f_bgc_hum      = 'm'   <br>
    f_bgc_PON      = 'm'<br>
    f_bgc_ml       = 'm'<br>
    f_upNO         = 'm'<br>
    f_upNH         = 'm'<br>
    f_bTin         = 'm'<br>
    f_bphi         = 'm' <br>
    f_iDi          = 'm'<br>
    f_iki          = 'm'<br>
    f_fbri         = 'm'  <br>
    f_hbri         = 'm'<br>
    f_zfswin       = 'm'<br>
    f_bionet       = 'm'<br>
    f_biosnow      = 'm'<br>
    f_grownet      = 'm'<br>
    f_PPnet        = 'm'<br>
    f_algalpeak    = 'm'<br>
    f_zbgc_frac    = 'm'<br>
/

&icefields_drag_nml<br>
    f_drag         = 'x'<br>
    f_Cdn_atm      = 'x'<br>
    f_Cdn_ocn      = 'x'<br>
/