Version 2 (modified by 12 years ago) ( diff ) | ,
---|
Clump Objects
The Clump object uses the sphere shape object but with a non-uniform setup
TYPE ClumpDef REAL(KIND=qPREC) :: density=10 !Clump peak density REAL(KIND=qPREC) :: pressure=0 !Clump Pressure REAL(KIND=qPREC) :: temperature=.1 !Clump temperature REAL(KIND=qPREC), DIMENSION(3) :: velocity= (/0,0,0/) !Clump velocity REAL(KIND=qPREC), DIMENSION(3) :: position = (/0,0,0/) !Clump velocity REAL(KIND=qPREC) :: radius=1 REAL(KIND=qPREC) :: phi=0 REAL(KIND=qPREC) :: theta=0 REAL(KIND=qPREC) :: thickness=.1 !thickness of clump smoothing region ! REAL(KIND=qPREC) :: B_tor=0 !Maximum Bfield for toroidal configuration ! REAL(KIND=qPREC) :: B_pol=0 !Maximum Bfield for poloidal configuration REAL(KIND=qPREC) :: B_phi=0 REAL(KIND=qPREC) :: B_theta=0 REAL(KIND=qPREC) :: B_toroidal=0 !Maximum Bfield for toroidal configuration REAL(KIND=qPREC) :: B_poloidal=0 !Maximum Bfield for poloidal configuration REAL(KIND=qPREC) :: omega=0 !Angular velocity for solid body rotation INTEGER :: RefineLevel=MaxDepth INTEGER :: SubSample=1 INTEGER :: iTracer=0 !Clump Tracer INTEGER :: density_profile=0 REAL(KIND=qPREC) :: m2A=0 !Amplitude of an m=2 density perturbation [ rho = rho_0*(1+m2A*cos(2*phi)) ] LOGICAL :: PersistInBoundaries = .false. INTEGER :: id TYPE(ShapeDef), POINTER :: Shape => Null() TYPE(PerturbationDef), POINTER :: DensityPerturbation => Null() TYPE(VectorPerturbationDef), POINTER :: VelocityPerturbation => Null() TYPE(PerturbationDef), DIMENSION(:), POINTER :: MagneticPerturbation => Null() INTEGER :: ObjId END TYPE ClumpDef
After specifying the clump properties you only have to run
CALL InitClump(Clump)
Note:
See TracWiki
for help on using the wiki.