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

Timestamp:

Dec 21, 2022, 10:55:05 AM (2 years ago)

Author:

streel

Message:

LMDZ.VENUS - chemical reaction rates update with JPL 2019

• Update of the main 1D subroutine with nitrogen species + generalisation of VMR output every 1/5 venusian day
• Added the H2 photolysis
• Added the traceur.def for all neutral species only
• Added the VMR def files for 1D model

NS

File:

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

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

@@ -89 +89 @@
       real :: dummy
 
-      character*8  specname(36)                  
-      real  mmol(36)
+      character*8  specname(40)                  
+      real  mmol(40)
 
 c=======================================================================
@@ -103 +103 @@
      &       99., 32., 48., 64., 80.,   
      &       96., 60., 81., 98., 64.,   
-     &       99., 83., 28., 4.,  18.,   
-     &       98./)
+     &       99., 83., 28.,  4., 46.,
+     &       30., 14., 14., 18., 98./)
 
       specname = (/"co2", "co", "h2", "h2o", "o1d",            
@@ -112 +112 @@
      &       "cocl2", "s", "so", "so2", "so3",           
      &       "s2o2", "ocs", "hso3", "h2so4", "s2",       
-     &       "clso2", "oscl", "n2", "he", "h2oliq",
-     &       "h2so4liq"/)
+     &       "clso2", "oscl", "n2", "he", "no2", "no",
+     &       "n", "n2d", "h2oliq", "h2so4liq"/)
 
 c ------------------------------------------------------
@@ -514 +514 @@
 c       ------------------------------------------------------------ 
 
-         if ((idt==1) .OR. (idt==9600) .OR. (idt==19200) .OR.
-     $      (idt==28800) .OR. (idt==38400) .OR. (idt==48000)) then
+         if (mod(idt,9600)==0) then   
               DO ilayer=1,nlayer 
                  write (5,'(100e12.4)')play(ilayer)/100., q(ilayer,:)
@@ -523 +522 @@
       ENDDO   ! fin de la boucle temporelle
 
-c      close(5)
+      close(5)
 c    ========================================================
 c    GESTION DES SORTIE