Changeset 460 for trunk

Timestamp:

Dec 8, 2011, 11:22:59 AM (13 years ago)

Author:

emillour

Message:

Mars GCM: more updates for photochemistry:
-aeronomars/inichim_newstart.F : add methane + bug correction in "oldnames" case
-aeronomars/read_phototable.F : adapt routine so that photolysis table file

name may be set in the def file with
phototable = filename
(default is "jmars.20111014")

EM

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/aeronomars/inichim_newstart.F (modified) (6 diffs)
• libf/aeronomars/read_phototable.F (modified) (2 diffs)

trunk/LMDZ.MARS/README

@@ -1317 +1317 @@
                      and rice.
  
+== 08/12/2011 == EM
+>> more updates for photochemistry:
+-aeronomars/inichim_newstart.F : add methane + bug correction in "oldnames" case
+-aeronomars/read_phototable.F : adapt routine so that photolysis table file
+                                name may be set in the def file with
+                                phototable = filename
+                                (default is "jmars.20111014")
+

trunk/LMDZ.MARS/libf/aeronomars/inichim_newstart.F

@@ -15 +15 @@
 c   Author:   Sebastien Lebonnois (08/11/2002)
 c   -------
-c Revised 07/2003 by Monica Angelats-i-Coll to use input files
-c Modified 10/2008 Identify tracers by their names Ehouarn Millour
+c   Revised 07/2003 by Monica Angelats-i-Coll to use input files
+c   Modified 10/2008 Identify tracers by their names Ehouarn Millour
+c   Modified 11/2011 Addition of methane (Franck Lefevre)
 c
 c   Arguments:
@@ -150 +151 @@
           end if
         endif ! of if (photochem .or. callthermos)
-        do iq=nqchem_start,nqchem_start+ncomp-1
+        do iq=nqchem_start,nqchem_start+ncomp-2
           noms(iq)=nomchem(iq-nqchem_start+1)
           mmol(iq)=mmolchem(iq-nqchem_start+1)
@@ -260 +261 @@
       igcm_ho2=0
       igcm_h2o2=0
+      igcm_ch4=0
       igcm_n2=0
       igcm_ar=0
@@ -358 +360 @@
           nbqchem=nbqchem+1
         endif
+        if (noms(iq).eq."ch4") then
+          igcm_ch4=iq
+          mmol(igcm_ch4)=16.
+          count=count+1
+          nbqchem=nbqchem+1
+        endif
         if (noms(iq).eq."n2") then
           igcm_n2=iq
@@ -408 +416 @@
       ! as in the old days, water vapour is last and water ice,
       ! if present, is just before water vapour
-      do iq=1,15 ! loop so as to avoid out-of-bounds on array
+      do iq=1,16 ! loop so as to avoid out-of-bounds on array
       if (iq==1) nqchem(i)=igcm_co2
       if (iq==2) nqchem(i)=igcm_co
@@ -420 +428 @@
       if (iq==10) nqchem(i)=igcm_ho2
       if (iq==11) nqchem(i)=igcm_h2o2
-      if (iq==12) nqchem(i)=igcm_n2
-      if (iq==13) nqchem(i)=igcm_ar
-      if (iq==14) nqchem(i)=igcm_h2o_ice
-      if (iq==15) nqchem(i)=igcm_h2o_vap
+      if (iq==12) nqchem(i)=igcm_ch4
+      if (iq==13) nqchem(i)=igcm_n2
+      if (iq==14) nqchem(i)=igcm_ar
+      if (iq==15) nqchem(i)=igcm_h2o_ice
+      if (iq==16) nqchem(i)=igcm_h2o_vap
       enddo
 

trunk/LMDZ.MARS/libf/aeronomars/read_phototable.F

@@ -19 +19 @@
 cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
 c
+! to use  'getin'
+      USE ioipsl_getincom 
+
       implicit none
 c
@@ -32 +35 @@
 c
       real xsza
+      character(len=128) :: phototable ! photolysis table file name
 c
 cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
 c
+! set photolysis table input file name
+      phototable="jmars.20111014" ! default
+      ! look for a " phototable= ..." option in def files
+      call getin("phototable",phototable)
+
       fic = 81
 
       open(fic, form = 'formatted', status = 'old',
-     &file =trim(datafile)//'/jmars.20101220',iostat=ierr)
+     &file =trim(datafile)//trim(phototable),iostat=ierr)
 
       IF (ierr.NE.0) THEN
-        write(*,*)'Error : cannot open photolysis lookup table'
+        write(*,*)'Error : cannot open photolysis lookup table ',
+     &             trim(phototable)
         write(*,*)'(in aeronomars/init_chemistry.F)'
         write(*,*)'It should be in :',trim(datafile),'/'
-        write(*,*)'1) You can change this directory address in '
-        write(*,*)'   file phymars/datafile.h'
+        write(*,*)'1) You can change this directory path in physiq.def'
+        write(*,*)'   with:'
+        write(*,*)'   datadir=/path/to/the/directory'
+        write(*,*)'2) You can change the input phototable file name in'
+        write(*,*)'   physiq.def with:'
+        write(*,*)'   phototable=filename'
         STOP
       ENDIF
 
 c
-      print*, 'read photolysis lookup table...'
+      print*, 'read photolysis lookup table ',trim(phototable)
 c
       do itau = 1,ntau