Changeset 4842 for LMDZ6

Timestamp:

Mar 3, 2024, 2:34:25 AM (9 months ago)

Author:

fhourdin

Message:

Suite poches

Location:

LMDZ6/trunk/libf/phylmd

Files:

• lmdz_wake.F90 (modified) (9 diffs)
• lmdz_wake_ini.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/lmdz_wake.F90

@@ -319 +319 @@
   REAL, DIMENSION (klon)                                :: h_zzz
 
-!print*,'WAKE LJYFo'
+!print*,'WAKE LJYFz'
 
   ! -------------------------------------------------------------------------
@@ -2422 +2422 @@
 USE lmdz_wake_ini , ONLY : RG
 USE lmdz_wake_ini , ONLY : hwmin
+USE lmdz_wake_ini , ONLY : iflag_wk_new_ptop
 
 IMPLICIT NONE
@@ -2457 +2458 @@
 REAL                                              :: ddd!!
 
-LOGICAL :: new_ptop
 
 INTEGER, SAVE :: ipas=0
@@ -2471 +2471 @@
 !print*,'compte=',compte
 
-new_ptop=.false.
-
- 
     ! Determine Ptop from buoyancy integral
     ! ---------------------------------------
@@ -2589 +2586 @@
     ! nouveau calcul de hw et ptop
     ! -----------------------------------------------------------------------
-if (new_ptop) then
+if (iflag_wk_new_ptop>0) then
     
     epsil = 0.05                            ! 5 pour cent
@@ -2611 +2608 @@
              endif
              int_dth(i,k) = int_dth(i,k-1) + ddd
+          !ELSE
+          !   int_dth(i,k) = 0.
           END IF
        END DO
@@ -2621 +2620 @@
        IF (wk_adv(i)) THEN
         frac_int_dth(i)=(1.-epsil)*int_dth(i,k_ptop_provis(i))
+       !ELSE
+       ! frac_int_dth(i)=0.
        ENDIF
     ENDDO
@@ -2626 +2627 @@
        DO i =1, klon
 !         print*,ipas,'yyy ',k,int_dth(i,k),frac_int_dth(i)
-          IF (wk_adv(i) .AND. int_dth(i,k)>=frac_int_dth(i)) THEN
-            ktop1(i) = min(k, k_ptop_provis(i))
-            !print*,ipas,'yyy ktop1= ',ktop1
-          END if
+          IF (wk_adv(i)) THEN
+            IF (int_dth(i,k)>=frac_int_dth(i)) THEN
+              ktop1(i) = min(k, k_ptop_provis(i))
+              !print*,ipas,'yyy ktop1= ',ktop1
+            ENDIF
+          ENDIF
        END DO
     END DO
@@ -2679 +2682 @@
 
 endif
+
+if (iflag_wk_new_ptop==2) then
+   do i=1,klon
+      ptop(i)=ptop1(i)
+      ktop(i)=ktop1(i)
+      hw_(i)=h_zzz(i)
+   enddo
+endif
+
 
  kupper = 0

LMDZ6/trunk/libf/phylmd/lmdz_wake_ini.F90

@@ -48 +48 @@
   INTEGER, SAVE, PROTECTED                                    :: iflag_wk_profile
   !$OMP THREADPRIVATE(iflag_wk_profile)
+
+  INTEGER, SAVE, PROTECTED                                    :: iflag_wk_new_ptop
+  !$OMP THREADPRIVATE(iflag_wk_new_ptop)
 
   REAL, SAVE, PROTECTED                                       :: wdensinit ! Minimum wake density used to restart wakes from a wake-free state
@@ -212 +215 @@
   CALL getin_p('iflag_wk_profile',iflag_wk_profile) ! switch between wdens prescribed 
                                                     ! and wdens prognostic
+  iflag_wk_profile = 0
+  CALL getin_p('iflag_wk_new_ptop',iflag_wk_new_ptop)
+
   tau_cv = 4000.
   CALL getin_p('tau_cv', tau_cv)