From b565ad89f97e3f48c0d9ee59cd9b32bdc5e876bf Mon Sep 17 00:00:00 2001
From: Cheng Gong <enigne@gmail.com>
Date: Fri, 12 Apr 2024 12:39:10 -0400
Subject: [PATCH] Add scale output depend on the result enum

---
 src/core/modules.jl | 17 +++++++++++++++--
 1 file changed, 15 insertions(+), 2 deletions(-)

diff --git a/src/core/modules.jl b/src/core/modules.jl
index 7067c76..7b9a447 100644
--- a/src/core/modules.jl
+++ b/src/core/modules.jl
@@ -153,13 +153,15 @@ function OutputResultsx(femmodel::FemModel, md::model, solutionkey::Symbol)# {{{
 		for i in 1:length(femmodel.results)
 			result = femmodel.results[i]
 			step   = femmodel.results[i].step
-			(output[step])[EnumToString(result.enum)] = result.value
+			scale = OutputEnumToScale(md, result.enum)
+			(output[step])[EnumToString(result.enum)] = result.value .* scale
 		end
 	else
 		output = Dict()
 		for i in 1:length(femmodel.results)
 			result = femmodel.results[i]
-			output[EnumToString(result.enum)] = result.value
+			scale = OutputEnumToScale(md, result.enum)
+			output[EnumToString(result.enum)] = result.value .* scale
 		end
 	end
 
@@ -421,3 +423,14 @@ function GetMaskOfIceVerticesLSMx0(femmodel::FemModel) #{{{
 
 	return nothing
 end#}}}
+function OutputEnumToScale(md::model, result::IssmEnum) #{{{
+	if (result == VxEnum) scale = md.constants.yts
+	elseif (result == VyEnum) scale = md.constants.yts
+	elseif (result == VzEnum) scale = md.constants.yts
+	elseif (result == VelEnum) scale = md.constants.yts
+	elseif (result == CalvingCalvingrateEnum) scale = md.constants.yts
+	else scale = 1.0
+	end
+
+	return scale
+end#}}}