reset : true  w 50 scriptmode : direct
emitmenu _ stabconc
baseyear : 1850 stabconc&inyear : 2150 stabconc&atlevel : 500 

pan0 : blankplot 
blankplot  : title~CO2Emissions 
CO2Emissions : clear : scale~giga&ton&carbon~1000~1000
CO2Emissions.name2 : (fossil+landuse) 
CO2Emissions.legw : 120
CO2Emissions.xmin : 1900 CO2Emissions.xmax : 2250
CO2Emissions.ymin : 0 CO2Emissions.ymax : 10000

pan1 : blankplot 
blankplot  : title~CO2Concentration 
CO2Concentration : clear : scale~ppm~50~1
CO2Concentration.legw : 120
CO2Concentration.xmin : 1900 CO2Concentration.xmax : 2250
CO2Concentration.ymin : 300 CO2Concentration.ymax : 600

pan2 : blankplot 
blankplot  : title~Temperature
Temperature : clear : scale~degc~0.2001~1.0
Temperature.legw : 120
Temperature.xmin : 1900 Temperature.xmax : 2250
Temperature.ymin : -0.2 Temperature.ymax : 6.0

pan3 : blankplot 
blankplot  : title~RadiativeForcing
RadiativeForcing : clear : scale~w&per&m2~0.2001~1.0
RadiativeForcing.legw : 120
RadiativeForcing.xmin : 1900 RadiativeForcing.xmax : 2250
RadiativeForcing.ymin : 0.0 RadiativeForcing.ymax : 6.5

othgasemit : sresscale  $r : 0 
sub sresloop

othgasemit : 2000scale  $r : 50 $b : 0 $s : ' '
sub gcmloop

othgasemit : sresfix  $r : 200 $s : 
sub sresloop

othgasemit : 2000fix   $r : 150 $b : 0  $s : ' '
sub gcmloop

end

sresloop (
	$b : 0 
	sresmenu : A1B 
	nextsres ( 
	$b : ( $b + 40 ) $s : sresmenu
	sub gcmloop
	sresmenu : next  
	if ( sresmenu !=  IS92a )  ( sub nextsres )
) )

gcmloop (
	climodmenu : gfdl 
	$g : 0
	nextgcm ( 
	$g :  $g + 35 
	sub main
	climodmenu : next  
	if ( climodmenu != ipccsar ) ( sub nextgcm )  
) )

main (
	$n : climodmenu ~ othgasemit ~ $s ~ ~true~
	$c : $r ~ ' , ' ~ $g ~ ' , ' ~ $b 
	CO2Emissions   : ( ' curve~carboncycle.totemit~ ' ~ $n  ~ $c )
	CO2Concentration  : ( ' curve~carboncycle.atppm~ ' ~ $n ~ $c  )
	RadiativeForcing  : ( ' curve~radfor.totalrf~ ' ~ $n  ~ $c )
	Temperature   : ( ' curve~heatflux.glotemp~ ' ~ $n ~ $c )
) 