Changeset 2044 for trunk/LMDZ.MARS/libf/aeronomars

Timestamp:

Nov 18, 2018, 8:09:06 PM (6 years ago)

Author:

flefevre

Message:

Photolyse on-line : desormais par defaut

Location:

trunk/LMDZ.MARS/libf/aeronomars

Files:

• calchim.F90 (modified) (1 diff)
• photochemistry.F90 (modified) (3 diffs)
• photolysis_mod.F90 (modified) (3 diffs)

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

@@ -193 +193 @@
 
          if (photochem) then
-            print*,'calchim: Read photolysis lookup table'
-            call read_phototable     ! off-line photolysis
-            print*,'calchim: Read UV absorption cross-sections'
-            call init_photolysis     ! on-line photolysis
+            if (jonline) then
+               print*,'calchim: Read UV absorption cross-sections'
+               call init_photolysis     ! on-line photolysis
+            else
+               print*,'calchim: Read photolysis lookup table'
+               call read_phototable     ! off-line photolysis
+            end if
          end if
          ! find index of chemical tracers to use

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

@@ -20 +20 @@
 use photolysis_mod, only : nb_phot_max,       &
                            nb_reaction_3_max, &
-                           nb_reaction_4_max
+                           nb_reaction_4_max, &
+                           jonline
 
 use param_v4_h, only: jdistot, jdistot_b
@@ -73 +74 @@
 integer :: lswitch
 logical, save :: firstcall = .true.
-logical :: jonline               ! switch for online photolysis
 logical :: jparam                ! switch for J parameterization
 
@@ -151 +151 @@
 !===================================================================
 
-jonline = .false.
 jparam= .false.
 

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

@@ -3 +3 @@
   implicit none
 
+! photolysis
+
+  logical, parameter :: jonline = .true.  ! true: on-line ! false: lookup table
   integer, parameter :: nphot = 13        ! number of photolysis
   integer, parameter :: nabs  = 10        ! number of absorbing gases
@@ -397 +400 @@
          fil = trim(datadir)//'/solar_fluxes/atlas3_thuillier_tuv.txt'
          print*, 'solar flux : ', fil
-         OPEN(UNIT=kin,FILE=fil,STATUS='old')
+         open(kin, file=fil, status='old', iostat=ierr)
+
+         if (ierr /= 0) THEN 
+            write(*,*)'cant find solar flux : ', fil
+            write(*,*)'It should be in :', trim(datadir),'/solar_fluxes'
+            write(*,*)'1) You can change this directory address in '
+            write(*,*)'   callphys.def with datadir=/path/to/dir'
+            write(*,*)'2) If necessary, /solar fluxes (and other datafiles)'
+            write(*,*)'   can be obtained online on:'
+            write(*,*)' http://www.lmd.jussieu.fr/~lmdz/planets/mars/datadir'
+            stop
+         end if
+
          nhead = 9
          n = 19193
@@ -928 +943 @@
       fil = trim(datadir)//'/cross_sections/o2_composite_2018_150K.txt'
       print*, 'section efficace O2 150K: ', fil
-      OPEN(UNIT=kin,FILE=fil,STATUS='old')
+      open(kin, file=fil, status='old', iostat=ierr)
+
+      if (ierr /= 0) THEN 
+         write(*,*)'cant find O2 cross-sections : ', fil
+         write(*,*)'It should be in :', trim(datadir),'/cross_sections'
+         write(*,*)'1) You can change this directory address in '
+         write(*,*)'   callphys.def with datadir=/path/to/dir'
+         write(*,*)'2) If necessary, /cross_sections (and other datafiles)'
+         write(*,*)'   can be obtained online on:'
+         write(*,*)' http://www.lmd.jussieu.fr/~lmdz/planets/mars/datadir'
+         stop
+      end if
 
       DO i = 1,nhead