Changeset 1036 for trunk/LMDZ.MARS/libf/phymars/initracer.F

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/phymars/initracer.F (modified) (19 diffs)

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

@@ -1 @@
-      SUBROUTINE initracer(qsurf,co2ice)
-
+      SUBROUTINE initracer(ngrid,nq,qsurf,co2ice)
+
+       use infotrac, only: tnom
+       use tracer_mod
        IMPLICIT NONE
 c=======================================================================
@@ -11 +13 @@
 c
 c   Test of dimension :
-c   Initialize COMMON tracer in tracer.h, using tracer names provided
-c   by the dynamics in "advtrac.h"
+c   Initialize tracer related data in tracer_mod, using tracer names provided
+c   by the dynamics in "infotrac"
 c
 c
@@ -26 +28 @@
 #include "comcstfi.h"
 #include "callkeys.h"
-#include "tracer.h"
-#include "advtrac.h"
+!#include "tracer.h"
+!#include "advtrac.h"
 #include "comgeomfi.h"
 
 #include "surfdat.h"
 
-      real qsurf(ngridmx,nqmx)       ! tracer on surface (e.g.  kg.m-2)
-      real co2ice(ngridmx)           ! co2 ice mass on surface (e.g.  kg.m-2)
+      integer,intent(in) :: ngrid ! number of atmospheric columns
+      integer,intent(in) :: nq ! number of tracers
+      real,intent(out) :: qsurf(ngrid,nq) ! tracer on surface (e.g.  kg.m-2)
+      real,intent(out) :: co2ice(ngrid) ! co2 ice mass on surface (e.g.  kg.m-2)
+
       integer iq,ig,count
       real r0_lift , reff_lift, nueff_lift
@@ -43 +48 @@
 
 c-----------------------------------------------------------------------
-c  radius(nqmx)      ! aerosol particle radius (m)
-c  rho_q(nqmx)       ! tracer densities (kg.m-3)
-c  alpha_lift(nqmx)  ! saltation vertical flux/horiz flux ratio (m-1)
-c  alpha_devil(nqmx) ! lifting coeeficient by dust devil
+c  radius(nq)      ! aerosol particle radius (m)
+c  rho_q(nq)       ! tracer densities (kg.m-3)
+c  alpha_lift(nq)  ! saltation vertical flux/horiz flux ratio (m-1)
+c  alpha_devil(nq) ! lifting coeeficient by dust devil
 c  rho_dust          ! Mars dust density
 c  rho_ice           ! Water ice density
@@ -55 +60 @@
 c-----------------------------------------------------------------------
 
+! Initialization: allocate arrays in tracer_mod
+      allocate(noms(nq))
+      allocate(mmol(nq))
+      allocate(radius(nq))
+      allocate(rho_q(nq))
+      allocate(alpha_lift(nq))
+      allocate(alpha_devil(nq))
+      allocate(dryness(ngridmx))
+      allocate(igcm_dustbin(nq))
+
 ! Initialization: get tracer names from the dynamics and check if we are
 !                 using 'old' tracer convention ('q01',q02',...)
@@ -60 +75 @@
       ! check if tracers have 'old' names
       count=0
-      do iq=1,nqmx
+      do iq=1,nq
         txt=" "
         write(txt,'(a1,i2.2)') 'q',iq
@@ -66 +81 @@
           count=count+1
         endif
-      enddo ! of do iq=1,nqmx
+      enddo ! of do iq=1,nq
       
-      if (count.eq.nqmx) then
+      if (count.eq.nq) then
         write(*,*) "initracer: tracers seem to follow old naming ",
      &             "convention (q01,q02,...)"
@@ -76 +91 @@
 
       ! copy tracer names from dynamics
-      do iq=1,nqmx
+      do iq=1,nq
         noms(iq)=tnom(iq)
       enddo
@@ -86 +101 @@
 ! Identify tracers by their names: (and set corresponding values of mmol)
       ! 0. initialize tracer indexes to zero:
-      do iq=1,nqmx
-        igcm_dustbin(iq)=0
-      enddo
+      igcm_dustbin(1:nq)=0
+
       igcm_dust_mass=0
       igcm_dust_number=0
@@ -130 +144 @@
       count=0
       if (dustbin.gt.0) then
-        do iq=1,nqmx
+        do iq=1,nq
           txt=" "
           write(txt,'(a4,i2.2)')'dust',count+1
@@ -138 +152 @@
             mmol(iq)=100.
           endif
-        enddo !do iq=1,nqmx
+        enddo !do iq=1,nq
       endif ! of if (dustbin.gt.0)
       if (doubleq) then
-        do iq=1,nqmx
+        do iq=1,nq
           if (noms(iq).eq."dust_mass") then
             igcm_dust_mass=iq
@@ -153 +167 @@
       endif ! of if (doubleq)
       if (microphys) then
-        do iq=1,nqmx
+        do iq=1,nq
           if (noms(iq).eq."ccn_mass") then
             igcm_ccn_mass=iq
@@ -165 +179 @@
       endif ! of if (microphys)
       if (submicron) then
-        do iq=1,nqmx
+        do iq=1,nq
           if (noms(iq).eq."dust_submicron") then
             igcm_dust_submicron=iq
@@ -174 +188 @@
       endif ! of if (submicron)
       ! 2. find chemistry and water tracers
-      do iq=1,nqmx
+      do iq=1,nq
         if (noms(iq).eq."co2") then
           igcm_co2=iq
@@ -337 +351 @@
         endif
 
-      enddo ! of do iq=1,nqmx
+      enddo ! of do iq=1,nq
       
       ! check that we identified all tracers:
-      if (count.ne.nqmx) then
+      if (count.ne.nq) then
         write(*,*) "initracer: found only ",count," tracers"
-        write(*,*) "               expected ",nqmx
+        write(*,*) "               expected ",nq
         do iq=1,count
           write(*,*)'      ',iq,' ',trim(noms(iq))
@@ -349 +363 @@
       else
         write(*,*) "initracer: found all expected tracers, namely:"
-        do iq=1,nqmx
+        do iq=1,nq
           write(*,*)'      ',iq,' ',trim(noms(iq))
         enddo
@@ -365 +379 @@
        ! as qsurf(i_h2o_vap) & as qsurf(i_h2o_ice), so to be clean:
        if (igcm_h2o_vap.ne.0) then
-         qsurf(1:ngridmx,igcm_h2o_vap)=0
+         qsurf(1:ngrid,igcm_h2o_vap)=0
        endif
       endif
 
 c------------------------------------------------------------
-c     Initialize tracers .... (in tracer.h)
+c     Initialize tracers .... (in tracer_mod)
 c------------------------------------------------------------
       ! start by setting everything to (default) zero
-      rho_q(1:nqmx)=0     ! tracer density (kg.m-3)
-      radius(1:nqmx)=0.   ! tracer particle radius (m)
-      alpha_lift(1:nqmx) =0.  ! tracer saltation vertical flux/horiz flux ratio (m-1)
-      alpha_devil(1:nqmx)=0.  ! tracer lifting coefficient by dust devils
+      rho_q(1:nq)=0     ! tracer density (kg.m-3)
+      radius(1:nq)=0.   ! tracer particle radius (m)
+      alpha_lift(1:nq) =0.  ! tracer saltation vertical flux/horiz flux ratio (m-1)
+      alpha_devil(1:nq)=0.  ! tracer lifting coefficient by dust devils
 
 
@@ -393 +407 @@
 c       iq=1: Q mass mixing ratio, iq=2: N number mixing ratio
 
-        if( (nqmx.lt.2).or.(water.and.(nqmx.lt.4)) ) then
-          write(*,*)'initracer: nqmx is too low : nqmx=', nqmx
+        if( (nq.lt.2).or.(water.and.(nq.lt.4)) ) then
+          write(*,*)'initracer: nq is too low : nq=', nq
           write(*,*)'water= ',water,' doubleq= ',doubleq   
         end if
@@ -499 +513 @@
          alpha_lift(igcm_h2o_vap) =0.
          alpha_devil(igcm_h2o_vap)=0.
-         if(water.and.(nqmx.ge.2)) then
+         if(water.and.(nq.ge.2)) then
            radius(igcm_h2o_ice)=3.e-6
            rho_q(igcm_h2o_ice)=rho_ice
            alpha_lift(igcm_h2o_ice) =0.
            alpha_devil(igcm_h2o_ice)=0.
-         elseif(water.and.(nqmx.lt.2)) then
-            write(*,*) 'nqmx is too low : nqmx=', nqmx
+         elseif(water.and.(nq.lt.2)) then
+            write(*,*) 'nq is too low : nq=', nq
             write(*,*) 'water= ',water
          endif