Index: trunk/libf/phyvenus/cpdet.F
===================================================================
--- trunk/libf/phyvenus/cpdet.F	(revision 3)
+++ 	(revision )
@@ -1,117 +1,0 @@
-! ADAPTATION GCM POUR CP(T)
-c======================================================================
-c S. Lebonnois, 10/2007:
-c
-c VENUS: Cp(T) = cpp*(T/T0)^nu 
-c avec T0=460. et nu=0.35
-c cpp=RCPD=cp0 = 1000.
-c R/RCPD = RKAPPA
-c
-c La fonction d'Exner reste pk = RCPD*(play/pref)**RKAPPA
-c 
-c T et teta (temperature potentielle) sont liees par:
-c 
-c   integrale[teta a T](cp/T dT) = integrale[pref a p](R/p dp)
-c
-c Dans le cas de l'expression pour Venus, ca donne:
-c
-c   teta**nu = T**nu - nu * T0**nu * ln[ (p/pref)**RKAPPA ]
-c ou
-c   teta**nu = T**nu - nu * T0**nu * ln[pk/RCPD]
-c
-c On passe de T a teta par t2tpot(t,teta,pk)
-c On passe de teta a T par tpot2t(teta,t,pk)
-c
-c Pour DT <-> Dteta, on utilise: dteta = dT *(T/teta)**(nu-1)
-c -> routine dt2dtpot(dt,dteta,t,teta) 
-c (utilisee seulement pour le contregradient)
-c
-c======================================================================
-
-      FUNCTION cpdet(t)
-      IMPLICIT none
-#include "cpdet.h"
-
-      real cpdet,t
-
-      cpdet = cp0*(t/t0)**nu
-
-      return
-      end
-      
-c======================================================================
-c======================================================================
-
-      SUBROUTINE t2tpot(npoints,yt, yteta, ypk)
-      IMPLICIT none
-c======================================================================
-c Arguments:
-c
-c yt   --------input-R- Temperature
-c yteta-------output-R- Temperature potentielle
-c ypk  --------input-R- Fonction d'Exner: RCPD*(pplay/pref)**RKAPPA
-c
-c======================================================================
-#include "cpdet.h"
-
-      integer npoints
-      REAL    yt(npoints), yteta(npoints), ypk(npoints)
-      
-      yteta = yt**nu - nu * t0**nu * log(ypk/cp0)
-      yteta = yteta**(1./nu)
-        
-      return
-      end
-
-c======================================================================
-c======================================================================
-
-      SUBROUTINE tpot2t(npoints,yteta, yt, ypk)
-      IMPLICIT none
-c======================================================================
-c Arguments:
-c
-c yteta--------input-R- Temperature potentielle
-c yt   -------output-R- Temperature
-c ypk  --------input-R- Fonction d'Exner: RCPD*(pplay/pref)**RKAPPA
-c
-c======================================================================
-#include "cpdet.h"
-
-      integer npoints
-      REAL yt(npoints), yteta(npoints), ypk(npoints)
-
-      yt = yteta**nu + nu * t0**nu * log(ypk/cp0)
-      yt = yt**(1./nu)
-      
-      return
-      end
-
-c======================================================================
-c======================================================================
-
-      SUBROUTINE dt2dtpot(npoints,ydt, ydteta, yt, yteta, ypk, ydpk)
-      IMPLICIT none
-c======================================================================
-c Arguments:
-c
-c ydt  --------input-R- D Temperature
-c ydteta------output-R- D Temperature potentielle
-c yt   --------input-R- Temperature
-c yteta--------input-R- Temperature potentielle
-c ypk  --------input-R- Fonction d'Exner: RCPD*(pplay/pref)**RKAPPA
-c ydpk  -------input-R- D Fonction d'Exner
-c
-c======================================================================
-#include "cpdet.h"
-
-      integer npoints
-      REAL    ydt(npoints), ydteta(npoints), ydpk(npoints)
-      REAL    yt(npoints), yteta(npoints), ypk(npoints)
-      
-      ydteta =  yteta**(1.-nu)* (ydt * yt**(nu-1.) 
-     .                           - t0**nu * ydpk / ypk)
-        
-      return
-      end
-
Index: trunk/libf/phyvenus/cpdet.h
===================================================================
--- trunk/libf/phyvenus/cpdet.h	(revision 3)
+++ 	(revision )
@@ -1,8 +1,0 @@
-c-----------------------------------------------------------------------
-c INCLUDE cpdet.h
-
-      real cp0,t0,nu
-      parameter (cp0 = 1000.) !doit etre egal a cpp (dyn) et RCPD (phy)
-      parameter (t0  = 460.)
-      parameter (nu  = 0.35)
-c-----------------------------------------------------------------------
Index: trunk/libf/phyvenus/cpdet.h.terre
===================================================================
--- trunk/libf/phyvenus/cpdet.h.terre	(revision 3)
+++ 	(revision )
@@ -1,6 +1,0 @@
-c-----------------------------------------------------------------------
-c INCLUDE cpdet.h
-
-      real cp0
-      parameter (cp0 = 900.) !doit etre egal a cpp (dyn) et RCPD (phy)
-c-----------------------------------------------------------------------
Index: trunk/libf/phyvenus/cpdet.h.venus
===================================================================
--- trunk/libf/phyvenus/cpdet.h.venus	(revision 3)
+++ 	(revision )
@@ -1,8 +1,0 @@
-c-----------------------------------------------------------------------
-c INCLUDE cpdet.h
-
-      real cp0,t0,nu
-      parameter (cp0 = 1000.) !doit etre egal a cpp (dyn) et RCPD (phy)
-      parameter (t0  = 460.)
-      parameter (nu  = 0.35)
-c-----------------------------------------------------------------------
Index: trunk/libf/phyvenus/cpdet.terre
===================================================================
--- trunk/libf/phyvenus/cpdet.terre	(revision 3)
+++ 	(revision )
@@ -1,114 +1,0 @@
-! ADAPTATION GCM POUR CP(T)
-c======================================================================
-c S. Lebonnois, 10/2007:
-c
-c TERRE: Cp(T) = cpp  
-c cpp=RCPD=cp0 = 1000.
-c R/RCPD = RKAPPA
-c
-c La fonction d'Exner reste pk = RCPD*(play/pref)**RKAPPA
-c 
-c T et teta (temperature potentielle) sont liees par:
-c 
-c   integrale[teta a T](cp/T dT) = integrale[pref a p](R/p dp)
-c
-c Dans le cas de la Terre, ca donne:
-c
-c   teta = T * (pref/p)**RKAPPA
-c ou
-c   teta = T * RCPD/pk
-c
-c On passe de T a teta par t2tpot(t,teta,pk)
-c On passe de teta a T par tpot2t(teta,t,pk)
-c
-c
-c Pour DT <-> Dteta, on utilise: dteta = dT * teta/T
-c -> routine dt2dtpot(dt,dteta,t,teta) 
-c (utilisee seulement pour le contregradient)
-c
-c======================================================================
-
-      FUNCTION cpdet(t)
-      IMPLICIT none
-#include "cpdet.h"
-
-      real cpdet,t
-
-      cpdet = cp0
-
-      return
-      end
-      
-c======================================================================
-c======================================================================
-
-      SUBROUTINE t2tpot(npoints, yt, yteta, ypk)
-      IMPLICIT none
-c======================================================================
-c Arguments:
-c
-c yt   --------input-R- Temperature
-c yteta-------output-R- Temperature potentielle
-c ypk  --------input-R- Fonction d'Exner: RCPD*(pplay/pref)**RKAPPA
-c
-c======================================================================
-#include "cpdet.h"
-
-      integer npoints
-      REAL yt(npoints), yteta(npoints), ypk(npoints)
-
-      yteta = yt * cp0/ypk
-      
-      return
-      end
-
-c======================================================================
-c======================================================================
-
-      SUBROUTINE tpot2t(npoints, yteta, yt, ypk)
-      IMPLICIT none
-c======================================================================
-c Arguments:
-c
-c yteta--------input-R- Temperature potentielle
-c yt   -------output-R- Temperature
-c ypk  --------input-R- Fonction d'Exner: RCPD*(pplay/pref)**RKAPPA
-c
-c======================================================================
-#include "cpdet.h"
-
-      integer npoints
-      REAL yt(npoints), yteta(npoints), ypk(npoints)
-
-      yt = yteta * ypk/cp0
-      
-      return
-      end
-
-c======================================================================
-c======================================================================
-
-      SUBROUTINE dt2dtpot(npoints, ydt, ydteta, yt, yteta, ypk, ydpk)
-      IMPLICIT none
-c======================================================================
-c Arguments:
-c
-c ydt  --------input-R- D Temperature
-c ydteta------output-R- D Temperature potentielle
-c yt   --------input-R- Temperature
-c yteta--------input-R- Temperature potentielle
-c ypk  --------input-R- Fonction d'Exner: RCPD*(pplay/pref)**RKAPPA
-c ydpk  -------input-R- D Fonction d'Exner
-c
-c======================================================================
-#include "cpdet.h"
-
-      integer npoints
-      REAL    ydt(npoints), ydteta(npoints), ydpk(npoints)
-      REAL    yt(npoints), yteta(npoints), ypk(npoints)
-      
-      ydteta = yteta* (ydt/yt - ydpk/ypk)
-        
-      return
-      end
-
Index: trunk/libf/phyvenus/cpdet.terre.old
===================================================================
--- trunk/libf/phyvenus/cpdet.terre.old	(revision 3)
+++ 	(revision )
@@ -1,115 +1,0 @@
-! ADAPTATION GCM POUR CP(T)
-c======================================================================
-c S. Lebonnois, 10/2007:
-c
-c TERRE: Cp(T) = cpp  
-c cpp=RCPD=cp0 = 1000.
-c R/RCPD = RKAPPA
-c
-c La fonction d'Exner reste pk = RCPD*(play/pref)**RKAPPA
-c 
-c T et teta (temperature potentielle) sont liees par:
-c 
-c   integrale[teta a T](cp/T dT) = integrale[pref a p](R/p dp)
-c
-c Dans le cas de la Terre, ca donne:
-c
-c   teta = T * (pref/p)**RKAPPA
-c ou
-c   teta = T * RCPD/pk
-c
-c On passe de T a teta par t2tpot(t,teta,pk)
-c On passe de teta a T par tpot2t(teta,t,pk)
-c
-c
-c Pour DT <-> Dteta, on utilise: dteta = dT * teta/T
-c -> routine dt2dtpot(dt,dteta,t,teta) 
-c (utilisee seulement pour le contregradient)
-c
-c======================================================================
-
-      FUNCTION cpdet(t)
-      IMPLICIT none
-#include "cpdet.h"
-
-      real cpdet,t
-
-      cpdet = cp0
-
-      return
-      end
-      
-c======================================================================
-c======================================================================
-
-      SUBROUTINE t2tpot(npoints, yt, yteta, ypk)
-      IMPLICIT none
-c======================================================================
-c Arguments:
-c
-c yt   --------input-R- Temperature
-c yteta-------output-R- Temperature potentielle
-c ypk  --------input-R- Fonction d'Exner: RCPD*(pplay/pref)**RKAPPA
-c
-c======================================================================
-#include "dimensions.h"
-#include "cpdet.h"
-
-      integer npoints
-      REAL yt(npoints), yteta(npoints), ypk(npoints)
-
-      yteta = yt * cp0/ypk
-      
-      return
-      end
-
-c======================================================================
-c======================================================================
-
-      SUBROUTINE tpot2t(npoints, yteta, yt, ypk)
-      IMPLICIT none
-c======================================================================
-c Arguments:
-c
-c yteta--------input-R- Temperature potentielle
-c yt   -------output-R- Temperature
-c ypk  --------input-R- Fonction d'Exner: RCPD*(pplay/pref)**RKAPPA
-c
-c======================================================================
-#include "dimensions.h"
-#include "cpdet.h"
-
-      integer npoints
-      REAL yt(npoints), yteta(npoints), ypk(npoints)
-
-      yt = yteta * ypk/cp0
-      
-      return
-      end
-
-c======================================================================
-c======================================================================
-
-      SUBROUTINE dt2dtpot(npoints, ydt, ydteta, yt, yteta)
-      IMPLICIT none
-c======================================================================
-c Arguments:
-c
-c ydt  --------input-R- D Temperature
-c ydteta------output-R- D Temperature potentielle
-c yt   --------input-R- Temperature
-c yteta--------input-R- Temperature potentielle
-c
-c======================================================================
-#include "dimensions.h"
-#include "cpdet.h"
-
-      integer npoints
-      REAL    ydt(npoints), ydteta(npoints)
-      REAL    yt(npoints), yteta(npoints)
-      
-      ydteta = ydt * (yteta/yt)
-        
-      return
-      end
-
Index: trunk/libf/phyvenus/cpdet.venus
===================================================================
--- trunk/libf/phyvenus/cpdet.venus	(revision 3)
+++ 	(revision )
@@ -1,117 +1,0 @@
-! ADAPTATION GCM POUR CP(T)
-c======================================================================
-c S. Lebonnois, 10/2007:
-c
-c VENUS: Cp(T) = cpp*(T/T0)^nu 
-c avec T0=460. et nu=0.35
-c cpp=RCPD=cp0 = 1000.
-c R/RCPD = RKAPPA
-c
-c La fonction d'Exner reste pk = RCPD*(play/pref)**RKAPPA
-c 
-c T et teta (temperature potentielle) sont liees par:
-c 
-c   integrale[teta a T](cp/T dT) = integrale[pref a p](R/p dp)
-c
-c Dans le cas de l'expression pour Venus, ca donne:
-c
-c   teta**nu = T**nu - nu * T0**nu * ln[ (p/pref)**RKAPPA ]
-c ou
-c   teta**nu = T**nu - nu * T0**nu * ln[pk/RCPD]
-c
-c On passe de T a teta par t2tpot(t,teta,pk)
-c On passe de teta a T par tpot2t(teta,t,pk)
-c
-c Pour DT <-> Dteta, on utilise: dteta = dT *(T/teta)**(nu-1)
-c -> routine dt2dtpot(dt,dteta,t,teta) 
-c (utilisee seulement pour le contregradient)
-c
-c======================================================================
-
-      FUNCTION cpdet(t)
-      IMPLICIT none
-#include "cpdet.h"
-
-      real cpdet,t
-
-      cpdet = cp0*(t/t0)**nu
-
-      return
-      end
-      
-c======================================================================
-c======================================================================
-
-      SUBROUTINE t2tpot(npoints,yt, yteta, ypk)
-      IMPLICIT none
-c======================================================================
-c Arguments:
-c
-c yt   --------input-R- Temperature
-c yteta-------output-R- Temperature potentielle
-c ypk  --------input-R- Fonction d'Exner: RCPD*(pplay/pref)**RKAPPA
-c
-c======================================================================
-#include "cpdet.h"
-
-      integer npoints
-      REAL    yt(npoints), yteta(npoints), ypk(npoints)
-      
-      yteta = yt**nu - nu * t0**nu * log(ypk/cp0)
-      yteta = yteta**(1./nu)
-        
-      return
-      end
-
-c======================================================================
-c======================================================================
-
-      SUBROUTINE tpot2t(npoints,yteta, yt, ypk)
-      IMPLICIT none
-c======================================================================
-c Arguments:
-c
-c yteta--------input-R- Temperature potentielle
-c yt   -------output-R- Temperature
-c ypk  --------input-R- Fonction d'Exner: RCPD*(pplay/pref)**RKAPPA
-c
-c======================================================================
-#include "cpdet.h"
-
-      integer npoints
-      REAL yt(npoints), yteta(npoints), ypk(npoints)
-
-      yt = yteta**nu + nu * t0**nu * log(ypk/cp0)
-      yt = yt**(1./nu)
-      
-      return
-      end
-
-c======================================================================
-c======================================================================
-
-      SUBROUTINE dt2dtpot(npoints,ydt, ydteta, yt, yteta, ypk, ydpk)
-      IMPLICIT none
-c======================================================================
-c Arguments:
-c
-c ydt  --------input-R- D Temperature
-c ydteta------output-R- D Temperature potentielle
-c yt   --------input-R- Temperature
-c yteta--------input-R- Temperature potentielle
-c ypk  --------input-R- Fonction d'Exner: RCPD*(pplay/pref)**RKAPPA
-c ydpk  -------input-R- D Fonction d'Exner
-c
-c======================================================================
-#include "cpdet.h"
-
-      integer npoints
-      REAL    ydt(npoints), ydteta(npoints), ydpk(npoints)
-      REAL    yt(npoints), yteta(npoints), ypk(npoints)
-      
-      ydteta =  yteta**(1.-nu)* (ydt * yt**(nu-1.) 
-     .                           - t0**nu * ydpk / ypk)
-        
-      return
-      end
-
Index: trunk/libf/phyvenus/cpdet.venus.old
===================================================================
--- trunk/libf/phyvenus/cpdet.venus.old	(revision 3)
+++ 	(revision )
@@ -1,117 +1,0 @@
-! ADAPTATION GCM POUR CP(T)
-c======================================================================
-c S. Lebonnois, 10/2007:
-c
-c VENUS: Cp(T) = cpp*(T/T0)^nu 
-c avec T0=460. et nu=0.35
-c cpp=RCPD=cp0 = 1000.
-c R/RCPD = RKAPPA
-c
-c La fonction d'Exner reste pk = RCPD*(play/pref)**RKAPPA
-c 
-c T et teta (temperature potentielle) sont liees par:
-c 
-c   integrale[teta a T](cp/T dT) = integrale[pref a p](R/p dp)
-c
-c Dans le cas de l'expression pour Venus, ca donne:
-c
-c   teta**nu = T**nu - nu * T0**nu * ln[ (p/pref)**RKAPPA ]
-c ou
-c   teta**nu = T**nu - nu * T0**nu * ln[pk/RCPD]
-c
-c On passe de T a teta par t2tpot(t,teta,pk)
-c On passe de teta a T par tpot2t(teta,t,pk)
-c
-c Pour DT <-> Dteta, on utilise: dteta = dT *(T/teta)**(nu-1)
-c -> routine dt2dtpot(dt,dteta,t,teta) 
-c (utilisee seulement pour le contregradient)
-c
-c======================================================================
-
-      FUNCTION cpdet(t)
-      IMPLICIT none
-#include "cpdet.h"
-
-      real cpdet,t
-
-      cpdet = cp0*(t/t0)**nu
-
-      return
-      end
-      
-c======================================================================
-c======================================================================
-
-      SUBROUTINE t2tpot(npoints,yt, yteta, ypk)
-      IMPLICIT none
-c======================================================================
-c Arguments:
-c
-c yt   --------input-R- Temperature
-c yteta-------output-R- Temperature potentielle
-c ypk  --------input-R- Fonction d'Exner: RCPD*(pplay/pref)**RKAPPA
-c
-c======================================================================
-#include "dimensions.h"
-#include "cpdet.h"
-
-      integer npoints
-      REAL    yt(npoints), yteta(npoints), ypk(npoints)
-      
-      yteta = yt**nu - nu * t0**nu * log(ypk/cp0)
-      yteta = yteta**(1./nu)
-        
-      return
-      end
-
-c======================================================================
-c======================================================================
-
-      SUBROUTINE tpot2t(npoints,yteta, yt, ypk)
-      IMPLICIT none
-c======================================================================
-c Arguments:
-c
-c yteta--------input-R- Temperature potentielle
-c yt   -------output-R- Temperature
-c ypk  --------input-R- Fonction d'Exner: RCPD*(pplay/pref)**RKAPPA
-c
-c======================================================================
-#include "dimensions.h"
-#include "cpdet.h"
-
-      integer npoints
-      REAL yt(npoints), yteta(npoints), ypk(npoints)
-
-      yt = yteta**nu + nu * t0**nu * log(ypk/cp0)
-      yt = yt**(1./nu)
-      
-      return
-      end
-
-c======================================================================
-c======================================================================
-
-      SUBROUTINE dt2dtpot(npoints,ydt, ydteta, yt, yteta)
-      IMPLICIT none
-c======================================================================
-c Arguments:
-c
-c ydt  --------input-R- D Temperature
-c ydteta------output-R- D Temperature potentielle
-c yt   --------input-R- Temperature
-c yteta--------input-R- Temperature potentielle
-c
-c======================================================================
-#include "dimensions.h"
-#include "cpdet.h"
-
-      integer npoints
-      REAL    ydt(npoints), ydteta(npoints)
-      REAL    yt(npoints), yteta(npoints)
-      
-      ydteta = ydt * (yt/yteta)**(1.-nu)
-        
-      return
-      end
-
