Changeset 1036 for trunk/LMDZ.MARS/libf/aeronomars/calchim.F90

Timestamp:

Sep 11, 2013, 2:34:44 PM (11 years ago)

Author:

emillour

Message:

Mars GCM: (a first step towards using parallel dynamics)

• IMPORTANT CHANGE: Implemented dynamic tracers. It is no longer necessary to compile the model with the '-t #' option, number of tracers is simply read from tracer.def file (as before). Adapted makegcm_* scripts (and co.) accordingly. Technical aspects of the switch to dynamic tracers are:
  • advtrac.h (in dyn3d) removed and replaced by module infotrac.F
  • tracer.h (in phymars) removed and replaced by module tracer_mod.F90 (which contains nqmx, the number of tracers, etc. and can be used anywhere in the physics).
• Included some side cleanups: removed unused files (in dyn3d) anldoppler2.F, anl_mcdstats.F and anl_stats-diag.F, and all the unecessary dimensions.* files in grid/dimension.
• Checked that changes are clean and that GCM yields identical results (in debug mode) to previous svn version.

EM

File:

• trunk/LMDZ.MARS/libf/aeronomars/calchim.F90 (modified) (8 diffs)

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

@@ -1 @@
-      subroutine calchim(ptimestep,pplay,pplev,pt,pdt,dist_sol,mu0,         &
+      subroutine calchim(nq,                                                &
+                         ptimestep,pplay,pplev,pt,pdt,dist_sol,mu0,         &
                          zzlev,zzlay,zday,pq,pdq,dqchim,dqschim,dqcloud,    &
                          dqscloud,tauref,co2ice,                            &
                          pu,pdu,pv,pdv,surfdust,surfice)
 
+      use tracer_mod, only: igcm_co2, igcm_co, igcm_o, igcm_o1d, igcm_o2, &
+                            igcm_o3, igcm_h, igcm_h2, igcm_oh, igcm_ho2,  &
+                            igcm_h2o2, igcm_ch4, igcm_n2, igcm_h2o_vap,   &
+                            igcm_no, igcm_n, igcm_no2, igcm_n2d,          &
+                            igcm_o2plus, igcm_co2plus, igcm_oplus,        &
+                            igcm_coplus, igcm_cplus, igcm_nplus,          &
+                            igcm_noplus, igcm_n2plus, igcm_hplus,         &
+                            igcm_hco2plus, igcm_elec, mmol
       implicit none
 
@@ -56 +65 @@
 #include "dimphys.h"
 #include "chimiedata.h"
-#include "tracer.h"
+!#include "tracer.h"
 #include "comcstfi.h"
 #include "callkeys.h"
@@ -63 +72 @@
 !     input:
 
+      integer,intent(in) :: nq ! number of tracers
       real :: ptimestep
       real :: pplay(ngridmx,nlayermx)    ! pressure at the middle of the layers
@@ -76 +86 @@
       real :: dist_sol                   ! distance of the sun (AU)
       real :: mu0(ngridmx)               ! cos of solar zenith angle (=1 when sun at zenith)
-      real :: pq(ngridmx,nlayermx,nqmx)  ! tracers mass mixing ratio
-      real :: pdq(ngridmx,nlayermx,nqmx) ! previous tendencies
+      real :: pq(ngridmx,nlayermx,nq)  ! tracers mass mixing ratio
+      real :: pdq(ngridmx,nlayermx,nq) ! previous tendencies
       real :: zday                       ! date (time since Ls=0, in martian days)
       real :: tauref(ngridmx)            ! optical depth at 7 hPa
@@ -86 +96 @@
 !     output:
 
-      real :: dqchim(ngridmx,nlayermx,nqmx) ! tendencies on pq due to chemistry
-      real :: dqschim(ngridmx,nqmx)         ! tendencies on qsurf 
-      real :: dqcloud(ngridmx,nlayermx,nqmx)! tendencies on pq due to condensation
-      real :: dqscloud(ngridmx,nqmx)        ! tendencies on qsurf 
+      real :: dqchim(ngridmx,nlayermx,nq) ! tendencies on pq due to chemistry
+      real :: dqschim(ngridmx,nq)         ! tendencies on qsurf 
+      real :: dqcloud(ngridmx,nlayermx,nq)! tendencies on pq due to condensation
+      real :: dqscloud(ngridmx,nq)        ! tendencies on qsurf 
 
 !     local variables:
 
       integer,save :: nbq                   ! number of tracers used in the chemistry
-      integer,save :: niq(nqmx)             ! array storing the indexes of the tracers
+      integer,allocatable,save :: niq(:)    ! array storing the indexes of the tracers
       integer :: iloc(1)            ! index of major species
       integer :: ig,l,i,iq,iqmax
@@ -133 +143 @@
 
       real    :: latvl1, lonvl1
-      real    :: zq(ngridmx,nlayermx,nqmx) ! pq+pdq*ptimestep before chemistry
+      real    :: zq(ngridmx,nlayermx,nq) ! pq+pdq*ptimestep before chemistry
                                            ! new mole fraction after
       real    :: zt(ngridmx,nlayermx)      ! temperature
@@ -149 +159 @@
       real :: ztemp(nlayermx)        !  Temperature (K)
       real :: zlocal(nlayermx)       !  Altitude (km)
-      real :: zycol(nlayermx,nqmx)   !  Composition (mole fractions)
+      real :: zycol(nlayermx,nq)   !  Composition (mole fractions)
       real :: szacol                 !  Solar zenith angle
       real :: surfice1d(nlayermx)    !  Ice surface area (cm2/cm3)
@@ -172 +182 @@
          end if
          ! find index of chemical tracers to use
+         allocate(niq(nq))
          ! Listed here are all tracers that can go into photochemistry
          nbq = 0 ! to count number of tracers