From 61caabcf2d895b6b708bc6bab75325f45a7b8963 Mon Sep 17 00:00:00 2001
From: ESCOBAR Juan <escj@nuwa>
Date: Thu, 10 Mar 2022 14:52:18 +0100
Subject: [PATCH] Juan 10/03/2022:tensor/discretisation.f90, upgrade zt1d size
 & public zt1d_discretisation_allocate3d

---
 .../tensorproductmultigrid_Source/discretisation.f90       | 7 ++++---
 1 file changed, 4 insertions(+), 3 deletions(-)

diff --git a/src/ZSOLVER/tensorproductmultigrid_Source/discretisation.f90 b/src/ZSOLVER/tensorproductmultigrid_Source/discretisation.f90
index fa4092f54..295e13e96 100644
--- a/src/ZSOLVER/tensorproductmultigrid_Source/discretisation.f90
+++ b/src/ZSOLVER/tensorproductmultigrid_Source/discretisation.f90
@@ -148,7 +148,7 @@ private
 
   real(kind=rl), pointer , contiguous , dimension(:) :: zt1d_discretisation
   INTEGER                                            :: zt1d_discretisation_size = 10000000
-  INTEGER , parameter                                :: zt1d_discretisation_size_factor = 2 * 4
+  INTEGER , parameter                                :: zt1d_discretisation_size_factor = 4 * 4
   INTEGER                                            :: nt1d_discretisation_top = 0 , nt1d_discretisation_max = 0
   
   INTEGER , ALLOCATABLE, DIMENSION (:)               :: nt1d_discretisation_pointer , nt1d_discretisation_pointer_ksize
@@ -182,6 +182,7 @@ public::DIRECTION_BACKWARD
 public::ORDERING_LEX
 public::ORDERING_RB
 public::zt1d_discretisation_init
+public::zt1d_discretisation_allocate3d
 
 contains
 
@@ -226,8 +227,8 @@ contains
     integer :: ksize, kbeg, kend
 
     if ( nt1d_discretisation_pointer_top == zt1d_discretisation_pointer_size ) then
-       print*,"ERROR zt1d_discretisation_init:: zt1d_discretisation_pointer_size to small=",zt1d_discretisation_pointer_size
-       print*,"ERROR zt1d_discretisation_init:: Augment it "
+       print*,"ERROR zt1d_discretisation_allocate3d:: zt1d_discretisation_pointer_size to small=",zt1d_discretisation_pointer_size
+       print*,"ERROR zt1d_discretisation_allocate3d:: Augment it "
        stop 
     end if
 
-- 
GitLab