simchan_yook@colostate_edu
Member
Hi all,
I'm trying to run a simulation with CAM5 coupled to a slab ocean in RCE setting.
I'm using QSC5 compset in cesm2_1_2 with resolution of f09_f09_mg17.
To run a preindustrial like RCE simulation following Popke et al. (2013), I've made changes on the CAM namelist:
omega=0
seasalt_emis_scale = 0.0
ch4vmr = 650.0e-9
co2vmr = 278.0e-6
f11vmr = 0.0
f12vmr = 0.0
flbc_list = ' '
n2ovmr = 270.0e-9
I also set the solar zenith angle to be spatially uniform but temporally varying with repeated diurnal cycle:
src/share/util/shr_orb_mod.F90
shr_orb_cosz = -cos((jday-floor(jday))*2.0_SHR_KIND_R8*pi)
I set the insolation to 1069.3W/m2 to make the global average of daily insolation to be 340 W/m2.
I did that via making alteration in the variable SOLIN from /src/physics/rrtmg/radsw.F90
solin(i) = sum(sfac(:)*solar_band_irrad(:)) * eccf * coszrs(i) * 1069.3_r8/1365.0_r8
And I set the global Q-flux to be zero.
Then, my model crashed or showed unrealistically high surface temperature values (>315 K) probably due to a runaway greenhouse effect.
(I've checked global average of FSNT and it was around 340 W/m^2, although it varied throughout the simulation).
Is this "too warm" atmosphere caused by the zero Q-flux? If so, how should I set the Q-flux to prevent the model from the "crash"?
Or are there any other factors that I should modify?
Thank you in advance
Shim
I'm trying to run a simulation with CAM5 coupled to a slab ocean in RCE setting.
I'm using QSC5 compset in cesm2_1_2 with resolution of f09_f09_mg17.
To run a preindustrial like RCE simulation following Popke et al. (2013), I've made changes on the CAM namelist:
omega=0
seasalt_emis_scale = 0.0
ch4vmr = 650.0e-9
co2vmr = 278.0e-6
f11vmr = 0.0
f12vmr = 0.0
flbc_list = ' '
n2ovmr = 270.0e-9
I also set the solar zenith angle to be spatially uniform but temporally varying with repeated diurnal cycle:
src/share/util/shr_orb_mod.F90
shr_orb_cosz = -cos((jday-floor(jday))*2.0_SHR_KIND_R8*pi)
I set the insolation to 1069.3W/m2 to make the global average of daily insolation to be 340 W/m2.
I did that via making alteration in the variable SOLIN from /src/physics/rrtmg/radsw.F90
solin(i) = sum(sfac(:)*solar_band_irrad(:)) * eccf * coszrs(i) * 1069.3_r8/1365.0_r8
And I set the global Q-flux to be zero.
Then, my model crashed or showed unrealistically high surface temperature values (>315 K) probably due to a runaway greenhouse effect.
(I've checked global average of FSNT and it was around 340 W/m^2, although it varied throughout the simulation).
Is this "too warm" atmosphere caused by the zero Q-flux? If so, how should I set the Q-flux to prevent the model from the "crash"?
Or are there any other factors that I should modify?
Thank you in advance
Shim