Changeset 894 for trunk/LMDZ.TITAN/libf

Timestamp:

Mar 6, 2013, 1:50:54 PM (12 years ago)

Author:

slebonnois

Message:

SL: small modifications on rcm1d (Venus, Titan)

Location:

trunk/LMDZ.TITAN/libf/phytitan

Files:

• profile.F (modified) (5 diffs)
• rcm1d.F (modified) (3 diffs)

trunk/LMDZ.TITAN/libf/phytitan/profile.F

@@ -1 @@
-      SUBROUTINE profile(unit,nlev,dzst,temp)
+      SUBROUTINE profile(unit,nlev,dzst,pres,temp)
       IMPLICIT NONE
 c=======================================================================
@@ -10 +10 @@
 c   differents profils d'atmospheres. T=f(z)
 c   entree:
-c     unit    unite de lecture de "testphys1d.def"
+c     unit    unite de lecture de "rcm1d.def"
 c     nlev    nombre de niveaux (nlev=llm+1, surf + couches 1 a llm)
 c     dzst    dz/T (avec dz = epaisseur de la couche en m)
+c     pres    pressure profile
 c     ichoice choix de l'atmosphere:
 c             1 Temperature constante
@@ -40 +41 @@
 
        INTEGER nlev, unit
-       REAL dzst(nlev),temp(nlev)
+       REAL dzst(nlev),pres(nlev),temp(nlev)
 
 c   local:
@@ -51 +52 @@
       REAL hauteur,tmp
       REAL zkm(nlev)    ! altitude en km
+      real a1,b1,c1,a2,b2,c2
 
       isin = 0
@@ -86 +88 @@
 
       ELSE IF(ichoice.EQ.2) THEN
-       print*,"Profil T a faire..."
-       stop
-c        temp(1) = 735.
-c        zkm(1)  = 0.0
-c        DO il=2,nlev
-c           zkm(il) = zkm(il-1)+temp(il-1)*dzst(il)/1000. ! approx avec T(l-1)
-c           if(zkm(il).lt.60.) then
-c             temp(il)=735.-7.95*zkm(il)
-c           else
-c             temp(il)=AMAX1(258.-3.*(zkm(il)-60.),168.)
-c           endif
-c           zkm(il) = zkm(il-1)+(temp(il-1)+temp(il))/2.*dzst(il)/1000.
-c        ENDDO
+       a1 =       142.1 
+       b1 =      -21.45 
+       c1 =       40.11 
+       a2 =       106.3 
+       b2 =       3183. 
+       c2 =       4737. 
+       DO il=1,nlev
+         temp(il)=a1*exp(-((pres(il)-b1)/c1)^2) 
+     .          + a2*exp(-((pres(il)-b2)/c2)^2)
+       ENDDO
+       zkm(1)  = 0.0
+       DO il=2,nlev
+          zkm(il) = zkm(il-1)+(temp(il-1)+temp(il))/2.*dzst(il)/1000.
+       ENDDO
 
 c-----------------------------------------------------------------------

trunk/LMDZ.TITAN/libf/phytitan/rcm1d.F

@@ -78 +78 @@
       REAL pk(llm),pks, w(llm)
       INTEGER l, ierr, aslun
-      REAL tmp1(0:llm),tmp2(0:llm)                        
+      REAL tmp1(0:llm),tmp2(0:llm),tmp3(0:llm)                        
 
       character*2 str2
@@ -85 +85 @@
       COMMON/cpdetvenus/cppdyn,nu_venus,t0_venus
       REAL cppdyn,nu_venus,t0_venus
+      real pi
 
 c=======================================================================
@@ -279 +280 @@
         tmp1(ilayer)=-log(play(ilayer)/play(ilayer-1))*r/g
       ENDDO
-      call profile(unit,nlayer+1,tmp1,tmp2)
+      DO ilayer=0,nlayer
+        tmp2(ilayer)=plev(ilayer+1)
+      ENDDO
+      call profile(unit,nlayer+1,tmp1,tmp2,tmp3)
       CLOSE(unit)
 
       print*,"               Pression        Altitude     Temperature"
       ilayer=1
-      tsurf(1)=tmp2(0)
-       temp(1)=tmp2(1)
-       zlay(1)=tmp2(1)*tmp1(1)
+      tsurf(1)=tmp3(0)
+       temp(1)=tmp3(1)
+       zlay(1)=tmp3(1)*tmp1(1)
       print*,"           0",tsurf(1)
       print*,ilayer,play(ilayer),zlay(ilayer),temp(ilayer)
       DO ilayer=2,nlayer
-        temp(ilayer)=tmp2(ilayer)
-        zlay(ilayer)=zlay(ilayer-1)+tmp2(ilayer)*tmp1(ilayer)
+        temp(ilayer)=tmp3(ilayer)
+        zlay(ilayer)=zlay(ilayer-1)+tmp3(ilayer)*tmp1(ilayer)
         print*,ilayer,play(ilayer),zlay(ilayer),temp(ilayer)
       ENDDO