Changeset 935 for codes/icosagcm/devel/src/kernels_hex/caldyn_fast.k90

Timestamp:

07/03/19 17:15:11 (5 years ago)

Author:

dubos

Message:

devel : interfaces for caldyn_fast and caldyn_slow_hydro

File:

• codes/icosagcm/devel/src/kernels_hex/caldyn_fast.k90 (modified) (2 diffs)

codes/icosagcm/devel/src/kernels_hex/caldyn_fast.k90

@@ -16 +16 @@
          !DIR$ SIMD
          DO ij=ij_begin, ij_end
-            berni(ij,l) = .5*(geopot(ij,l)+geopot(ij,l+1))
+            Phi_ik = .5*(geopot(ij,l)+geopot(ij,l+1))
+            berni(ij,l) = Phi_ik
          END DO
       END DO
@@ -23 +24 @@
          !DIR$ SIMD
          DO ij=ij_begin, ij_end
-            berni(ij,l) = .5*(geopot(ij,l)+geopot(ij,l+1))
-            berni(ij,l) = berni(ij,l) + pk(ij,l)*(cpp-theta(ij,l,1)) ! Gibbs = Cp.T-Ts = T(Cp-s)
+            Phi_ik = .5*(geopot(ij,l)+geopot(ij,l+1))
+            berni(ij,l) = Phi_ik + pk(ij,l)*(cpp-theta(ij,l,1)) ! Gibbs = Cp.T-Ts = T(Cp-s)
          END DO
       END DO
    CASE(thermo_variable_Cp)
-      ! thermodynamics with variable Cp
-      ! Cp(T) = Cp0 * (T/T0)^nu
-      ! => h = Cp(T).T/(nu+1)
       DO l = ll_begin, ll_end
          !DIR$ SIMD
          DO ij=ij_begin, ij_end
-            berni(ij,l) = .5*(geopot(ij,l)+geopot(ij,l+1))
+            Phi_ik = .5*(geopot(ij,l)+geopot(ij,l+1))
+            ! thermodynamics with variable Cp
+            ! Cp(T) = Cp0 * (T/T0)^nu
+            ! => h = Cp(T).T/(nu+1)
             cp_ik = cpp*(pk(ij,l)/Treff)**nu
-            berni(ij,l) = berni(ij,l) + pk(ij,l)*(cp_ik/(nu+1.)-theta(ij,l,1)) ! Gibbs = h-Ts = T(Cp/(nu+1)-s)
+            berni(ij,l) = Phi_ik + pk(ij,l)*(cp_ik/(nu+1.)-theta(ij,l,1)) ! Gibbs = h-Ts = T(Cp/(nu+1)-s)
          END DO
       END DO