RASM simulations: Quality Control for Community Based Sea Ice Model Development

RASM simulations: Quality Control for Community Based Sea Ice Model Development
-----------------------------------------------------------------------------------

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