Changeset 807 for codes/icosagcm/devel/Python/dynamico/LAM.py

Timestamp:

02/08/19 18:53:34 (5 years ago)

Author:

jisesh

Message:

devel/Python : polar projection + Williamson test case 2 on LAM

File:

• codes/icosagcm/devel/Python/dynamico/LAM.py (modified) (1 diff)

codes/icosagcm/devel/Python/dynamico/LAM.py

@@ -17 +17 @@
             if ci > c1 : m[i]=0.
         return m
+    def relax_RSW(self, llm, step, flow):
+        beta_i, beta_e = self.beta_i, self.beta_e
+        m,u=flow
+        step.mass[:]        = beta_i*step.mass[:]        + (1.-beta_i)*m[:]
+        step.theta_rhodz[:] = beta_i*step.theta_rhodz[:] + (1.-beta_i)*m[:]
+        step.u[:]           = beta_e*step.u[:]           + (1.-beta_e)*u[:]
+    def relax_hydro(self, llm, step, flow):
+        beta_i, beta_e = self.beta_i, self.beta_e
+        m,S,u=flow
+        for l in range(llm):
+            step.mass[:,l]        = beta_i*step.mass[:,l]        + (1.-beta_i)*m[:,l]
+            step.theta_rhodz[:,l] = beta_i*step.theta_rhodz[:,l] + (1.-beta_i)*S[:,l]
+            step.u[:,l]           = beta_e*step.u[:,l]           + (1.-beta_e)*u[:,l]
     def relax(self, llm, step, flow):
         beta_i, beta_e = self.beta_i, self.beta_e