[524] | 1 | ! |
---|
| 2 | ! $Header$ |
---|
| 3 | ! |
---|
| 4 | c |
---|
| 5 | c |
---|
| 6 | SUBROUTINE fxyhyper ( yzoom, grossy, dzoomy,tauy , |
---|
| 7 | , xzoom, grossx, dzoomx,taux , |
---|
| 8 | , rlatu,yprimu,rlatv,yprimv,rlatu1, yprimu1, rlatu2, yprimu2 , |
---|
| 9 | , rlonu,xprimu,rlonv,xprimv,rlonm025,xprimm025,rlonp025,xprimp025) |
---|
| 10 | |
---|
| 11 | IMPLICIT NONE |
---|
| 12 | c |
---|
| 13 | c Auteur : P. Le Van . |
---|
| 14 | c |
---|
| 15 | c d'apres formulations de R. Sadourny . |
---|
| 16 | c |
---|
| 17 | c |
---|
| 18 | c Ce spg calcule les latitudes( routine fyhyp ) et longitudes( fxhyp ) |
---|
| 19 | c par des fonctions a tangente hyperbolique . |
---|
| 20 | c |
---|
| 21 | c Il y a 3 parametres ,en plus des coordonnees du centre du zoom (xzoom |
---|
| 22 | c et yzoom ) : |
---|
| 23 | c |
---|
| 24 | c a) le grossissement du zoom : grossy ( en y ) et grossx ( en x ) |
---|
| 25 | c b) l' extension du zoom : dzoomy ( en y ) et dzoomx ( en x ) |
---|
| 26 | c c) la raideur de la transition du zoom : taux et tauy |
---|
| 27 | c |
---|
| 28 | c N.B : Il vaut mieux avoir : grossx * dzoomx < pi ( radians ) |
---|
| 29 | c ****** |
---|
| 30 | c et grossy * dzoomy < pi/2 ( radians ) |
---|
| 31 | c |
---|
| 32 | #include "dimensions.h" |
---|
| 33 | #include "paramet.h" |
---|
| 34 | |
---|
| 35 | |
---|
| 36 | c ..... Arguments ... |
---|
| 37 | c |
---|
| 38 | REAL xzoom,yzoom,grossx,grossy,dzoomx,dzoomy,taux,tauy |
---|
| 39 | REAL rlatu(jjp1), yprimu(jjp1),rlatv(jjm), yprimv(jjm), |
---|
| 40 | , rlatu1(jjm), yprimu1(jjm), rlatu2(jjm), yprimu2(jjm) |
---|
| 41 | REAL rlonu(iip1),xprimu(iip1),rlonv(iip1),xprimv(iip1), |
---|
| 42 | , rlonm025(iip1),xprimm025(iip1), rlonp025(iip1),xprimp025(iip1) |
---|
| 43 | REAL*8 dxmin, dxmax , dymin, dymax |
---|
| 44 | |
---|
| 45 | c .... var. locales ..... |
---|
| 46 | c |
---|
| 47 | INTEGER i,j |
---|
| 48 | c |
---|
| 49 | |
---|
| 50 | CALL fyhyp ( yzoom, grossy, dzoomy,tauy , |
---|
| 51 | , rlatu, yprimu,rlatv,yprimv,rlatu2,yprimu2,rlatu1,yprimu1 , |
---|
| 52 | , dymin,dymax ) |
---|
| 53 | |
---|
| 54 | CALL fxhyp(xzoom,grossx,dzoomx,taux,rlonm025,xprimm025,rlonv, |
---|
| 55 | , xprimv,rlonu,xprimu,rlonp025,xprimp025 , dxmin,dxmax ) |
---|
| 56 | |
---|
| 57 | |
---|
| 58 | DO i = 1, iip1 |
---|
| 59 | IF(rlonp025(i).LT.rlonv(i)) THEN |
---|
| 60 | WRITE(6,*) ' Attention ! rlonp025 < rlonv',i |
---|
| 61 | STOP |
---|
| 62 | ENDIF |
---|
| 63 | |
---|
| 64 | IF(rlonv(i).LT.rlonm025(i)) THEN |
---|
| 65 | WRITE(6,*) ' Attention ! rlonm025 > rlonv',i |
---|
| 66 | STOP |
---|
| 67 | ENDIF |
---|
| 68 | |
---|
| 69 | IF(rlonp025(i).GT.rlonu(i)) THEN |
---|
| 70 | WRITE(6,*) ' Attention ! rlonp025 > rlonu',i |
---|
| 71 | STOP |
---|
| 72 | ENDIF |
---|
| 73 | ENDDO |
---|
| 74 | |
---|
| 75 | WRITE(6,*) ' *** TEST DE COHERENCE OK POUR FX **** ' |
---|
| 76 | |
---|
| 77 | c |
---|
| 78 | DO j = 1, jjm |
---|
| 79 | c |
---|
| 80 | IF(rlatu1(j).LE.rlatu2(j)) THEN |
---|
| 81 | WRITE(6,*)'Attention ! rlatu1 < rlatu2 ',rlatu1(j), rlatu2(j),j |
---|
| 82 | STOP 13 |
---|
| 83 | ENDIF |
---|
| 84 | c |
---|
| 85 | IF(rlatu2(j).LE.rlatu(j+1)) THEN |
---|
| 86 | WRITE(6,*)'Attention ! rlatu2 < rlatup1 ',rlatu2(j),rlatu(j+1),j |
---|
| 87 | STOP 14 |
---|
| 88 | ENDIF |
---|
| 89 | c |
---|
| 90 | IF(rlatu(j).LE.rlatu1(j)) THEN |
---|
| 91 | WRITE(6,*)' Attention ! rlatu < rlatu1 ',rlatu(j),rlatu1(j),j |
---|
| 92 | STOP 15 |
---|
| 93 | ENDIF |
---|
| 94 | c |
---|
| 95 | IF(rlatv(j).LE.rlatu2(j)) THEN |
---|
| 96 | WRITE(6,*)' Attention ! rlatv < rlatu2 ',rlatv(j),rlatu2(j),j |
---|
| 97 | STOP 16 |
---|
| 98 | ENDIF |
---|
| 99 | c |
---|
| 100 | IF(rlatv(j).ge.rlatu1(j)) THEN |
---|
| 101 | WRITE(6,*)' Attention ! rlatv > rlatu1 ',rlatv(j),rlatu1(j),j |
---|
| 102 | STOP 17 |
---|
| 103 | ENDIF |
---|
| 104 | c |
---|
| 105 | IF(rlatv(j).ge.rlatu(j)) THEN |
---|
| 106 | WRITE(6,*) ' Attention ! rlatv > rlatu ',rlatv(j),rlatu(j),j |
---|
| 107 | STOP 18 |
---|
| 108 | ENDIF |
---|
| 109 | c |
---|
| 110 | ENDDO |
---|
| 111 | c |
---|
| 112 | WRITE(6,*) ' *** TEST DE COHERENCE OK POUR FY **** ' |
---|
| 113 | c |
---|
| 114 | WRITE(6,18) |
---|
| 115 | WRITE(6,*) ' Latitudes ' |
---|
| 116 | WRITE(6,*) ' *********** ' |
---|
| 117 | WRITE(6,18) |
---|
| 118 | WRITE(6,3) dymin, dymax |
---|
| 119 | WRITE(6,*) ' Si cette derniere est trop lache , modifiez les par |
---|
| 120 | ,ametres grossism , tau , dzoom pour Y et repasser ! ' |
---|
| 121 | c |
---|
| 122 | WRITE(6,18) |
---|
| 123 | WRITE(6,*) ' Longitudes ' |
---|
| 124 | WRITE(6,*) ' ************ ' |
---|
| 125 | WRITE(6,18) |
---|
| 126 | WRITE(6,3) dxmin, dxmax |
---|
| 127 | WRITE(6,*) ' Si cette derniere est trop lache , modifiez les par |
---|
| 128 | ,ametres grossism , tau , dzoom pour Y et repasser ! ' |
---|
| 129 | WRITE(6,18) |
---|
| 130 | c |
---|
| 131 | 3 Format(1x, ' Au centre du zoom , la longueur de la maille est', |
---|
| 132 | , ' d environ ',f8.2 ,' degres ', |
---|
| 133 | , ' alors que la maille en dehors de la zone du zoom est d environ |
---|
| 134 | , ', f8.2,' degres ' ) |
---|
| 135 | 18 FORMAT(/) |
---|
| 136 | |
---|
| 137 | RETURN |
---|
| 138 | END |
---|
| 139 | |
---|