CO2 ramp experiment is not ramping

[peter_m_kalmus@jpl_nasa_gov]

Hi,I ran a control experiment and two CO2 ramp experiments at low resolution. However I'm missing something: the ramp experiments are not ramping. Surface temperatures do not increase over decades of model time, and co2vmr stays constant in the atm.log.* file. I've pasted in the atm_in file below (for 20% ramp!). This is CESM 1_0_4, -compset B_2000 -res f45_g37.I made a user_nl_cam with just three entries as reflected in atm_in:&camexp
scenario_ghg = 'RAMP_CO2_ONLY'
ramp_co2_annual_rate = 20.0
ramp_co2_start_ymd = 20000101
/ Thanks,Peter &aerodep_flx_nlaerodep_flx_datapath = '/data2/pkalmus/cesmtest/input//atm/cam/chem/trop_mozart_aero/aero'
aerodep_flx_file = 'aerosoldep_monthly_1849-2006_1.9x2.5_c090803.nc'
aerodep_flx_type = 'CYCLICAL'
aerodep_flx_ymd = 20000101
/
&cam_inparm
absems_data = '/data2/pkalmus/cesmtest/input//atm/cam/rad/abs_ems_factors_fastvx.c030508.nc'
bnd_topo = '/data2/pkalmus/cesmtest/input//atm/cam/topo/USGS-gtopo30_4x5_remap_c050520.nc'
cam_branch_file = ' '
del2coef = 3.e+5
div24del2flag = 2
dtime = 1800
fft_flt = 1
ncdata = '/data2/pkalmus/cesmtest/input//atm/cam/inic/fv/cami_0001-01-01_4x5_L26_c060608.nc'
npr_yz = 10,2,2,10
phys_loadbalance = 2
/
&chem_surfvals_nl
ch4vmr = 1760.0e-9
co2vmr = 368.9e-6
f11vmr = 653.45e-12
f12vmr = 535.0e-12
n2ovmr = 316.0e-9
ramp_co2_annual_rate = 20.0
ramp_co2_start_ymd = 20000101
scenario_ghg = 'RAMP_CO2_ONLY'
/
&cldfrc_nl
cldfrc_freeze_dry = .true.
cldfrc_ice = .false.
cldfrc_rhminh = 0.800D0
cldfrc_rhminl = 0.900D0
/
&cldsed_nl
cldsed_ice_stokes_fac = 1.0D0
/
&cldwat_nl
cldwat_conke = 5.0e-6
cldwat_icritc = 9.5e-6
cldwat_icritw = 2.0e-4
/
&gw_drag_nl
fcrit2 = 1.0
/
&hkconv_nl
hkconv_c0 = 2.0e-4
hkconv_cmftau = 1800.0D0
/
&phys_ctl_nl
cam_chempkg = 'none'
cam_physpkg = 'cam4'
conv_water_in_rad = 0
deep_scheme = 'ZM'
do_iss = .false.
do_tms = .false.
eddy_scheme = 'HB'
history_microphysics = .false.
microp_scheme = 'RK'
shallow_scheme = 'Hack'
srf_flux_avg = 0
/
&prescribed_aero_nl
prescribed_aero_datapath = '/data2/pkalmus/cesmtest/input//atm/cam/chem/trop_mozart_aero/aero'
prescribed_aero_file = 'aero_1.9x2.5_L26_2000clim_c090803.nc'
prescribed_aero_specifier = 'sulf:SO4', 'bcar1:CB1', 'bcar2:CB2', 'ocar1:OC1', 'ocar2:OC2', 'sslt1:SSLT01', 'sslt2:SSLT02',
'sslt3:SSLT03', 'sslt4:SSLT04', 'dust1:DST01', 'dust2:DST02', 'dust3:DST03', 'dust4:DST04'
prescribed_aero_type = 'CYCLICAL'
prescribed_aero_ymd = 20000101
/
&prescribed_ozone_nl
prescribed_ozone_datapath = '/data2/pkalmus/cesmtest/input//atm/cam/ozone'
prescribed_ozone_file = 'ozone_1.9x2.5_L26_2000clim_c090803.nc'
prescribed_ozone_name = 'O3'
prescribed_ozone_type = 'CYCLICAL'
prescribed_ozone_ymd = 20000101
/
&rad_cnst_nl
rad_climate = 'P_Q:H2O', 'D_O2:O2', 'D_CO2:CO2', 'D_ozone:O3', 'D_N2O:N2O', 'D_CH4:CH4', 'D_CFC11:CFC11',
'D_CFC12:CFC12', 'D_sulf:/data2/pkalmus/cesmtest/input//atm/cam/physprops/sulfate_camrt_c080918.nc',
'D_dust1:/data2/pkalmus/cesmtest/input//atm/cam/physprops/dust1_camrt_c080918.nc',
'D_dust2:/data2/pkalmus/cesmtest/input//atm/cam/physprops/dust2_camrt_c080918.nc',
'D_dust3:/data2/pkalmus/cesmtest/input//atm/cam/physprops/dust3_camrt_c080918.nc',
'D_dust4:/data2/pkalmus/cesmtest/input//atm/cam/physprops/dust4_camrt_c080918.nc',
'D_bcar1:/data2/pkalmus/cesmtest/input//atm/cam/physprops/bcpho_camrt_c080918.nc',
'D_bcar2:/data2/pkalmus/cesmtest/input//atm/cam/physprops/bcphi_camrt_c080918.nc',
'D_ocar1:/data2/pkalmus/cesmtest/input//atm/cam/physprops/ocpho_camrt_c080918.nc',
'D_ocar2:/data2/pkalmus/cesmtest/input//atm/cam/physprops/ocphi_camrt_c080918.nc',
'D_SSLTA:/data2/pkalmus/cesmtest/input//atm/cam/physprops/ssam_camrt_c080918.nc',
'D_SSLTC:/data2/pkalmus/cesmtest/input//atm/cam/physprops/sscm_camrt_c080918.nc'
/&solar_inparm
solar_const = 1361.27
solar_htng_spctrl_scl = .false.
/
&tropopause_nl
tropopause_climo_file = '/data2/pkalmus/cesmtest/input//atm/cam/chem/trop_mozart/ub/clim_p_trop.nc'
/
&zmconv_nl
zmconv_c0_lnd = 0.0035D0
zmconv_c0_ocn = 0.0035D0
zmconv_ke = 1.0E-6
/ 

 

[mai]

Is this a hybrid run? If so, did you set the RUN_STARTDATE in the env_run.xml file?

 

[peter_m_kalmus@jpl_nasa_gov]

After looking in the code I got the runs to ramp. Maybe taking a minute to write what I did can help someone else starting out. I set ramp_co2_start_ymd = 0and ramps happen. In this case, the code determines the first day of the run and starts the ramping then. Apparently, by default the first day of my B_2000 run was not 20000101 as I had assumed. I don't yet know what it means that aerodep, prescribed_aero, and prescribed_ozone are all set to 20000101... but do beware of possible inconsistencies with the date handling in the default B_2000 case.In the cpl.log.* file, there is the line (and many others like it, counting up by days):tStamp_write: model date =    10107However, the "date" field in the *cam2.h0.*.nc history files starts counting at year 2000 as expected. So it is obviously important to understand how time formatting is handled in the model.  

 

[peter_m_kalmus@jpl_nasa_gov]

Thanks for your response. It's a startup run, and the start date (returned by the subroutine get_start_date()) is 10101 at the start of the run, which is the default for RUN_STARTDATE in env_conf.xml. I did not know to set it -- maybe it would be good to validate unreasonably large discrepancies between YMD namelist variables and the model time. RUN_STARTDATE is used for both hybrid and startup runs.