Updated to the current version of MKL.
[ptslat.git] / pzodrv.f90
blob496e523fde9c207e6b4cf85fbdfac3b30fc6dc5d
1 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
2 ! Driver for stain0.f subroutine
3 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
5 PROGRAM PZODRV
7 Use Dot_Geometry
8 Use Auxiliar_Procedures
9 Use Input_Data
10 Use NCPACK
12 IMPLICIT NONE
14 !! GRID RESULTS OF STRAIN0
16 REAL, ALLOCATABLE, DIMENSION(:,:,:) :: P_SPONT,P_PIEZO
18 INTERFACE
19 SUBROUTINE PIEZO(P_SPONT,P_PIEZO)
20 REAL,DIMENSION(:,:,:) :: P_SPONT,P_PIEZO
21 END SUBROUTINE
22 END INTERFACE
24 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
26 CALL READ_INPUT()
28 ALLOCATE (P_SPONT(1:XDim,1:YDim,1:ZDim),&
29 P_PIEZO(1:XDim,1:YDim,1:ZDim) )
32 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
33 ! Quantum Dot Dimensions (A)
34 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
37 RC=Rqd_Base+23.E0
38 ZC=Hqd+23.E0
41 RD = Rqd_Base/RC ! RD = Rqd_Base normalized to RC
42 HD = Hqd/ZC ! HD = HQD normalized to ZC
44 D=DWL/ZC ! Wetting layer thick. normalized to ZC
46 CALL PIEZO(P_SPONT,P_PIEZO)
47 CALL NCPACK_PZ(P_SPONT,P_PIEZO)
50 END PROGRAM PZODRV