blob: 0f30a3a30111a75695d0756899ba50f296fbf26d (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
|
* Eugeniy is grumpy on 2011/07/10:
Gosh, this is July 2011 and we are modifying Fortran 77 version to be Fortran 95.
Is there a sense?
* Eugeniy on 2011/08/23:
Navy chose very bad algorithm.
It seems to be a good idea to put normilezed dt=dz to avoid some unnesesary calculations
on each step. But results are devastating since code counter intuitively runs longer for
a shorter cell, and not surprisingly proportionally to time interval.
Execution time proportional to the total number of grid nodes which can be calculated as
Ngrid = tptz*zpts = (tmax-tmin)*Gamma_super / ((zmax*Gamma_super/(c*zpts)) *zpts
Ngrid = (tmax-tmin)*c/zmax * zpts^2
where 'c' is speed of light, 'zmax' is length of the cell, and 'zpts' is number
of desired points along z axis.
|
