-
Notifications
You must be signed in to change notification settings - Fork 19
/
WriteGridInfo.F90
40 lines (30 loc) · 1.16 KB
/
WriteGridInfo.F90
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
38
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! !
! FILE: WriteGridInfo.F90 !
! CONTAINS: subroutine WriteGridInfo !
! !
! PURPOSE: Write the grid information in !
! cordin_info.h5 !
! !
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
subroutine WriteGridInfo
use mpih
use param
use hdf5
IMPLICIT none
character*70 namfile
character*30 :: dsetname
if (ismaster) then
namfile='cordin_info.h5'
call HdfCreateBlankFile(namfile)
dsetname = trim('xm')
call HdfSerialWriteReal1D(dsetname,namfile,xm,nxm)
dsetname = trim('xc')
call HdfSerialWriteReal1D(dsetname,namfile,xc,nx)
dsetname = trim('ym')
call HdfSerialWriteReal1D(dsetname,namfile,ym,nym)
dsetname = trim('zm')
call HdfSerialWriteReal1D(dsetname,namfile,zm,nzm)
endif
return
end