Changeset 1000 for LMDZ4/trunk/libf/dyn3dpar/initfluxsto_p.F

Timestamp:

Oct 6, 2008, 10:43:22 AM (16 years ago)

Author:

Laurent Fairhead

Message:

• Modifs sur le parallelisme: masquage dans la physique
• Inclusion strato
• mise en coherence etat0
• le mode offline fonctionne maintenant en parallele,
• les fichiers de la dynamiques sont correctement sortis et peuvent etre reconstruit avec rebuild
• la version parallele de la dynamique peut s'executer sans MPI (sur 1 proc)
• L'OPENMP fonctionne maintenant sans la parallelisation MPI.

YM
LF

File:

• LMDZ4/trunk/libf/dyn3dpar/initfluxsto_p.F (modified) (3 diffs)

LMDZ4/trunk/libf/dyn3dpar/initfluxsto_p.F

@@ -74 +74 @@
       logical ok_sync
       integer :: jjb,jje,jjn
+
+! definition du domaine d'ecriture pour le rebuild
+
+      INTEGER,DIMENSION(2) :: ddid
+      INTEGER,DIMENSION(2) :: dsg
+      INTEGER,DIMENSION(2) :: dsl
+      INTEGER,DIMENSION(2) :: dpf
+      INTEGER,DIMENSION(2) :: dpl
+      INTEGER,DIMENSION(2) :: dhs
+      INTEGER,DIMENSION(2) :: dhe 
+      
+      INTEGER :: dynu_domain_id
+      INTEGER :: dynv_domain_id
+
+
 C
 C  Initialisations
@@ -100 +115 @@
       jje=jj_end
       jjn=jj_nb
+
+      ddid=(/ 1,2 /)
+      dsg=(/ iip1,jjp1 /)
+      dsl=(/ iip1,jjn /)
+      dpf=(/ 1,jjb /)
+      dpl=(/ iip1,jje /)
+      dhs=(/ 0,0 /)
+      dhe=(/ 0,0 /)
+
+      call flio_dom_set(mpi_size,mpi_rank,ddid,dsg,dsl,dpf,dpl,dhs,dhe, 
+     .                 'box',dynu_domain_id)
        
-      call histbeg(trim(infile)//'_'//trim(int2str(mpi_rank))//'.nc',
-     .             iip1, rlong(:,1), jjp1, rlat(1,jjb:jje),
-     .             1, iip1, 1, jjn,
-     .             tau0, zjulian, tstep, uhoriid, fileid)
+      call histbeg(trim(infile),iip1, rlong(:,1), jjn, rlat(1,jjb:jje),
+     .             1, iip1, 1, jjn, tau0, zjulian, tstep, uhoriid,
+     .             fileid,dynu_domain_id)
 C
 C  Creation du fichier histoire pour la grille en V (oblige pour l'instant,
@@ -124 +149 @@
       if (pole_sud) jjn=jj_nb-1
 
-      call histbeg('fluxstokev_'//trim(int2str(mpi_rank))//'.nc',
-     .             iip1, rlong(:,1), jjm, rlat(1,jjb:jje),
-     .             1, iip1, 1, jjn,
-     .             tau0, zjulian, tstep, vhoriid, filevid)
+      ddid=(/ 1,2 /)
+      dsg=(/ iip1,jjm /)
+      dsl=(/ iip1,jjn /)
+      dpf=(/ 1,jjb /)
+      dpl=(/ iip1,jje /)
+      dhs=(/ 0,0 /)
+      dhe=(/ 0,0 /)
+
+      call flio_dom_set(mpi_size,mpi_rank,ddid,dsg,dsl,dpf,dpl,dhs,dhe, 
+     .                 'box',dynv_domain_id)
+     
+      call histbeg('fluxstokev',iip1, rlong(:,1), jjn, rlat(1,jjb:jje),
+     .             1, iip1, 1, jjn,tau0, zjulian, tstep, vhoriid,
+     .             filevid,dynv_domain_id)
         
       rl(1,1) = 1.