strdrv.f90

   1 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
   2 ! Driver for stain0.f subroutine
   3 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
   4
   5   PROGRAM STRAINDRV
   6
   7         Use Dot_Geometry
   8         Use Auxiliar_Procedures
   9         Use Input_Data
  10         Use NCPACK
  11
  12       IMPLICIT NONE
  13
  14 !! GRID RESULTS OF STRAIN0
  15
  16       REAL, ALLOCATABLE, DIMENSION(:,:,:) :: EXX,EYY,EZZ,EXZ,EYZ,EXY
  17
  18       INTERFACE
  19        SUBROUTINE STRAIN0(EXX,EYY,EZZ,EXY,EXZ,EYZ)
  20          REAL,DIMENSION(:,:,:) ::  EXX,EYY,EZZ,EXY,EXZ,EYZ
  21        END SUBROUTINE
  22       END INTERFACE
  23
  24 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
  25
  26       CALL READ_INPUT()
  27
  28       ALLOCATE (EXX(1:XDim,1:YDim,1:ZDim),EYY(1:XDim,1:YDim,1:ZDim), &
  29                 EZZ(1:XDim,1:YDim,1:ZDim),EXZ(1:XDim,1:YDim,1:ZDim), &
  30                 EXY(1:XDim,1:YDim,1:ZDim),EYZ(1:XDim,1:YDim,1:ZDim) )
  31
  32
  33 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
  34 !  Quantum Dot Dimensions (A)
  35 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
  36
  37
  38        RC=Rqd_Base
  39        ZC=Hqd
  40
  41
  42        RD  = Rqd_Base/RC    ! RD = Rqd_Base normalized to RC
  43        HD  = Hqd/ZC         ! HD = HQD normalized to ZC
  44
  45        D=DWL/ZC ! Wetting layer thick. normalized to ZC
  46
  47        CALL STRAIN0(EXX,EYY,EZZ,EXY,EXZ,EYZ)
  48
  49
  50        IF (KCOOR.EQ.0) THEN
  51            CALL NCPACK_CART(EXX,EYY,EZZ,EXY,EXZ,EYZ)
  52        ELSE
  53            CALL NCPACK_CYL(EXX,EYY,EZZ,EXZ)
  54        END IF
  55
  56
  57   END PROGRAM STRAINDRV