&setup_nml
    days_per_year     = 365
    use_leap_years    = .true.
    year_init         = 1993
    month_init        = 1
    day_init          = 1
    sec_init          = 0
    istep0            = 0
    dt                = 1800
    ndtd              = 1
    npt_unit = 'm'
    npt = 12
    runtype = 'initial'
    ice_ic            = 'none'
    restart_ext       = .false.
    use_restart_time = .false.
    restart_format    = 'hdf5'
    lcdf64            = .true.
    numin             = 21
    numax             = 89
    restart_dir       = './restart/'
    restart_file      = 'iced'
    pointer_file      = './ice.restart_file'
    dumpfreq          = 'y'
    dumpfreq_n        = 3
    dumpfreq_base     = 'init'
    dump_last         = .true.
    bfbflag           = 'off'
    diagfreq          = 96
    diag_type         = 'file'
    diag_file         = 'ice_diag.d'
    debug_model       = .false.
    debug_model_step  = 0
    debug_model_i     = 100
    debug_model_j     = 100
    debug_model_iblk  = 10
    debug_model_task  = 1
    debug_forcing     = .false.
    print_global      = .false.
    print_points      = .false.
    timer_stats       = .false.
    memory_stats      = .false.
    conserv_check     = .false.
    latpnt(1)         =  -65.
    lonpnt(1)         =  0.
    latpnt(2)         = -65.
    lonpnt(2)         = -179.
    histfreq          = 'd','x','x','x','x'
    histfreq_n        =  1 , 1 , 1, 1 , 1
    histfreq_base     = 'zero'
    hist_avg          = .true.
    history_dir       = './history/'
    history_file      = 'iceh'
    history_precision = 4
    history_format    = 'default'
    write_ic          = .false.
    incond_dir        = './history/'
    incond_file       = 'iceh_ic'
    version_name      = 'CICE_6.4.0'
/

&grid_nml
    grid_format       = 'nc'
    grid_type         = 'tripole'
    grid_ice          = 'B'
    grid_atm          = 'A'
    grid_ocn          = 'B'
    grid_file         = '/g/data/gv90/da1339/grids/ACCESS-OM2-025_grid.nc'
    kmt_type          = 'file'
    kmt_file = '/g/data/gv90/da1339/grounded_icebergs/MyThin_db/kmt_mod_thinned-0p75_v1p50_iter0.nc'
    bathymetry_file   = '/g/data/gv90/da1339/grids/ACCESS-OM2-025_topog.nc'
    bathymetry_format = 'default'
    use_bathymetry    = .false.
    gridcpl_file      = 'unknown_gridcpl_file'
    kcatbound         = 0
    dxrect            = 30.e5
    dyrect            = 30.e5
    lonrefrect        = -156.50d0
    latrefrect        =   71.35d0
    scale_dxdy        = .false.
    dxscale           = 1.d0
    dyscale           = 1.d0
    close_boundaries  = .false.
    ncat              = 5
    nfsd              = 1
    nilyr             = 4
    nslyr             = 1
    nblyr             = 1
    orca_halogrid     = .false.
/

&tracer_nml
    n_iso             = 0
    n_aero            = 0
    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             = .false.
    restart_FY        = .false.
    tr_lvl            = .false.
    restart_lvl       = .false.
    tr_pond_topo      = .false.
    restart_pond_topo = .false.
    tr_pond_lvl       = .false.
    restart_pond_lvl  = .false.
    tr_snow           = .false.
    restart_snow      = .false.
    tr_iso            = .false.
    restart_iso       = .false.
    tr_aero           = .false.
    restart_aero      = .false.
    tr_fsd            = .false.
    restart_fsd       = .false.
/

&thermo_nml
    kitd              = 1
    ktherm            = 2
    conduct           = 'bubbly'
    ksno              = 0.3d0
    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
    hfrazilmin        = 0.05d0
    floediam          = 300.0d0
/

&dynamics_nml
    kdyn = 1
    ndte = 240
    revised_evp = .false.
    evp_algorithm        = 'standard_2d'
    brlx                 = 300.0
    arlx                 = 300.0
    advection            = 'remap'
    kstrength            = 0
    krdg_partic          = 1
    krdg_redist          = 1
    mu_rdg               = 3
    Pstar = 27500
    Cstar = 20
    Cf                   = 17.
    Ktens = 0.2
    e_yieldcurve = 1.2
    e_plasticpot = 1.2
    visc_method          = 'avg_zeta'
    elasticDamp          = 0.36d0
    deltaminEVP          = 1e-11
    deltaminVP           = 2e-9
    capping_method       = 'max'
    seabed_stress        = .false.
    seabed_stress_method = 'LKD'
    k1                   = 7.5
    k2                   = 15.
    alphab               = 20.  
    threshold_hw         = 30.
    coriolis             = 'latitude'
    kridge               = 1
    ktransport           = 1
    ssh_stress           = 'geostrophic'
    maxits_nonlin        = 10
    precond              = 'pgmres'
    dim_fgmres           = 50
    dim_pgmres           = 5
    maxits_fgmres        = 1
    maxits_pgmres        = 1
    monitor_nonlin       = .false.
    monitor_fgmres       = .false.
    monitor_pgmres       = .false.
    ortho_type           = 'mgs'
    reltol_nonlin        = 1e-8
    reltol_fgmres        = 1e-2
    reltol_pgmres        = 1e-6
    algo_nonlin          = 'picard'
    use_mean_vrel        = .true.
/

&shortwave_nml
    shortwave       = 'dEdd'
    albedo_type     = 'ccsm3'
    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.6
    sw_redist       = .false.
    sw_frac         = 0.9d0
    sw_dtemp        = 0.02d0
/

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

&snow_nml
    snwredist       = 'none'
    snwgrain        = .false.
    use_smliq_pnd   = .false.
    rsnw_fall       =  100.0
    rsnw_tmax       = 1500.0
    rhosnew         =  100.0
    rhosmin         =  100.0
    rhosmax         =  450.0
    windmin         =   10.0
    drhosdwind      =   27.3
    snwlvlfac       =    0.3
    snw_aging_table = 'test'
    snw_filename    = 'unknown'
    snw_rhos_fname  = 'unknown'
    snw_Tgrd_fname  = 'unknown'
    snw_T_fname     = 'unknown'
    snw_tau_fname   = 'unknown'
    snw_kappa_fname = 'unknown'
    snw_drdt0_fname = 'unknown'
/

&forcing_nml
    formdrag        = .false.
    atmbndy         = 'default'
    rotate_wind     = .true.
    calc_strair     = .true.
    calc_Tsfc       = .true.
    highfreq        = .false.
    natmiter        = 5
    atmiter_conv    = 0.0d0
    ustar_min       = 0.0005
    iceruf          = 0.0005
    calc_dragio     = .true.
    iceruf_ocn      = 0.03
    emissivity      = 0.985
    fbot_xfer_type  = 'constant'
    update_ocn_f    = .false.
    l_mpond_fresh   = .false.
    tfrz_option     = 'mushy'
    oceanmixed_ice  = .true.
    wave_spec_type  = 'none'
    wave_spec_file  = 'unknown_wave_spec_file'
    nfreq           = 25
    restart_coszen  = .false.
    restore_ice     = .false.
    restore_ocn     = .true.
    trestore        =  1
    precip_units    = 'mks'
    default_season  = 'winter'
    atm_data_type   = 'ERA5'
    ocn_data_type   = 'AFIM'
    ocn_data_freq   = 'daily'
    bgc_data_type   = 'default'
    fe_data_type    = 'default'
    ice_data_type   = 'default'
    ice_data_conc   = 'default'
    ice_data_dist   = 'uniform'
    fyear_init = 1993
    ycycle          = 1
    atm_data_format = 'nc'
    atm_data_dir    = '/g/data/gv90/da1339/afim_input/ERA5/0p25/'
    bgc_data_dir    = ''
    ocn_data_format = 'nc'
    ocn_data_dir    = '/g/data/gv90/da1339/afim_input/ORAS/daily/sfc/fin/'
    oceanmixed_file = 'ORAS_.nc'
/

&domain_nml
    nprocs = 1232
    nx_global         = 1440
    ny_global         = 1080
    block_size_x = 20
    block_size_y = 20
      max_blocks = 10
    processor_shape = 'square-ice'
    distribution_type = 'rake'
    distribution_wght = 'latitude'
    distribution_wght_file = 'unknown'
    ew_boundary_type  = 'cyclic'
    ns_boundary_type  = 'tripoleT'
    maskhalo_dyn      = .false.
    maskhalo_remap    = .false.
    maskhalo_bound    = .false.
    add_mpi_barriers  = .false.
    debug_blocks      = .false.
/

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

&icefields_mechred_nml
    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       = 'x'
    f_apeffn       = 'x'
    f_hpondn       = 'x'
    f_apond        = 'x'
    f_hpond        = 'x'
    f_ipond        = 'x'
    f_apeff        = 'x'
    f_apond_ai     = 'x'
    f_hpond_ai     = 'x'
    f_ipond_ai     = 'x'
    f_apeff_ai     = 'x'
/

&icefields_snow_nml
    f_smassicen    = 'x'
    f_smassliqn    = 'x'
    f_rhos_cmpn    = 'x'
    f_rhos_cntn    = 'x'
    f_rsnwn        = 'x'
    f_smassice     = 'x'
    f_smassliq     = 'x'
    f_rhos_cmp     = 'x'
    f_rhos_cnt     = 'x'
    f_rsnw         = 'x'
    f_meltsliq     = 'x'
    f_fsloss       = 'x'
/

&icefields_bgc_nml
    f_fiso_atm     = 'x'
    f_fiso_ocn     = 'x'
    f_iso          = '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         = 'md'
    f_Cdn_atm      = 'md'
    f_Cdn_ocn      = 'md'
/

&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'
/
