Changes between Version 7 and Version 8 of u/bliu/pnMapping

Timestamp:

11/19/21 13:42:05 (4 years ago)

Author:

Baowei Liu

Comment:

—

u/bliu/pnMapping

@@ -173 +173 @@
 
 {{{
+SUBROUTINE ProblemModuleInit()
+    INTEGER :: n,i,j,k,edge
+    real(KIND=qPREC) :: rhoambient,pambient,mass,xloc(3),vel(3),alpha,radiusw
+    integer :: grav_soft_rad
+    TYPE(ProjectionDef), POINTER :: Projection
+    NAMELIST/AmbientData/ rhoambient,pambient
+    NAMELIST /ParticleData/ mass,xloc,vel,alpha,radiusw,grav_soft_rad
+    NAMELIST/ProblemData/ outflowstarttime, outflowendtime, outflowdensity, mintracer
+    open(unit=problem_data_handle,file='problem.data',status='old')
+    read(problem_data_handle,nml=ProblemData)
+
+    ! Make sure interpopts starts at zero for initial grid setup
+    if (.not. lRestart) then
+        InterpOpts=0
+        MinDensity=4d-6  ! << set mindensity higher to overwrite the odd spherical region, change to 1d-10 later in ProblemBeforeGlobalStep 
+        MinTemp=20       ! << set mintemp higher for frame 0, change to 10 later in ProblemBeforeGlobalStep 
+    end if
+
+    if (lRestart) then
+        particle=>SinkParticles%self                ! handle to the particle 
+        Call addtracer(particle%outflowobj%iTracer,'Outflow_Tracer')  ! <<< tracers need to be added in the same order when the code is restarted, but don't re-creat the particle - or the outflow object   
+    end if
+
+    ! Physical bounds of boxes
+    box_0=64000*reshape((/-1,-1,-1,1,1,1/),(/3,2/))
+    box_4=4000*reshape((/-1,-1,-1,1,1,1/),(/3,2/))
+    box_7=750*reshape((/-1,-1,-1,1,1,1/),(/3,2/))
+
+    ! index of boxes in each level's index space
+    ibox_0=nint(box_0-spread(GxBounds(:,1),2,2))/levels(0)%dx + spread((/1,0/),1,3)
+    if (maxlevel >= 4) then
+       ibox_4=nint(box_4-spread(GxBounds(:,1),2,2))/levels(4)%dx + spread((/1,0/),1,3)
+       if (maxlevel >= 7) then
+          ibox_7=nint(box_7-spread(GxBounds(:,1),2,2))/levels(7)%dx + spread((/1,0/),1,3)
+       end if
+    end if
+
+    CALL CreateRefinement(RefineJet)
+    CALL CreateShape(RefineJet%shape, RECTANGULAR_PRISM, (/750d0,750d0,0d0/)) !position defaults to 0,0,0
+    RefineJet%type=REFINE_INSIDE
+
     nx=Info%mx(1)
     ny=Info%mx(2)