no NetCDF.close anymore

Author: milankl

src/DefaultParameters.jl

@@ -185,7 +185,7 @@ Creates a Parameter struct with following options and default values
 
     # SURFACE FORCING
     surface_forcing::Bool=false         # yes?
-    ωFη::Real=1.0                       # (annual) frequency [1/year]
+    ωFη::Real=1.0                       # frequency [1/year] for surfance forcing
     A::Real=3e-5                        # Amplitude [m/s]
     ϕk::Real=ϕ                          # Central latitude of Kelvin wave pumping
     wk::Real=10e3                       # width [m] in y of Gaussian used for surface forcing

src/InitialConditions.jl

@@ -53,15 +53,15 @@ function initial_conditions(::Type{T},S::ModelSetup) where {T<:AbstractFloat}
         end
 
         u = ncu.vars["u"][:,:,init_starti]
-        NetCDF.close(ncu)
+        # NetCDF.close(ncu)
 
         ncv = NetCDF.open(joinpath(inirunpath,"v.nc"))
         v = ncv.vars["v"][:,:,init_starti]
-        NetCDF.close(ncv)
+        # NetCDF.close(ncv)
 
         ncη = NetCDF.open(joinpath(inirunpath,"eta.nc"))
         η = ncη.vars["eta"][:,:,init_starti]
-        NetCDF.close(ncη)
+        # NetCDF.close(ncη)
 
         # remove singleton time dimension
         u = reshape(u,size(u)[1:2])
@@ -147,7 +147,7 @@ function initial_conditions(::Type{T},S::ModelSetup) where {T<:AbstractFloat}
     if initial_cond == "ncfile" && sst_initial == "restart"
         ncsst = NetCDF.open(joinpath(inirunpath,"sst.nc"))
         sst = ncsst.vars["sst"][:,:,init_starti]
-        NetCDF.close(ncsst)
+        # NetCDF.close(ncsst)
 
         sst = reshape(sst,size(sst)[1:2])
     end
@@ -164,11 +164,3 @@ function initial_conditions(::Type{T},S::ModelSetup) where {T<:AbstractFloat}
 
     return PrognosticVars{T}(u,v,η,sst)
 end
-
-# if injection_region == "south"
-#     sst_inj_region = sst_south
-# elseif injection_region == "west"
-#     sst_inj_region = sst_west
-# elseif injection_region == "rect"
-#     sst_inj_region = sst_rect
-# end

src/Output.jl

@@ -147,11 +147,11 @@ function output_close!(ncs::NcFiles,feedback::Feedback,S::ModelSetup)
     @unpack output = S.parameters
 
     if output
-        for nc in (ncs.u,ncs.v,ncs.η,ncs.sst,ncs.q,ncs.ζ,ncs.du,ncs.dv,ncs.dη)
-            if nc != nothing
-                NetCDF.close(nc)
-            end
-        end
+        # for nc in (ncs.u,ncs.v,ncs.η,ncs.sst,ncs.q,ncs.ζ,ncs.du,ncs.dv,ncs.dη)
+        #     if nc != nothing
+        #         NetCDF.close(nc)
+        #     end
+        # end
         println("All data stored.")
         write(feedback.progress_txt,"All data stored.")
         close(feedback.progress_txt)