Changeset 1001 for LMDZ4/trunk/libf/phylmd/iophy.F90

Timestamp:

Oct 6, 2008, 11:11:53 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/phylmd/iophy.F90 (modified) (2 diffs)

LMDZ4/trunk/libf/phylmd/iophy.F90

@@ -16 +16 @@
   INTERFACE histwrite_phy
     MODULE PROCEDURE histwrite2d_phy,histwrite3d_phy
-  END INTERFACE 
+  END INTERFACE
+
 
 contains
+
+  subroutine init_iophy_new(rlat,rlon)
+  USE dimphy
+  USE mod_phys_lmdz_para
+  USE mod_grid_phy_lmdz
+  USE ioipsl
+  implicit none
+  include 'dimensions.h'   
+    real,dimension(klon),intent(in) :: rlon
+    real,dimension(klon),intent(in) :: rlat
+
+    REAL,dimension(klon_glo)        :: rlat_glo
+    REAL,dimension(klon_glo)        :: rlon_glo
+    
+    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 :: i    
+
+    CALL gather(rlat,rlat_glo)
+    CALL bcast(rlat_glo)
+    CALL gather(rlon,rlon_glo)
+    CALL bcast(rlon_glo)
+    
+!$OMP MASTER  
+    ALLOCATE(io_lat(jjm+1-1/iim))
+    io_lat(1)=rlat_glo(1)
+    io_lat(jjm+1-1/iim)=rlat_glo(klon_glo)
+    IF (iim > 1) then
+      DO i=2,jjm
+        io_lat(i)=rlat_glo(2+(i-2)*iim)
+      ENDDO
+    ENDIF
+
+    ALLOCATE(io_lon(iim))
+    io_lon(:)=rlon_glo(2-1/iim:iim+1-1/iim)
+
+
+    allocate(tmp_tab2d(iim,jj_nb))
+    allocate(tmp_tab3d(iim,jj_nb,klev))
+    allocate(ndex2d(iim*jj_nb))
+    allocate(ndex3d(iim*jj_nb*klev))
+    ndex2d(:)=0
+    ndex3d(:)=0
+    
+    ddid=(/ 1,2 /)
+    dsg=(/ iim, jjm+1-1/iim /)
+    dsl=(/ iim, jj_nb /)
+    dpf=(/ 1,jj_begin /)
+    dpl=(/ iim, jj_end /)
+    dhs=(/ ii_begin-1,0 /)
+    if (mpi_rank==mpi_size-1) then
+      dhe=(/0,0/)
+    else
+      dhe=(/ iim-ii_end,0 /)  
+    endif
+    
+    call flio_dom_set(mpi_size,mpi_rank,ddid,dsg,dsl,dpf,dpl,dhs,dhe, &
+                      'APPLE',phys_domain_id)
+
+!$OMP END MASTER
+      
+  end subroutine init_iophy_new
 
   subroutine init_iophy(lat,lon)
@@ -139 +207 @@
   
   
-!  subroutine phy2dyn(field_phy,field_dyn,nlev)
-!  USE dimphy_old
-!  implicit none
-!  include 'dimensions.h'
-!  
-!    real,dimension(klon_mpi,nlev),intent(in) :: field_phy
-!    real,dimension(iim,jjphy_nb,nlev),intent(out) :: field_dyn
-!    integer,intent(in) :: nlev
-!    
-!    integer :: next
-!    integer :: j,l
-!    
-!      do l=1,nlev
-!               
-!       if (jjphy_begin==jjphy_end) then
-!         field_dyn(:,1,l)=0.
-!         field_dyn(iiphy_begin:iiphy_end,1,l)=field_phy(1:klon_mpi,l)
-!       else
-!       
-!        if (jjphy_begin==1) then
-!           field_dyn(:,1,l)=field_phy(1,l)
-!           next=2
-!        else
-!          field_dyn(:,1,l)=0.
-!          next=iim-iiphy_begin+2
-!          field_dyn(iiphy_begin:iim,1,l)=field_phy(1:next-1,l)   
-!        endif
-!       
-!         do j=2,jjphy_nb-1
-!           field_dyn(:,j,l)=field_phy(next:next+iim-1,l)
-!           next=next+iim
-!         enddo
-!         
-!         if (jjphy_end==jjm+1-1/iim) then
-!             field_dyn(:,jjphy_nb,l)=field_phy(klon_mpi,l)
-!         else
-!          field_dyn(:,jjphy_nb,l)=0.
-!          field_dyn(1:iiphy_end,jjphy_nb,l)=field_phy(next:next+iiphy_end-1,l)   
-!         endif
-!         
-!       endif
-!     
-!     enddo
-!        
-!    end subroutine phy2dyn         
-  
-          
+
 end module iophy