import numpy as np import pdb import os from gfdl.experiment import Experiment, DiagTable num_cores = 1 run_idb=True #Step 1: Do a single experiment to compile code then run remaining experiments off this compile. compiled_exp = Experiment('debug',run_idb=run_idb) compiled_exp.inputfiles = [os.path.join(os.getcwd(),'/scratch/pm366/GFDL_2013_FMS/GFDLmoistModel/input/ozone_1990.nc')] #Step 3: Set up the diag table diag = DiagTable() #specify time periods diag.add_file('atmos_daily', 1, 'days', time_units='days') diag.add_file('atmos_monthly', -1, 'days', time_units='days') #Both daily and monthly data - don't specify files and all will be used diag.add_field('dynamics', 'slp', time_avg=True) diag.add_field('dynamics', 'height', time_avg=True) diag.add_field('mixed_layer', 't_surf', time_avg=True) diag.add_field('dynamics', 'ucomp', time_avg=True) diag.add_field('dynamics', 'vcomp', time_avg=True) diag.add_field('dynamics', 'omega', time_avg=True) diag.add_field('dynamics', 'temp', time_avg=True)
import numpy as np import os from gfdl.experiment import Experiment, DiagTable import f90nml #Define our base experiment to compile base_dir = os.getcwd() GFDL_BASE = os.environ['GFDL_BASE'] baseexp = Experiment('mima_test_experiment', overwrite_data=False) baseexp.inputfiles = [ os.path.join(GFDL_BASE, 'input/rrtm_input_files/ozone_1990.nc') ] #Tell model how to write diagnostics diag = DiagTable() diag.add_file('atmos_monthly', 30, 'days', time_units='days') #Tell model which diagnostics to write diag.add_field('dynamics', 'ps', time_avg=True) diag.add_field('dynamics', 'bk') diag.add_field('dynamics', 'pk') diag.add_field('atmosphere', 'precipitation', time_avg=True) diag.add_field('mixed_layer', 't_surf', time_avg=True) diag.add_field('dynamics', 'sphum', time_avg=True) 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)
import numpy as np import os from gfdl.experiment import Experiment, DiagTable import f90nml #Define our base experiment to compile base_dir=os.getcwd() GFDL_BASE = os.environ['GFDL_BASE'] baseexp = Experiment('realistic_continents_test_experiment', overwrite_data=False) baseexp.inputfiles = [os.path.join(GFDL_BASE,'input/land_masks/era_land_t42.nc'),os.path.join(GFDL_BASE,'input/rrtm_input_files/ozone_1990.nc'), os.path.join(base_dir,'input/ami_qflux_ctrl_ice_4320.nc'), os.path.join(base_dir,'input/siconc_clim_amip.nc')] #Tell model how to write diagnostics diag = DiagTable() diag.add_file('atmos_monthly', 30, 'days', time_units='days') #Tell model which diagnostics to write diag.add_field('dynamics', 'ps', time_avg=True) diag.add_field('dynamics', 'bk') diag.add_field('dynamics', 'pk') diag.add_field('atmosphere', 'precipitation', time_avg=True) diag.add_field('mixed_layer', 't_surf', time_avg=True) diag.add_field('dynamics', 'sphum', time_avg=True) 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)
import numpy as np import os from gfdl.experiment import Experiment, DiagTable import f90nml #Define our base experiment to compile base_dir = os.getcwd() GFDL_BASE = os.environ['GFDL_BASE'] baseexp = Experiment('giant_planet_test_experiment', overwrite_data=False) #Tell model how to write diagnostics diag = DiagTable() diag.add_file('atmos_monthly', 30, 'days', time_units='days') #Tell model which diagnostics to write 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) baseexp.use_diag_table(diag) #Turn off the full, slow radiation scheme compilation baseexp.disable_rrtm()
import numpy as np import os from gfdl.experiment import Experiment, DiagTable import f90nml #Define our base experiment to compile base_dir=os.getcwd() GFDL_BASE = os.environ['GFDL_BASE'] baseexp = Experiment('bucket_test_experiment_3', overwrite_data=False) #Add any input files that are necessary for a particular experiment. baseexp.inputfiles = [os.path.join(base_dir,'input/land.nc'),os.path.join(GFDL_BASE,'input/rrtm_input_files/ozone_1990.nc')] #Tell model how to write diagnostics diag = DiagTable() diag.add_file('atmos_monthly', 30, 'days', time_units='days') #Tell model which diagnostics to write diag.add_field('dynamics', 'ps', time_avg=True) diag.add_field('dynamics', 'bk') diag.add_field('dynamics', 'pk') diag.add_field('atmosphere', 'precipitation', time_avg=True) diag.add_field('atmosphere', 'bucket_depth', time_avg=True) diag.add_field('mixed_layer', 't_surf', time_avg=True) diag.add_field('dynamics', 'sphum', time_avg=True) 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)
import numpy as np import os from gfdl.experiment import Experiment, DiagTable import f90nml #Define our base experiment to compile base_dir=os.getcwd() GFDL_BASE = os.environ['GFDL_BASE'] baseexp = Experiment('held_suarez_test_case', overwrite_data=False) #Tell model how to write diagnostics diag = DiagTable() diag.add_file('atmos_monthly', 30, 'days', time_units='days') #Tell model which diagnostics to write 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) baseexp.use_diag_table(diag) #Turn off the full, slow radiation scheme compilation baseexp.disable_rrtm()
import numpy as np import os from gfdl.experiment import Experiment, DiagTable import f90nml #Define our base experiment to compile base_dir = os.getcwd() GFDL_BASE = os.environ['GFDL_BASE'] baseexp = Experiment('variable_co2_grey', overwrite_data=False) baseexp.inputfiles = [os.path.join(base_dir, 'input/co2.nc')] #Tell model how to write diagnostics diag = DiagTable() diag.add_file('atmos_monthly', 30, 'days', time_units='days') #Tell model which diagnostics to write diag.add_field('dynamics', 'ps', time_avg=True) diag.add_field('dynamics', 'bk') diag.add_field('dynamics', 'pk') diag.add_field('atmosphere', 'precipitation', time_avg=True) diag.add_field('mixed_layer', 't_surf', time_avg=True) diag.add_field('dynamics', 'sphum', time_avg=True) 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) diag.add_field('two_stream', 'co2', time_avg=True)
import numpy as np import os from gfdl.experiment import Experiment, DiagTable import f90nml #Define our base experiment to compile base_dir = os.getcwd() GFDL_BASE = os.environ['GFDL_BASE'] baseexp = Experiment('frierson_test_experiment', overwrite_data=False) #Tell model how to write diagnostics diag = DiagTable() diag.add_file('atmos_monthly', 30, 'days', time_units='days') #Tell model which diagnostics to write diag.add_field('dynamics', 'ps', time_avg=True) diag.add_field('dynamics', 'bk') diag.add_field('dynamics', 'pk') diag.add_field('atmosphere', 'precipitation', time_avg=True) diag.add_field('mixed_layer', 't_surf', time_avg=True) diag.add_field('dynamics', 'sphum', time_avg=True) 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) baseexp.use_diag_table(diag) #Turn off the full, slow radiation scheme compilation
import numpy as np import os from gfdl.experiment import Experiment, DiagTable baseexp = Experiment('top_down_test', overwrite_data=True) #s Define input files for experiment - by default they are found in exp_dir/input/ #s Define srcmods - by default they are found in exp_dir/srcmods/ diag = DiagTable() #diag.add_file('atmos_hourly', 1, 'hours', time_units='days') #diag.add_file('atmos_hourly_mk2', 1, 'hours', time_units='days') diag.add_file('atmos_daily', 1, 'days', time_units='days') #diag.add_file('atmos_daily_mk2', 1, 'days', time_units='days') # Define diag table entries diag.add_field('dynamics', 'ps', time_avg=True, files=['atmos_daily']) diag.add_field('dynamics', 'bk', time_avg=True, files=['atmos_daily']) diag.add_field('dynamics', 'pk', time_avg=True, files=['atmos_daily']) diag.add_field('dynamics', 'vor', time_avg=True, files=['atmos_daily']) diag.add_field('dynamics', 'div', time_avg=True, files=['atmos_daily']) diag.add_field('dynamics', 'ucomp', time_avg=True, files=['atmos_daily']) diag.add_field('dynamics', 'vcomp', time_avg=True, files=['atmos_daily']) diag.add_field('dynamics', 'temp', time_avg=True, files=['atmos_daily']) diag.add_field('atmosphere', 'rh', time_avg=True, files=['atmos_daily']) diag.add_field('dynamics', 'slp', time_avg=True, files=['atmos_daily']) diag.add_field('dynamics', 'omega', time_avg=True, files=['atmos_daily'])
import numpy as np import os from gfdl.experiment import Experiment, DiagTable import f90nml #Define our base experiment to compile base_dir=os.getcwd() GFDL_BASE = os.environ['GFDL_BASE'] baseexp = Experiment('variable_co2_rrtm', overwrite_data=False) baseexp.inputfiles = [os.path.join(GFDL_BASE,'input/rrtm_input_files/ozone_1990.nc'),os.path.join(base_dir,'input/co2.nc')] #Tell model how to write diagnostics diag = DiagTable() diag.add_file('atmos_monthly', 30, 'days', time_units='days') #Tell model which diagnostics to write diag.add_field('dynamics', 'ps', time_avg=True) diag.add_field('dynamics', 'bk') diag.add_field('dynamics', 'pk') diag.add_field('atmosphere', 'precipitation', time_avg=True) diag.add_field('mixed_layer', 't_surf', time_avg=True) diag.add_field('dynamics', 'sphum', time_avg=True) 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) diag.add_field('rrtm_radiation', 'co2', time_avg=True)