72 |
|
|
73 |
rlamda = iim / pi / colat0 * grossismx / sqrt(- eignvl(2: iim)) |
rlamda = iim / pi / colat0 * grossismx / sqrt(- eignvl(2: iim)) |
74 |
print *, "1 / rlamda(iim) = ", 1. / rlamda(iim) |
print *, "1 / rlamda(iim) = ", 1. / rlamda(iim) |
|
|
|
75 |
! This is demonstrated in the notes but just to be sure: |
! This is demonstrated in the notes but just to be sure: |
76 |
call assert(rlamda(iim) * colat0 >= 1. - epsilon(0.), & |
call assert(rlamda(iim) * colat0 >= 1. - 2. * epsilon(0.), & |
77 |
"inifilr rlamda(iim) * colat0") |
"inifilr rlamda(iim) * colat0") |
78 |
|
|
79 |
call new_unit(unit) |
call new_unit(unit) |