Changeset 1130 for trunk/LMDZ.MARS/libf/phymars/start2archive.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/phymars/start2archive.F (moved) (moved from trunk/LMDZ.MARS/libf/dyn3d/start2archive.F) (6 diffs)

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

@@ -19 +19 @@
 c=======================================================================
 
-      use infotrac, only: iniadvtrac, nqtot, tnom
+      use infotrac, only: iniadvtrac, nqtot, tname
       use comsoil_h, only: nsoilmx, inertiedat
       use surfdat_h, only: ini_surfdat_h
       use comsoil_h, only: ini_comsoil_h
+      use comgeomphy, only: initcomgeomphy
       implicit none
 
@@ -34 +35 @@
 #include "logic.h"
 #include "temps.h"
-#include "control.h"
+!#include "control.h"
 #include "ener.h"
 #include "description.h"
@@ -119 +120 @@
       CALL defrun_new(99, .TRUE. )
       grireg   = .TRUE.
+! initialize "serial/parallel" related stuff
+      CALL init_phys_lmdz(iim,jjp1,llm,1,(/(jjm-1)*iim+2/))
+      call initcomgeomphy
 
 c=======================================================================
@@ -355 +359 @@
 c        call write_archive(nid,ntime,'q'//str2,'tracer','kg/kg',
 c     .         3,q(1,1,iq))
-        call write_archive(nid,ntime,tnom(iq),'tracer','kg/kg',
+        call write_archive(nid,ntime,tname(iq),'tracer','kg/kg',
      &         3,q(1,1,iq))
       end do
@@ -365 +369 @@
 c       call write_archive(nid,ntime,'qsurf'//str2,'Tracer on surface',
 c     $  'kg.m-2',2,qsurfS(1,iq))
-        txt=trim(tnom(iq))//"_surf"
+        txt=trim(tname(iq))//"_surf"
         call write_archive(nid,ntime,txt,'Tracer on surface',
      &  'kg.m-2',2,qsurfS(1,iq))
@@ -404 +408 @@
 c-----------------------------------------------------------------------
 
+      write(*,*) "startarchive: all is well that ends well"
+      
       end