166 |
use conf_gcm_m, ONLY : fxyhypb, ysinus |
use conf_gcm_m, ONLY : fxyhypb, ysinus |
167 |
use fxy_m, only: fxy |
use fxy_m, only: fxy |
168 |
use fxyhyper_m, only: fxyhyper |
use fxyhyper_m, only: fxyhyper |
169 |
|
use fxysinus_m, only: fxysinus |
170 |
use jumble, only: new_unit |
use jumble, only: new_unit |
171 |
use nr_util, only: pi |
use nr_util, only: pi |
172 |
USE paramet_m, ONLY : iip1, jjp1 |
USE paramet_m, ONLY : iip1, jjp1 |
173 |
USE serre, ONLY : alphax, alphay, clat, clon, dzoomx, dzoomy, grossismx, & |
USE serre, ONLY : alphax, alphay, clat, clon, dzoomx, dzoomy, grossismx, & |
174 |
grossismy, pxo, pyo, taux, tauy, transx, transy |
grossismy, pxo, pyo, taux, tauy, transx, transy |
175 |
! Modifies pxo, pyo, transx, transy |
! Modifiés pxo, pyo, transx, transy |
|
|
|
|
! Variables locales |
|
176 |
|
|
177 |
|
! Local: |
178 |
INTEGER i, j, itmax, itmay, iter, unit |
INTEGER i, j, itmax, itmay, iter, unit |
179 |
REAL cvu(iip1, jjp1), cuv(iip1, jjm) |
REAL cvu(iip1, jjp1), cuv(iip1, jjm) |
180 |
REAL ai14, ai23, airez, un4rad2 |
REAL ai14, ai23, airez, un4rad2 |
216 |
print *, "gamdi_grot = ", gamdi_grot |
print *, "gamdi_grot = ", gamdi_grot |
217 |
print *, "gamdi_h = ", gamdi_h |
print *, "gamdi_h = ", gamdi_h |
218 |
|
|
219 |
IF (.NOT. fxyhypb) THEN |
IF (fxyhypb) THEN |
220 |
|
! Utilisation de fxyhyper, f(x, y) à dérivée tangente hyperbolique |
221 |
|
print *, 'inigeom: Y = latitude, dérivée tangente hyperbolique' |
222 |
|
CALL fxyhyper(clat, grossismy, dzoomy, tauy, clon, grossismx, dzoomx, & |
223 |
|
taux, rlatu, yprimu, rlatv, yprimv, rlatu1, yprimu1, rlatu2, & |
224 |
|
yprimu2, rlonu, xprimu, rlonv, xprimv, rlonm025, xprimm025, & |
225 |
|
rlonp025, xprimp025) |
226 |
|
ELSE |
227 |
IF (ysinus) THEN |
IF (ysinus) THEN |
228 |
print *, ' Inigeom, Y = Sinus (Latitude) ' |
print *, 'inigeom: Y = sin(latitude)' |
229 |
! utilisation de f(x, y) avec y = sinus de la latitude |
! Utilisation de f(x, y) avec y = sinus de la latitude |
230 |
CALL fxysinus(rlatu, yprimu, rlatv, yprimv, rlatu1, yprimu1, & |
CALL fxysinus(rlatu, yprimu, rlatv, yprimv, rlatu1, yprimu1, & |
231 |
rlatu2, yprimu2, rlonu, xprimu, rlonv, xprimv, rlonm025, & |
rlatu2, yprimu2, rlonu, xprimu, rlonv, xprimv, rlonm025, & |
232 |
xprimm025, rlonp025, xprimp025) |
xprimm025, rlonp025, xprimp025) |
275 |
yprimu2, rlonu, xprimu, rlonv, xprimv, rlonm025, xprimm025, & |
yprimu2, rlonu, xprimu, rlonv, xprimv, rlonm025, xprimm025, & |
276 |
rlonp025, xprimp025) |
rlonp025, xprimp025) |
277 |
END IF |
END IF |
|
ELSE |
|
|
! Utilisation de fxyhyper, f(x, y) à dérivée tangente hyperbolique |
|
|
print *, 'Inigeom, Y = Latitude, dérivée tangente hyperbolique' |
|
|
CALL fxyhyper(clat, grossismy, dzoomy, tauy, clon, grossismx, dzoomx, & |
|
|
taux, rlatu, yprimu, rlatv, yprimv, rlatu1, yprimu1, rlatu2, & |
|
|
yprimu2, rlonu, xprimu, rlonv, xprimv, rlonm025, xprimm025, & |
|
|
rlonp025, xprimp025) |
|
278 |
END IF |
END IF |
279 |
|
|
280 |
rlatu(1) = pi / 2. |
rlatu(1) = pi / 2. |