Scheduled Downtime
On Tuesday 24 October 2023 @ 5pm MT the forums will be in read only mode in preparation for the downtime. On Wednesday 25 October 2023 @ 5am MT, this website will be down for maintenance and expected to return online later in the morning.
Normal Operations
The forums are back online with normal operations. If you notice any issues or errors related to the forums, please reach out to help@ucar.edu

CO2 run with prescribed SST

Hello All,

I'm trying to run a simulation with increasing CO2 and prescribed SST using CAM3. I've tried two ways: 1. with observed CO2 increase during end-1949 to 2005 and 2. with 1% increase per year.
I'm comparing these runs with another forced by prescribed SST and fixed CO2 concentration. However, it seems that the model is not using the CO2 values, because the output of global averaged air temperature at high and low levels are the very similar in all these runs. I would expect to see a significant increase in temperature at least for the 1%/yr CO2 run and I don't see it. The log file archive/CO2RUN/atm/atm.log.091207-230417 shows that
"ghg_surfvals_set: ncdate= 20060813 co2vmr= XXX" is increasing though.

I wonder if I'm missing any parameter when setting up the run. For example, in the 1%/yr run I set up:

- env_conf:
"setenv IPCC_MODE RAMP_CO2_ONLY"
"setenv RAMP_CO2_START_YMD 19490916"

- Buildnml_Prestage/cam.buildnml_prestage.csh
scenario_ghg = 'RAMP_CO2_ONLY'
ramp_co2_start_ymd = 19490916
co2vmr = 3.107e-4
ramp_co2_annual_rate = 1


Does anyone know of any other parameter that I need to set up? Or maybe the CO2 runs are only valid in the coupled model (I'm using the atmospheric component only)?

Thanks in advance,
Andrea.
 
Hi, I'm giving this a bump as I'm having a similar problem -- no observable increases. Only, in my atm.log.* file co2vmr stays constant throught the run at the initial co2vmr setpoint, 3.689000000000000E-004. (Also, I'm not prescribing the SST.)I've attached the atm_in file (CESM 1_0_4, -compset B_2000 -res f45_g37). Maybe I'm missing something in the namelist config, as I'm new at this.Thanks,Peter &aerodep_flx_nl aerodep_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           = 1.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 / 
 
Top