0001 function conf = get_glimmer_default_config_new()
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015
0016
0017
0018
0019
0020 conf.sections=struct('names',{'Gen_pdd_climate','GLIMMER annual pdd','grid','options','isostasy','time','parameters','CF default','CF input','CF output','CF output'}...
0021 );
0022
0023 conf.Gen_pdd_climate = struct( ...
0024 'global_mean_precipitation' , 0.0, ...
0025 'global_precipitation_range' , 0.0, ...
0026 'global_mean_sl_temperature' , -12.0,...
0027 'global_sl_temperature_range' , 9.0,...
0028 'global_lapse_rate' , -0.0065, ...
0029 'snow_threshold' , 1.0, ...
0030 'pptm_file_name' , 'pptm_s1.txt', ...
0031 'sstm_file_name' , 'sstm_m12.txt');
0032 conf.GLIMMER_annual_pdd = struct( ...
0033 'dx' , 1.0, ...
0034 'dy' , 1.0, ...
0035 'ix' , 0.0, ...
0036 'iy' , -80, ...
0037 'nx' , 51, ...
0038 'ny' , 101, ...
0039 'wmax' , 0.6, ...
0040 'pddfac_ice' , 0.008, ...
0041 'pddfac_snow' , 0.003);
0042 conf.grid = struct( ...
0043 'ewn' , 258, ...
0044 'nsn' , 237, ...
0045 'upn' , 11, ...
0046 'dew' , 10000, ...
0047 'dns' , 10000);
0048 conf.time = struct( ...
0049 'tstart' , 0., ...
0050 'tend' , 30000., ...
0051 'dt' , 0.5, ...
0052 'ntem' , 1., ...
0053 'nvel' , 1., ...
0054 'niso' , 1.);
0055 conf.options = struct( ...
0056 'temperature' , 1, ...
0057 'flow_law' , 0, ...
0058 'marine_margin' , 3, ...
0059 'stress_calc' , 2, ...
0060 'evolution' , 0, ...
0061 'basal_water' , 0, ...
0062 'slip_coeff' , 4, ...
0063 'sliding_law' , 0, ...
0064 'vertical_integration' , 0, ...
0065 'hotstart' , 0, ...
0066 'topo_is_relaxed' , 1);
0067 conf.parameters = struct( ...
0068 'flow_factor' , 1, ...
0069 'geothermal' , -42e-3, ...
0070 'basal_tract_const' , 1.e-5, ...
0071 'basal_tract_max', 0.001,...
0072 'basal_tract_slope', 0.02,...
0073 'ice_limit' , 500, ...
0074 'calving_fraction' , 0.8,...
0075 'marine_limit' , -500);
0076 conf.isostasy = struct( ...
0077 'lithosphere' , 1, ...
0078 'asthenosphere' , 1);
0079 conf.CF_default = struct( ...
0080 'title' , 'Parameteric uncertainty test file', ...
0081 'comment', 'Comment', ...
0082 'name', 'filename.nc');
0083 conf.CF_input = struct( ...
0084 'name' , 'scand_aea_10k.nc', ...
0085 'time', 1);
0086 conf.CF_output = struct( ...
0087 'name' , '', ...
0088 'frequency', 1000, ...
0089 'variables' , {'thk','usurf','ablt','acab','bmlt','sstm','btrc','pptm','snow','artm','btemp','uvel','vvel','wvel','slc','topg','calving'});
0090 conf.CF_output2 = struct( ...
0091 'name' , '', ...
0092 'frequency' , 100, ...
0093 'variables' , {'thk'});