Changeset 1036 for trunk/LMDZ.MARS/libf/phymars/simpleclouds.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/simpleclouds.F (modified) (4 diffs)

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

@@ -4 +4 @@
      &             nq,tau,rice)
       USE updaterad
+      use tracer_mod, only: igcm_h2o_vap, igcm_h2o_ice
       implicit none
 c------------------------------------------------------------------
@@ -32 +33 @@
 #include "comcstfi.h"
 #include "callkeys.h"
-#include "tracer.h"
+!#include "tracer.h"
 #include "comgeomfi.h"
 #include "dimradmars.h"
@@ -73 +74 @@
       INTEGER ig,l
 
-      REAL zq(ngridmx,nlayermx,nqmx)  ! local value of tracers
-      REAL zq0(ngridmx,nlayermx,nqmx) ! local initial value of tracers
+      REAL zq(ngridmx,nlayermx,nq)    ! local value of tracers
+      REAL zq0(ngridmx,nlayermx,nq)   ! local initial value of tracers
       REAL zt(ngridmx,nlayermx)       ! local value of temperature
       REAL zqsat(ngridmx,nlayermx)    ! saturation
@@ -91 +92 @@
 c    -----------------
 
-c    On "update" la valeur de q(nqmx) (water vapor) et temperature.
+c    On "update" la valeur de q(nq) (water vapor) et temperature.
 c    On effectue qqes calculs preliminaires sur les couches :