13 |
c klevel et teta sont des arguments d'entree pour le s-prog |
c klevel et teta sont des arguments d'entree pour le s-prog |
14 |
c divgra est un argument de sortie pour le s-prog |
c divgra est un argument de sortie pour le s-prog |
15 |
c |
c |
16 |
|
use grad_m, only: grad |
17 |
use dimens_m |
use dimens_m |
18 |
use paramet_m |
use paramet_m |
19 |
use comgeom |
use comgeom |
23 |
c |
c |
24 |
c ............ variables en arguments .......... |
c ............ variables en arguments .......... |
25 |
c |
c |
26 |
INTEGER klevel |
INTEGER, intent(in):: klevel |
27 |
REAL teta( ip1jmp1,klevel ), divgra( ip1jmp1,klevel ) |
REAL teta( ip1jmp1,klevel ), divgra( ip1jmp1,klevel ) |
28 |
REAL cuvsga(ip1jm) , cvusga( ip1jmp1 ),unsaigam(ip1jmp1), |
REAL cuvsga(ip1jm) , cvusga( ip1jmp1 ),unsaigam(ip1jmp1), |
29 |
* unsapolnga, unsapolsga |
* unsapolnga, unsapolsga |
33 |
REAL ghy(ip1jm,llm), ghx(ip1jmp1,llm) |
REAL ghy(ip1jm,llm), ghx(ip1jmp1,llm) |
34 |
c ...................................................... |
c ...................................................... |
35 |
|
|
|
c |
|
|
c |
|
|
c ... cvuscugam = ( cvu/ cu ) ** (- gamdissip ) |
|
|
c ... cuvscvgam = ( cuv/ cv ) ** (- gamdissip ) calcules dans inigeom .. |
|
|
c ... unsairegam = 1. / aire ** (- gamdissip ) |
|
|
c |
|
|
|
|
36 |
CALL SCOPY ( ip1jmp1 * klevel, teta, 1, divgra, 1 ) |
CALL SCOPY ( ip1jmp1 * klevel, teta, 1, divgra, 1 ) |
37 |
c |
c |
38 |
CALL grad ( klevel, divgra, ghx, ghy ) |
CALL grad ( klevel, divgra, ghx, ghy ) |