Merge pull request #466 from matthewhoffman/landice/exo_convert_stiff…

…_log Apply log transform to stiffnessFactor conversion The Albany optimization now calculates stiffening_factor as a log, so conversion_exodus_init_to_mpasli_mesh.py has been updated accordingly.
MPAS-Dev · Jun 20, 2022 · 9510ed6 · 9510ed6
2 parents 1ae1309 + f09e724
commit 9510ed6
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/landice/mesh_tools_li/conversion_exodus_init_to_mpasli_mesh.py b/landice/mesh_tools_li/conversion_exodus_init_to_mpasli_mesh.py
@@ -217,6 +217,8 @@
                 dataset.variables[var_name][0,cellID_array-1] = np.exp(data_exo_layer) * 1000.0
             elif var_name == "muFriction":
                 dataset.variables[var_name][0,cellID_array-1] = np.exp(data_exo_layer)
+            elif var_name == "stiffnessFactor":
+                dataset.variables[var_name][0,cellID_array-1] = np.exp(data_exo_layer)
             elif var_name == "uReconstructX" or var_name == "uReconstructY":
                 dataset.variables[var_name][0,cellID_array-1, nVert] = data_exo_layer / (60. * 60. * 24 * 365)
             elif var_name == "thickness":