20 |
! dlonid(idat): abscisses des interfaces des mailles donnees |
! dlonid(idat): abscisses des interfaces des mailles donnees |
21 |
! rlonimod(imod): abscisses des interfaces des mailles modele |
! rlonimod(imod): abscisses des interfaces des mailles modele |
22 |
! (L'indice 1 correspond a l'interface mailLE 1 / maille 2) |
! (L'indice 1 correspond a l'interface mailLE 1 / maille 2) |
23 |
! (Les abscisses sont exprimées en degres) |
! (Les abscisses sont exprim\'ees en degres) |
24 |
|
|
25 |
use nr_util, only: assert_eq, pi |
use nr_util, only: assert_eq, pi |
26 |
|
|
36 |
REAL xxid(size(dlonid)+1), xxd(size(dlonid)+1), fdd(size(dlonid)+1) |
REAL xxid(size(dlonid)+1), xxd(size(dlonid)+1), fdd(size(dlonid)+1) |
37 |
REAL fxd(size(dlonid)+1), xchan(size(dlonid)+1), fdchan(size(dlonid)+1) |
REAL fxd(size(dlonid)+1), xchan(size(dlonid)+1), fdchan(size(dlonid)+1) |
38 |
REAL xxim(size(rlonimod)) |
REAL xxim(size(rlonimod)) |
|
|
|
39 |
REAL x0, xim0, dx, dxm |
REAL x0, xim0, dx, dxm |
40 |
REAL chmin, chmax |
REAL chmax |
|
|
|
41 |
INTEGER imod, idat, i, ichang, id0, id1, nid, idatmax1 |
INTEGER imod, idat, i, ichang, id0, id1, nid, idatmax1 |
42 |
|
|
43 |
!----------------------------------------------------- |
!----------------------------------------------------- |
51 |
xxim(imod) = rlonimod(imod) |
xxim(imod) = rlonimod(imod) |
52 |
ENDDO |
ENDDO |
53 |
|
|
54 |
CALL minmax(imodmax, xxim, chmin, chmax) |
chmax = maxval(xxim) |
55 |
IF(chmax < 6.50) THEN |
IF(chmax < 6.50) THEN |
56 |
DO imod = 1, imodmax |
DO imod = 1, imodmax |
57 |
xxim(imod) = xxim(imod) * 180./pi |
xxim(imod) = xxim(imod) * 180./pi |
70 |
xxd(idat) = dlonid(idat) |
xxd(idat) = dlonid(idat) |
71 |
ENDDO |
ENDDO |
72 |
|
|
73 |
CALL minmax(idatmax, xxd, chmin, chmax) |
chmax = maxval(xxd(:idatmax)) |
74 |
IF(chmax < 6.50) THEN |
IF(chmax < 6.50) THEN |
75 |
DO idat = 1, idatmax |
DO idat = 1, idatmax |
76 |
xxd(idat) = xxd(idat) * 180./pi |
xxd(idat) = xxd(idat) * 180./pi |