Changeset 894 for trunk/LMDZ.VENUS/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.VENUS/libf/phyvenus

Files:

• profile.F (modified) (1 diff)
• rcm1d.F (modified) (4 diffs)
• testphys1d.F (deleted)

trunk/LMDZ.VENUS/libf/phyvenus/profile.F

@@ -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)

trunk/LMDZ.VENUS/libf/phyvenus/rcm1d.F

@@ -52 +52 @@
       REAL gru,grv   ! prescribed "geostrophic" background wind
       REAL temp(llm)   ! temperature at the middle of the layers
-      REAL q(llm,nqtot) ! tracer mixing ratio (e.g. kg/kg)
+      REAL,allocatable :: q(:,:) ! tracer mixing ratio (e.g. kg/kg)
       REAL tsoil(nsoilmx)   ! subsurface soik temperature (K)
       REAL zlay(llm)   ! altitude estimee dans les couches (km)
@@ -67 +67 @@
       REAL dudyn(llm),dvdyn(llm),dtempdyn(llm)
       REAL dpsurf    
-      REAL dq(llm,nqtot)
-      REAL dqdyn(llm,nqtot)
+      REAL,allocatable :: dq(:,:)
 
 c   Various intermediate variables
@@ -82 +81 @@
       COMMON/cpdetvenus/cppdyn,nu_venus,t0_venus
       REAL cppdyn,nu_venus,t0_venus
+      real pi
 
 c=======================================================================
@@ -113 +113 @@
       r= 8.314511E+0 *1000.E+0/mugaz
       rcp= r/cpp
+
+c-----------------------------------------------------------------------
+c   Initialisation des traceurs
+c   ---------------------------
+c  Choix du nombre de traceurs et du schema pour l'advection
+c  dans fichier traceur.def
+      call infotrac_init
+
+c Allocation de la tableau q : champs advectes   
+      allocate(q(llm,nqtot))
+      allocate(dq(llm,nqtot))
 
 c ------------------------------------------------------