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

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

@@ -1 +1 @@
       SUBROUTINE soil_tifeedback(ngrid,nsoil,icecover,newtherm_i)
+
+      use tracer_mod, only: nqmx, igcm_h2o_ice, rho_ice
       IMPLICIT NONE
 
@@ -24 +26 @@
 #include "dimphys.h"
 #include "comsoil.h"
-#include "tracer.h"
+!#include "tracer.h"
 #include "surfdat.h"
 
@@ -41 +43 @@
 
       REAL icecover(ngrid,nqmx)         ! tracer on the surface (kg.m-2)
-                                        ! last one (iq=nqmx) is surface
+                                        ! (iq=igcm_h2o_ice) is surface
                                         ! water ice
 c Outputs