Changeset 203 for LMDZ.3.3/trunk/libf/dyn3d/fxyhyper.F
- Timestamp:
- Apr 12, 2001, 11:30:31 AM (23 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ.3.3/trunk/libf/dyn3d/fxyhyper.F
r2 r203 1 SUBROUTINE fxyhyper ( yzoom, grossy, dzoomy,tauy,deltay , 2 , xzoom, grossx, dzoomx,taux, 1 c 2 c $Header 3 c 4 SUBROUTINE fxyhyper ( yzoom, grossy, dzoomy,tauy , 5 , xzoom, grossx, dzoomx,taux , 3 6 , rlatu,yprimu,rlatv,yprimv,rlatu1, yprimu1, rlatu2, yprimu2 , 4 7 , rlonu,xprimu,rlonv,xprimv,rlonm025,xprimm025,rlonp025,xprimp025) … … 19 22 c a) le grossissement du zoom : grossy ( en y ) et grossx ( en x ) 20 23 c b) l' extension du zoom : dzoomy ( en y ) et dzoomx ( en x ) 21 c c) la raideur du zoom : tau , ici = 1.24 c c) la raideur de la transition du zoom : taux et tauy 22 25 c 23 c N.B : le produit ( grossissement * extension ) doit etre < 1. 24 c ******* 25 c En plus , il y a un autre parametre , moins important mais quand 26 c meme utile , c'est deltay , deplacement de la zone du zoom en 27 c latitude : Si on deplace de 10. degres vers le nord ( deltay 28 c = 10. dans inigeom ) . 29 c 26 c N.B : Il vaut mieux avoir : grossx * dzoomx < pi ( radians ) 27 c ****** 28 c et grossy * dzoomy < pi/2 ( radians ) 30 29 c 31 30 #include "dimensions.h" … … 40 39 REAL rlonu(iip1),xprimu(iip1),rlonv(iip1),xprimv(iip1), 41 40 , rlonm025(iip1),xprimm025(iip1), rlonp025(iip1),xprimp025(iip1) 42 REAL deltay43 41 44 42 c .... var. locales ..... … … 47 45 c 48 46 49 CALL fyhyp ( yzoom, grossy, dzoomy,tauy , deltay,47 CALL fyhyp ( yzoom, grossy, dzoomy,tauy , 50 48 , rlatu, yprimu,rlatv,yprimv,rlatu2,yprimu2,rlatu1,yprimu1 ) 51 52 49 53 50 CALL fxhyp(xzoom,grossx,dzoomx,taux,rlonm025,xprimm025,rlonv, … … 55 52 56 53 57 DO i = 1, ii m54 DO i = 1, iip1 58 55 IF(rlonp025(i).LT.rlonv(i)) THEN 59 PRINT *,' Attention ! rlonp025 < rlonv',i56 WRITE(6,*) ' Attention ! rlonp025 < rlonv',i 60 57 STOP 61 58 ENDIF 62 59 63 60 IF(rlonv(i).LT.rlonm025(i)) THEN 64 PRINT *,' Attention ! rlonm025 > rlonv',i61 WRITE(6,*) ' Attention ! rlonm025 > rlonv',i 65 62 STOP 66 63 ENDIF 67 64 68 65 IF(rlonp025(i).GT.rlonu(i)) THEN 69 PRINT *,' Attention ! rlonp025 > rlonu',i66 WRITE(6,*) ' Attention ! rlonp025 > rlonu',i 70 67 STOP 71 68 ENDIF 72 69 ENDDO 73 70 74 PRINT *,' *** TEST DE COHERENCE OK POUR FX **** '71 WRITE(6,*) ' *** TEST DE COHERENCE OK POUR FX **** ' 75 72 76 73 c … … 78 75 c 79 76 IF(rlatu1(j).LE.rlatu2(j)) THEN 80 PRINT *,'Attention ! rlatu1 < rlatu2 ',rlatu1(j), rlatu2(j),j77 WRITE(6,*)'Attention ! rlatu1 < rlatu2 ',rlatu1(j), rlatu2(j),j 81 78 STOP 13 82 79 ENDIF 83 80 c 84 81 IF(rlatu2(j).LE.rlatu(j+1)) THEN 85 PRINT *,'Attention ! rlatu2 < rlatup1 ',rlatu2(j),rlatu(j+1),j82 WRITE(6,*)'Attention ! rlatu2 < rlatup1 ',rlatu2(j),rlatu(j+1),j 86 83 STOP 14 87 84 ENDIF 88 85 c 89 86 IF(rlatu(j).LE.rlatu1(j)) THEN 90 PRINT *,' Attention ! rlatu < rlatu1 ',rlatu(j),rlatu1(j),j87 WRITE(6,*)' Attention ! rlatu < rlatu1 ',rlatu(j),rlatu1(j),j 91 88 STOP 15 92 89 ENDIF 93 90 c 94 91 IF(rlatv(j).LE.rlatu2(j)) THEN 95 PRINT *,' Attention ! rlatv < rlatu2 ',rlatv(j),rlatu2(j),j92 WRITE(6,*)' Attention ! rlatv < rlatu2 ',rlatv(j),rlatu2(j),j 96 93 STOP 16 97 94 ENDIF 98 95 c 99 96 IF(rlatv(j).ge.rlatu1(j)) THEN 100 PRINT *,' Attention ! rlatv > rlatu1 ',rlatv(j),rlatu1(j),j97 WRITE(6,*)' Attention ! rlatv > rlatu1 ',rlatv(j),rlatu1(j),j 101 98 STOP 17 102 99 ENDIF 103 100 c 104 101 IF(rlatv(j).ge.rlatu(j)) THEN 105 PRINT *,'Attention ! rlatv > rlatu ',rlatv(j),rlatu(j),j102 WRITE(6,*) ' Attention ! rlatv > rlatu ',rlatv(j),rlatu(j),j 106 103 STOP 18 107 104 ENDIF … … 109 106 ENDDO 110 107 c 111 PRINT *,' *** TEST DE COHERENCE OK POUR FY **** ' 112 108 WRITE(6,*) ' *** TEST DE COHERENCE OK POUR FY **** ' 109 WRITE(6,25) 110 25 FORMAT(//) 113 111 c 114 112
Note: See TracChangeset
for help on using the changeset viewer.