From b58ab983ee60662cdd3a1605049ac1b51b136b14 Mon Sep 17 00:00:00 2001 From: Kento Kaneko Date: Fri, 27 Sep 2024 12:56:48 -0400 Subject: [PATCH] use lcs for coefficient arrays --- code/MOR | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/code/MOR b/code/MOR index 44ee0ec4..5362560c 100644 --- a/code/MOR +++ b/code/MOR @@ -20,8 +20,8 @@ c parameters and common blocks for MOR $ pic(lx2*ly2*lz2*lelm), $ tic(lx1*ly1*lz1*lelm) - common /morcoef/ uk((lub+1)*(lsu+1)),tk((ltb+1)*(lst+1)), - $ timek(lcs),edk((lub+1)*(lsu+1)) + common /morcoef/ uk((lub+1)*(lcs+1)),tk((ltb+1)*(lcs+1)), + $ timek(lcs),edk((lub+1)*(lcs+1)) common /morcoef2/ ukp((lub+1)*lsu),tkp((ltb+1)*lst) common /morusnap/ us0(lx1*ly1*lz1*lelm,ldim,lsu) @@ -105,7 +105,7 @@ c parameters and common blocks for MOR common /moritmp/ itmp1(lub),itmp2(lub) common /morrtmp/ rtmp1(lub,lub),rtmp2(lb,lb),rtmp3(lb,lb) - common /morstmp/ stmp(lb+1,ls) + common /morstmp/ stmp(lb+1,max(ls,lcs)) common /morfld/ wm1(lx1*ly1*lz1*lelm) common /morqoi/ rnus,rdgx(0:lub),rdgy(0:lub),rdgz(0:lub)