fix default temperature gradient

Author: zhenwu0728

Project.toml

@@ -1,7 +1,7 @@
 name = "PlanktonIndividuals"
 uuid = "dfeb4ede-c75e-11e9-050d-41ab7a730e9e"
 authors = ["zhenwu0728 <zhenwu@mit.edu>", "gaelforget <gforget@mit.edu>"]
-version = "0.6.4"
+version = "0.6.5"
 
 [deps]
 Adapt = "79e6a3ab-5dfb-504d-930d-738a2a938a0e"

src/Parameters/param_default.jl

@@ -31,7 +31,7 @@ function default_temperature(grid, ΔT, iterations)
     end
     # vertical temperature gradient
     for j in grid.Nz-1:-1:1
-        CUDA.@allowscalar temp_day[:,:,j,:] .= temp_day[:,:,j+1,:] .+ 0.04 * grid.zC[j]
+        CUDA.@allowscalar temp_day[:,:,j,:] .= temp_day[:,:,j+1,:] .+ 4.0e-4 * grid.zC[j]
     end
     temp_domain = repeat(temp_day, outer = (1,1,1,total_days))
     return temp_domain