Tek-Tips is the largest IT community on the Internet today!
Members share and learn making Tek-Tips Forums the best source of peer-reviewed technical information on the Internet!
-
Congratulations IamaSherpa on being selected by the Tek-Tips community for having the most helpful posts in the forums last week. Way to Go!
-
heat equation using crank-nicolsan scheme in fortran
The error statement is pretty clear. These lines: Real, allocatable :: x(:),u(:),a(:),b(:),c(:),d(:) Real:: m,dx,dt,Tmax Integer:: j,NI,JI have to appear before this line: call thomas(a,b,c,d,JI) I'm not sure why you're calling the Thomas solver before you even assemble a, b, c, and d... -
Help in coding - sub-diving area into circular cells.
If you want to have your domain covered by circles that don't overlap, but each touch the circles next to them, then stop thinking about circles for a minute. Instead, think of a hexagonal tessellation, e.g. Now, inscribe each of these circles with a circle. Out of all possible ways to... -
Test for intentionally undefined variables
If xwb's tip doesn't work with your compiler, you may have another option. Do you know anything about the possible values that X could be initialized to? For example, the logic of your code might only allow X to be a positive number. If this is the case, you could initialize X to a negative... -
Modifying a matriz of dynamic product
If you want to allow for the possibility that the matrix is not square, then you will need to change this: ! allocate space enogh of matriz ALLOCATE ( matriz(siz_matriz, & siz_matriz), & stat=situation_of_allocation ) to something like this: ALLOCATE ( matrix(number_of_rows... -
Problem with FORTRAN 90 code
If I remember correctly, a double pendulum is described by a system of two second-order non-linear ordinary differential equations. Because of the non-linear part, you'll need to resort to numerical techniques. There are many algorithms for solving ODEs. The easiest to code would probably be...
