) # compile the source code to working directory $GFDL_WORK/codebase exp_name = 'held_suarez_default' exp = Experiment(exp_name, codebase=cb) diag = DiagTable() diag.add_file('atmos_monthly', 30, 'days', time_units='days') diag.add_field('dynamics', 'ps', time_avg=True) diag.add_field('dynamics', 'bk') diag.add_field('dynamics', 'pk') diag.add_field('dynamics', 'ucomp', time_avg=True) diag.add_field('dynamics', 'vcomp', time_avg=True) diag.add_field('dynamics', 'temp', time_avg=True) diag.add_field('dynamics', 'vor', time_avg=True) diag.add_field('dynamics', 'div', time_avg=True) exp.diag_table = diag namelist = Namelist({ 'main_nml': { 'dt_atmos': 600, 'days': 30, 'calendar': 'thirty_day', 'current_date': [2000, 1, 1, 0, 0, 0] }, 'atmosphere_nml': { 'idealized_moist_model': False # False for Newtonian Cooling. True for Isca/Frierson }, 'spectral_dynamics_nml': { 'damping_order': 4, # default: 2 'water_correction_limit': 200.e2, # default: 0
experiments = [('qflux_on', { 'mixed_layer_nml': { 'do_qflux': True }, 'qflux_nml': { 'qflux_amp': 30.0, 'qflux_width': 16.0 } }), ('qflux_off', { 'mixed_layer_nml': { 'do_qflux': False } })] exps = [] for exp_name, parameters in experiments: exp = Experiment(exp_name, codebase=codebase) # Get some basic paramters from config.py. We'll copy the Namelist, # DiagTable and input files list so for multiple Experiments in the # same file, each has their own version that can be edited independently. exp.diag_table = diag_table.copy() exp.namelist = namelist.copy() exp.inputfiles = inputfiles[::] exp.set_resolution('T42', 40) exp.update_namelist(parameters) exps.append(exp) run_cli(exps)