Changeset 5 for trunk/libf/dyn3d/vlspltqs.F

Timestamp:

Oct 14, 2010, 3:18:14 PM (14 years ago)

Author:

slebonnois

Message:

• Ajout des fichiers .def Venus et Titan (tels qu'ils sont utilisés

actuellement) dans les deftanks.

• Ajout d'une doc sur Cp(T).
• Modifications dans dyn3d concernant Cp(T), cf le log (v5) dans chantiers
• Premières modifs de l'appel à la physique dans dyn3d/calfis, cf log (v5)
• Elimination des cpdet.* dans phytitan et phyvenus (remplacés par cpdet.F

dans dyn3d).

File:

• trunk/libf/dyn3d/vlspltqs.F (modified) (2 diffs)

trunk/libf/dyn3d/vlspltqs.F

@@ -65 +65 @@
       REAL rtt,retv,r2es,r3les,r3ies,r4les,r4ies,play
       REAL ptarg,pdelarg,foeew,zdelta
-      REAL tempe(ip1jmp1)
+! ADAPTATION GCM POUR CP(T)
+      REAL tempe(ip1jmp1,llm)
 
 c    fonction psat(T)
@@ -84 +85 @@
 c-- approximation: au milieu des couches play(l)=(p(l)+p(l+1))/2
 c   pour eviter une exponentielle.
+
+! ADAPTATION GCM POUR CP(T)
+        call tpot2t(ip1jmp1*llm,teta,tempe,pk)
         DO l = 1, llm
          DO ij = 1, ip1jmp1
-          tempe(ij) = teta(ij,l) * pk(ij,l) /cpp
-         ENDDO
-         DO ij = 1, ip1jmp1
-          zdelta = MAX( 0., SIGN(1., rtt - tempe(ij)) )
+          zdelta = MAX( 0., SIGN(1., rtt - tempe(ij,l)) )
           play   = 0.5*(p(ij,l)+p(ij,l+1))
-          qsat(ij,l) = MIN(0.5, r2es* FOEEW(tempe(ij),zdelta) / play )
+          qsat(ij,l) = MIN(0.5, r2es* FOEEW(tempe(ij,l),zdelta) / play )
           qsat(ij,l) = qsat(ij,l) / ( 1. - retv * qsat(ij,l) )
          ENDDO