Changeset 212 for LMDZ.3.3/trunk/libf/dyn3d
- Timestamp:
- Apr 24, 2001, 3:12:45 PM (24 years ago)
- Location:
- LMDZ.3.3/trunk/libf/dyn3d
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ.3.3/trunk/libf/dyn3d/fxhyp.F
r207 r212 3 3 c 4 4 SUBROUTINE fxhyp ( xzoomdeg,grossism,dzoom,tau , 5 , rlonm025,xprimm025,rlonv,xprimv,rlonu,xprimu,rlonp025,xprimp025) 5 , rlonm025,xprimm025,rlonv,xprimv,rlonu,xprimu,rlonp025,xprimp025, 6 , champmin,champmax ) 6 7 7 8 c Auteur : P. Le Van … … 49 50 INTEGER is2 50 51 SAVE is2 52 51 53 REAL*8 heavyside 52 54 EXTERNAL coefpoly,heavyside … … 57 59 xzoom = xzoomdeg * pi/180. 58 60 59 WRITE(6,24) xzoomdeg,grossism,tau,dzoom60 61 IF( dzoom.LT.1.) THEN 61 62 dzoom = dzoom * depi … … 68 69 ENDIF 69 70 71 WRITE(6,*) ' xzoom( rad.),grossism,tau,dzoom (radians)' 72 WRITE(6,24) xzoom,grossism,tau,dzoom 73 70 74 DO i = 0, nmax2 71 75 xtild(i) = - pi + FLOAT(i) * depi /nmax2 … … 117 121 118 122 beta = ( grossism * ffdx - pi ) / ( ffdx - pi ) 119 ccc WRITE(6,*) ' ** X beta **',beta,grossism120 123 121 124 IF( 2.*beta - grossism.LE. 0.) THEN … … 178 181 c ..... xuv = 0.5 si calcul aux pts U ........ 179 182 c 183 WRITE(6,18) 180 184 c 181 185 DO 5000 ik = 1, 4 … … 271 275 GO TO 1600 272 276 ELSE 273 WRITE(6,18)274 277 WRITE(6,*) 'Reorganisation des longitudes pour avoir entre - pi', 275 278 , ' et pi ' … … 331 334 332 335 IF( ik.EQ.1 ) THEN 333 WRITE(6,*) ' XLON aux pts. V-0.25 apres ( en deg. ) '334 WRITE(6,18)335 WRITE(6,68) xvrai336 c WRITE(6,*) ' XLON aux pts. V-0.25 apres ( en deg. ) ' 337 c WRITE(6,18) 338 c WRITE(6,68) xvrai 336 339 ccc WRITE(6,*) ' XPRIM k ',ik 337 340 ccc WRITE(6,566) xprimm … … 343 346 344 347 ELSE IF( ik.EQ.2 ) THEN 345 WRITE(6,18)346 WRITE(6,*) ' XLON aux pts. V apres ( en deg. ) '347 WRITE(6,68) xvrai348 c WRITE(6,18) 349 c WRITE(6,*) ' XLON aux pts. V apres ( en deg. ) ' 350 c WRITE(6,68) xvrai 348 351 ccc WRITE(6,*) ' XPRIM k ',ik 349 352 ccc WRITE(6,566) xprimm … … 355 358 356 359 ELSE IF( ik.EQ.3) THEN 357 WRITE(6,18)358 WRITE(6,*) ' XLON aux pts. U apres ( en deg. ) '359 WRITE(6,68) xvrai360 c WRITE(6,18) 361 c WRITE(6,*) ' XLON aux pts. U apres ( en deg. ) ' 362 c WRITE(6,68) xvrai 360 363 ccc WRITE(6,*) ' XPRIM ik ',ik 361 364 ccc WRITE(6,566) xprimm … … 367 370 368 371 ELSE IF( ik.EQ.4 ) THEN 369 WRITE(6,18)370 WRITE(6,*) ' XLON aux pts. V+0.25 apres ( en deg. ) '371 WRITE(6,68) xvrai372 c WRITE(6,18) 373 c WRITE(6,*) ' XLON aux pts. V+0.25 apres ( en deg. ) ' 374 c WRITE(6,68) xvrai 372 375 ccc WRITE(6,*) ' XPRIM ik ',ik 373 376 ccc WRITE(6,566) xprimm … … 381 384 382 385 5000 CONTINUE 386 c 387 WRITE(6,18) 383 388 c 384 389 c ........... fin de la boucle do 5000 ............ … … 396 401 champmax = champmax * 180./pi 397 402 398 WRITE(6,18)399 WRITE(6,*) ' Longitudes '400 WRITE(6,18)401 WRITE(6,3) champmin, champmax402 WRITE(6,*) ' Si cette derniere est trop lache , modifiez les par403 ,ametres grossism , tau , dzoom pour X et repasser ! '404 WRITE(6,18)405 406 3 Format(1x, ' Au centre du zoom , la longueur de la maille est',407 , ' d environ ',f8.2 ,' degres ',408 , ' alors que la maille en dehors de la zone du zoom est d environ409 , ', f8.2,' degres ' )410 403 18 FORMAT(/) 411 24 FORMAT(2x,' Parametres yzoom,gross,tau ,dzoom pour fyhyp ',4f8.3)404 24 FORMAT(2x,'Parametres xzoom,gross,tau ,dzoom pour fxhyp ',4f8.3) 412 405 68 FORMAT(1x,7f9.2) 413 406 -
LMDZ.3.3/trunk/libf/dyn3d/fxyhyper.F
r207 r212 39 39 REAL rlonu(iip1),xprimu(iip1),rlonv(iip1),xprimv(iip1), 40 40 , rlonm025(iip1),xprimm025(iip1), rlonp025(iip1),xprimp025(iip1) 41 REAL*8 dxmin, dxmax , dymin, dymax 41 42 42 43 c .... var. locales ..... … … 46 47 47 48 CALL fyhyp ( yzoom, grossy, dzoomy,tauy , 48 , rlatu, yprimu,rlatv,yprimv,rlatu2,yprimu2,rlatu1,yprimu1 ) 49 , rlatu, yprimu,rlatv,yprimv,rlatu2,yprimu2,rlatu1,yprimu1 , 50 , dymin,dymax ) 49 51 50 52 CALL fxhyp(xzoom,grossx,dzoomx,taux,rlonm025,xprimm025,rlonv, 51 , xprimv,rlonu,xprimu,rlonp025,xprimp025 )53 , xprimv,rlonu,xprimu,rlonp025,xprimp025 , dxmin,dxmax ) 52 54 53 55 … … 107 109 c 108 110 WRITE(6,*) ' *** TEST DE COHERENCE OK POUR FY **** ' 109 WRITE(6,25)110 25 FORMAT(//)111 111 c 112 WRITE(6,18) 113 WRITE(6,*) ' Latitudes ' 114 WRITE(6,*) ' *********** ' 115 WRITE(6,18) 116 WRITE(6,3) dymin, dymax 117 WRITE(6,*) ' Si cette derniere est trop lache , modifiez les par 118 ,ametres grossism , tau , dzoom pour Y et repasser ! ' 119 c 120 WRITE(6,18) 121 WRITE(6,*) ' Longitudes ' 122 WRITE(6,*) ' ************ ' 123 WRITE(6,18) 124 WRITE(6,3) dxmin, dxmax 125 WRITE(6,*) ' Si cette derniere est trop lache , modifiez les par 126 ,ametres grossism , tau , dzoom pour Y et repasser ! ' 127 WRITE(6,18) 128 c 129 3 Format(1x, ' Au centre du zoom , la longueur de la maille est', 130 , ' d environ ',f8.2 ,' degres ', 131 , ' alors que la maille en dehors de la zone du zoom est d environ 132 , ', f8.2,' degres ' ) 133 18 FORMAT(/) 112 134 113 135 RETURN -
LMDZ.3.3/trunk/libf/dyn3d/fyhyp.F
r207 r212 3 3 c 4 4 SUBROUTINE fyhyp ( yzoomdeg, grossism, dzoom,tau , 5 , rrlatu,yyprimu,rrlatv,yyprimv,rlatu2,yprimu2,rlatu1,yprimu1 ) 5 , rrlatu,yyprimu,rrlatv,yyprimv,rlatu2,yprimu2,rlatu1,yprimu1 , 6 , champmin,champmax ) 6 7 7 8 cc ... Version du 01/04/2001 .... … … 86 87 ENDIF 87 88 88 WRITE(6,*) ' yzoom ,dzoomy (radians),tau',y0,dzoom,tau 89 WRITE(6,18) 90 WRITE(6,*) ' yzoom( rad.),grossism,tau,dzoom (radians)' 91 WRITE(6,24) y0,grossism,tau,dzoom 89 92 90 93 DO i = 0, nmax2 … … 181 184 c ..... yuv = 0.5 si calcul des latitudes aux pts. V ..... 182 185 c 186 WRITE(6,18) 183 187 c 184 188 DO 5000 ik = 1,4 … … 253 257 ENDDO 254 258 255 WRITE(6,18)256 259 WRITE(6,*) 'Reorganisation des latitudes pour avoir entre - pi/2' 257 260 , ,' et pi/2 ' … … 296 299 297 300 IF( ik.EQ.1 ) THEN 298 WRITE(6,18)299 WRITE(6,*) ' YLAT en U apres ( en deg. ) '300 WRITE(6,68) (yvrai(j),j=1,jlat)301 c WRITE(6,18) 302 c WRITE(6,*) ' YLAT en U apres ( en deg. ) ' 303 c WRITE(6,68) (yvrai(j),j=1,jlat) 301 304 cc WRITE(6,*) ' YPRIM ' 302 305 cc WRITE(6,445) ( yprim(j),j=1,jlat) … … 308 311 309 312 ELSE IF ( ik.EQ. 2 ) THEN 310 WRITE(6,18)311 WRITE(6,*) ' YLAT en V apres ( en deg. ) '312 WRITE(6,68) (yvrai(j),j=1,jlat)313 c WRITE(6,18) 314 c WRITE(6,*) ' YLAT en V apres ( en deg. ) ' 315 c WRITE(6,68) (yvrai(j),j=1,jlat) 313 316 cc WRITE(6,*)' YPRIM ' 314 317 cc WRITE(6,445) ( yprim(j),j=1,jlat) … … 320 323 321 324 ELSE IF ( ik.EQ. 3 ) THEN 322 WRITE(6,18)323 WRITE(6,*) ' YLAT en U + 0.75 apres ( en deg. ) '324 WRITE(6,68) (yvrai(j),j=1,jlat)325 c WRITE(6,18) 326 c WRITE(6,*) ' YLAT en U + 0.75 apres ( en deg. ) ' 327 c WRITE(6,68) (yvrai(j),j=1,jlat) 325 328 cc WRITE(6,*) ' YPRIM ' 326 329 cc WRITE(6,445) ( yprim(j),j=1,jlat) … … 332 335 333 336 ELSE IF ( ik.EQ. 4 ) THEN 334 WRITE(6,18)335 WRITE(6,*) ' YLAT en U + 0.25 apres ( en deg. ) '336 WRITE(6,68)(yvrai(j),j=1,jlat)337 c WRITE(6,18) 338 c WRITE(6,*) ' YLAT en U + 0.25 apres ( en deg. ) ' 339 c WRITE(6,68)(yvrai(j),j=1,jlat) 337 340 cc WRITE(6,*) ' YPRIM ' 338 341 cc WRITE(6,68) ( yprim(j),j=1,jlat) … … 346 349 347 350 5000 CONTINUE 351 c 352 WRITE(6,18) 348 353 c 349 354 c ..... fin de la boucle do 5000 ..... … … 361 366 champmax = champmax * 180./pi 362 367 363 WRITE(6,18) 364 WRITE(6,*) ' Latitudes ' 365 WRITE(6,18) 366 WRITE(6,3) champmin, champmax 367 WRITE(6,*) ' Si cette derniere est trop lache , modifiez les par 368 ,ametres grossism , tau , dzoom pour Y et repasser ! ' 369 WRITE(6,18) 370 c 371 3 Format(1x, ' Au centre du zoom , la longueur de la maille est', 372 , ' d environ ',f8.2 ,' degres ', 373 , ' alors que la maille en dehors de la zone du zoom est d environ 374 , ', f8.2,' degres ' ) 368 24 FORMAT(2x,'Parametres yzoom,gross,tau ,dzoom pour fyhyp ',4f8.3) 375 369 18 FORMAT(/) 376 370 68 FORMAT(1x,7f9.2)
Note: See TracChangeset
for help on using the changeset viewer.
