print_starwall_m_ee.f90 918 Bytes
Newer Older
Serhiy Mochalskyy's avatar
Serhiy Mochalskyy committed
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
subroutine print_starwall_m_ee

use mpi_v
use resistive
use solv
use sca
!-----------------------------------------------------------------------
implicit none
include "mpif.h"

integer  :: i_loc,j_loc
real     :: num
!-----------------------------------------------------------------------
   
   CALL DESCINIT(DESCA, nd_bez,nd_bez, NB, NB, 0, 0, CONTEXT, LDA_ee, INFO_A )
   if(INFO_A .NE. 0) then
          write(6,*) "Something is wrong in output  CALL DESCINIT DESCA, INFO_A=",INFO_A
          stop
   endif


   if(rank==0) then
      open(60,file='starwall_m_ee',form="formatted",iostat=ier)
      write(60,'(2i8)') nd_bez,nd_bez
 
   endif

   DO i_loc = 1,nd_bez
       DO j_loc = 1,nd_bez
            CALL pdelget('A','D',num, a_ee_loc,i_loc,j_loc,DESCA)
            if(rank==0)     write(60,'(2i8,1pe16.8)') i_loc,j_loc,num
       END DO
   END DO

if(rank==0) close(60)

end subroutine print_starwall_m_ee