|
From: David J. P. M. <djp...@an...> - 2013-06-28 22:43:10
|
I am using DISCUS 3.6.2 on Linux (Scientific Linux 6.3). I tried to run mod.3.mac from the Discus Cookbook macro/mod folder but keep getting the error "Too many points in reciprocal layer". This example wants to use 351*351 points. It would be useful in the future if I can use more pixels than this so looked into how to change it. I have tried to change the maximum number of pixels in the reciprocal space layer (MAXQXY) in config.inc and then recompile. I have even removed all previously compiled discus files and then recompiled. But still get the same error. Previously MAXQXY=301*301 I've set this to MAXQXY=501*501 When I use the Linux command grep "MAXQXY=" * to find if there is another file that could be overwriting my setting of MAXQXY the only files that define this value are config.inc, config,reinhard and config.thomas. Neither of the last two are called in fourier.f. I have tried entering 501*501 directly into line 421 of fourier.f if(inc(1)*inc(2).le.MAXQXY) then then recompiling but that still only gave 90601 when I typed "show config". I compile by creating the DiffuseBuild folder then running INST_unix which I have not changed from the downloaded version. My fortran compiler is gfortran. Is there another trick to getting my value of MAXQXY to be accepted as the maximum number of pixels in the Q-layer? Best wishes, Jonathan |