Changeset 212 for LMDZ.3.3/trunk/libf/dyn3d/fyhyp.F
- Timestamp:
- Apr 24, 2001, 3:12:45 PM (23 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
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.