There are two ways, a quick and a flexible one. For both first determine the equivalent CO2 concentration in ppm and convert it into mixing ratio by dividing by 658157.16
The quick one:
In src_radiation.f90 after the ico2_rad select case assign your mixing ratio to the variable zqco2
CASE (10) ! RCP85 scenario (for 1950 <= zyear <= 2150) ! eff. CO2 (all GreeHouseGases (GHG) considered) zqco2 = ( – 4.0501145412E+07_ireals & + 7.9386473439E+04_ireals * zyear & – 5.8292720579E+01_ireals * zyear*zyear & + 1.9002921793E-02_ireals * zyear**3 & – 2.3202412328E-06_ireals * zyear**4) * 1.519E-06_irealsEND SELECT zqco2 = …. ! insert your mixing ratio here
The flexible one:
in organize_physics.f90 change
IF (ico2_rad.GT.10 ) ico2_rad = ico2_rad_dto
IF (ico2_rad.GT.11 ) ico2_rad = ico2_rad_d
in src_radiation change
zyearmin = 1950to
IF (ico2_rad < 11) THEN zyearmin = 1950 ELSE zyearmin = ! .... put a year lower than the start year of your simulation period hereand change
CASE (10) ! RCP85 scenario (for 1950 <= zyear <= 2150) ! eff. CO2 (all GreeHouseGases (GHG) considered) zqco2 = ( – 4.0501145412E+07_ireals & + 7.9386473439E+04_ireals * zyear & – 5.8292720579E+01_ireals * zyear*zyear & + 1.9002921793E-02_ireals * zyear**3 & – 2.3202412328E-06_ireals * zyear**4) * 1.519E-06_irealsEND SELECT to
CASE (10) ! RCP85 scenario (for 1950 <= zyear <= 2150) ! eff. CO2 (all GreeHouseGases (GHG) considered) zqco2 = ( – 4.0501145412E+07_ireals & + 7.9386473439E+04_ireals * zyear & – 5.8292720579E+01_ireals * zyear*zyear & + 1.9002921793E-02_ireals * zyear**3 & – 2.3202412328E-06_ireals * zyear**4) * 1.519E-06_ireals CASE (11) !… last millenium zqco2 = !… insert your mixing ratio here ENDSELECTIn the CCLM namelist PHYCTL set
ico2_rad=11In the same way you can define several different mixing ratios with additional numbers 12, 13, 14 etc.
Push notifications in your browser are not yet configured.
There are two ways, a quick and a flexible one. For both first determine the equivalent CO2 concentration in ppm and convert it into mixing ratio by dividing by 658157.16
The quick one:
END SELECT zqco2 = …. ! insert your mixing ratio hereIn src_radiation.f90 after the ico2_rad select case assign your mixing ratio to the variable zqco2
The flexible one:
toin organize_physics.f90 change
in src_radiation change
to
and change
END SELECT to In the CCLM namelist PHYCTL set In the same way you can define several different mixing ratios with additional numbers 12, 13, 14 etc.