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

description: 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:

&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>
&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>
&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>
&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>
&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>
&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>
&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>
&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>
&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>
&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>
&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>
&icefields_drag_nml<br>
 f_cdn_atm = 'hxxxx'<br>
 f_cdn_ocn = 'hxxxx'<br>
 f_drag = 'hxxxx'<br>
/<br>
&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>
&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>
&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>
&ice_prescribed_nml<br>
 prescribed_ice = .false.<br>
/