diag.add_field('atmosphere', 'rh', time_avg=True)
diag.add_field('atmosphere', 'convection_rain', time_avg=True)
diag.add_field('atmosphere', 'condensation_rain', time_avg=True)
diag.add_field('atmosphere', 'precipitation', time_avg=True)
#for pressure interpolations
diag.add_field('dynamics', 'bk', time_avg=True)
diag.add_field('dynamics', 'pk', time_avg=True)
diag.add_field('dynamics', 'ps', time_avg=True)
#Output data to test if model is stable after spin-up
diag.add_field('dynamics', 'sphum', time_avg=True, files=['atmos_monthly'])
compiled_exp.use_diag_table(diag)
compiled_exp.compile()
compiled_exp.clear_rundir()
#namelists are defined here: /scratch/pm366/GFDL_2013_FMS/GFDLmoistModel/src/extra/python/gfdl/templates/
compiled_exp.namelist['main_nml'] = {
'days' : 30,
'hours' : 0,
'minutes': 0,
'seconds': 0,
'dt_atmos':720,
'current_date' : [0001,1,1,0,0,0],
'calendar' : 'thirty_day'
}
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)
#Compile model if not already compiled
baseexp.compile()
#Empty the run directory ready to run
baseexp.clear_rundir()
#Define values for the 'core' namelist
baseexp.namelist['main_nml'] = f90nml.Namelist({
'days' : 30,
'hours' : 0,
'minutes': 0,
'seconds': 0,
'dt_atmos':720,
'current_date' : [0001,1,1,0,0,0],
'calendar' : 'thirty_day'
})