Ignore:
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:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.MARS/libf/grid/dimension/makdim

    r38 r1036  
    1 nqmx=$1
    2 shift
     1
    32for i in $* ; do
    43   list=$list.$i
    54done
    6 fichdim=dimensions${list}.t${nqmx}
     5fichdim=dimensions${list}
    76
    87
     
    1514   lm=$3
    1615   n2=$1
    17    ndm=1
     16#   ndm=1
    1817
    1918#  while [ "$n2" -gt 2 ]; do
     
    3029       jm=$1
    3130      lm=$2
    32       ndm=1
     31#      ndm=1
    3332 else if [ $# -ge 1 ] ; then
    3433         im=1
    3534         jm=1
    3635         lm=$1
    37          ndm=1
     36#         ndm=1
    3837 else
    3938         echo il faut au moins une dimension
     
    4948!   dimensions.h contient les dimensions du modele
    5049!   ndm est tel que iim=2**ndm
    51 !   nqmx est la dimension de la variable traceur q
    5250!-----------------------------------------------------------------------
    5351
    54       INTEGER, parameter :: iim= $im
     52      INTEGER, parameter :: iim=$im
    5553      INTEGER, parameter :: jjm=$jm
    5654      INTEGER, parameter :: llm=$lm
    57       INTEGER, parameter :: ndm=$ndm
    58 
    59       integer, parameter :: nqmx=$nqmx
    6055
    6156!-----------------------------------------------------------------------
Note: See TracChangeset for help on using the changeset viewer.