Open HenryDane opened 1 year ago
state = climlab.column_state(num_lev=20, num_lat=20)
state['q'] = climlab.radiation.water_vapor.ManabeWaterVapor(state=state).q
model = climlab.TimeDependentProcess(state=state, timestep=climlab.utils.constants.seconds_per_day)
model.add_subprocess('RRTMG', climlab.radiation.RRTMG(state=state))
model.add_subprocess('EC', climlab.convection.EmanuelConvection(state=state))
model.add_subprocess('Sun', climlab.radiation.insolation.FixedInsolation(state=state,
domains=state.Ts.domain))
#model.add_subprocess('cloud', climlab.convection.LargeScaleCloud(state=state))
model.compute_diagnostics()
model.integrate_years()
plt.pcolormesh(model.lat, model.lev, model.relative_humidity.T)
plt.gca().invert_yaxis();
plt.colorbar(label='RH')
plt.ylabel('Lev (hPa)'); plt.xlabel('Latitude (deg)')
The above produces a near crash (super odd values).
try looking at #153?