Changeset 4825 for LMDZ6

Timestamp:

Feb 16, 2024, 11:56:24 AM (10 months ago)

Author:

fhourdin

Message:

Correction plantage debug wake

Location:

LMDZ6/trunk/libf/phylmd

Files:

• lmdz_wake.F90 (modified) (4 diffs)
• yamada4.F90 (modified) (1 diff)

LMDZ6/trunk/libf/phylmd/lmdz_wake.F90

@@ -2476 +2476 @@
 REAL,     DIMENSION (klon)         :: sum_dth
 
+!INTEGER, SAVE :: compte=0
+
 ! LJYF : a priori z, dz sum_dth sont aussi des variables internes
 ! Les eliminer apres verification convergence numerique
+
+!compte=compte+1
+!print*,'compte=',compte
 
   delta_t_min = 0.2
@@ -2488 +2493 @@
 
     DO i = 1, klon
-      IF (wk_adv(i)) THEN
         ptop_provis(i) = ph(i, 1)
-      END IF
     END DO
 
     DO k = 2, klev
       DO i = 1, klon
+!       if (compte==92) then
+!            print*,'debug xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx'
+!            print*,'debug i =',i
+!            print*,'debug k =',k
+!            print*,'debug wk_adv(i) =',wk_adv(i)
+!            print*,'debug ptop_provis(i) =',ptop_provis(i)
+!            print*,'debug ph(i,1) =',ph(i,1)
+!            print*,'debug dth(i,k) =',dth(i,k)
+!            print*,'debug delta_t_min =',delta_t_min
+!            print*,'debug p(i,k-1) =',p(i,k-1)
+!            print*,'debug dth(i,k-1) =',dth(i,k-1)
+!            print*,'debug p(i,k) =',p(i,k)
+!       endif
         IF (wk_adv(i) .AND. ptop_provis(i)==ph(i,1) .AND. &
 ! LJYF changer :           dth(i,k)>=-delta_t_min .AND. dth(i,k-1)<-delta_t_min) THEN
@@ -2561 +2577 @@
 
     DO i = 1, klon
-      IF (wk_adv(i)) THEN
         ptop_new(i) = ptop(i)
-      END IF
     END DO
 
@@ -2580 +2594 @@
 
     DO i = 1, klon
-      IF (wk_adv(i)) THEN
         ptop(i) = ptop_new(i)
-      END IF
     END DO
 

LMDZ6/trunk/libf/phylmd/yamada4.F90

@@ -532 +532 @@
     !.......................................
 
-print*,'iflag_vdif_q2 =',iflag_vdif_q2
   IF (iflag_vdif_q2==1) THEN
     q2(1:ngrid, 1) = q2(1:ngrid, 2)