Changeset 232 for LMDZ.3.3/branches/rel-LF/libf/dyn3d/inigeom.F
- Timestamp:
- Jun 20, 2001, 3:53:15 PM (23 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ.3.3/branches/rel-LF/libf/dyn3d/inigeom.F
r2 r232 1 c 2 c $Header$ 3 c 1 4 SUBROUTINE inigeom 2 5 c 3 6 c Auteur : P. Le Van 4 c ..................... 5 c 6 c ............ Version du 20/12/98 ........................ 7 c 8 c ............ Version du 01/04/2001 ........................ 7 9 c 8 10 c Calcul des elongations cuij1,.cuij4 , cvij1,..cvij4 aux memes en- 9 11 c endroits que les aires aireij1,..aireij4 . 12 10 13 c Choix entre f(y) a derivee sinusoid. ou a derivee tangente hyperbol. 11 14 c 12 15 c 13 16 IMPLICIT NONE 14 c15 REAL deltay, tauy, taux16 PARAMETER ( deltay = 0. , tauy = 1. , taux = 1. )17 c18 c deltay est ( en degres ) le deplacement eventuel en Y du zoom19 cc taux et tauy sont les raideurs du zoom20 c21 17 c 22 18 #include "dimensions.h" … … 50 46 SAVE rlonm025,xprimm025,rlonp025,xprimp025 51 47 52 53 54 c----------------------------------------------------------------------55 48 REAL SSUM 56 49 EXTERNAL SSUM 57 c58 50 c 59 51 c … … 62 54 c - calcul des coeff. ( cu, cv , 1./cu**2, 1./cv**2 ) - 63 55 c - - 64 c ------------------------------------------------------------------65 56 c ------------------------------------------------------------------ 66 57 c … … 168 159 c 169 160 c 170 PRINT 3161 WRITE(6,3) 171 162 3 FORMAT( // 10x,' .... INIGEOM date du 01/06/98 ..... ', 172 163 * //5x,' Calcul des elongations cu et cv comme sommes des 4 ' / … … 191 182 ENDIF 192 183 193 PRINT *,' gamdi_gd ',gamdi_gdiv,gamdi_grot,gamdi_h,coefdis,184 WRITE(6,*) ' gamdi_gd ',gamdi_gdiv,gamdi_grot,gamdi_h,coefdis, 194 185 * nitergdiv,nitergrot,niterh 195 186 c 196 187 pi = 2.* ASIN(1.) 197 188 c 198 PRINT 990189 WRITE(6,990) 199 190 200 191 c ---------------------------------------------------------------- … … 205 196 IF( ysinus ) THEN 206 197 c 207 PRINT *,' *** Inigeom , Y = Sinus ( Latitude ) *** '198 WRITE(6,*) ' *** Inigeom , Y = Sinus ( Latitude ) *** ' 208 199 c 209 200 c .... utilisation de f(x,y ) avec y = sinus de la latitude ..... … … 215 206 ELSE 216 207 c 217 PRINT *,'*** Inigeom , Y = Latitude , der. sinusoid . ***'208 WRITE(6,*) '*** Inigeom , Y = Latitude , der. sinusoid . ***' 218 209 219 210 c .... utilisation de f(x,y) a tangente sinusoidale , y etant la latit. ... … … 267 258 ELSE 268 259 c 269 c .... utilisation de fxyhyper , f(x,y) a derivee tangente hyperbol.260 c .... Utilisation de fxyhyper , f(x,y) a derivee tangente hyperbol. 270 261 c ..................................................................... 271 262 272 PRINT *,'*** Inigeom , Y = Latitude , der.tg. hyperbolique ***'273 274 CALL fxyhyper( clat, grossismy, dzoomy, tauy , deltay ,275 , clon, grossismx, dzoomx, taux,263 WRITE(6,*)'*** Inigeom , Y = Latitude , der.tg. hyperbolique ***' 264 265 CALL fxyhyper( clat, grossismy, dzoomy, tauy , 266 , clon, grossismx, dzoomx, taux , 276 267 , rlatu,yprimu,rlatv, yprimv,rlatu1, yprimu1,rlatu2,yprimu2 , 277 268 , rlonu,xprimu,rlonv,xprimv,rlonm025,xprimm025,rlonp025,xprimp025 ) … … 387 378 c 388 379 IF ( j. eq. jjp1 ) THEN 389 390 380 yprp = yprimu2(j-1) 391 381 rlatp = rlatu2 (j-1) 392 cc yprp = fyprim( FLOAT(j) - 0.25 )393 cc rlatp = fy ( FLOAT(j) - 0.25 )382 ccc yprp = fyprim( FLOAT(j) - 0.25 ) 383 ccc rlatp = fy ( FLOAT(j) - 0.25 ) 394 384 c 395 385 coslatp = COS( rlatp ) … … 425 415 rlatm = rlatu1 ( j ) 426 416 yprm = yprimu1( j ) 427 c rlatp = fy ( FLOAT(j) - 0.25 )428 c yprp = fyprim( FLOAT(j) - 0.25 )429 c rlatm = fy ( FLOAT(j) + 0.25 )430 c yprm = fyprim( FLOAT(j) + 0.25 )417 cc rlatp = fy ( FLOAT(j) - 0.25 ) 418 cc yprp = fyprim( FLOAT(j) - 0.25 ) 419 cc rlatm = fy ( FLOAT(j) + 0.25 ) 420 cc yprm = fyprim( FLOAT(j) + 0.25 ) 431 421 432 422 coslatm = COS( rlatm ) … … 490 480 36 CONTINUE 491 481 c 492 c .... Modif P. Le Van ( 4/07/96 ) .....493 482 c 494 483 aire (iip1,j) = aire (1,j) … … 620 609 aiuscv2gam(iip1,j) = aiuscv2gam(1,j) 621 610 ENDDO 622 c 611 623 612 c 624 613 c calcul des aires aux poles : … … 666 655 c----------------------------------------------------------------------- 667 656 c 668 PRINT *,' INIGEOM RLONV ' 657 WRITE(6,*) ' *** Coordonnees de la grille *** ' 658 WRITE(6,995) 659 c 660 WRITE(6,*) ' LONGITUDES aux pts. V ( degres ) ' 661 WRITE(6,995) 669 662 DO i=1,iip1 670 663 rlonvv(i) = rlonv(i)*180./pi 671 664 ENDDO 672 PRINT 400,rlonvv 673 c 674 PRINT *,' RLATV ' 665 WRITE(6,400) rlonvv 666 c 667 WRITE(6,995) 668 WRITE(6,*) ' LATITUDES aux pts. V ( degres ) ' 669 WRITE(6,995) 675 670 DO i=1,jjm 676 671 rlatuu(i)=rlatv(i)*180./pi 677 672 ENDDO 678 PRINT 400,(rlatuu(i),i=1,jjm)673 WRITE(6,400) (rlatuu(i),i=1,jjm) 679 674 c 680 675 DO i=1,iip1 681 676 rlonvv(i)=rlonu(i)*180./pi 682 677 ENDDO 683 PRINT *,' RLONU ' 684 PRINT 400,rlonvv 685 c 686 PRINT *,' RLATU ' 678 WRITE(6,995) 679 WRITE(6,*) ' LONGITUDES aux pts. U ( degres ) ' 680 WRITE(6,995) 681 WRITE(6,400) rlonvv 682 WRITE(6,995) 683 684 WRITE(6,*) ' LATITUDES aux pts. U ( degres ) ' 685 WRITE(6,995) 687 686 DO i=1,jjp1 688 687 rlatuu(i)=rlatu(i)*180./pi 689 688 ENDDO 690 PRINT 400,(rlatuu(i),i=1,jjp1) 691 c 689 WRITE(6,400) (rlatuu(i),i=1,jjp1) 690 WRITE(6,995) 691 c 692 444 format(f10.3,f6.0) 692 693 400 FORMAT(1x,8f8.2) 693 694 990 FORMAT(//) 695 995 FORMAT(/) 694 696 c 695 697 RETURN
Note: See TracChangeset
for help on using the changeset viewer.