From josef.stoeckl at uibk.ac.at Mon Jun 2 03:22:00 2008 From: josef.stoeckl at uibk.ac.at (=?ISO-8859-1?Q?Josef_St=F6ckl?=) Date: Mon, 02 Jun 2008 10:22:00 +0200 Subject: [FLASH-BUGS] Restart Bug with HDF5 Message-ID: <4843ADA8.6070306@uibk.ac.at> Hello, I have noticed and tracked down a restart bug (segmentation fault) in Flash 3 with the serial HDF5 IO Unit, when a 3D problem gets restarted. Basically the wrong buffer variable (unk instead of faceZBuf) gets used, which at that time is not allocated. Below you find the fix to this in a patch-like manner. Best regards, Josef IO/IOMain/hdf5/serial/PM/io_readData.F90: if(NDIM .gt. 2) then allocate(faceZBuf(NUNK_VARS, NXB, NYB, NZB+1, localNumBlocks)) - call MPI_RECV(unk(i,:,:,:,1:localNumBlocks), & + call MPI_RECV(faceZBuf(i,:,:,:,1:localNumBlocks), & NXB*NYB*(NZB+1)*localNumBlocks, & FLASH_REAL, MASTER_PE, & 9+i+NUNK_VARS+(NFACE_VARS*2), & MPI_COMM_WORLD, status, ierr) facevarz(i,io_ilo:io_ihi, io_jlo:io_jhi, io_klo:io_khi+1,1:localNumBlocks) faceZBuf(i,1:NXB,1:NYB,1:NZB+1,1:localNumBlocks) deallocate(faceZBuf) end if From josef.stoeckl at uibk.ac.at Mon Jun 2 08:41:08 2008 From: josef.stoeckl at uibk.ac.at (=?ISO-8859-1?Q?Josef_St=F6ckl?=) Date: Mon, 02 Jun 2008 15:41:08 +0200 Subject: [FLASH-BUGS] A Cosmology restart Bug Message-ID: <4843F874.2060304@uibk.ac.at> Hi again, I've fixed another restart bug which had to do with Cosmology not being initialized before Cosmology_getRedshift was called in IO_init. This led to a division by zero that caused io_nextCheckpointZ to be infinity (with PGI 7.1), which in turn led to the IO_output routine being stuck in a while loop. To quickly fix this, I have added a "call Cosmology_init(myPE, io_restart)" to IO_init.F90 around line 318 (before the Cosmology_getRedshift call). Additionally the calls to Cosmology_init could be omitted from Driver_initFlash, although a better solution would probably be to let the Cosmology module handly the io_*Z variables. Best regards, Josef From absefko at sandia.gov Thu Jun 19 15:04:20 2008 From: absefko at sandia.gov (A. B. Sefkow) Date: Thu, 19 Jun 2008 14:04:20 -0600 Subject: [FLASH-BUGS] HDF5 issue at compile Message-ID: <485ABBC4.8060001@sandia.gov> Dear Flash Bugs, I possibly uncovered a pseudo bug in the Flash3 distribution while trying to install it on Mac OS X 10.4. In the Prototypes/Darwin/Makefile.h.gfortran file, which I used as my template, there is an extra space in LIB_HDF5 between -L and ${HDF5_PATH}/lib. That error and the -static option in LFLAGS_OPT caused the creation of the Flash3 executable to not be made because of errors in executing /usr/bin/ld after "Linking in Units and Libraries" near the end of the make output. At first, I tried different versions of HDF5 (binaries and trying to compile from source code) to no avail, and only now (apparently) succeeded by disabling static linking (and enabling fortran support) when configuring/compiling HDF5-1.6.7. After a bit of Googling, Apple says Mac OS X doesn't support static linking of binaries: http://developer.apple.com/qa/qa2001/qa1118.html And only after these errors could I get the flash3 executable, although I haven't yet run it with the Sedov test problem (that's next). I do still have this warning, but get a "success": > ate.o tmr_etime.o tmr_findTimerIndex.o tmr_getMaxCallStackDepth.o tmr_getMaxTimerParents.o tmr_init.o tmr_lookupIndex.o tmr_stackLib.o tree.o umap.o user_coord_transfm.o workspace.o -L/usr/local/hdf5/hdf5-1.6.7/lib -lhdf5 -lz > /usr/bin/ld: warning can't open dynamic library: /libgcc_s.1.dylib referenced from: /usr/local/lib/gcc/i386-apple-darwin8.10.1/4.3.0/../../../libgfortran.dylib (checking for undefined symbols may be affected) (No such file or directory, errno = 2) > SUCCESS Since I'm a relative newbie to both coding and Mac OS X, it's been frustrating trying to figure out the configure/compile for C, Fortran, MPI, and HDF5 required to get flash3 running, but I must compliment the Flash people for making the manual quite a bit more coherent and helpful than most others. Thanks, Adam Sefkow Sandia National Labs From absefko at sandia.gov Thu Jun 26 15:52:29 2008 From: absefko at sandia.gov (A. B. Sefkow) Date: Thu, 26 Jun 2008 14:52:29 -0600 Subject: [FLASH-BUGS] [FLASH-USERS] MRT question In-Reply-To: <5b942a610806201652r1c145c7ex56f893c67daacdf7@mail.gmail.com> References: <485C2866.4060807@sandia.gov> <5b942a610806201652r1c145c7ex56f893c67daacdf7@mail.gmail.com> Message-ID: <4864018D.5060501@sandia.gov> Hi Flash-bugs, It is possible I found another bug. I'm simply going through the provided test problems and, after a "./setup Cellular -2d -auto", "cd object", and "make", the code compiles pretty far until the bn_burner.F90 file fails: ------------------ ... ... ... mpif90 -I../lib/ma28/include -c -O2 -fdefault-real-8 -fdefault-double-8 -ffree-line-length-none -DMAXBLOCKS=1000 -DNXB=8 -DNYB=8 -DNZB=1 -DN_DIM=2 avisco.F90 mpif90 -I../lib/ma28/include -c -O2 -fdefault-real-8 -fdefault-double-8 -ffree-line-length-none -DMAXBLOCKS=1000 -DNXB=8 -DNYB=8 -DNZB=1 -DN_DIM=2 bn_azbar.F90 mpif90 -I../lib/ma28/include -c -O2 -fdefault-real-8 -fdefault-double-8 -ffree-line-length-none -DMAXBLOCKS=1000 -DNXB=8 -DNYB=8 -DNZB=1 -DN_DIM=2 bn_burner.F90 bn_burner.F90:164.86: & bn_network,bn_networkSparseJakob,bn_networkSparsePointers,bn_b 1 Error: Type/rank mismatch in argument 'steper' at (1) bn_burner.F90:173.85: & bn_network,bn_networkDenseJakob,bn_networkSparsePointers,bn_b 1 Error: Type/rank mismatch in argument 'steper' at (1) bn_burner.F90:187.86: & bn_network,bn_networkSparseJakob,bn_networkSparsePointers,bn_r 1 Error: Type/rank mismatch in argument 'steper' at (1) bn_burner.F90:195.85: & bn_network,bn_networkDenseJakob,bn_networkSparsePointers,bn_r 1 Error: Type/rank mismatch in argument 'steper' at (1) make: *** [bn_burner.o] Error 1 ------------------------ I tried looking for an easy fix by myself, but with no luck. I did notice that the parameters for the odestepper and algebra are missing the prefix "bn_" in the Config file (while "bn_" is used in the code), but that didn't fix the problem. Has anyone else reported that the Cellular test problem did work for them? Thanks, Adam Anshu Dubey wrote: > Dear Adam, > > First, thanks for reporting the bug in the Makefile. We have made a note of > it. About RT setup, I believe Dongwook(dongwook at flash.uchicago.edu > ) > will be able to help you, but currently > he is on vacation, so it may be a week or two before you get a response > from him > > Anshu > > > On Fri, Jun 20, 2008 at 5:00 PM, A. B. Sefkow > wrote: > > Dear Flash Users, > > I am a new Flash user and became interested in the code in order to > attempt a Magnetic Rayleigh Taylor (MRT) simulation. I'm still > going through the provided examples and the manual in order to learn > how to use the code, but so far it isn't obvious to me how one can > specify magnetic fields and pressure *beyond* the initialization. > In other words, would I need a circuit model to do this in Flash, > or is there some way to tell the code how I want I(t) > [Bx(t),By(t),Bz(t)] to evolve (say, at a boundary, and let it > diffuse in)? Along those lines, can nonideal MHD parameters be set > for different fluids? I essentially want to start really simple, > say by just exchanging gravity for magnetics in a straight forward > RT simulation. > > If anybody has a simple RT or MRT setup they would be willing to > share in order to help a new user learn, I would be much indebted. > > Thanks, > Adam > > From cdaley at flash.uchicago.edu Thu Jun 26 16:14:02 2008 From: cdaley at flash.uchicago.edu (Chris Daley) Date: Thu, 26 Jun 2008 16:14:02 -0500 Subject: [FLASH-BUGS] [FLASH-USERS] MRT question In-Reply-To: <4864018D.5060501@sandia.gov> References: <485C2866.4060807@sandia.gov> <5b942a610806201652r1c145c7ex56f893c67daacdf7@mail.gmail.com> <4864018D.5060501@sandia.gov> Message-ID: <4864069A.3030101@flash.uchicago.edu> A. B. Sefkow wrote: > Hi Flash-bugs, > > It is possible I found another bug. I'm simply going through the > provided test problems and, after a "./setup Cellular -2d -auto", "cd > object", and "make", the code compiles pretty far until the > bn_burner.F90 file fails: > > > ------------------ > ... > ... > ... > mpif90 -I../lib/ma28/include -c -O2 -fdefault-real-8 > -fdefault-double-8 -ffree-line-length-none -DMAXBLOCKS=1000 -DNXB=8 > -DNYB=8 -DNZB=1 -DN_DIM=2 avisco.F90 > mpif90 -I../lib/ma28/include -c -O2 -fdefault-real-8 > -fdefault-double-8 -ffree-line-length-none -DMAXBLOCKS=1000 -DNXB=8 > -DNYB=8 -DNZB=1 -DN_DIM=2 bn_azbar.F90 > mpif90 -I../lib/ma28/include -c -O2 -fdefault-real-8 > -fdefault-double-8 -ffree-line-length-none -DMAXBLOCKS=1000 -DNXB=8 > -DNYB=8 -DNZB=1 -DN_DIM=2 bn_burner.F90 > bn_burner.F90:164.86: > > & bn_network,bn_networkSparseJakob,bn_networkSparsePointers,bn_b > > 1 > Error: Type/rank mismatch in argument 'steper' at (1) > bn_burner.F90:173.85: > > & bn_network,bn_networkDenseJakob,bn_networkSparsePointers,bn_b > > 1 > Error: Type/rank mismatch in argument 'steper' at (1) > bn_burner.F90:187.86: > > & bn_network,bn_networkSparseJakob,bn_networkSparsePointers,bn_r > > 1 > Error: Type/rank mismatch in argument 'steper' at (1) > bn_burner.F90:195.85: > > & bn_network,bn_networkDenseJakob,bn_networkSparsePointers,bn_r > > 1 > Error: Type/rank mismatch in argument 'steper' at (1) > make: *** [bn_burner.o] Error 1 > > ------------------------ > > > I tried looking for an easy fix by myself, but with no luck. I did > notice that the parameters for the odestepper and algebra are missing > the prefix "bn_" in the Config file (while "bn_" is used in the code), > but that didn't fix the problem. > > Has anyone else reported that the Cellular test problem did work for > them? > > Thanks, > Adam > > > > > > Anshu Dubey wrote: >> Dear Adam, >> >> First, thanks for reporting the bug in the Makefile. We have made a >> note of >> it. About RT setup, I believe Dongwook(dongwook at flash.uchicago.edu >> ) >> will be able to help you, but currently >> he is on vacation, so it may be a week or two before you get a response >> from him >> >> Anshu >> >> >> On Fri, Jun 20, 2008 at 5:00 PM, A. B. Sefkow > > wrote: >> >> Dear Flash Users, >> >> I am a new Flash user and became interested in the code in order to >> attempt a Magnetic Rayleigh Taylor (MRT) simulation. I'm still >> going through the provided examples and the manual in order to learn >> how to use the code, but so far it isn't obvious to me how one can >> specify magnetic fields and pressure *beyond* the initialization. >> In other words, would I need a circuit model to do this in Flash, >> or is there some way to tell the code how I want I(t) >> [Bx(t),By(t),Bz(t)] to evolve (say, at a boundary, and let it >> diffuse in)? Along those lines, can nonideal MHD parameters be set >> for different fluids? I essentially want to start really simple, >> say by just exchanging gravity for magnetics in a straight forward >> RT simulation. >> >> If anybody has a simple RT or MRT setup they would be willing to >> share in order to help a new user learn, I would be much indebted. >> >> Thanks, >> Adam >> >> Hi Adam, I was unable to reproduce this error. However, your compilation error message seems to imply that a single line character limit is exceeded: e.g. line 164 in bn_burner.f90 & bn_network,bn_networkSparseJakob,bn_networkSparsePointers,bn_b should read: & bn_network,bn_networkSparseJakob,bn_networkSparsePointers,bn_baderMa28) I suspect the problem can be cured by simply splitting each failing line into two separate lines. Let me know if this helps. Regards, Chris From absefko at sandia.gov Thu Jun 26 16:34:00 2008 From: absefko at sandia.gov (A. B. Sefkow) Date: Thu, 26 Jun 2008 15:34:00 -0600 Subject: [FLASH-BUGS] [FLASH-USERS] MRT question In-Reply-To: <4864069A.3030101@flash.uchicago.edu> References: <485C2866.4060807@sandia.gov> <5b942a610806201652r1c145c7ex56f893c67daacdf7@mail.gmail.com> <4864018D.5060501@sandia.gov> <4864069A.3030101@flash.uchicago.edu> Message-ID: <48640B48.6010105@sandia.gov> Chris Daley wrote: > A. B. Sefkow wrote: >> Hi Flash-bugs, >> >> It is possible I found another bug. I'm simply going through the >> provided test problems and, after a "./setup Cellular -2d -auto", "cd >> object", and "make", the code compiles pretty far until the >> bn_burner.F90 file fails: >> >> >> ------------------ >> ... >> ... >> ... >> mpif90 -I../lib/ma28/include -c -O2 -fdefault-real-8 >> -fdefault-double-8 -ffree-line-length-none -DMAXBLOCKS=1000 -DNXB=8 >> -DNYB=8 -DNZB=1 -DN_DIM=2 avisco.F90 >> mpif90 -I../lib/ma28/include -c -O2 -fdefault-real-8 >> -fdefault-double-8 -ffree-line-length-none -DMAXBLOCKS=1000 -DNXB=8 >> -DNYB=8 -DNZB=1 -DN_DIM=2 bn_azbar.F90 >> mpif90 -I../lib/ma28/include -c -O2 -fdefault-real-8 >> -fdefault-double-8 -ffree-line-length-none -DMAXBLOCKS=1000 -DNXB=8 >> -DNYB=8 -DNZB=1 -DN_DIM=2 bn_burner.F90 >> bn_burner.F90:164.86: >> >> & bn_network,bn_networkSparseJakob,bn_networkSparsePointers,bn_b >> >> 1 >> Error: Type/rank mismatch in argument 'steper' at (1) >> bn_burner.F90:173.85: >> >> & bn_network,bn_networkDenseJakob,bn_networkSparsePointers,bn_b >> >> 1 >> Error: Type/rank mismatch in argument 'steper' at (1) >> bn_burner.F90:187.86: >> >> & bn_network,bn_networkSparseJakob,bn_networkSparsePointers,bn_r >> >> 1 >> Error: Type/rank mismatch in argument 'steper' at (1) >> bn_burner.F90:195.85: >> >> & bn_network,bn_networkDenseJakob,bn_networkSparsePointers,bn_r >> >> 1 >> Error: Type/rank mismatch in argument 'steper' at (1) >> make: *** [bn_burner.o] Error 1 >> >> ------------------------ >> >> >> I tried looking for an easy fix by myself, but with no luck. I did >> notice that the parameters for the odestepper and algebra are missing >> the prefix "bn_" in the Config file (while "bn_" is used in the code), >> but that didn't fix the problem. >> >> Has anyone else reported that the Cellular test problem did work for >> them? >> >> Thanks, >> Adam >> >> >> >> >> >> Anshu Dubey wrote: >>> Dear Adam, >>> >>> First, thanks for reporting the bug in the Makefile. We have made a >>> note of >>> it. About RT setup, I believe Dongwook(dongwook at flash.uchicago.edu >>> ) >>> will be able to help you, but currently >>> he is on vacation, so it may be a week or two before you get a response >>> from him >>> >>> Anshu >>> >>> >>> On Fri, Jun 20, 2008 at 5:00 PM, A. B. Sefkow >> > wrote: >>> >>> Dear Flash Users, >>> >>> I am a new Flash user and became interested in the code in order to >>> attempt a Magnetic Rayleigh Taylor (MRT) simulation. I'm still >>> going through the provided examples and the manual in order to learn >>> how to use the code, but so far it isn't obvious to me how one can >>> specify magnetic fields and pressure *beyond* the initialization. >>> In other words, would I need a circuit model to do this in Flash, >>> or is there some way to tell the code how I want I(t) >>> [Bx(t),By(t),Bz(t)] to evolve (say, at a boundary, and let it >>> diffuse in)? Along those lines, can nonideal MHD parameters be set >>> for different fluids? I essentially want to start really simple, >>> say by just exchanging gravity for magnetics in a straight forward >>> RT simulation. >>> >>> If anybody has a simple RT or MRT setup they would be willing to >>> share in order to help a new user learn, I would be much indebted. >>> >>> Thanks, >>> Adam >>> >>> > Hi Adam, > > I was unable to reproduce this error. However, your compilation error > message seems to imply that a single line character limit is exceeded: > > e.g. line 164 in bn_burner.f90 > & bn_network,bn_networkSparseJakob,bn_networkSparsePointers,bn_b > should read: > & bn_network,bn_networkSparseJakob,bn_networkSparsePointers,bn_baderMa28) > > I suspect the problem can be cured by simply splitting each failing line > into two > separate lines. Let me know if this helps. > > Regards, > Chris > > Hi Chris, Thanks for the note. Unfortunately I still get the same error, but now with the full line and closing ) ------------------------- bn_burner.F90:165.31: & bn_networkSparsePointers,bn_baderMa28) 1 Error: Type/rank mismatch in argument 'steper' at (1) bn_burner.F90:175.44: & bn_networkSparsePointers,bn_baderGift) 1 Error: Type/rank mismatch in argument 'steper' at (1) bn_burner.F90:190.44: & bn_networkSparsePointers,bn_rosenMa28) 1 Error: Type/rank mismatch in argument 'steper' at (1) bn_burner.F90:199.44: & bn_networkSparsePointers,bn_rosenGift) 1 Error: Type/rank mismatch in argument 'steper' at (1) make: *** [bn_burner.o] Error 1 ----------------------------------- Any other ideas? I don't suppose it would work if I got a fresh and functional copy of that file from someone else? Also, why wouldn't it matter that bn_burner.F90 is using bn_algebra and bn_odeStepper in its code, but the flash.par file specifies odeStepper and algebra (without the bn_)... a tiny bug? I guess there is a default setting that gets used probably. Thanks, Adam From cdaley at flash.uchicago.edu Fri Jun 27 11:53:58 2008 From: cdaley at flash.uchicago.edu (Chris Daley) Date: Fri, 27 Jun 2008 11:53:58 -0500 Subject: [FLASH-BUGS] [FLASH-USERS] MRT question In-Reply-To: <48640B48.6010105@sandia.gov> References: <485C2866.4060807@sandia.gov> <5b942a610806201652r1c145c7ex56f893c67daacdf7@mail.gmail.com> <4864018D.5060501@sandia.gov> <4864069A.3030101@flash.uchicago.edu> <48640B48.6010105@sandia.gov> Message-ID: <48651B26.3070402@flash.uchicago.edu> A. B. Sefkow wrote: > > > Chris Daley wrote: >> A. B. Sefkow wrote: >>> Hi Flash-bugs, >>> >>> It is possible I found another bug. I'm simply going through the >>> provided test problems and, after a "./setup Cellular -2d -auto", "cd >>> object", and "make", the code compiles pretty far until the >>> bn_burner.F90 file fails: >>> >>> >>> ------------------ >>> ... >>> ... >>> ... >>> mpif90 -I../lib/ma28/include -c -O2 -fdefault-real-8 >>> -fdefault-double-8 -ffree-line-length-none -DMAXBLOCKS=1000 -DNXB=8 >>> -DNYB=8 -DNZB=1 -DN_DIM=2 avisco.F90 >>> mpif90 -I../lib/ma28/include -c -O2 -fdefault-real-8 >>> -fdefault-double-8 -ffree-line-length-none -DMAXBLOCKS=1000 -DNXB=8 >>> -DNYB=8 -DNZB=1 -DN_DIM=2 bn_azbar.F90 >>> mpif90 -I../lib/ma28/include -c -O2 -fdefault-real-8 >>> -fdefault-double-8 -ffree-line-length-none -DMAXBLOCKS=1000 -DNXB=8 >>> -DNYB=8 -DNZB=1 -DN_DIM=2 bn_burner.F90 >>> bn_burner.F90:164.86: >>> >>> & bn_network,bn_networkSparseJakob,bn_networkSparsePointers,bn_b >>> >>> 1 >>> Error: Type/rank mismatch in argument 'steper' at (1) >>> bn_burner.F90:173.85: >>> >>> & bn_network,bn_networkDenseJakob,bn_networkSparsePointers,bn_b >>> >>> 1 >>> Error: Type/rank mismatch in argument 'steper' at (1) >>> bn_burner.F90:187.86: >>> >>> & bn_network,bn_networkSparseJakob,bn_networkSparsePointers,bn_r >>> >>> 1 >>> Error: Type/rank mismatch in argument 'steper' at (1) >>> bn_burner.F90:195.85: >>> >>> & bn_network,bn_networkDenseJakob,bn_networkSparsePointers,bn_r >>> >>> 1 >>> Error: Type/rank mismatch in argument 'steper' at (1) >>> make: *** [bn_burner.o] Error 1 >>> >>> ------------------------ >>> >>> >>> I tried looking for an easy fix by myself, but with no luck. I did >>> notice that the parameters for the odestepper and algebra are missing >>> the prefix "bn_" in the Config file (while "bn_" is used in the code), >>> but that didn't fix the problem. >>> >>> Has anyone else reported that the Cellular test problem did work for >>> them? >>> >>> Thanks, >>> Adam >>> >>> >>> >>> >>> >>> Anshu Dubey wrote: >>>> Dear Adam, >>>> >>>> First, thanks for reporting the bug in the Makefile. We have made a >>>> note of >>>> it. About RT setup, I believe Dongwook(dongwook at flash.uchicago.edu >>>> ) >>>> will be able to help you, but currently >>>> he is on vacation, so it may be a week or two before you get a >>>> response >>>> from him >>>> >>>> Anshu >>>> >>>> >>>> On Fri, Jun 20, 2008 at 5:00 PM, A. B. Sefkow >>> > wrote: >>>> >>>> Dear Flash Users, >>>> >>>> I am a new Flash user and became interested in the code in >>>> order to >>>> attempt a Magnetic Rayleigh Taylor (MRT) simulation. I'm still >>>> going through the provided examples and the manual in order to >>>> learn >>>> how to use the code, but so far it isn't obvious to me how one can >>>> specify magnetic fields and pressure *beyond* the initialization. >>>> In other words, would I need a circuit model to do this in Flash, >>>> or is there some way to tell the code how I want I(t) >>>> [Bx(t),By(t),Bz(t)] to evolve (say, at a boundary, and let it >>>> diffuse in)? Along those lines, can nonideal MHD parameters be >>>> set >>>> for different fluids? I essentially want to start really simple, >>>> say by just exchanging gravity for magnetics in a straight forward >>>> RT simulation. >>>> >>>> If anybody has a simple RT or MRT setup they would be willing to >>>> share in order to help a new user learn, I would be much indebted. >>>> >>>> Thanks, >>>> Adam >>>> >>>> >> Hi Adam, >> >> I was unable to reproduce this error. However, your compilation error >> message seems to imply that a single line character limit is exceeded: >> >> e.g. line 164 in bn_burner.f90 >> & bn_network,bn_networkSparseJakob,bn_networkSparsePointers,bn_b >> should read: >> & >> bn_network,bn_networkSparseJakob,bn_networkSparsePointers,bn_baderMa28) >> >> I suspect the problem can be cured by simply splitting each failing line >> into two >> separate lines. Let me know if this helps. >> >> Regards, >> Chris >> >> > > Hi Chris, > > Thanks for the note. Unfortunately I still get the same error, but > now with the full line and closing ) > > ------------------------- > bn_burner.F90:165.31: > > & bn_networkSparsePointers,bn_baderMa28) > 1 > Error: Type/rank mismatch in argument 'steper' at (1) > bn_burner.F90:175.44: > > & bn_networkSparsePointers,bn_baderGift) > 1 > Error: Type/rank mismatch in argument 'steper' at (1) > bn_burner.F90:190.44: > > & bn_networkSparsePointers,bn_rosenMa28) > 1 > Error: Type/rank mismatch in argument 'steper' at (1) > bn_burner.F90:199.44: > > & bn_networkSparsePointers,bn_rosenGift) > 1 > Error: Type/rank mismatch in argument 'steper' at (1) > make: *** [bn_burner.o] Error 1 > ----------------------------------- > > Any other ideas? I don't suppose it would work if I got a fresh and > functional copy of that file from someone else? > > Also, why wouldn't it matter that bn_burner.F90 is using bn_algebra > and bn_odeStepper in its code, but the flash.par file specifies > odeStepper and algebra (without the bn_)... a tiny bug? I guess there > is a default setting that gets used probably. > > Thanks, > Adam > > > > > Hi Adam, I get the same compile time error when using the gfortran Fortran compiler. It does not appear straightforward to resolve cleanly. I will mention this to my work colleague who knows more about the burn unit. If you wish to compile and run the simulation, one workaround is to remove the type/rank checking: In bn_burner.F90: 1. Comment out - use bnIntegrate_interface, ONLY: bn_baderMa28, bn_baderGift, & bn_rosenMa28, bn_rosenGift 2. Uncomment - external bn_baderMa28, bn_baderGift, bn_rosenMa28, bn_rosenGift This is not an ideal solution, but the program should work the same in this particular case, as there are no assumed shape array arguments or optional arguments in the subroutine calls. > Also, why wouldn't it matter that bn_burner.F90 is using bn_algebra and bn_odeStepper in its code, but the flash.par file specifies odeStepper and algebra (without the bn_)... a tiny bug? I guess there is a default setting that gets used probably. These parameters are stored in FLASH variables in Burn_init.F90: call RuntimeParameters_get("odeStepper", bn_odeStepper ) call RuntimeParameters_get("algebra", bn_algebra ) Regards, Chris