-
Notifications
You must be signed in to change notification settings - Fork 1
/
param_em3p8_Ca.h
70 lines (70 loc) · 2.99 KB
/
param_em3p8_Ca.h
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
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
! param_em3p8_Ca.h
!
! restart: numbr2, numbr1
! EM parallel: num_proc
!
! Give Lengthx in Angstrom, then is converted to cm !!
!
integer(C_INT) &
ns0,np0,nq0,kstart,kgrp,isizeX,isizeY,isizeZ, &
num_proc,lxy3,mx,my,mz,mza,mxyza,mxh,myh,mzh, &
npq0,n00,n10,nbxc,nbxs,nbx2,nc3
real(C_DOUBLE) &
sconv,Lenx3,Leny3,Lenz3,xmax3,ymax3,zmax3, &
xmin3,ymin3,zmin3
logical iflinx
character :: sname*6,cname*6,numbr2*2,numbr1*2,numbr0*1
! ----------------------------------------------------------
! e(-6) ns0/2 + e(-70/4) 4*ns0/2
parameter (ns0=110600,np0=10120,nq0=np0+5*ns0/2)
!
parameter (sname='Cntemp',cname='cntemp') ! cntemp
parameter (kstart=0,kgrp=1)
! parameter (kstart=1,kgrp=1)
parameter (numbr2='Ca',numbr1='Ca',numbr0='C') ! .Ca
!
parameter (iflinx=.true.,num_proc=12) ! 600/12=50
! parameter (iflinx=.false.,num_proc=96) ! 768/96= 8
!
parameter (isizeX=50,isizeY=50,isizeZ=150) ! isizeX 10 Ang !
parameter (mx=201,my=201,mz=600) ! Grid: 2.5 Ang, mz=600
parameter (Lenx3=500.d0,Leny3=500.d0,Lenz3=1500.d0)
parameter (mza=50) ! mz=600, 12 ranks
parameter (mxyza=mx*my*mza)
parameter (mxh=101,myh=101,mzh=300)
!
! parameter (isizeX=50,isizeY=50,isizeZ=192) ! isizeX 10 Ang !
! parameter (Lenx3=500.d0,Leny3=500.d0,Lenz3=1920.d0)
! parameter (mza=8) ! mz=768, 96 ranks
! parameter (mx=201,my=201,mz=768) ! Grid: 2.5 Ang, mz=768
! parameter (mxyza=mx*my*mza)
! parameter (mxh=101,myh=101,mzh=384)
! ----------------------------------------------------------
parameter (lxy3= 3*mx*my)
parameter (sconv=1.0d-8)
!
parameter (npq0=ns0+np0+nq0)
parameter (n00=npq0/num_proc+1,n10=ns0/num_proc+1)
! parameter (nbxs=7000,nbxc=5000,nbx2=5000)
parameter (nbxs=3500,nbxc=2000,nbx2=2000)
parameter (nc3=isizeX*isizeY*isizeZ)
!
! zero at ymin3< 0 <ymax3
! E*B direction is -y
parameter (xmax3= (100/200.d0)*sconv*Lenx3, & ! EM field 200*400 meshes
xmin3=-(100/200.d0)*sconv*Lenx3, & !
ymax3= (100/200.d0)*sconv*Leny3, & ! E*H field
ymin3=-(100/200.d0)*sconv*Leny3, & !
zmax3= (300/600.d0)*sconv*Lenz3, & ! axial direction
zmin3=-(300/600.d0)*sconv*Lenz3) !
! zmax3= (384/768.d0)*sconv*Lenz3, & ! axial direction
! zmin3=-(384/768.d0)*sconv*Lenz3) !
!
character praefixs*33,praefixc*24,praefixe*24, &
praefixi*24,suffix2*2,suffix1*2,suffix0*1
! character praefixs*33,praefixc*32,praefixe*32, &
! praefixi*32,suffix2*2,suffix1*2,suffix0*1
common/filname/ praefixs,praefixc,praefixe, &
suffix2,suffix1,suffix0
common/filnam2/ praefixi
!