forums.silverfrost.com

pashupati · Joined: 16 Aug 2012 Posts: 8

I can only get first data point after that it shows error.

Please help

PROGRAM TDGL

IMPLICIT NONE

!c **************** Initialization of variables *********************
INTEGER, PARAMETER :: IMAX=100,JMAX=100
INTEGER Nx ! number of cells in x
INTEGER Ny ! number of cells in x

REAL Ax ! mesh spacing in x
REAL Ay ! mesh spacing in y

REAL dt ! time step
REAL k ! kappa ( = lambda / xi )
REAL Eta ! eta
REAL T ! temperature ( scaled with Tc )
REAL He ! applied external magnetic field
REAL Ho ! initial applied external magnetic field
REAL Hi ! magnetic field in the interior hole
REAL dHe ! applied external magnetic field increment
REAL Ruido ! Standard deviation of random noise
REAL Eo ! Parameter for random noise
INTEGER TotSteps! number of time steps
INTEGER Step ! step counter
INTEGER OutSteps! number of steps between outputs
INTEGER OutStp2 ! number of steps between run-time information
INTEGER HeSteps ! number of steps between applied field increments

INTEGER idum ! auxiliary
INTEGER i,j ! auxiliary
CHARACTER*36 Label ! auxiliary
INTEGER iHe ! auxiliary for naming restart files

COMPLEX F(IMAX+1,JMAX+1) ! order parameter
COMPLEX Ux(IMAX,JMAX+1) ! link variable in x
COMPLEX Uy(IMAX+1,JMAX) ! link variable in y
COMPLEX bloop(IMAX,JMAX) ! exp(-i a_x a_y B_z)
!c The following three arrays are auxiliary to avoid copies
COMPLEX G(IMAX+1,JMAX+1) ! order parameter
COMPLEX Vx(IMAX,JMAX+1) ! link variable in x
COMPLEX Vy(IMAX+1,JMAX) ! link variable in y
!c
LOGICAL Bulk(IMAX,JMAX) ! bulk cell indicator
INTEGER iswnod (IMAX+1,JMAX+1) ! vertex indicator
INTEGER iswlx (IMAX,JMAX+1) ! x-link indicator
INTEGER iswly (IMAX+1,JMAX) ! y-link indicator

REAL rAx2 ! auxiliar
REAL rAy2 ! auxiliar
REAL rAxAy ! auxiliar
REAL k2 ! auxiliar
REAL rEta ! auxiliar
REAL T1 ! auxiliar
REAL rT1 ! auxiliar
REAL k2rT12 ! auxiliar
REAL AxAy ! auxiliar
REAL AxAyHe ! auxiliar
REAL AxAyHi ! auxiliar
REAL AxAydHi ! auxiliar
COMPLEX EiAxAyHe! auxiliar, approx. exp(-i*ax*ay*He)
COMPLEX EiAxAyHi! auxiliar, approx. exp(-i*ax*ay*Hi)
REAL rotA ! auxiliar
COMPLEX EiMAxAydHi ! var. aux., approx. exp(-i*M*ax*ay*dHi)
COMPLEX EiMAxAyHi ! var. aux., approx. exp(-i*M*ax*ay*dHi)
REAL FFF ! auxiliar

REAL FreeEnergy ! !
REAL Bz ! mean field
REAL Magnetization ! !

COMPLEX GaussDev! complex numbers generator (Gaussian distrib.)

INTEGER d ! auxiliar
REAL Flux ! total fluxoid
REAL Flux1 ! internal fluxoid
REAL PATHx1 !
REAL PATHx2 !
REAL PATHy1 !
REAL PATHy2 ! circulations
COMPLEX PTHx1 !
COMPLEX PTHx2 !
COMPLEX PTHy1 !
COMPLEX PTHy2 ! circulations
REAL MinF !
REAL MinF1 !

INTEGER N ! Hole N boundary
INTEGER S ! Hole S boundary
INTEGER E ! Hole E boundary
INTEGER W ! Hole W boundary
INTEGER M ! auxiliar
REAL rM ! auxiliar

integer iaux

!c ******************************************************************

iHe = 9000000

!c ********************** Input data **************************

OPEN (8,FILE='tdgl.in

jjgermis · Posted: Wed Aug 22, 2012 9:52 am Post subject:

Thanks for posting some code - this is alway helps to get down to the problem. However, since the number of code lines posted in the forum is limited, you should try to reduce the code so that one can at least reproduce it.

I always like cut, paste and compile examples Smile

Give it try - and you will get the help you are looking for.

LitusSaxonicum · Posted: Wed Aug 22, 2012 11:36 am Post subject:

If the code is long, you could always upload it to one of the (legal) file sharing websites and provide a link to it. I often use DropBox (www.dropbox.com) to share files that are too long even for e mail attachments or are of types like .EXE that many e mail systems won't accept.

Eddie

jjgermis · Posted: Wed Aug 22, 2012 11:45 am Post subject:

Hi Eddie,

thanks for the DropBox tipp. It seems like one need to install some software before one can use it. We for example are not allowed to install software that is not approved - difficult situation. But on the otherhand our IT always try to avoid the worst case scenario.

pashupati · Joined: 16 Aug 2012 Posts: 8

Anyone willing to help me to take a look at the code?
Please help

my email :dklpashupati at yahoo.com

jjgermis · Posted: Thu Aug 23, 2012 3:16 pm Post subject:

We are willing to help if you post a working example of your code. The part of the program given above is only variable declaration. It is not fair to expect from anybody to figure out what is happening with such little information.

The only "real" code is the line where a file is opened. So for a start you must at least specify:
1.) are you able to read the file?
2.) occur the error before/after any calls to subroutines or functions
3.) what sort of files are you reading.

You need to reduce you problem. Once you have done that we can help you. Give it a try Wink

pashupati · Joined: 16 Aug 2012 Posts: 8

Sorry about the code, I didn't realize that only a part of the code is posted Sad

.
I am reading .ini files with the input parameters. After I start "run", it reads the parameter and generate the first data point. After that it gives "floating point overflow". floating point co-processor fault at address 0040a7e5 in the line 917". I am running this in window 7 32 bit laptop.

I used a gfortan in mac computer the same code works fine. I am not planning to buy Mac. Crying or Very sad

jjgermis · Posted: Thu Aug 23, 2012 3:49 pm Post subject:

You have to post something like this:

pashupati · Joined: 16 Aug 2012 Posts: 8

pashupati · Joined: 16 Aug 2012 Posts: 8

LitusSaxonicum · Posted: Thu Aug 23, 2012 6:42 pm Post subject:

Jaques,

Drop Box is only one of several legal file sharing sites. You only need to instal software to post to it, as reading from it is a simple internet download. You get 2 Gb free, but if you buy a licence, you can have much more space, and then your IT section could add it to the list of accepted applications (or whatever system you chose would probably work the same). In my experience, IT organisations distrust free software, but they are happy to pay for stuff you can get free - and then they won't pay for the updates!

I found it useful for sharing EXE files.

Eddie

pashupati · Joined: 16 Aug 2012 Posts: 8

Code is very long
May I email to someone willing to help me?

email: dklpashupati at yahoo.com

Thank you !

jjgermis · Posted: Fri Aug 24, 2012 7:51 am Post subject:

Seeing some code is a help - we are getting there Smile

However, where does the error occur. Can you isolate it:
1.) are you able to read the file?
2.) occur the error before/after any calls to subroutines or functions
3.) what sort of files are you reading.

pashupati · Joined: 16 Aug 2012 Posts: 8

I can read the input file and it gives the first data point (valid point) and it terminate due floating point error

DanRRight · Posted: Sat Aug 25, 2012 6:34 pm Post subject:

In the DOS prompt window (better use for that the Total Commander - to me it's #1 program on all computers and the best and the most useful app of all times) compile it like this

FTN95 yourfile.for /checkmate /link

Yourfile.for is the name of your file of course. If this is free form source add
/free to this line

Then run the debugger
>SDBG yourfile.exe

You will find all your errors in seconds. If not ask here again what's specifically wrong. Post here how many errors Apple compiler missed