normalization of Gn

martalmeida · Dec 12, 2019 · f19cebd · f19cebd
1 parent 84f8a04
commit f19cebd
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 0 deletions.
diff --git a/nmf3d/vertical_structure.py b/nmf3d/vertical_structure.py
@@ -106,6 +106,9 @@ def vse(Tprof,Plev,**kargs):
   nkMax=Plev.size
   Gn=Gn_all[::-1][:nkMax]
 
+  # normalize:
+  Gn*=2**.5
+
   # The equivalent heights re-ordered
   cond=S==0
   S[cond]=1

diff --git a/nmf3d_mat/vertical_structure.m b/nmf3d_mat/vertical_structure.m
@@ -105,6 +105,9 @@
 Gn=flipud(Gn_all);
 Gn=Gn(1:nkMax,:);
 
+% normalize:
+Gn=Gn*sqrt(2);
+
 % The equivalent heights re-ordered
 cond=S==0;
 S(cond)=1;