Pre-Update Solubility in Atmosphere module

Octave/make_movie.m

@@ -6,7 +6,7 @@ filenames = sprintf('%s%s', fpath, '/data_snap*')
 files = glob(filenames);
 for i=1:numel(files)
   [~, name] = fileparts(files{i});
-  eval(sprintf('plot_snap(30,4000,-1,1,0,"%s", "-ascii");', files{i}));
+  eval(sprintf('plot_snap(30,3000,-1,1,0,"%s", "-ascii");', files{i}));
 %  eval(sprintf('%s = load("%s", "-ascii");', name, files{i}));
 endfor
 

src/atmosphere.f90

@@ -56,7 +56,7 @@ contains
     endif
 
     pi__ = acos(-1.0_dp)
-    FD=(pF%D**3)*pF%rho ! FD correspond to Final Dimensionalization; values in atmos_outgassing scaled up to 3D sphere, hence "*D^3"
+    FD = pF%D**3 * pF%rho ! FD correspond to Final Dimensionalization; values in atmos_outgassing scaled up to 3D sphere, hence "*D^3" -> kg
 
     M_w_loss = 0.0_dp
     M_c_loss = 0.0_dp

src/input.f90

@@ -829,7 +829,7 @@ contains
           pI%TinclTref = int(val)
           if (rank.eq.0) write(*,*) var,"=",pI%TinclTref
         case("Psurf")
-          pI%Psurf = val*1.0e+9_dp*pF%D**2/(pF%kappa*pF%eta_ref) 
+          pI%Psurf = val*1.0e+9_dp*pF%D**2/(pF%kappa*pF%eta_ref)  ! only for stress induced by overlying pressure
           if (rank.eq.0) write(*,*) var,"=",val," -> ",pI%Psurf," (pI)"
           pE%Psurf = val*1.0e+9_dp/(pF%D*pF%rho*pF%g)
           pM%Psurf = pE%Psurf

src/parameters.f90

@@ -140,7 +140,7 @@ module parameters
     real(dp)::a1=1.0_dp,a2=0.0_dp,a3=0.0_dp,a4=0.0_dp ! compressibility parameters fro depth-dependent alpha-function
     real(dp)::Gr=1.0_dp ! Grueneisen Parameter
     integer::TinclTref=1 ! TinclTref = 1: simulation T is T' + Tref, TinclTref = 0: simulation T is only T', -1: like 0 but Work for T instead of T', 2: like 1 but with Work with T' instead of T
-    real(dp)::Psurf=0.0_dp,mu=0.6_dp
+    real(dp)::Psurf=0.0_dp,mu=0.6_dp ! only for stress induced by overlying pressure
     integer::read_profs=0 ! read-in of profiles for density, alpha etc.
     real(dp)::X_FeM=0.1_dp ! 10 wt-% iron / Mg# 0.9
     real(dp)::X_M0=0.0_dp ! depression of liquidus melt temp by other minerals to solidus after Stixrude et al., 2014 (e.g. 0.21 for MgO and SiO2)