Avoid zeros in jacobian when reaching tmp limits

Author: nichollsh

res/config/planets/toi561b.toml

@@ -15,7 +15,7 @@ title = "TOI-561 b"
     wind_speed      = 2.0
 
 [files]
-    input_sf        = "res/spectral_files/Dayspring/48/Dayspring.sf"
+    input_sf        = "greygas"
     input_star      = "res/stellar_spectra/sun.txt"
     output_dir      = "out/"
 
@@ -33,7 +33,7 @@ title = "TOI-561 b"
     max_runtime     = 400
     num_levels      = 45
     continua        = true
-    rayleigh        = true
+    rayleigh        = false
     cloud           = false
     overlap_method  = "ee"
     real_gas        = false
@@ -46,10 +46,10 @@ title = "TOI-561 b"
     rainout         = false
     solution_type   = 3
     solver          = "newton"
-    dx_max          = 400.0
-    initial_state   = ["iso","2000"] #["ncdf", "out/atm.nc"]
-    linesearch      = 0
-    easy_start      = false
+    dx_max          = 300.0
+    initial_state   = ["iso","3000"] #["ncdf", "out/atm.nc"]
+    linesearch      = 2
+    easy_start      = true
     converge_atol   = 1.0e-2
     converge_rtol   = 1.0e-3
     perturb_all     = false

src/solver.jl

@@ -178,7 +178,7 @@ module solver
         # Execution parameters
         # --------------------
         #    padding
-        tmp_pad::Float64 =  10.0        # do not allow the solver to get closer than this to tmp_floor
+        tmp_pad::Float64 =  5.0        # do not allow the solver to get closer than this to tmp_floor
 
         #    easy_start
         easy_incr::Float64 = 2.0        # Factor by which to increase easy_sf at each step
@@ -251,7 +251,6 @@ module solver
         # Calculate the (remaining) temperatures from known temperatures
         function _set_tmps!(_x::Array{Float64,1})
             # Read new guess
-            clamp!(_x, atmos.tmp_floor+1.0, atmos.tmp_ceiling-1.0)
             for i in 1:atmos.nlev_c
                 atmos.tmp[i] = _x[i]
             end