Changeset 5098 for LMDZ6/branches/Amaury_dev/libf/filtrez
- Timestamp:
- Jul 22, 2024, 6:53:44 PM (2 months ago)
- Location:
- LMDZ6/branches/Amaury_dev/libf/filtrez
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/filtrez/filtreg_mod.F90
r5082 r5098 43 43 LOGICAL, SAVE :: first_call_inifilr = .TRUE. 44 44 45 #ifdef CRAY46 45 INTEGER ISMIN 47 46 EXTERNAL ISMIN 48 47 INTEGER iymin 49 48 INTEGER ixmineq 50 #endif51 49 ! 52 50 ! ------------------------------------------------------------ … … 104 102 ENDDO 105 103 ! 106 #ifdef CRAY 107 iymin = ISMIN( jjm, dlatu, 1 ) 108 ixmineq = ISMIN( iim, dlonu, 1 ) 109 dymin = dlatu( iymin ) 110 dxmin = dlonu( ixmineq ) 111 #else 104 112 105 dxmin = dlonu(1) 113 106 DO i = 2, iim … … 118 111 dymin = MIN( dymin,dlatu(j) ) 119 112 ENDDO 120 #endif121 113 ! 122 114 ! For a regular grid, we want the filter to start at latitudes … … 372 364 ENDDO 373 365 ENDDO ! of DO i=1,iim 374 #ifdef CRAY 375 CALL MXM( eignfnv,iim,eignft,iim,matriceun(1,1,j),iim ) 376 #else 366 377 367 #ifdef BLAS 378 368 CALL SGEMM ('N', 'N', iim, iim, iim, 1.0, & … … 389 379 ENDDO ! of DO k = 1, iim 390 380 #endif 391 #endif392 381 393 382 ENDDO ! of DO j = 2, jfiltnu … … 402 391 ENDDO 403 392 ENDDO ! of DO i=1,iim 404 #ifdef CRAY405 CALL MXM(eignfnv,iim,eignft,iim,matriceus(1,1,j-jfiltsu+1),iim)406 #else407 393 #ifdef BLAS 408 394 CALL SGEMM ('N', 'N', iim, iim, iim, 1.0, & … … 420 406 ENDDO ! of DO k = 1, iim 421 407 #endif 422 #endif423 408 424 409 ENDDO ! of DO j = jfiltsu, jjm … … 439 424 ENDDO 440 425 ENDDO 441 #ifdef CRAY 442 CALL MXM( eignfnu,iim,eignft,iim,matricevn(1,1,j),iim ) 443 #else 426 444 427 #ifdef BLAS 445 428 CALL SGEMM ('N', 'N', iim, iim, iim, 1.0, & … … 456 439 ENDDO 457 440 #endif 458 #endif459 441 460 442 ENDDO ! of DO j = 1, jfiltnv … … 469 451 ENDDO 470 452 ENDDO 471 #ifdef CRAY 472 CALL MXM(eignfnu,iim,eignft,iim,matricevs(1,1,j-jfiltsv+1),iim) 473 #else 453 474 454 #ifdef BLAS 475 455 CALL SGEMM ('N', 'N', iim, iim, iim, 1.0, & … … 487 467 ENDDO 488 468 #endif 489 #endif490 469 491 470 ENDDO ! of DO j = jfiltsv, jjm … … 506 485 ENDDO 507 486 ENDDO 508 #ifdef CRAY 509 CALL MXM( eignfnv,iim,eignft,iim,matrinvn(1,1,j),iim ) 510 #else 487 511 488 #ifdef BLAS 512 489 CALL SGEMM ('N', 'N', iim, iim, iim, 1.0, & … … 523 500 ENDDO 524 501 #endif 525 #endif526 502 527 503 ENDDO ! of DO j = 2, jfiltnu … … 536 512 ENDDO 537 513 ENDDO 538 #ifdef CRAY539 CALL MXM(eignfnv,iim,eignft,iim,matrinvs(1,1,j-jfiltsu+1),iim)540 #else541 514 #ifdef BLAS 542 515 CALL SGEMM ('N', 'N', iim, iim, iim, 1.0, & … … 553 526 ENDDO 554 527 #endif 555 #endif556 528 557 529 ENDDO ! of DO j = jfiltsu, jjm … … 569 541 ! 570 542 334 FORMAT(1x,24i3) 571 755 FORMAT(1x,6f10.3,i3) 572 573 RETURN 543 574 544 END SUBROUTINE inifilr 575 545 -
LMDZ6/branches/Amaury_dev/libf/filtrez/inifgn.F
r5079 r5098 67 67 END DO 68 68 c 69 #ifdef CRAY70 CALL MXM(eignfnu,iim,eignfnv,iim,vec ,iim)71 CALL MXM(eignfnv,iim,eignfnu,iim,vec1,iim)72 #else73 69 DO j = 1, iim 74 70 DO i = 1, iim … … 81 77 ENDDO 82 78 ENDDO 83 #endif84 79 85 80 c
Note: See TracChangeset
for help on using the changeset viewer.