37 |
REAL ghy(ip1jm, klevel), ghx(ip1jmp1, klevel) |
REAL ghy(ip1jm, klevel), ghx(ip1jmp1, klevel) |
38 |
! ...................................................... |
! ...................................................... |
39 |
|
|
40 |
CALL scopy(ip1jmp1*klevel, teta, 1, divgra, 1) |
divgra = teta |
41 |
CALL grad(klevel, divgra, ghx, ghy) |
CALL grad(klevel, divgra, ghx, ghy) |
42 |
CALL diverg_gam(klevel, cuvsga, cvusga, unsaigam, unsapolnga, unsapolsga, & |
CALL diverg_gam(klevel, cuvsga, cvusga, unsaigam, unsapolnga, unsapolsga, & |
43 |
ghx, ghy, divgra) |
ghx, ghy, divgra) |