Index: LMDZ6/trunk/libf/phylmd/wake.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/wake.F90	(revision 4452)
+++ LMDZ6/trunk/libf/phylmd/wake.F90	(revision 4453)
@@ -27,5 +27,5 @@
   USE wake_ini_mod , ONLY : wake_ini
   USE wake_ini_mod , ONLY : prt_level,epsim1,RG,RD
-  USE wake_ini_mod , ONLY : stark, wdens_ref, coefgw, alpk, pupperbyphs
+  USE wake_ini_mod , ONLY : stark, wdens_ref, coefgw, alpk, wk_pupper
   USE wake_ini_mod , ONLY : crep_upper, crep_sol, tau_cv, rzero, aa0, flag_wk_check_trgl
   USE wake_ini_mod , ONLY : iflag_wk_act, iflag_wk_check_trgl, iflag_wk_pop_dyn, wdensmin
@@ -2313,5 +2313,5 @@
 SUBROUTINE pkupper (klon, klev, ptop, ph, pupper, kupper)
 
-USE wake_ini_mod , ONLY : pupperbyphs
+USE wake_ini_mod , ONLY : wk_pupper
 IMPLICIT NONE
 
@@ -2326,5 +2326,5 @@
  kupper = 0
  
-IF (pupperbyphs<1.) THEN
+IF (wk_pupper<1.) THEN
  ! Choose an integration bound well above wake top
   ! -----------------------------------------------------------------
@@ -2335,5 +2335,5 @@
   DO i = 1, klon
   !  pupper(i) = 0.6*ph(i, 1)
-    pupper(i) = pupperbyphs*ph(i, 1)
+    pupper(i) = wk_pupper*ph(i, 1)
     pupper(i) = max(pupper(i), 45000.)
     ! cc        Pupper(i) = 60000.
@@ -2343,6 +2343,6 @@
   
   DO i=1, klon
-     ! pupper(i) = pupperbyphs*ptop(i)+(1.-pupperbyphs)*ph(i, 1) 
-     pupper(i) = min( pupperbyphs*ptop(i)+(1.-pupperbyphs)*ph(i, 1) , ptop(i)-5000.) 
+     ! 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)-5000.) 
   END DO
 END IF
Index: LMDZ6/trunk/libf/phylmd/wake_ini_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/wake_ini_mod.F90	(revision 4452)
+++ LMDZ6/trunk/libf/phylmd/wake_ini_mod.F90	(revision 4453)
@@ -25,8 +25,8 @@
   INTEGER, SAVE, PROTECTED                                         :: prt_level
   REAL, SAVE, PROTECTED, DIMENSION(2)                              :: wdens_ref
-  REAL, SAVE, PROTECTED                                            :: stark, coefgw, alpk, pupperbyphs
+  REAL, SAVE, PROTECTED                                            :: stark, coefgw, alpk, wk_pupper
 !>jyg
   REAL, SAVE, PROTECTED                                            :: crep_upper, crep_sol  
-  !$OMP THREADPRIVATE(stark, wdens_ref, coefgw, alpk, pupperbyphs, crep_upper, crep_sol)
+  !$OMP THREADPRIVATE(stark, wdens_ref, coefgw, alpk, wk_pupper, crep_upper, crep_sol)
 
   REAL, SAVE, PROTECTED                                            :: tau_cv
@@ -153,6 +153,6 @@
   CALL getin_p('alpk',alpk)
   
-  pupperbyphs=0.6
-  CALL getin_p('pupperbyphs',pupperbyphs)
+  wk_pupper=0.6
+  CALL getin_p('wk_pupper',wk_pupper)
 
 
@@ -195,5 +195,5 @@
   WRITE(*,*) 'stark=', stark
   WRITE(*,*) 'alpk=', alpk
-  WRITE(*,*) 'pupperbyphs=', pupperbyphs
+  WRITE(*,*) 'wk_pupper=', wk_pupper
 !jyg<
 !!  WRITE(*,*) 'wdens_ref=', wdens_ref