units of swflux (issue #10)

Alexander-Barth · Oct 26, 2023 · bfece00 · bfece00
1 parent 5a050d7
commit bfece00
Show file tree

Hide file tree

Showing 3 changed files with 7 additions and 2 deletions.
diff --git a/src/Atmosphere/prepare_ecmwf.jl b/src/Atmosphere/prepare_ecmwf.jl
@@ -81,7 +81,7 @@ function prepare_ecmwf(
             ECMWFname = "",
             accumulation = true,
             output = "swflux",
-            scale  = 100.0/Δt_seconds*(24*3600.0),
+            scale  = 1.0/Δt_seconds,
         ),
         (
             Vname  = "swrad",

diff --git a/src/Atmosphere/prepare_gfs.jl b/src/Atmosphere/prepare_gfs.jl
@@ -454,6 +454,11 @@ function prepare_gfs(
 
                 # check sign
                 field = -(field - evaporation_rate)
+
+                density_fresh_water = 998 #  kg m⁻³
+
+                # kg m⁻² s⁻¹ / (kg m⁻³) = m s⁻¹
+                field = field / density_fresh_water
             elseif Vname == "swrad"
                 field = swrad_down - field
             else

diff --git a/src/metadata.jl b/src/metadata.jl
@@ -143,7 +143,7 @@ metadata = Dict(
         Tname = "swf_time",
         ncattrib = (
             long_name = "surface net freshwater flux (E-P)",
-            units = "centimeter day-1",
+            units = "m s-1",
             positive_value = "net evaporation",
             negative_value = "net precipitation",
         ),