6 |
|
|
7 |
private llm, hybrid, funcd, y, ya, compute_ab |
private llm, hybrid, funcd, y, ya, compute_ab |
8 |
|
|
9 |
real, save:: ap(llm+1), pa ! in Pa |
real, save:: ap(llm+1) ! in Pa |
10 |
real, save:: bp(llm+1) |
real, save:: bp(llm+1) |
11 |
|
|
12 |
REAL s(llm+1) |
REAL s(llm+1) |
30 |
! variables "ap", "bp", "presnivs". "pa" should be defined before |
! variables "ap", "bp", "presnivs". "pa" should be defined before |
31 |
! this procedure is called. |
! this procedure is called. |
32 |
|
|
33 |
|
! Libraries: |
34 |
use jumble, only: read_column, new_unit |
use jumble, only: read_column, new_unit |
35 |
use nr_util, only: pi, assert |
use nr_util, only: pi, assert |
36 |
|
|
37 |
|
use dynetat0_chosen_m, only: pa |
38 |
use unit_nml_m, only: unit_nml |
use unit_nml_m, only: unit_nml |
39 |
|
|
40 |
! Local: |
! Local: |
176 |
|
|
177 |
subroutine compute_ab |
subroutine compute_ab |
178 |
|
|
179 |
|
use dynetat0_chosen_m, only: pa |
180 |
|
|
181 |
! Calcul de "ap" et "bp". |
! Calcul de "ap" et "bp". |
182 |
|
|
183 |
where (s >= 1. / sqrt(1. - log(tiny(0.)))) |
where (s >= 1. / sqrt(1. - log(tiny(0.)))) |