Changeset 3519 for trunk/LMDZ.VENUS/libf

Timestamp:

Nov 15, 2024, 3:39:29 PM (5 weeks ago)

Author:

slebonnois

Message:

SL: some cleaning of rcm1d Venus

File:

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

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

@@ -96 +96 @@
       real :: dummy
 
-      character*8  specname(40)                  
-      real  mmol(40)
-
 c=======================================================================
 c INITIALISATION
@@ -108 +105 @@
       call init_parallel
 #endif
-
-      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., 46.,
-     &       30., 14., 14., 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", "no2", "no",
-     &       "n", "n2d", "h2oliq", "h2so4liq"/)
 
 c ------------------------------------------------------
@@ -155 +134 @@
 c  dans fichier traceur.def
       call infotrac_init
+      iflag_trac=0
       if (nqtot.gt.1) iflag_trac=1
 
@@ -360 +340 @@
       ENDDO
 
+      if (iflag_trac.eq.1) then
       print*,"rcm1d: Loading chemistry profiles from init_1D.txt"
       ! check if the file is indeed there
@@ -376 +357 @@
         stop
       endif ! of if(file_is_present)
+      endif ! iflag_trac
 
 c    Initialisation des parametres d'oro
@@ -422 +404 @@
 c=======================================================================
 c
-!TEMPOAIRE
+!TEMPORAIRE
 
       firstcall=.true.
       lastcall=.false.
-
-c     Ouverture du fichier d'écriture des VMR
-
-      OPEN(15,file='chem.txt',form = 'formatted')
-      WRITE(15,'(2x,100a12)')'hpa',(specname(i), i=1,nqtot)
 
 !     debut de boucle temporelle
@@ -471 +448 @@
      ,     presnivs,
      ,     u,v,temp,q,
-c     ,     plev,temp, ! "planetary mean" plev and temperature
      ,     w,
 C - sorties
@@ -535 +511 @@
              play(ilayer)=aps(ilayer)+psurf*bps(ilayer)
            ENDDO 
-
-c       ------------------------------------------------------------ 
-c       sortie des VMR tous les 20% d'une journée vénusienne dans le
-c       fichier chem.txt
-c       ------------------------------------------------------------ 
-
-         if (mod(idt,9600)==0) then   
-              DO ilayer=1,nlayer 
-                 write (15,'(100e12.4)')play(ilayer)/100., q(ilayer,:)
-     $                 *mmean(1,ilayer)/mmol(:)
-              ENDDO   
-         end if
+        
       ENDDO   ! fin de la boucle temporelle