Skip to content

Commit

Permalink
Added Amery ice shelf example.
Browse files Browse the repository at this point in the history
The datasets are not yet in place so this example
can not be run yet.
  • Loading branch information
MJHarrison-GFDL committed Sep 25, 2014
1 parent 1aa03f4 commit d1d22af
Show file tree
Hide file tree
Showing 7 changed files with 1,174 additions and 0 deletions.
1 change: 1 addition & 0 deletions examples/ocean_SIS/Amery/INPUT/.datasets
895 changes: 895 additions & 0 deletions examples/ocean_SIS/Amery/MOM_input

Large diffs are not rendered by default.

Empty file.
60 changes: 60 additions & 0 deletions examples/ocean_SIS/Amery/data_table
Original file line number Diff line number Diff line change
@@ -0,0 +1,60 @@
"ATM" , "p_surf" , "" , "" , .false. , 1.0e5
"ATM" , "p_bot" , "" , "" , .false. , 1.0e5
"ATM" , "t_bot" , "" , "" , .false. , 273.0
"ATM" , "sphum_bot" , "" , "" , .false. , 0.01
"ATM" , "u_bot" , "" , "" , .false. , 0.0
"ATM" , "u_flux" , "" , "" , .false. , 0.0
"ATM" , "v_flux" , "" , "" , .false. , 0.0
"ATM" , "v_bot" , "" , "" , .false. , 0.0
"ATM" , "z_bot" , "" , "" , .false. , 10.0
"ATM" , "gust" , "" , "" , .false. , 1.0e-4
"ATM" , "flux_sw" , "" , "" , .false. , 0.0
"ATM" , "flux_lw" , "" , "" , .false. , 0.0
"ATM" , "lprec" , "" , "" , .false. , 0.0
"ATM" , "fprec" , "" , "" , .false. , 0.0
"LND" , "t_surf" , "" , "" , .false. , 273.0
"LND" , "t_ca" , "" , "" , .false. , 273.0
"LND" , "q_ca" , "" , "" , .false. , 0.01
"LND" , "rough_mom" , "" , "" , .false. , 1.0e-3
"LND" , "rough_heat" , "" , "" , .false. , 1.0e-3
"LND" , "albedo" , "" , "" , .false. , 0.0
"LND" , "sphum_surf" , "" , "" , .false. , 0.01
"LND" , "sphum_ca" , "" , "" , .false. , 0.01
"LND" , "t_flux" , "" , "" , .false. , 0.0
"LND" , "sphum_flux" , "" , "" , .false. , 0.0
"LND" , "lw_flux" , "" , "" , .false. , 0.0
"LND" , "sw_flux" , "" , "" , .false. , 0.0
"LND" , "lprec" , "" , "" , .false. , 0.0
"LND" , "fprec" , "" , "" , .false. , 0.0
"LND" , "dhdt" , "" , "" , .false. , 1.0e-5
"LND" , "dedt" , "" , "" , .false. , 1.0e-5
"LND" , "dedq" , "" , "" , .false. , 1.0e-5
"LND" , "drdt" , "" , "" , .false. , 1.0e-5
"LND" , "drag_q" , "" , "" , .false. , 1.0e-3
"LND" , "p_surf" , "" , "" , .false. , 1.0e5
"ICE" , "lw_flux_dn" , "" , "" , .false. , 0.0
"ICE" , "sw_flux_vis_dir_dn" , "" , "" , .false. , 0.0
"ICE" , "sw_flux_vis_dif_dn" , "" , "" , .false. , 0.0
"ICE" , "sw_flux_nir_dir_dn" , "" , "" , .false. , 0.0
"ICE" , "sw_flux_nir_dif_dn" , "" , "" , .false. , 0.0
"ICE" , "runoff" , "" , "" , .false. , 0.0
"ICE" , "calving" , "" , "" , .false. , 0.0
"ICE" , "u_flux" , "" , "" , .false. , 0.0
"ICE" , "v_flux" , "" , "" , .false. , 0.0
"ICE" , "t_flux" , "" , "" , .false. , 0.0
"ICE" , "q_flux" , "" , "" , .false. , 0.0
"ICE" , "lw_flux_dn" , "" , "" , .false. , 0.0
"ICE" , "lw_flux" , "" , "" , .false. , 0.0
"ICE" , "lprec" , "" , "" , .false. , 0.0
"ICE" , "fprec" , "" , "" , .false. , 0.0
"ICE" , "dhdt" , "" , "" , .false. , 80.0
"ICE" , "dedt" , "" , "" , .false. , 2.0e-6
"ICE" , "drdt" , "" , "" , .false. , 10.0
"ICE" , "sic_obs" , "" , "" , .true. , 0.00
"ICE" , "sit_obs" , "" , "" , .true. , 0.00
"ICE" , "sst_obs" , "SST" , "./INPUT/sst_ice_clim.nc" , .false. , 1.0
"OCN" , "u_flux" , "" , "" , .true. , 0.00
"OCN" , "v_flux" , "" , "" , .true. , 0.00
"OCN" , "t_flux" , "" , "" , .true. , 0.00
"OCN" , "q_flux" , "" , "" , .true. , 0.00
"OCN" , "salt_flux" , "" , "" , .true. , 0.00
78 changes: 78 additions & 0 deletions examples/ocean_SIS/Amery/diag_table
Original file line number Diff line number Diff line change
@@ -0,0 +1,78 @@
Amery_baseline
1992 1 1 0 0 0
## output files
"ocean_month", 1, "months", 1, "days", "time"
"ocean_month_snap", 1, "months", 1, "days", "time"
"ocean_daily", 1, "days", 1, "days", "time"
## output variables
#
#==================
# OCEAN DIAGNOSTICS
#==================
# Monthly snapshots of tracers, interface positions and Kd
"ocean_model","e","e","ocean_month_snap","all",.false.,"none",2
"ocean_model","temp","temp","ocean_month_snap","all",.false.,"none",2
"ocean_model","salt","salt","ocean_month_snap","all",.false.,"none",2
"ocean_model","u","u","ocean_month_snap","all",.false.,"none",2
"ocean_model","v","v","ocean_month_snap","all",.false.,"none",2
"ocean_model","salt_int","salt_int","ocean_month_snap","all",.false.,"none",2
"ocean_model","temp_int","temp_int","ocean_month_snap","all",.false.,"none",2
#"ocean_model","Kd_shear","Kd_shear","ocean_month_snap","all",.false.,"none",2
"ocean_model","Kd_layer","Kd_layer","ocean_month_snap","all",.false.,"none",2
#"ocean_model","Kd_work","Kd_work","ocean_month_snap","all",.false.,"none",2
# Monthly Avg transports in sigma2
"ocean_model","uh","uh","ocean_month","all",.true.,"none",2
"ocean_model","vh","vh","ocean_month","all",.true.,"none",2
"ocean_model","wd","wd","ocean_month","all",.true.,"none",2
#"ocean_model","Res_fn","Res_fn","ocean_month","all",.true.,"none",2
#"ocean_model","col_mass","col_mass","ocean_month","all",.true.,"none",2
#"ocean_model","Tflx_dia_diff","Tflx_dia_diff","ocean_month","all",.true.,"none",2
#"ocean_model","Tflx_dia_adv","Tflx_dia_adv","ocean_month","all",.true.,"none",2
#"ocean_model","Sflx_dia_diff","Sflx_dia_diff","ocean_month","all",.true.,"none",2
#"ocean_model","Sflx_dia_adv","Sflx_dia_adv","ocean_month","all",.true.,"none",2
#"ocean_model","T_adx_2d","T_adx_2d","ocean_month","all",.true.,"none",2
#"ocean_model","T_ady_2d","T_ady_2d","ocean_month","all",.true.,"none",2
#"ocean_model","S_adx_2d","S_adx_2d","ocean_month","all",.true.,"none",2
#"ocean_model","S_ady_2d","S_ady_2d","ocean_month","all",.true.,"none",2
#"ocean_model","T_diffx_2d","T_diffx_2d","ocean_month","all",.true.,"none",2
#"ocean_model","T_diffy_2d","T_diffy_2d","ocean_month","all",.true.,"none",2
#"ocean_model","S_diffx_2d","S_diffx_2d","ocean_month","all",.true.,"none",2
#"ocean_model","S_diffy_2d","S_diffy_2d","ocean_month","all",.true.,"none",2
"ocean_model","age","age","ocean_month","all",.true.,"none",2
"ocean_model","temp_h","temp_h","ocean_month","all",.true.,"none",2
"ocean_model","salt_h","salt_h","ocean_month","all",.true.,"none",2
"ocean_model","h","h","ocean_month","all",.true.,"none",2
"ocean_model","Kd_work","Kd_work","ocean_month","all",.true.,"none",2
# Monthly Avg Stress
#"ocean_model","taux", "taux", "ocean_month","all",.true.,"none",2
#"ocean_model","tauy", "tauy", "ocean_month","all",.true.,"none",2
#===========================
# Mixed Layer and TKE Budget Terms:
#===========================
"ocean_model","h_ML","h_ML","ocean_daily","all",.true.,"none",2
"ocean_model","ustar_shelf","ustar_shelf","ocean_daily","all",.true.,"none",2
"ocean_model","SST","SST","ocean_daily","all",.true.,"none",2
"ocean_model","SST","SST_max","ocean_daily","all",max,"none",2
"ocean_model","SST","SST_min","ocean_daily","all",min,"none",2
"ocean_model","SSS","SSS","ocean_daily","all",.true.,"none",2
"ocean_model","SSH","SSH","ocean_daily","all",.true.,"none",2
"ocean_model","SSU","SSU","ocean_daily","all",.true.,"none",2
"ocean_model","SSV","SSV","ocean_daily","all",.true.,"none",2
"ocean_model","melt","melt","ocean_daily","all",.true.,"none",2
"ocean_model","uhbt","uhbt","ocean_daily","all",.true.,"none",2
"ocean_model","vhbt","vhbt","ocean_daily","all",.true.,"none",2
#==================
# Surface Forcing:
#=================
"ocean_model","taux", "taux", "ocean_daily","all",.true.,"none",2
"ocean_model","tauy", "tauy", "ocean_daily","all",.true.,"none",2
"ocean_model","ustar", "ustar", "ocean_daily","all",.true.,"none",2
"ocean_model","PRCmE", "PmE", "ocean_daily","all",.true.,"none",2
#"ocean_model","SW", "SW", "ocean_daily","all",.true.,"none",2
#"ocean_model","LwLatSens","LwLatSens","ocean_daily","all",.true.,"none",2
#"ocean_model","salt_flux","salt_flux","ocean_daily","all",.true.,"none",2
"ocean_model","exch_vel_t", "exch_vel_t", "ocean_daily","all",.true.,"none",2
"ocean_model","exch_vel_s", "exch_vel_s", "ocean_daily","all",.true.,"none",2
"ocean_model","frazil", "frazil", "ocean_daily","all",.true.,"none",2
#"ocean_model","net_heat_surface","net_heat_surface","ocean_daily","all",.true.,"none",2
#
19 changes: 19 additions & 0 deletions examples/ocean_SIS/Amery/field_table
Original file line number Diff line number Diff line change
@@ -0,0 +1,19 @@
# specific humidity for moist runs
## check if "profile_type", "fixed", "surface_value=3.e-6" is needed for coupled model
"TRACER", "atmos_mod", "sphum"
"longname", "specific humidity"
"units", "kg/kg"
"profile_type", "fixed", "surface_value=3.e-6" /
# prognostic cloud scheme tracers
"TRACER", "atmos_mod", "liq_wat"
"longname", "cloud liquid specific humidity"
"units", "kg/kg" /
"TRACER", "atmos_mod", "ice_wat"
"longname", "cloud ice water specific humidity"
"units", "kg/kg" /
"TRACER", "atmos_mod", "cld_amt"
"longname", "cloud fraction"
"units", "none" /
"TRACER", "land_mod", "sphum"
"longname", "specific humidity"
"units", "kg/kg" /
121 changes: 121 additions & 0 deletions examples/ocean_SIS/Amery/input.nml
Original file line number Diff line number Diff line change
@@ -0,0 +1,121 @@
&MOM_input_nml
output_directory = 'RESTART/',
input_filename = 'n'
restart_input_dir = 'INPUT/',
restart_output_dir = 'RESTART/',
parameter_filename = 'MOM_input',
'MOM_override'
/

&coupler_nml
months = 12,
days = 0,
current_date = 1992,1,1,0,0,0,
calendar = 'JULIAN',
dt_cpld = 3600,
dt_atmos = 3600,
do_atmos = .true.,
do_land = .false.,
do_ice = .true.,
do_ocean = .true.,
do_flux = .true.,
atmos_npes = 32,
ocean_npes = 240,
concurrent = .true.
use_lag_fluxes=.true.
check_stocks = 0
/

&diag_manager_nml
max_axes = 100,
max_num_axis_sets = 100,
max_input_fields = 699
max_output_fields = 699
mix_snapshot_average_fields=.false.
conserve_water=.true.
/

&flux_exchange_nml
debug_stocks = .FALSE.
divert_stocks_report = .TRUE.
do_area_weighted_flux = .FALSE.
/

&fms_io_nml
fms_netcdf_restart=.true.
threading_read='multi'
threading_write='single'
fileset_write='single'
max_files_r = 200
max_files_w = 200
/

&fms_nml
clock_grain='ROUTINE'
clock_flags='NONE'
domains_stack_size = 5000000
stack_size =0
/

&ice_albedo_nml
t_range = 10.
/

&ice_model_nml
layout=8,4
nsteps_dyn=144
nsteps_adv=2
num_part = 6
wd_turn = 0.0
spec_ice=.false.
ice_bulk_salin = 0.010
alb_sno = 0.85 ! keep CM2 setting
alb_ice = 0.65 ! keep CM2 setting
t_range_melt = 1.0 ! NOTE: CM2 uses 1.0
heat_rough_ice = 5.0e-4
cm2_bugs = .false.
do_icebergs = .false.
atmos_winds=.true.
add_diurnal_sw = .false.
do_ice_limit=.false.
max_ice_limit=10.0
channel_viscosity=5.e-5
chan_cfl_limit=0.001
h_lo_lim = 1.e-10
io_layout=1,1
/

&icebergs_nml
verbose=.FALSE.
traj_sample_hrs=0
verbose_hrs=2400
/

&monin_obukhov_nml
stable_option = 2,
rich_crit = 10.0,
zeta_trans = 0.5
/

&ocean_albedo_nml
ocean_albedo_option = 2
/

&ocean_rough_nml
rough_scheme = 'beljaars'
/

&sat_vapor_pres_nml
construct_table_wrt_liq = .true.
construct_table_wrt_liq_and_ice = .true.
/

&surface_flux_nml
gust_min = 1.e-10
/

&xgrid_nml
make_exchange_reproduce = .false.
interp_method = 'second_order'
/

0 comments on commit d1d22af

Please sign in to comment.