Changeset 1130 for trunk/LMDZ.MARS/libf/dyn3d/lect_start_archive.F

Timestamp:

Dec 20, 2013, 4:04:56 PM (11 years ago)

Author:

emillour

Message:

Mars GCM:
Series of changes to enable running in parallel (using LMDZ.COMMON dynamics);
Current LMDZ.MARS can still notheless be compiled and run in serial mode
"as previously".
Summary of main changes:

• Main programs (newstart, start2archive, xvik) that used to be in dyn3d have been moved to phymars.
• dyn3d/control.h is now module control_mod.F90
• rearanged input/outputs routines everywhere to handle serial/MPI cases. physdem.F => phyredem.F90 , phyetat0.F => phyetat0.F90 ; all read/write routines for startfi files are gathered in module iostart.F90
• added parallelism related routines init_phys_lmdz.F90, comgeomphy.F90, dimphy.F90, iniphysiq.F90, mod_grid_phy_lmdz.F90, mod_phys_lmdz_mpi_data.F90, mod_phys_lmdz_mpi_transfert.F90, mod_phys_lmdz_omp_data.F90, mod_phys_lmdz_omp_transfert.F90, mod_phys_lmdz_para.F90, mod_phys_lmdz_transfert_para.F90 in phymars and mod_const_mpi.F90 in dyn3d (for compliance with parallel case)
• created generic routines 'planetwide_maxval' and 'planetwide_minval', in module "planetwide_mod", that enable obtaining the min and max of a field over the whole planet.

EM

File:

• trunk/LMDZ.MARS/libf/dyn3d/lect_start_archive.F (modified) (4 diffs)

trunk/LMDZ.MARS/libf/dyn3d/lect_start_archive.F

@@ -17 +17 @@
 c
 c=======================================================================
-      use infotrac, only: tnom
+      use infotrac, only: tname
       use comsoil_h, only: nsoilmx, layer, mlayer, volcapa, inertiedat
       implicit none
@@ -32 +32 @@
 #include "comvert.h"
 #include "comgeom2.h"
-#include "control.h"
+!#include "control.h"
 #include "logic.h"
 #include "description.h"
 #include "ener.h"
 #include "temps.h"
-#include "lmdstd.h"
+!#include "lmdstd.h"
 #include "netcdf.inc"
 !#include "tracer.h"
@@ -760 +760 @@
           write(txt,'(a5,i2.2)')'qsurf',iq
         ELSE
-          txt=trim(tnom(iq))//"_surf"
+          txt=trim(tname(iq))//"_surf"
           if (txt.eq."h2o_vap") then
             ! There is no surface tracer for h2o_vap;
@@ -948 +948 @@
           write(txt,'(a1,i2.2)')'q',iq
         ELSE
-          txt=tnom(iq)
+          txt=tname(iq)
         ENDIF
         write(*,*)"lect_start_archive: loading tracer ",trim(txt)