undo scaling and Tini introduced

Author: milankl

src/default_parameters.jl

@@ -4,6 +4,7 @@
 
     Tprog=T                   # number format for prognostic variables
     Tcomm=Tprog               # number format for ghost-point copies
+    Tini=Tprog                # number format to reduce precision for initial conditions
 
     # DOMAIN RESOLUTION AND RATIO
     nx::Int=100                         # number of grid cells in x-direction
@@ -157,6 +158,7 @@ Creates a Parameter struct with following options and default values
 
     Tprog=T                   # number format for prognostic variables
     Tcomm=Tprog               # number format for ghost-point copies
+    Tini=Tprog                # number format to reduce precision for initial conditions
 
     # DOMAIN RESOLUTION AND RATIO
     nx::Int=100                         # number of grid cells in x-direction

src/ghost_points.jl

@@ -42,10 +42,11 @@ function remove_halo(   u::Array{T,2},
     @unpack halo,haloη,halosstx,halossty = S.grid
     @unpack scale_inv = S.constants
 
+    # undo scaling as well
     ucut = scale_inv*u[halo+1:end-halo,halo+1:end-halo]
     vcut = scale_inv*v[halo+1:end-halo,halo+1:end-halo]
     ηcut = η[haloη+1:end-haloη,haloη+1:end-haloη]
-    sstcut = sst[halosstx+1:end-halosstx,halossty+1:end-halossty]
+    sstcut = scale_inv*sst[halosstx+1:end-halosstx,halossty+1:end-halossty]
 
     return ucut,vcut,ηcut,sstcut
 end

src/initial_conditions.jl

@@ -28,6 +28,7 @@ function initial_conditions(::Type{T},S::ModelSetup) where {T<:AbstractFloat}
     ## PROGNOSTIC VARIABLES U,V,η
     @unpack nux,nuy,nvx,nvy,nx,ny = S.grid
     @unpack initial_cond = S.parameters
+    @unpack Tini = S.parameters
 
     if initial_cond == "rest"
 
@@ -146,10 +147,11 @@ function initial_conditions(::Type{T},S::ModelSetup) where {T<:AbstractFloat}
     end
 
     # Convert to number format T
-    sst = T.(sst)
-    u = T.(u)
-    v = T.(v)
-    η = T.(η)
+    # allow for comparable initial conditions via Tini
+    sst = T.(Tini.(sst))
+    u = T.(Tini.(u))
+    v = T.(Tini.(v))
+    η = T.(Tini.(η))
 
     #TODO SST INTERPOLATION
     u,v,η,sst = add_halo(u,v,η,sst,S)