Changes between Version 35 and Version 36 of ModulesOnAstroBear

Timestamp:

01/14/13 22:05:40 (13 years ago)

Author:

Jonathan

Comment:

—

ModulesOnAstroBear

@@ -276 +276 @@
     END DO
   END SUBROUTINE
-}}}
-
-
+ ...
+}}}
+
+There are a lot of other ways we could modify this simple example to have the clump be located anywhere, to have a density profile that is smoothed at the edge, to be a different temperature, or move with a particular velocity, etc...  Fortunately, clumps are a commonly used astrophysical object and there is a clump object module to assist users in easily creating clumps.  For example the above module could be rewritten as:
+
+{{{
+MODULE Problem
+  USE GlobalDeclarations
+  USE DataDeclarations
+  USE Clumps
+  USE Ambients
+  IMPLICIT NONE
+  SAVE
+  PUBLIC ProblemModuleInit, ProblemGridInit, &
+       ProblemBeforeStep, ProblemAfterStep, ProblemSetErrFlag, ProblemBeforeGlobalStep
+  PRIVATE REAL(KIND=qPREC) :: rho, radius
+
+CONTAINS
+
+  SUBROUTINE ProblemModuleInit()
+    TYPE(AmbientDef), POINTER :: Ambient
+    TYPE(ClumpDef), POINTER :: Clump
+    NAMELIST/ProblemData/ rho, radius
+    OPEN(UNIT=PROBLEM_DATA_HANDLE, FILE='problem.data', STATUS="OLD")
+    READ(PROBLEM_DATA_HANDLE,NML=ProblemData)
+    CLOSE(PROBLEM_DATA_HANDLE)
+    CALL CreateAmbient(Ambient)
+    CALL CreateClump(Clump)
+    Clump%density=rho
+    Clump%radius=radius
+    CALL UpdateClump(Clump)
+  END SUBROUTINE
+  
+  SUBROUTINE ProblemGridInit(Info)
+    TYPE(InfoDef), POINTER :: Info
+  END SUBROUTINE
+ 
+ ...
+}}}
+
+So what is going on here?  First we've added USE statements for the Clumps and Ambients modules.  Then in our ProblemModuleInit() routine we've declared an ambientdef and a clumpdef pointer.  Then in CreateAmbient() we create an ambient object with the default properties - density = 1, pressure = 1, velocity = 0, etc... which sets the entire grid to be uniform.  Then we create a clump object and modify its properties (density, and radius) before updated the clump object.  And we're done.
+
+
+And if we want to get more complicated - we can modify other clump attributes.