Changeset 5808

Timestamp:

Sep 12, 2025, 6:09:44 PM (7 weeks ago)

Author:

fhourdin

Message:

Sommet des poches froides

File:

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

LMDZ6/trunk/libf/phylmd/lmdz_wake_pkupper.f90

@@ -50 +50 @@
 REAL                                              :: ddd!!
 
+REAL :: www
 
 INTEGER, SAVE :: ipas=0
@@ -344 +345 @@
  kupper = 0
  
-IF (wk_pupper<1.) THEN
+IF (0.<wk_pupper .and. wk_pupper<1.) THEN
  ! Choose an integration bound well above wake top
   ! -----------------------------------------------------------------
@@ -358 +359 @@
   END DO
 
-ELSE
+ELSE IF (1.<=wk_pupper) THEN
   DO i=1, klon
      ! pupper(i) = wk_pupper*ptop(i)+(1.-wk_pupper)*ph(i, 1) 
      !  pupper(i) = min( wk_pupper*ptop(i)+(1.-wk_pupper)*ph(i, 1) , ptop(i)-50.) 
       pupper(i) = min( wk_pupper*ptop(i)+(1.-wk_pupper)*ph(i, 1) , ptop(i)-5000.) 
+  END DO
+ELSE
+  www=abs(wk_pupper)
+  DO i=1, klon
+     ! pupper(i) = wk_pupper*ptop(i)+(1.-wk_pupper)*ph(i, 1) 
+     !  pupper(i) = min( wk_pupper*ptop(i)+(1.-wk_pupper)*ph(i, 1) , ptop(i)-50.) 
+      pupper(i) = min(   exp(www*log(ptop(i))+(1.-www)*log(ph(i, 1)) ) , ptop(i)-5000.) 
   END DO
 END IF