Rasm Simulations: Quality Control For Community Based Sea Ice Model Development

RASM simulations: Quality Control for Community Based Sea Ice Model Development<br> ----------------------------------------------------------------------------------- Data from the Regional Arctic System Model (RASM) simulations by Andrew F. Roberts at Naval Postgraduate School, 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 RASM with CICE sea ice model version 5.2 with the namelist: &amp;setup_nml<br> bfbflag = .false. <br> diagfreq = 24 <br> hist_avg = .true. <br> histfreq = 'h','x','x','x','x'<br> histfreq_n = 1,0,0,0,0<br> ice_ic = '/work4/projects/apcraig_dir/RACM/inputdata/ccsm4_init/R1008Gcaaa01a_loop2/1979-09-01/R1008Gcaaa01a_loop2.cice.r.1979-09-01-00000.nc'<br> latpnt = 90.0,-65.0 <br> lcdf64 = .true. <br> lonpnt = 0.0,-45.0 <br> ndtd = 1 <br> pointer_file = './rpointer.ice'<br> print_global = .true. <br> print_points = .false. <br> restart_ext = .false. <br> restart_format = 'pio'<br> year_init = 1 <br> /<br> &amp;grid_nml<br> grid_file = '/work4/projects/apcraig_dir/RACM/inputdata/ice/cice/ar9v3_20101118.grid'<br> grid_format = 'bin'<br> grid_type = 'regional'<br> gridcpl_file = '/work4/projects/apcraig_dir/RACM/inputdata/ice/cice/domain.ocn.ar9v4.100920.nc'<br> kcatbound = 0 <br> kmt_file = '/work4/projects/apcraig_dir/RACM/inputdata/ice/cice/ar9v3_20110106.kmt'<br> /<br> &amp;tracer_nml<br> restart_aero = .false. <br> restart_age = .false. <br> restart_fy = .false. <br> restart_lvl = .false. <br> restart_pond_cesm = .false. <br> restart_pond_lvl = .false. <br> restart_pond_topo = .false. <br> tr_aero = .false.<br> tr_fy = .true. <br> tr_iage = .true. <br> tr_lvl = .true. <br> tr_pond_cesm = .false. <br> tr_pond_lvl = .true. <br> tr_pond_topo = .false. <br> /<br> &amp;thermo_nml<br> a_rapid_mode = 0.5e-03 <br> aspect_rapid_mode = 1.0 <br> conduct = 'bubbly'<br> dsdt_slow_mode = -1.5e-07 <br> kitd = 1 <br> ktherm = 2 <br> phi_c_slow_mode = 0.05 <br> phi_i_mushy = 0.85 <br> rac_rapid_mode = 10 <br> /<br> &amp;dynamics_nml<br> advection = 'remap'<br> cf = 21.3<br> kdyn = 2 (EAP) or kdyn = 2 (EVP)<br> krdg_partic = 1<br> krdg_redist = 1<br> kstrength = 1 <br> mu_rdg = 3.0<br> ndte = 600 <br> revised_evp = .false. <br> /<br> &amp;shortwave_nml<br> ahmax = 0.3 <br> albedo_type = 'default'<br> albicei = 0.45<br> albicev = 0.75<br> albsnowi = 0.73<br> albsnowv = 0.98<br> dt_mlt = 1.00 <br> kalg = 0.0 <br> r_ice = 0.0 <br> r_pnd = 0.0 <br> r_snw = 0.50 <br> rsnw_mlt = 1000. <br> shortwave = 'dEdd'<br> /<br> &amp;ponds_nml<br> dpscale = 1.0e-3 <br> frzpnd = 'hlid'<br> hp1 = 0.01 <br> hs0 = 0.03 <br> hs1 = 0.03 <br> pndaspect = 0.8 <br> rfracmax = 0.85 <br> rfracmin = 0.15 <br> /<br> &amp;forcing_nml<br> fbot_xfer_type = 'constant'<br> formdrag = .false.<br> highfreq = .true.<br> l_mpond_fresh = .false.<br> natmiter = 25 <br> tfrz_option = 'mushy'<br> /<br> &amp;domain_nml<br> distribution_type = 'roundrobin'<br> ew_boundary_type = 'cyclic'<br> maskhalo_bound = .true.<br> maskhalo_dyn = .true.<br> maskhalo_remap = .true.<br> ns_boundary_type = 'open'<br> processor_shape = 'null'<br> /<br> &amp;zbgc_nml<br> bgc_data_dir = 'unknown_bgc_data_dir'<br> bgc_flux_type = 'Jin2006'<br> nit_data_type = 'unknown'<br> phi_snow = 0.5 <br> restart_bgc = .false. <br> restart_hbrine = .false. <br> restore_bgc = .false. <br> sil_data_type = 'unknown'<br> skl_bgc = .false. <br> tr_bgc_am_sk = .false. <br> tr_bgc_c_sk = .false. <br> tr_bgc_chl_sk = .false. <br> tr_bgc_dms_sk = .false. <br> tr_bgc_dmspd_sk = .false. <br> tr_bgc_dmspp_sk = .false. <br> tr_bgc_sil_sk = .false. <br> tr_brine = .false. <br> /<br> &amp;icefields_bgc_nml<br> f_aero = 'mxxxx'<br> f_aeron = 'xxxxx'<br> f_bgc_am_ml = 'xxxxx'<br> f_bgc_am_sk = 'xxxxx'<br> f_bgc_c = 'xxxxx'<br> f_bgc_c_sk = 'xxxxx'<br> f_bgc_chl = 'xxxxx'<br> f_bgc_chl_sk = 'xxxxx'<br> f_bgc_dms = 'xxxxx'<br> f_bgc_dms_ml = 'xxxxx'<br> f_bgc_dms_sk = 'xxxxx'<br> f_bgc_dmsp_ml = 'xxxxx'<br> f_bgc_dmspd = 'xxxxx'<br> f_bgc_dmspd_sk = 'xxxxx'<br> f_bgc_dmspp = 'xxxxx'<br> f_bgc_dmspp_sk = 'xxxxx'<br> f_bgc_n = 'xxxxx'<br> f_bgc_n_sk = 'xxxxx'<br> f_bgc_nh = 'xxxxx'<br> f_bgc_nit_ml = 'xxxxx'<br> f_bgc_nit_sk = 'xxxxx'<br> f_bgc_no = 'xxxxx'<br> f_bgc_s = 'xxxxx'<br> f_bgc_sil = 'xxxxx'<br> f_bgc_sil_ml = 'xxxxx'<br> f_bgc_sil_sk = 'xxxxx'<br> f_bphi = 'xxxxx'<br> f_btin = 'xxxxx'<br> f_faero_atm = 'mxxxx'<br> f_faero_ocn = 'mxxxx'<br> f_fbri = 'xxxxx'<br> f_fn = 'xxxxx'<br> f_fn_ai = 'xxxxx'<br> f_fnh = 'xxxxx'<br> f_fnh_ai = 'xxxxx'<br> f_fno = 'xxxxx'<br> f_fno_ai = 'xxxxx'<br> f_fsil = 'xxxxx'<br> f_fsil_ai = 'xxxxx'<br> f_grownet = 'xxxxx'<br> f_hbri = 'xxxxx'<br> f_ppnet = 'xxxxx'<br> /<br> &amp;icefields_drag_nml<br> f_cdn_atm = 'hxxxx'<br> f_cdn_ocn = 'hxxxx'<br> f_drag = 'hxxxx'<br> /<br> &amp;icefields_mechred_nml<br> f_alvl = 'hxxxx'<br> f_aparticn = 'xxxxx'<br> f_araftn = 'xxxxx'<br> f_ardg = 'hxxxx'<br> f_ardgn = 'hxxxx'<br> f_aredistn = 'xxxxx'<br> f_dardg1dt = 'xxxxx'<br> f_dardg1ndt = 'xxxxx'<br> f_dardg2dt = 'xxxxx'<br> f_dardg2ndt = 'xxxxx'<br> f_dvirdgdt = 'xxxxx'<br> f_dvirdgndt = 'xxxxx'<br> f_krdgn = 'xxxxx'<br> f_opening = 'hxxxx'<br> f_vlvl = 'hxxxx'<br> f_vraftn = 'xxxxx'<br> f_vrdg = 'hxxxx'<br> f_vrdgn = 'hxxxx'<br> f_vredistn = 'xxxxx'<br> /<br> &amp;icefields_pond_nml<br> f_apeff = 'hxxxx'<br> f_apeff_ai = 'hxxxx'<br> f_apeffn = 'hxxxx'<br> f_apond = 'hxxxx'<br> f_apond_ai = 'hxxxx'<br> f_apondn = 'hxxxx'<br> f_hpond = 'hxxxx'<br> f_hpond_ai = 'hxxxx'<br> f_hpondn = 'hxxxx'<br> f_ipond = 'hxxxx'<br> f_ipond_ai = 'hxxxx'<br> /<br> &amp;icefields_nml<br> f_a11 = 'hxxxx'<br> f_a12 = 'hxxxx'<br> f_aice = 'hxxxx'<br> f_aicen = 'hxxxx'<br> f_aisnap = 'hxxxx'<br> f_albice = 'hxxxx'<br> f_albpnd = 'hxxxx'<br> f_albsni = 'hxxxx'<br> f_albsno = 'hxxxx'<br> f_alidf = 'hxxxx'<br> f_alidr = 'hxxxx'<br> f_alvdf = 'hxxxx'<br> f_alvdr = 'hxxxx'<br> f_angle = .true.<br> f_anglet = .true.<br> f_blkmask = .true.<br> f_bounds = .false.<br> f_congel = 'hxxxx'<br> f_coszen = 'hxxxx'<br> f_daidtd = 'hxxxx'<br> f_daidtt = 'hxxxx'<br> f_divu = 'hxxxx'<br> f_dsnow = 'hxxxx'<br> f_dvidtd = 'hxxxx'<br> f_dvidtt = 'hxxxx'<br> f_dxt = .false.<br> f_dxu = .false.<br> f_dyt = .false.<br> f_dyu = .false.<br> f_e11 = 'xxxxx'<br> f_e12 = 'xxxxx'<br> f_e22 = 'xxxxx'<br> f_evap = 'hxxxx'<br> f_evap_ai = 'hxxxx'<br> f_fcondtop_ai = 'hxxxx'<br> f_fcondtopn_ai = 'hxxxx'<br> f_fhocn = 'hxxxx'<br> f_fhocn_ai = 'hxxxx'<br> f_flat = 'hxxxx'<br> f_flat_ai = 'hxxxx'<br> f_flatn_ai = 'hxxxx'<br> f_flwdn = 'hxxxx'<br> f_flwup = 'hxxxx'<br> f_flwup_ai = 'hxxxx'<br> f_fmeltt_ai = 'hxxxx'<br> f_fmelttn_ai = 'hxxxx'<br> f_frazil = 'hxxxx'<br> f_fresh = 'hxxxx'<br> f_fresh_ai = 'hxxxx'<br> f_frz_onset = 'hxxxx'<br> f_frzmlt = 'hxxxx'<br> f_fsalt = 'hxxxx'<br> f_fsalt_ai = 'hxxxx'<br> f_fsens = 'hxxxx'<br> f_fsens_ai = 'hxxxx'<br> f_fsensn_ai = 'hxxxx'<br> f_fsurf_ai = 'hxxxx'<br> f_fsurfn_ai = 'hxxxx'<br> f_fswabs = 'hxxxx'<br> f_fswabs_ai = 'hxxxx'<br> f_fswdn = 'hxxxx'<br> f_fswfac = 'hxxxx'<br> f_fswint_ai = 'xxxxx'<br> f_fswthru = 'hxxxx'<br> f_fswthru_ai = 'hxxxx'<br> f_fy = 'hxxxx'<br> f_hi = 'hxxxx'<br> f_hisnap = 'xxxxx'<br> f_hs = 'hxxxx'<br> f_hte = .false.<br> f_htn = .false.<br> f_iage = 'hxxxx'<br> f_icepresent = 'hxxxx'<br> f_keffn_top = 'xxxxx'<br> f_meltb = 'hxxxx'<br> f_meltl = 'hxxxx'<br> f_melts = 'hxxxx'<br> f_meltt = 'hxxxx'<br> f_mlt_onset = 'hxxxx'<br> f_ncat = .true.<br> f_qref = 'hxxxx'<br> f_rain = 'hxxxx'<br> f_rain_ai = 'hxxxx'<br> f_s11 = 'hxxxx'<br> f_s12 = 'hxxxx'<br> f_s22 = 'hxxxx'<br> f_shear = 'hxxxx'<br> f_sice = 'hxxxx'<br> f_sig1 = 'hxxxx'<br> f_sig2 = 'hxxxx'<br> f_sinz = 'xxxxx'<br> f_snoice = 'hxxxx'<br> f_snow = 'hxxxx'<br> f_snow_ai = 'hxxxx'<br> f_sss = 'hxxxx'<br> f_sst = 'hxxxx'<br> f_strairx = 'hxxxx'<br> f_strairy = 'hxxxx'<br> f_strcorx = 'hxxxx'<br> f_strcory = 'hxxxx'<br> f_strength = 'hxxxx'<br> f_strintx = 'hxxxx'<br> f_strinty = 'hxxxx'<br> f_strocnx = 'hxxxx'<br> f_strocny = 'hxxxx'<br> f_strtltx = 'hxxxx'<br> f_strtlty = 'hxxxx'<br> f_tair = 'hxxxx'<br> f_tarea = .true.<br> f_tinz = 'xxxxx'<br> f_tmask = .true.<br> f_tref = 'hxxxx'<br> f_trsig = 'hxxxx'<br> f_tsfc = 'hxxxx'<br> f_tsnz = 'xxxxx'<br> f_uarea = .true.<br> f_uatm = 'hxxxx'<br> f_uocn = 'hxxxx'<br> f_uvel = 'hxxxx'<br> f_vatm = 'hxxxx'<br> f_vgrdb = .true.<br> f_vgrdi = .true.<br> f_vgrds = .true.<br> f_vicen = 'hxxxx'<br> f_vocn = 'hxxxx'<br> f_vsnon = 'hxxxx'<br> f_vvel = 'hxxxx'<br> f_yieldstress11 = 'hxxxx'<br> f_yieldstress12 = 'hxxxx'<br> f_yieldstress22 = 'hxxxx'<br> /<br> &amp;ice_prescribed_nml<br> prescribed_ice = .false.<br> /