Changeset 855 for codes/icosagcm/devel

Timestamp:

05/05/19 23:11:40 (5 years ago)

Author:

dubos

Message:

devel : cleaner fix to duplicate nu in compute_fast.F90

File:

• codes/icosagcm/devel/src/dynamics/compute_caldyn_fast.F90 (modified) (2 diffs)

codes/icosagcm/devel/src/dynamics/compute_caldyn_fast.F90

@@ -24 +24 @@
 
     INTEGER :: i,j,ij,l
-    !REAL(rstd) :: cp_ik, qv, temp, chi, nu, due, due_right, due_lup, due_ldown
-    REAL(rstd) :: cp_ik, qv, temp, chi, due, due_right, due_lup, due_ldown
+    REAL(rstd) :: cp_ik, qv, temp, chi, log_p_preff, due, due_right, due_lup, due_ldown
 
     CALL trace_start("compute_caldyn_fast")
@@ -68 +67 @@
                 temp = pk(ij,l)
                 chi = log(temp/Treff)
-                nu = (chi*(cpp+qv*cppv)-theta(ij,l,1))/(Rd+qv*Rv) ! log(p/preff)
+                log_p_preff = (chi*(cpp+qv*cppv)-theta(ij,l,1))/(Rd+qv*Rv) ! log(p/preff)
                 berni(ij,l) = .5*(geopot(ij,l)+geopot(ij,l+1)) &
-                     + temp*(cpp*(1.-chi)+Rd*nu)
+                     + temp*(cpp*(1.-chi)+Rd*log_p_preff)
                 berniv(ij,l) = .5*(geopot(ij,l)+geopot(ij,l+1)) &
-                     + temp*(cppv*(1.-chi)+Rv*nu)
+                     + temp*(cppv*(1.-chi)+Rv*log_p_preff)
             END DO
           END SELECT