allocate in reconstruct

clawpack · Mar 18, 2016 · 4a23e24 · 4a23e24
1 parent 9d28e7b
commit 4a23e24
Showing 1 changed file with 11 additions and 1 deletion.
diff --git a/src/pyclaw/sharpclaw/reconstruct.f90 b/src/pyclaw/sharpclaw/reconstruct.f90
@@ -48,6 +48,9 @@ subroutine alloc_recon_workspace(maxnx,num_ghost,num_eqn,num_waves,lim_type,char
             case(3)
                 allocate(uu(2,maxnx+2*num_ghost))
                 allocate( dq1m(maxnx+2*num_ghost))
+            case(4)
+                allocate(uu(2,maxnx+2*num_ghost))
+                allocate( dq1m(maxnx+2*num_ghost))
         end select
         recon_alloc = .True.
 
@@ -83,8 +86,15 @@ subroutine dealloc_recon_workspace(lim_type,char_decomp)
                     deallocate(hh)
                     deallocate(uh)
             end select
-            recon_alloc = .False.
+            case(3)
+                deallocate(uu)
+                deallocate(dq1m)
+            case(4)
+                deallocate(uu)
+                deallocate(dq1m)
         end select
+        recon_alloc = .False.
+
     end subroutine dealloc_recon_workspace
 
 ! ===================================================================