allocate in reconstruct

Author: sanromd

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
 
 ! ===================================================================