Changeset 5711 for LMDZ6/trunk/libf/phylmd/cv_routines.f90

Timestamp:

Jun 16, 2025, 6:24:37 PM (3 weeks ago)

Author:

yann meurdesoif

Message:

Convection GPU porting : Bug fix in cv_yield, scalar variable dpinv, depending of the column indice, is reused in other loop without been reintialized. May change result (only for iflag_con=4)

YM

File:

• LMDZ6/trunk/libf/phylmd/cv_routines.f90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/cv_routines.f90

@@ -1523 +1523 @@
       DO ij = 1, ncum
         IF (i<=inb(ij)) THEN
+          dpinv = 0.01/(ph(ij,i)-ph(ij,i+1))
           awat = elij(ij, k, i) - (1.-ep(ij,i))*clw(ij, i)
           awat = max(awat, 0.0)
@@ -1540 +1541 @@
       DO ij = 1, ncum
         IF ((i<=inb(ij)) .AND. (k<=inb(ij))) THEN
+          dpinv = 0.01/(ph(ij,i)-ph(ij,i+1))
           fq(ij, i) = fq(ij, i) + g*dpinv*ment(ij, k, i)*(qent(ij,k,i)-q(ij,i &
             ))
@@ -1551 +1553 @@
     DO ij = 1, ncum
       IF (i<=inb(ij)) THEN
+        dpinv = 0.01/(ph(ij,i)-ph(ij,i+1))
         fq(ij, i) = fq(ij, i) + sigd*evap(ij, i) + g*(mp(ij,i+1)*(qp(ij, &
           i+1)-q(ij,i))-mp(ij,i)*(qp(ij,i)-q(ij,i-1)))*dpinv