Changeset 2684 for trunk/LMDZ.VENUS/libf/phyvenus/dyn1d

Timestamp:

May 23, 2022, 6:29:21 PM (3 years ago)

Author:

flefevre

Message:

Venus: Modifications pour tourner en 1D avec la chimie. Provisoire (cles a ajouter dans le futur pour gerer les cas avec/sans chimie) mais fonctionne.

File:

• trunk/LMDZ.VENUS/libf/phyvenus/dyn1d/rcm1d.F (modified) (8 diffs)

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

@@ -11 +11 @@
      .                       aps,bps,scaleheight,pseudoalt,
      .                       disvert_type,pressure_exner
-      use conc, only: rho
+      use conc, only: rho,mmean
       USE iniphysiq_mod, ONLY: iniphysiq
       USE mod_const_mpi, ONLY: comm_lmdz
@@ -49 +49 @@
       INTEGER unitstart      ! unite d'ecriture de "startphy.nc"
       INTEGER nlayer,nlevel,nsoil,ndt
-      INTEGER ilayer,ilevel,isoil,idt,iq
+      INTEGER ilayer,ilevel,isoil,idt,iq,i
       LOGICAl firstcall,lastcall
 c
@@ -84 +84 @@
       real pi
 
-c=======================================================================
+!     initialisation des traceurs
+
+      integer :: idummy
+      real :: dummy
+
+      character*8  specname(36)                  
+      real  mmol(36)
 
 c=======================================================================
 c INITIALISATION
-c=======================================================================
 
       lunout = 6
+
+      mmol =     (/44., 28.,  2., 18., 16.,   
+     &       16., 32., 32., 48.,  1.,   
+     &       17., 33., 34., 35., 51.,   
+     &       71., 36., 52., 63., 95.,   
+     &       99., 32., 48., 64., 80.,   
+     &       96., 60., 81., 98., 64.,   
+     &       99., 83., 28., 4.,  18.,   
+     &       98./)
+
+      specname = (/"co2", "co", "h2", "h2o", "o1d",            
+     &       "o"  , "o2", "o2dg", "o3", "h",             
+     &       "oh" , "ho2", "h2o2", "cl", "clo",          
+     &       "cl2", "hcl", "hocl", "clco", "clco3",      
+     &       "cocl2", "s", "so", "so2", "so3",           
+     &       "s2o2", "ocs", "hso3", "h2so4", "s2",       
+     &       "clso2", "oscl", "n2", "he", "h2oliq",
+     &       "h2so4liq"/)
 
 c ------------------------------------------------------
@@ -169 +192 @@
 c     print*,nbapp_rad
 c LU DANS PHYSIQ.DEF...
-      nbapp_rad = 1000.
+      nbapp_rad = 24000
 
       PRINT *,'nombre de jours simules ?'
@@ -223 +246 @@
      &            (/cufi,0.,0.,0./),
      &            (/cvfi,0./),
-     &            rad,g,r,cpp,1)
+     &             rad,g,r,cpp,1)
 
 c   le geopotentiel au sol est inutile en 1D car tout est controle
@@ -323 +346 @@
       ENDDO
 
-c FULL CHEMISTRY !! AJOUTER INIT AURELIEN...
-C Faudrait lire les cles avant pour mettre ca en option....
-c ou alors mettre ca dans physiq
+      print*,"lecture des profils chimiques"
+      open(21, form = 'formatted', file = 'init_1D.txt')
+      read(21,*)
+      do ilayer = nlayer,1,-1
+         read(21,*) idummy, dummy, dummy, (q(ilayer,iq), iq = 1,nqtot)
+!        print*, idummy, q(ilayer,1), q(ilayer,nqtot)
+      end do
+      close(21)
 
 c    Initialisation des parametres d'oro
@@ -367 +395 @@
 c=======================================================================
 c
+!TEMPOAIRE
+
       firstcall=.true.
       lastcall=.false.
 
-      DO idt=1,ndt
+c     Ouverture du fichier d'écriture des VMR
+
+      OPEN(5,file='chem.txt',form = 'formatted')
+      WRITE(5,'(2x,100a12)')'hpa',(specname(i), i=1,nqtot)
+
+!     debut de boucle temporelle
+
+      DO idt=1,48000
         IF (idt.eq.ndt) then
          lastcall=.true.
@@ -470 +507 @@
            DO ilayer=1,nlayer
              play(ilayer)=aps(ilayer)+psurf*bps(ilayer)
-           ENDDO
-
+           ENDDO 
+
+c       ------------------------------------------------------------ 
+c       sortie des VMR tous les 20% d'une journée vénusienne dans le
+c       fichier chem.txt
+c       ------------------------------------------------------------ 
+
+         if ((idt==1) .OR. (idt==9600) .OR. (idt==19200) .OR.
+     $      (idt==28800) .OR. (idt==38400) .OR. (idt==48000)) then
+              DO ilayer=1,nlayer 
+                 write (5,'(100e12.4)')play(ilayer)/100., q(ilayer,:)
+     $                 *mmean(1,ilayer)/mmol(:)
+              ENDDO   
+         end if
       ENDDO   ! fin de la boucle temporelle
 
+c      close(5)
 c    ========================================================
 c    GESTION DES SORTIE