Changeset 635 for trunk/LMDZ.MARS/libf/phymars

Timestamp:

Apr 26, 2012, 3:22:19 PM (13 years ago)

Author:

emillour

Message:

Mars GCM: Update of the chemistry package, including:

• 93 reactions are accounted for (instead of 22); tracking 28 species (instead of 11)
• computation of photoabsorption using raytracing
• improved time stepping in the photochemistry
• updated parameters (cross-sections); with this new version input files

are in 'EUV/param_v5' of "datafile" directory.

• transition between lower and upper atmosphere chemistry set to 0.1 Pa (calchim.F90)
• Lots of code clean-up: removed obsolete files column.F, param_v3.h, flujo.F, phdisrate.F, ch.F, interpfast.F, paramfoto.F, getch.F Converted chemtermos.F -> chemthermos.F90 and euvheat.F -> euvheat.F90. Added paramfoto_compact.F , param_v4.h and iono.h
• Upadted surfacearea.F
• Cleaned initracer.F and callkeys.h (removed use of obsolete "nqchem" and "oldnames" case when initializing tracers).
• Minor correction in "callsedim": compute "rdust" and/or "rice" only when it makes sense.

FGG+FL+EM

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• callkeys.h (modified) (2 diffs)
• callsedim.F (modified) (4 diffs)
• initracer.F (modified) (8 diffs)
• physiq.F (modified) (1 diff)
• tracer.h (modified) (2 diffs)

trunk/LMDZ.MARS/libf/phymars/callkeys.h

@@ -16 +16 @@
      
       COMMON/callkeys_i/iradia,iaervar,iddist,ilwd,ilwb,ilwn,ncouche    &
-     &   ,dustbin,nqchem_min,nltemodel,nircorr
+     &   ,dustbin,nltemodel,nircorr
      
       COMMON/callkeys_r/topdustref,solarcondate,semi,alphan,euveff,     &
@@ -54 +54 @@
       logical sedimentation,activice,water,microphys,caps
       logical photochem
-      integer nqchem_min
       integer nltemodel
       integer nircorr

trunk/LMDZ.MARS/libf/phymars/callsedim.F

@@ -434 +434 @@
 c     Update the dust particle size "rdust"
 c     -------------------------------------
-      DO l = 1, nlay
+      if (doubleq) then
+       DO l = 1, nlay
         DO ig=1,ngrid
           rdust(ig,l)=
@@ -441 +442 @@
           rdust(ig,l)=min(max(rdust(ig,l),1.e-10),500.e-6)
         ENDDO
-      ENDDO
+       ENDDO
+      endif ! of if (doubleq)
       
 c     Update the ice particle size "rice"
 c     -------------------------------------
-      IF(scavenging) THEN 
+      if (water) then
+       IF(scavenging) THEN 
         DO l = 1, nlay
           DO ig=1,ngrid
@@ -464 +467 @@
           ENDDO
         ENDDO
-      ELSE
+       ELSE
         DO l = 1, nlay
           DO ig=1,ngrid
@@ -475 +478 @@
           ENDDO
         ENDDO
-      ENDIF
+       ENDIF ! of IF(scavenging)
+      endif ! of if (water)
 
       RETURN

trunk/LMDZ.MARS/libf/phymars/initracer.F

@@ -24 +24 @@
 c                  q(iq=2) is the dust number mixing ratio 
 
-c   If (photochem.or.thermochem) there is "ncomp" chemical species (ncomp
-c   is set in aeronomars/chimiedata.h) using the ncomp iq values starting at
-c      iq=nqchem_min = dustbin+1   (nqchem_min is defined in inifis.F)
-c   
 c
 c   author: F.Forget
@@ -43 +39 @@
 #include "advtrac.h"
 #include "comgeomfi.h"
-#include "chimiedata.h"
 
 #include "surfdat.h"
@@ -55 +50 @@
 c       et al. (2002), a value of 25 has been deduced;
       real, parameter :: popratio = 25.
-      logical :: oldnames ! =.true. if old tracer naming convention (q01,...)
       character(len=20) :: txt ! to store some text
 
@@ -71 +65 @@
 c-----------------------------------------------------------------------
 
-      integer :: nqchem_start
-
 ! Initialization: get tracer names from the dynamics and check if we are
 !                 using 'old' tracer convention ('q01',q02',...)
 !                 or new convention (full tracer names)
       ! check if tracers have 'old' names
-      oldnames=.false.
       count=0
       do iq=1,nqmx
@@ -90 +81 @@
         write(*,*) "initracer: tracers seem to follow old naming ",
      &             "convention (q01,q02,...)"
-        write(*,*) "   => will work for now ... "
-        write(*,*) "      but you should run newstart to rename them"
-        oldnames=.true.
+        write(*,*) "you should run newstart to rename them"
+        stop
       endif
 
-      ! copy/set tracer names
-      if (oldnames) then ! old names (derived from iq & option values)
-        ! 1. dust:
-        if (dustbin.ne.0) then ! transported dust
-          do iq=1,dustbin
-            txt=" "
-            write(txt,'(a4,i2.2)') 'dust',iq
-            noms(iq)=txt
-            mmol(iq)=100.
-          enddo
-          ! special case if doubleq
-          if (doubleq) then
-            if (dustbin.ne.2) then
-              write(*,*) 'initracer: error expected dustbin=2'
-            else
-!              noms(1)='dust_mass'   ! dust mass mixing ratio
-!              noms(2)='dust_number' ! dust number mixing ratio
-! same as above, but avoid explict possible out-of-bounds on array
-               noms(1)='dust_mass'   ! dust mass mixing ratio
-               do iq=2,2
-               noms(iq)='dust_number' ! dust number mixing ratio
-               enddo
-            endif
-          endif
-        endif
-        ! 2. water & ice
-        if (water) then
-          noms(nqmx)='h2o_vap'
-          mmol(nqmx)=18.
-!            noms(nqmx-1)='h2o_ice'
-!            mmol(nqmx-1)=18.
-          !"loop" to avoid potential out-of-bounds in array
-          do iq=nqmx-1,nqmx-1
-            noms(iq)='h2o_ice'
-            mmol(iq)=18.
-          enddo
-        endif
-        ! 3. Chemistry
-        if (photochem .or. callthermos) then
-          if (doubleq) then
-            nqchem_start=3
-          else
-            nqchem_start=dustbin+1
-          end if
-          do iq=nqchem_start,nqchem_start+ncomp-2
-            noms(iq)=nomchem(iq-nqchem_start+1)
-            mmol(iq)=mmolchem(iq-nqchem_start+1)
-          enddo
-        endif ! of if (photochem .or. callthermos)
-        ! 4. Other tracers ????
-        if ((dustbin.eq.0).and.(.not.water)) then
-          noms(1)='co2'
-          mmol(1)=44
-          if (nqmx.eq.2) then
-            noms(nqmx)='Ar_N2'
-            mmol(nqmx)=30
-          endif
-        endif
-      else
-        ! copy tracer names from dynamics
-        do iq=1,nqmx
-          noms(iq)=tnom(iq)
-        enddo
-        ! mmol(:) array is initialized later (see below)
-      endif ! of if (oldnames)
-
-      ! special modification when using 'old' tracers:
-      ! qsurf(nqmx) was h2o ice whereas q(nqmx) was water vapour
-      ! and (if iceparty) q(nqmx-1) was null whereas q(nqmx-1) was water ice
-      if (oldnames.and.water) then
-        write(*,*)'initracer: moving surface water ice to index ',nqmx-1
-        ! "loop" to avoid potential out-of-bounds on arrays
-        do iq=nqmx-1,nqmx-1
-          qsurf(1:ngridmx,iq)=qsurf(1:ngridmx,iq+1)
-        enddo
-        qsurf(1:ngridmx,nqmx)=0
-      endif 
-      
-c------------------------------------------------------------
-c     Test Dimensions tracers
-c------------------------------------------------------------
-
-! Ehouarn: testing number of tracers is obsolete...
-!      if(photochem.or.thermochem) then
-!          if (water) then
-!              if ((dustbin+ncomp+2).ne.nqmx) then
-!                 print*,'initracer: tracer dimension problem:'
-!                 print*,"(dustbin+ncomp+2).ne.nqmx"
-!                 print*,"ncomp: ",ncomp
-!                 print*,"dustbin: ",dustbin
-!                 print*,"nqmx: ",nqmx
-!                 print*,'Change ncomp in chimiedata.h'
-!               endif
-!          else
-!              if ((dustbin+ncomp+1).ne.nqmx) then
-!                 print*,'initracer: tracer dimension problem:'
-!                 print*,"(dustbin+ncomp+1).ne.nqmx"
-!                 print*,"ncomp: ",ncomp
-!                 print*,"dustbin: ",dustbin
-!                 print*,"nqmx: ",nqmx
-!                 print*,'Change ncomp in chimiedata.h'
-!                 STOP
-!               endif
-!            endif
-!      endif
+      ! copy tracer names from dynamics
+      do iq=1,nqmx
+        noms(iq)=tnom(iq)
+      enddo
 
 c------------------------------------------------------------
 c         NAME and molar mass of the tracer 
 c------------------------------------------------------------
-
    
 ! Identify tracers by their names: (and set corresponding values of mmol)
@@ -246 +134 @@
       igcm_n2plus=0
       igcm_hplus=0
+      igcm_hco2plus=0
       igcm_elec=0
-
 
       ! 1. find dust tracers
@@ -432 +320 @@
           count=count+1
         endif
+        if (noms(iq).eq."hco2plus") then
+          igcm_hco2plus=iq
+          mmol(igcm_hco2plus)=45.
+          count=count+1
+        endif
         if (noms(iq).eq."elec") then
           igcm_elec=iq
@@ -454 +347 @@
         endif
 
-
       enddo ! of do iq=1,nqmx
-!      count=count+nbqchem
       
       ! check that we identified all tracers:

trunk/LMDZ.MARS/libf/phymars/physiq.F

@@ -131 +131 @@
 #include "chimiedata.h"
 #include "param.h"
-#include "param_v3.h"
+#include "param_v4.h"
 #include "conc.h"
 

trunk/LMDZ.MARS/libf/phymars/tracer.h

@@ -66 +66 @@
       integer :: igcm_n2plus 
       integer :: igcm_hplus
+      integer :: igcm_hco2plus
       integer :: igcm_elec
       ! other tracers
@@ -85 +86 @@
      & igcm_co2plus,igcm_oplus,igcm_o2plus,igcm_coplus,igcm_cplus,      &
      & igcm_nplus,igcm_noplus,igcm_n2plus,igcm_hplus,igcm_elec,         &
-     & igcm_ar_n2!,nbqchem,niqchem
+     & igcm_hco2plus,igcm_ar_n2!,nbqchem,niqchem
       COMMON/tracer3/noms
 !-----------------------------------------------------------------------