Changeset 2433 for trunk/LMDZ.MARS/libf/aeronomars/photolysis_mod.F90

Timestamp:

Nov 14, 2020, 4:51:23 PM (4 years ago)

Author:

flefevre

Message:

1) ajout de la photodissociation de HDO (commentee pour l'instant)

• necessite la presence des sections efficaces de HDO dans le datadir:

datadir/hdo_composite_295K.txt

2) CO + OH -> CO2 + H : cinetique JPL 2015 (attention ce changement affecte les resultats)
3) chimie heterogene : inactive par defaut.

File:

• trunk/LMDZ.MARS/libf/aeronomars/photolysis_mod.F90 (modified) (4 diffs)

trunk/LMDZ.MARS/libf/aeronomars/photolysis_mod.F90

@@ -6 +6 @@
 
   integer, parameter :: nphot = 13        ! number of photolysis
+! integer, parameter :: nphot = 14        ! number of photolysis (with hdo)
   integer, parameter :: nabs  = 10        ! number of absorbing gases
 
@@ -37 +38 @@
   real, dimension(nw), save :: xsn2                                   ! n2 absorption cross-section (cm2)
   real, dimension(nw), save :: yieldn2                                ! n2 photodissociation yield
+  real, dimension(nw), save :: xshdo                                  ! hdo absorption cross-section (cm2)
   real, dimension(nw), save :: albedo                                 ! surface albedo
 
@@ -97 +99 @@
 
   call rdxsn2(nw,wl,xsn2,yieldn2)
+
+! read and grid hdo cross-sections
+
+  call rdxshdo(nw,wl,xshdo)
 
 ! set surface albedo
@@ -1247 +1253 @@
 !==============================================================================
 
+      subroutine rdxshdo(nw, wl, yg)
+
+!-----------------------------------------------------------------------------*
+!=  PURPOSE:                                                                 =*
+!=  Read HDO molecular absorption cross section.  Re-grid data to match      =*
+!=  specified wavelength working grid.                                       =*
+!-----------------------------------------------------------------------------*
+!=  PARAMETERS:                                                              =*
+!=  NW     - INTEGER, number of specified intervals + 1 in working        (I)=*
+!=           wavelength grid                                                 =*
+!=  WL     - REAL, vector of lower limits of wavelength intervals in      (I)=*
+!=           working wavelength grid                                         =*
+!=  YG     - REAL, molecular absoprtion cross section (cm^2) of HDO at    (O)=*
+!=           each specified wavelength                                       =*
+!-----------------------------------------------------------------------------*
+
+      use datafile_mod, only: datadir
+
+      IMPLICIT NONE
+
+!     input
+
+      integer :: nw               ! number of wavelength grid points
+      real, dimension(nw) :: wl   ! lower and central wavelength for each interval
+
+!     output
+
+      real, dimension(nw) :: yg   ! hdo cross-sections (cm2)
+
+!     local
+
+      integer, parameter :: kdata = 900
+      real, parameter :: deltax = 1.e-4
+      REAL x1(kdata)
+      REAL y1(kdata)
+      INTEGER ierr
+      INTEGER i, n
+      CHARACTER*100 fil
+      integer :: kin, kout ! input/output logical units
+
+      kin = 10
+
+      fil = trim(datadir)//'/cross_sections/hdo_composite_295K.txt'
+      print*, 'section efficace HDO: ', fil
+
+      OPEN(UNIT=kin,FILE=fil,STATUS='old') 
+
+      DO i = 1,17 
+         read(kin,*)
+      END DO
+
+      n = 806 
+      DO i = 1, n
+         READ(kin,*) x1(i), y1(i)
+      END DO
+      CLOSE (kin)
+
+      CALL addpnt(x1,y1,kdata,n,x1(1)*(1.-deltax),0.)
+      CALL addpnt(x1,y1,kdata,n,          0.,0.)
+      CALL addpnt(x1,y1,kdata,n,x1(n)*(1.+deltax),0.)
+      CALL addpnt(x1,y1,kdata,n,      1.e+38,0.)
+      CALL inter2(nw,wl,yg,n,x1,y1,ierr)
+      IF (ierr .NE. 0) THEN
+         WRITE(*,*) ierr, fil
+         STOP
+      ENDIF
+      
+      end subroutine rdxshdo
+
+!==============================================================================
+
       subroutine rdxsh2o2(nw, wl, xsh2o2)