(* Content-type: application/mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 6.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       145,          7]
NotebookDataLength[     23925,        529]
NotebookOptionsPosition[     22905,        491]
NotebookOutlinePosition[     23240,        506]
CellTagsIndexPosition[     23197,        503]
WindowFrame->Normal
ContainsDynamic->False*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{"Model", " ", "from", " ", "Boeker"}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"(*", " ", "parameters", " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"{", 
     RowBox[{"als", ",", "ta", ",", "ala", ",", "ta2", ",", "ala2"}], "}"}], 
    " ", "=", 
    RowBox[{"{", " ", 
     RowBox[{
     "0.11", ",", " ", "0.53", ",", "0.30", ",", "0.06", ",", "0.31"}], 
     "}"}]}], "\[IndentingNewLine]", 
   RowBox[{"cvec", "=", "2.7"}]}]}]], "Input",
 CellChangeTimes->{{3.433760725981313*^9, 3.433760734358175*^9}, {
  3.433760773517944*^9, 3.43376078913193*^9}, {3.433760822710953*^9, 
  3.433760925054064*^9}, {3.433761039214822*^9, 3.433761040056542*^9}, {
  3.433762148773091*^9, 3.433762197951057*^9}, {3.433762279330711*^9, 
  3.433762279600665*^9}, {3.433762319850661*^9, 3.433762455569988*^9}, {
  3.43376251215958*^9, 3.433762573490525*^9}, {3.43376261532931*^9, 
  3.433762618048963*^9}, {3.433762666369145*^9, 3.433762669048928*^9}, {
  3.433762703678677*^9, 3.433762706168726*^9}, {3.433762829429775*^9, 
  3.43376285464806*^9}, {3.433762907168054*^9, 3.433763006277419*^9}, {
  3.433763239267047*^9, 3.433763257456043*^9}, {3.433763402635514*^9, 
  3.433763426985435*^9}, {3.433763480575048*^9, 3.433763482585093*^9}, {
  3.433763533405022*^9, 3.433763549803319*^9}, {3.433763582094673*^9, 
  3.433763617324424*^9}, {3.433763687274136*^9, 3.433763897952946*^9}, {
  3.433764044492373*^9, 3.433764051342828*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0.11`", ",", "0.53`", ",", "0.3`", ",", "0.06`", ",", "0.31`"}], 
  "}"}]], "Output",
 CellChangeTimes->{
  3.433761637616286*^9, {3.433762149951985*^9, 3.43376219824713*^9}, 
   3.433762280058683*^9, {3.433762325214707*^9, 3.433762455868886*^9}, 
   3.43376251287335*^9, {3.433762557763707*^9, 3.43376258140784*^9}, 
   3.433762618382325*^9, 3.433762669467396*^9, 3.433762707445916*^9, {
   3.433762830230106*^9, 3.433762859111004*^9}, {3.433762908182954*^9, 
   3.433763006635845*^9}, {3.43376324100246*^9, 3.433763257882562*^9}, {
   3.433763403372797*^9, 3.433763427356822*^9}, 3.43376348896404*^9, 
   3.433763550711402*^9, {3.433763603714973*^9, 3.433763618162343*^9}, {
   3.433763666401475*^9, 3.433763898907053*^9}, 3.433764051949919*^9, 
   3.433764211563368*^9, 3.433764509242845*^9, 3.434198774614166*^9, 
   3.434198874993128*^9, 3.434199219113088*^9, 3.43419962847497*^9}],

Cell[BoxData["2.7`"], "Output",
 CellChangeTimes->{
  3.433761637616286*^9, {3.433762149951985*^9, 3.43376219824713*^9}, 
   3.433762280058683*^9, {3.433762325214707*^9, 3.433762455868886*^9}, 
   3.43376251287335*^9, {3.433762557763707*^9, 3.43376258140784*^9}, 
   3.433762618382325*^9, 3.433762669467396*^9, 3.433762707445916*^9, {
   3.433762830230106*^9, 3.433762859111004*^9}, {3.433762908182954*^9, 
   3.433763006635845*^9}, {3.43376324100246*^9, 3.433763257882562*^9}, {
   3.433763403372797*^9, 3.433763427356822*^9}, 3.43376348896404*^9, 
   3.433763550711402*^9, {3.433763603714973*^9, 3.433763618162343*^9}, {
   3.433763666401475*^9, 3.433763898907053*^9}, 3.433764051949919*^9, 
   3.433764211563368*^9, 3.433764509242845*^9, 3.434198774614166*^9, 
   3.434198874993128*^9, 3.434199219113088*^9, 3.43419962856895*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"combi", "=", 
  RowBox[{
   RowBox[{"ta", "*", 
    RowBox[{"(", 
     RowBox[{"1", "-", "als"}], ")"}]}], "+", 
   RowBox[{"(", 
    RowBox[{"1", "-", "ala"}], ")"}]}]}]], "Input",
 CellChangeTimes->{{3.433761232594208*^9, 3.433761278996307*^9}, {
  3.434199110247169*^9, 3.434199119359761*^9}}],

Cell[BoxData["1.1717`"], "Output",
 CellChangeTimes->{
  3.433760927487608*^9, {3.433761042554125*^9, 3.433761050394458*^9}, 
   3.43376128661201*^9, {3.433761634825095*^9, 3.433761640014834*^9}, 
   3.433762151934177*^9, {3.433762185092364*^9, 3.433762201704271*^9}, 
   3.43376228300267*^9, {3.433762327283577*^9, 3.433762440277561*^9}, 
   3.433762515138015*^9, {3.433762559548645*^9, 3.433762584282128*^9}, 
   3.43376262044052*^9, 3.433762673235194*^9, {3.433762832081505*^9, 
   3.433762860411264*^9}, {3.433762909798429*^9, 3.433763009016544*^9}, {
   3.433763242825131*^9, 3.433763259679583*^9}, {3.433763406028495*^9, 
   3.433763429274337*^9}, 3.43376349060657*^9, 3.43376355275635*^9, {
   3.433763605376203*^9, 3.433763620093722*^9}, {3.433763668045092*^9, 
   3.433763901181102*^9}, 3.433764054464126*^9, 3.433764213423388*^9, 
   3.433764514499942*^9, 3.434198775054836*^9, 3.434198877878827*^9, 
   3.434199120151657*^9, {3.434199219508913*^9, 3.434199245584779*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", "constants", " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"sigma", "=", "5.67*^-8"}], "\[IndentingNewLine]", 
   RowBox[{"T0", "=", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"1370", "/", "4"}], "/", "sigma"}], ")"}], "^", "0.25"}]}], 
   "\[IndentingNewLine]", 
   RowBox[{"gamma", "=", 
    RowBox[{
     RowBox[{"cvec", "/", 
      RowBox[{"T0", "^", "3"}]}], "/", "sigma"}]}]}]}]], "Input",
 CellChangeTimes->{{3.433760930107686*^9, 3.433761014336613*^9}, {
  3.433761045854717*^9, 3.433761046550682*^9}, {3.433761558112498*^9, 
  3.43376155848237*^9}, {3.433762711188859*^9, 3.433762718058422*^9}, {
  3.433763956864729*^9, 3.433763979952701*^9}, {3.433764057320762*^9, 
  3.433764096382225*^9}}],

Cell[BoxData["5.67`*^-8"], "Output",
 CellChangeTimes->{{3.433761047601139*^9, 3.433761052017112*^9}, 
   3.43376129046794*^9, 3.433761567948202*^9, {3.433761642755586*^9, 
   3.433761657703389*^9}, 3.433762153392406*^9, {3.433762186586839*^9, 
   3.433762204116418*^9}, 3.43376228458596*^9, {3.433762328652349*^9, 
   3.433762442037402*^9}, 3.433762516760392*^9, {3.433762561390675*^9, 
   3.433762586310362*^9}, 3.433762622271688*^9, 3.433762675198543*^9, {
   3.433762712952498*^9, 3.43376271856325*^9}, {3.433762833703173*^9, 
   3.433762861516506*^9}, {3.433762911266669*^9, 3.433763011451442*^9}, {
   3.433763244302463*^9, 3.433763261134733*^9}, {3.433763407481496*^9, 
   3.433763430897374*^9}, 3.433763492079009*^9, 3.433763554415772*^9, {
   3.433763606855188*^9, 3.433763621330022*^9}, {3.43376366951418*^9, 
   3.433763902808086*^9}, 3.433763961085038*^9, {3.433764061302895*^9, 
   3.433764097599237*^9}, 3.433764215245058*^9, 3.433764516406924*^9, 
   3.434198775275832*^9, 3.434198880092783*^9, {3.434199219544953*^9, 
   3.434199249249781*^9}}],

Cell[BoxData["278.78498481970155`"], "Output",
 CellChangeTimes->{{3.433761047601139*^9, 3.433761052017112*^9}, 
   3.43376129046794*^9, 3.433761567948202*^9, {3.433761642755586*^9, 
   3.433761657703389*^9}, 3.433762153392406*^9, {3.433762186586839*^9, 
   3.433762204116418*^9}, 3.43376228458596*^9, {3.433762328652349*^9, 
   3.433762442037402*^9}, 3.433762516760392*^9, {3.433762561390675*^9, 
   3.433762586310362*^9}, 3.433762622271688*^9, 3.433762675198543*^9, {
   3.433762712952498*^9, 3.43376271856325*^9}, {3.433762833703173*^9, 
   3.433762861516506*^9}, {3.433762911266669*^9, 3.433763011451442*^9}, {
   3.433763244302463*^9, 3.433763261134733*^9}, {3.433763407481496*^9, 
   3.433763430897374*^9}, 3.433763492079009*^9, 3.433763554415772*^9, {
   3.433763606855188*^9, 3.433763621330022*^9}, {3.43376366951418*^9, 
   3.433763902808086*^9}, 3.433763961085038*^9, {3.433764061302895*^9, 
   3.433764097599237*^9}, 3.433764215245058*^9, 3.433764516406924*^9, 
   3.434198775275832*^9, 3.434198880092783*^9, {3.434199219544953*^9, 
   3.434199249264052*^9}}],

Cell[BoxData["2.1977210482137055`"], "Output",
 CellChangeTimes->{{3.433761047601139*^9, 3.433761052017112*^9}, 
   3.43376129046794*^9, 3.433761567948202*^9, {3.433761642755586*^9, 
   3.433761657703389*^9}, 3.433762153392406*^9, {3.433762186586839*^9, 
   3.433762204116418*^9}, 3.43376228458596*^9, {3.433762328652349*^9, 
   3.433762442037402*^9}, 3.433762516760392*^9, {3.433762561390675*^9, 
   3.433762586310362*^9}, 3.433762622271688*^9, 3.433762675198543*^9, {
   3.433762712952498*^9, 3.43376271856325*^9}, {3.433762833703173*^9, 
   3.433762861516506*^9}, {3.433762911266669*^9, 3.433763011451442*^9}, {
   3.433763244302463*^9, 3.433763261134733*^9}, {3.433763407481496*^9, 
   3.433763430897374*^9}, 3.433763492079009*^9, 3.433763554415772*^9, {
   3.433763606855188*^9, 3.433763621330022*^9}, {3.43376366951418*^9, 
   3.433763902808086*^9}, 3.433763961085038*^9, {3.433764061302895*^9, 
   3.433764097599237*^9}, 3.433764215245058*^9, 3.433764516406924*^9, 
   3.434198775275832*^9, 3.434198880092783*^9, {3.434199219544953*^9, 
   3.434199249318314*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
   "Solve", " ", "the", " ", "energy", " ", "flow", " ", "equations", " ", 
    "by", " ", "iteration"}], " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"{", 
     RowBox[{"Gtemp", ",", "Htemp"}], "}"}], "=", 
    RowBox[{"{", 
     RowBox[{"0.5", ",", "0.5"}], "}"}]}], "\[IndentingNewLine]", 
   RowBox[{"Do", "[", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      RowBox[{"G", "=", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"1", "-", "ala", "-", 
          RowBox[{"ta2", "*", 
           RowBox[{"Htemp", "^", "4"}]}]}], ")"}], "^", "0.25"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"Hp4", "=", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"combi", "+", 
          RowBox[{"gamma", "*", 
           RowBox[{"(", 
            RowBox[{"Gtemp", "-", "Htemp"}], ")"}]}]}], ")"}], "/", 
        RowBox[{"(", 
         RowBox[{"1", "-", "ala2", "+", "ta2"}], ")"}]}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"H", "=", 
       RowBox[{
        RowBox[{"(", "Hp4", ")"}], "^", "0.25"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"Gtemp", ",", "Htemp"}], "}"}], "=", 
       RowBox[{"{", 
        RowBox[{"G", ",", "H"}], "}"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{"Git", "[", "i", "]"}], ",", 
         RowBox[{"Hit", "[", "i", "]"}]}], "}"}], "=", 
       RowBox[{"{", 
        RowBox[{"G", ",", "H"}], "}"}]}]}], ",", "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "20"}], "}"}]}], "]"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"Tatmos", ",", "Tsurf"}], "}"}], "=", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"T0", "*", "G"}], ",", 
      RowBox[{"T0", "*", "H"}]}], "}"}]}]}]}]], "Input",
 CellChangeTimes->{{3.433761060623834*^9, 3.433761220685442*^9}, {
  3.433761300086512*^9, 3.43376137590465*^9}, {3.433761661674566*^9, 
  3.433761663154273*^9}, {3.433761717526026*^9, 3.4337617495032*^9}, {
  3.433761780480653*^9, 3.433761805282953*^9}, {3.43376206183175*^9, 
  3.433762063141904*^9}, {3.433762094832682*^9, 3.43376211550602*^9}, {
  3.433762294841458*^9, 3.433762295447172*^9}, {3.433763642074512*^9, 
  3.433763642499299*^9}, {3.433764303169037*^9, 3.433764346058909*^9}, {
  3.433764550737977*^9, 3.433764563264568*^9}, {3.434199144971214*^9, 
  3.434199207524476*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0.5`", ",", "0.5`"}], "}"}]], "Output",
 CellChangeTimes->{{3.433761196842268*^9, 3.433761221350867*^9}, {
   3.433761292413004*^9, 3.433761302318066*^9}, {3.433761369242177*^9, 
   3.4337613767705*^9}, 3.433761571649964*^9, 3.433761663533405*^9, {
   3.433761742236086*^9, 3.433761749865443*^9}, 3.433761806438009*^9, 
   3.433762063473076*^9, {3.433762096417957*^9, 3.433762115912978*^9}, 
   3.433762155556228*^9, {3.433762187704506*^9, 3.433762206230115*^9}, {
   3.433762286312324*^9, 3.433762296395434*^9}, {3.433762330309409*^9, 
   3.433762443546047*^9}, 3.433762518946389*^9, {3.433762563340222*^9, 
   3.433762588421018*^9}, 3.433762624507562*^9, 3.433762675871326*^9, 
   3.433762721708473*^9, {3.433762835361957*^9, 3.433762862764231*^9}, {
   3.433762912755718*^9, 3.43376301713779*^9}, {3.433763245867337*^9, 
   3.433763262846632*^9}, {3.433763409280536*^9, 3.433763432513398*^9}, 
   3.433763493723583*^9, 3.43376355581626*^9, {3.433763608493482*^9, 
   3.43376390513217*^9}, 3.433763964377705*^9, {3.433764063973965*^9, 
   3.433764085149058*^9}, 3.433764218409181*^9, {3.433764316548959*^9, 
   3.433764346630007*^9}, 3.433764519328882*^9, {3.433764551472335*^9, 
   3.43376456379604*^9}, 3.434198775401656*^9, 3.4341988826268*^9, {
   3.434199208631167*^9, 3.434199219958176*^9}, {3.434199252167634*^9, 
   3.434199265864002*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"248.50770121520316`", ",", "288.3127069858217`"}], "}"}]], "Output",\

 CellChangeTimes->{{3.433761196842268*^9, 3.433761221350867*^9}, {
   3.433761292413004*^9, 3.433761302318066*^9}, {3.433761369242177*^9, 
   3.4337613767705*^9}, 3.433761571649964*^9, 3.433761663533405*^9, {
   3.433761742236086*^9, 3.433761749865443*^9}, 3.433761806438009*^9, 
   3.433762063473076*^9, {3.433762096417957*^9, 3.433762115912978*^9}, 
   3.433762155556228*^9, {3.433762187704506*^9, 3.433762206230115*^9}, {
   3.433762286312324*^9, 3.433762296395434*^9}, {3.433762330309409*^9, 
   3.433762443546047*^9}, 3.433762518946389*^9, {3.433762563340222*^9, 
   3.433762588421018*^9}, 3.433762624507562*^9, 3.433762675871326*^9, 
   3.433762721708473*^9, {3.433762835361957*^9, 3.433762862764231*^9}, {
   3.433762912755718*^9, 3.43376301713779*^9}, {3.433763245867337*^9, 
   3.433763262846632*^9}, {3.433763409280536*^9, 3.433763432513398*^9}, 
   3.433763493723583*^9, 3.43376355581626*^9, {3.433763608493482*^9, 
   3.43376390513217*^9}, 3.433763964377705*^9, {3.433764063973965*^9, 
   3.433764085149058*^9}, 3.433764218409181*^9, {3.433764316548959*^9, 
   3.433764346630007*^9}, 3.433764519328882*^9, {3.433764551472335*^9, 
   3.43376456379604*^9}, 3.434198775401656*^9, 3.4341988826268*^9, {
   3.434199208631167*^9, 3.434199219958176*^9}, {3.434199252167634*^9, 
   3.434199265877616*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
   "Check", " ", "that", " ", "the", " ", "iteration", " ", "converged", " ", 
    "to", " ", "a", " ", "solution"}], " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"thline", "=", 
    RowBox[{"AbsoluteThickness", "[", "2", "]"}]}], "\[IndentingNewLine]", 
   RowBox[{"epi", "=", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"Blue", ",", 
        RowBox[{"Text", "[", 
         RowBox[{"\"\<H : surface\>\"", ",", 
          RowBox[{"{", 
           RowBox[{"10", ",", "1.00"}], "}"}]}], "]"}]}], "}"}], ",", 
      "\[IndentingNewLine]", 
      RowBox[{"{", 
       RowBox[{"Red", ",", 
        RowBox[{"Text", "[", 
         RowBox[{"\"\<G : atmosphere\>\"", ",", 
          RowBox[{"{", 
           RowBox[{"10", ",", "0.98"}], "}"}]}], "]"}]}], "}"}]}], "}"}]}], 
   "\[IndentingNewLine]", 
   RowBox[{"ListPlot", "[", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Table", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{"i", ",", 
           RowBox[{"Git", "[", "i", "]"}]}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{"i", ",", "1", ",", "20"}], "}"}]}], "]"}], ",", 
       "\[IndentingNewLine]", 
       RowBox[{"Table", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{"i", ",", 
           RowBox[{"Hit", "[", "i", "]"}]}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{"i", ",", "1", ",", "20"}], "}"}]}], "]"}]}], "}"}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{"PlotRange", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0", ",", "20"}], "}"}], ",", "Automatic"}], "}"}]}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"Red", ",", "thline"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"Blue", ",", "thline"}], "}"}]}], "}"}]}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{"Joined", "\[Rule]", "True"}], ",", "\[IndentingNewLine]", 
     RowBox[{"Epilog", "\[Rule]", "epi"}]}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.433761813131198*^9, 3.433762050221621*^9}, {
  3.433762220840897*^9, 3.433762249585194*^9}, {3.433764365158352*^9, 
  3.433764377770907*^9}, {3.433764411821616*^9, 3.43376448750042*^9}, {
  3.434198764935616*^9, 3.434198769106155*^9}, {3.43419891696136*^9, 
  3.434199068009307*^9}}],

Cell[BoxData[
 RowBox[{"AbsoluteThickness", "[", "2", "]"}]], "Output",
 CellChangeTimes->{{3.433761938831176*^9, 3.43376195995574*^9}, {
   3.433762002651796*^9, 3.433762015362087*^9}, {3.433762050673389*^9, 
   3.433762068961873*^9}, 3.433762120414708*^9, 3.433762161212147*^9, {
   3.433762192232417*^9, 3.433762250553584*^9}, {3.433762308301287*^9, 
   3.433762385684227*^9}, {3.433762428621056*^9, 3.433762446237253*^9}, 
   3.433762523144894*^9, {3.433762565920109*^9, 3.433762590500327*^9}, 
   3.433762627558958*^9, 3.433762678779986*^9, 3.433762724162212*^9, {
   3.433762838430146*^9, 3.433762864799591*^9}, {3.433762915218375*^9, 
   3.433763018957449*^9}, {3.433763248957221*^9, 3.43376326574717*^9}, {
   3.43376341165928*^9, 3.433763434869318*^9}, 3.433763497816434*^9, 
   3.433763567698607*^9, 3.433763611388249*^9, {3.433763648232731*^9, 
   3.433763674979517*^9}, 3.433764224752771*^9, 3.433764378889594*^9, {
   3.433764483481775*^9, 3.433764488249258*^9}, 3.433764523273517*^9, {
   3.434198770512607*^9, 3.434198775460965*^9}, 3.43419888617363*^9, {
   3.434199027768611*^9, 3.434199069207075*^9}, 3.434199220732354*^9, 
   3.434199272676988*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"RGBColor", "[", 
      RowBox[{"0", ",", "0", ",", "1"}], "]"}], ",", 
     RowBox[{"Text", "[", 
      RowBox[{"\<\"H : surface\"\>", ",", 
       RowBox[{"{", 
        RowBox[{"10", ",", "1.`"}], "}"}]}], "]"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"RGBColor", "[", 
      RowBox[{"1", ",", "0", ",", "0"}], "]"}], ",", 
     RowBox[{"Text", "[", 
      RowBox[{"\<\"G : atmosphere\"\>", ",", 
       RowBox[{"{", 
        RowBox[{"10", ",", "0.98`"}], "}"}]}], "]"}]}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{{3.433761938831176*^9, 3.43376195995574*^9}, {
   3.433762002651796*^9, 3.433762015362087*^9}, {3.433762050673389*^9, 
   3.433762068961873*^9}, 3.433762120414708*^9, 3.433762161212147*^9, {
   3.433762192232417*^9, 3.433762250553584*^9}, {3.433762308301287*^9, 
   3.433762385684227*^9}, {3.433762428621056*^9, 3.433762446237253*^9}, 
   3.433762523144894*^9, {3.433762565920109*^9, 3.433762590500327*^9}, 
   3.433762627558958*^9, 3.433762678779986*^9, 3.433762724162212*^9, {
   3.433762838430146*^9, 3.433762864799591*^9}, {3.433762915218375*^9, 
   3.433763018957449*^9}, {3.433763248957221*^9, 3.43376326574717*^9}, {
   3.43376341165928*^9, 3.433763434869318*^9}, 3.433763497816434*^9, 
   3.433763567698607*^9, 3.433763611388249*^9, {3.433763648232731*^9, 
   3.433763674979517*^9}, 3.433764224752771*^9, 3.433764378889594*^9, {
   3.433764483481775*^9, 3.433764488249258*^9}, 3.433764523273517*^9, {
   3.434198770512607*^9, 3.434198775460965*^9}, 3.43419888617363*^9, {
   3.434199027768611*^9, 3.434199069207075*^9}, 3.434199220732354*^9, 
   3.43419927282211*^9}],

Cell[BoxData[
 GraphicsBox[{{}, {}, 
   {RGBColor[1, 0, 0], AbsoluteThickness[2], 
    LineBox[{{1., 0.9134637176231641}, {2., 0.8823998896906985}, {3., 
     0.8952031207844309}, {4., 0.8892421119354257}, {5., 
     0.8925335632364138}, {6., 0.8907599187660983}, {7., 
     0.8917414811894989}, {8., 0.8912047340761261}, {9., 
     0.8915003199725885}, {10., 0.8913381576202649}, {11., 
     0.8914273085938943}, {12., 0.8913783529719996}, {13., 
     0.891405253035731}, {14., 0.8913904771584693}, {15., 
     0.8913985949150626}, {16., 0.8913941355473788}, {17., 
     0.8913965853749966}, {18., 0.8913952395640945}, {19., 
     0.8913959788971078}, {20., 0.8913955727419122}}]}, 
   {RGBColor[0, 0, 1], AbsoluteThickness[2], 
    LineBox[{{1., 1.1179922464766747`}, {2., 0.9906028163569839}, {3., 
     1.0563547646178382`}, {4., 1.0217886923292394`}, {5., 
     1.0408895347715006`}, {6., 1.0304643799112934`}, {7., 
     1.0362083838704415`}, {8., 1.033058311386857}, {9., 
     1.0347904321779833`}, {10., 1.0338393696809214`}, {11., 
     1.0343619882674522`}, {12., 1.0340749293093812`}, {13., 
     1.0342326401017834`}, {14., 1.0341460048726758`}, {15., 
     1.0341935996228557`}, {16., 1.0341674535590695`}, {17., 
     1.0341818171510868`}, {18., 1.0341739264680694`}, {19., 
     1.0341782612673718`}, {20., 1.0341758799251044`}}]}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  AxesOrigin->{0, 0.9},
  Epilog->{{
     RGBColor[0, 0, 1], 
     InsetBox[
      FormBox["\"H : surface\"", TraditionalForm], {10, 1.}]}, {
     RGBColor[1, 0, 0], 
     InsetBox[
      FormBox["\"G : atmosphere\"", TraditionalForm], {10, 0.98}]}},
  PlotRange->{{0, 20}, {0.8823998896906985, 1.1179922464766747`}},
  PlotRangeClipping->True,
  PlotRangePadding->{Automatic, 
    Scaled[0.02]}]], "Output",
 CellChangeTimes->{{3.433761938831176*^9, 3.43376195995574*^9}, {
   3.433762002651796*^9, 3.433762015362087*^9}, {3.433762050673389*^9, 
   3.433762068961873*^9}, 3.433762120414708*^9, 3.433762161212147*^9, {
   3.433762192232417*^9, 3.433762250553584*^9}, {3.433762308301287*^9, 
   3.433762385684227*^9}, {3.433762428621056*^9, 3.433762446237253*^9}, 
   3.433762523144894*^9, {3.433762565920109*^9, 3.433762590500327*^9}, 
   3.433762627558958*^9, 3.433762678779986*^9, 3.433762724162212*^9, {
   3.433762838430146*^9, 3.433762864799591*^9}, {3.433762915218375*^9, 
   3.433763018957449*^9}, {3.433763248957221*^9, 3.43376326574717*^9}, {
   3.43376341165928*^9, 3.433763434869318*^9}, 3.433763497816434*^9, 
   3.433763567698607*^9, 3.433763611388249*^9, {3.433763648232731*^9, 
   3.433763674979517*^9}, 3.433764224752771*^9, 3.433764378889594*^9, {
   3.433764483481775*^9, 3.433764488249258*^9}, 3.433764523273517*^9, {
   3.434198770512607*^9, 3.434198775460965*^9}, 3.43419888617363*^9, {
   3.434199027768611*^9, 3.434199069207075*^9}, 3.434199220732354*^9, 
   3.434199273083623*^9}]
}, Open  ]]
},
WindowSize->{640, 750},
WindowMargins->{{216, Automatic}, {Automatic, 98}},
FrontEndVersion->"6.0 for Linux x86 (32-bit) (April 20, 2007)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[590, 23, 1511, 29, 99, "Input"],
Cell[2104, 54, 926, 15, 31, "Output"],
Cell[3033, 71, 832, 12, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[3902, 88, 320, 9, 32, "Input"],
Cell[4225, 99, 982, 14, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[5244, 118, 780, 19, 99, "Input"],
Cell[6027, 139, 1060, 15, 33, "Output"],
Cell[7090, 156, 1070, 15, 31, "Output"],
Cell[8163, 173, 1070, 15, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[9270, 193, 2496, 66, 231, "Input"],
Cell[11769, 261, 1389, 21, 31, "Output"],
Cell[13161, 284, 1420, 22, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[14618, 311, 2465, 66, 253, "Input"],
Cell[17086, 379, 1167, 17, 31, "Output"],
Cell[18256, 398, 1702, 35, 52, "Output"],
Cell[19961, 435, 2928, 53, 220, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)