Try to make the solidus temperature can dependent on radius or pressure.


On Mon, May 19, 2014 at 5:29 PM, Wolfgang Bangerth <notifications@github.com
> wrote:

> In source/initial_conditions/solidus.cc:
>
> > +          	{
> > +                double T,p;
> > +                in>>T>>p;
> > +                if(!in.fail())
> > +                {
> > +                    //Unit switching
> > +                    if(T_Unit=="C")    T+=273.15; // Degree C to K
> > +                    if(P_Unit=="kbar") p*=1.e8;   // kbar to Pa
> > +                    if(P_Unit=="GPa")  p*=1.e9;   // GPa to Pa
> > +                    is_radius=false;              // Second column is pressure
> > +                    if(P_Unit=="km")
> > +                    {
> > +                        is_radius=true;           // Second column in radius instead of pressure
> > +                        p*=1.e3;                  // km to meters
> > +                    }
> > +                    T_array.push_back(T);
>
> Same here: error checking for P_Unit. That said, is there ever a case
> where you want P_unit to be kilometers? P is the pressure!?
>
> —
> Reply to this email directly or view it on GitHub<https://github.com/geodynamics/aspect/pull/19/files#r12819911>
> .
>



-- 
Siqi Zhang
Research Associate
Department of Earth and Planetary Sciences
Faculty of Science
Macquarie University

E7A, Rm 404
Macquarie University
Macquarie Park, NSW 2113