Ignore:
Timestamp:
Oct 21, 2024, 2:58:45 PM (23 hours ago)
Author:
abarral
Message:

Convert fixed-form to free-form sources .F -> .{f,F}90
(WIP: some .F remain, will be handled in subsequent commits)

File:
1 moved

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/filtrez/eigen.f90

    r5245 r5246  
    22! $Header$
    33!
    4       SUBROUTINE eigen( e,d)
    5       IMPLICIT NONE
    6       INCLUDE "dimensions.h"
    7       real :: e( iim,iim ), d( iim )
    8       real :: asm( iim )
    9       integer :: im,i,j
    10       im=iim
    11 c
    12       DO 48 i = 1,im
    13         asm( i ) = d( im-i+1 )
    14  48   CONTINUE
    15       DO 49 i = 1,iim
    16         d( i ) = asm( i )
    17  49   CONTINUE
    18 c
    19 c    PRINT 70,d
     4SUBROUTINE eigen( e,d)
     5  IMPLICIT NONE
     6  INCLUDE "dimensions.h"
     7  real :: e( iim,iim ), d( iim )
     8  real :: asm( iim )
     9  integer :: im,i,j
     10  im=iim
     11  !
     12  DO i = 1,im
     13    asm( i ) = d( im-i+1 )
     14  END DO
     15  DO i = 1,iim
     16    d( i ) = asm( i )
     17  END DO
     18  !
     19  ! PRINT 70,d
    2020 70   FORMAT(5x,'Valeurs propres',/,8(1x,8f10.4,/),/)
    21                 print *
    22 c
    23       DO 51 i = 1,im
    24          DO 52 j = 1,im
    25             asm( j ) = e( i , im-j+1 )
    26  52      CONTINUE
    27          DO 50 j = 1,im
    28             e( i,j ) = asm( j )
    29  50      CONTINUE
    30  51   CONTINUE
     21            print *
     22  !
     23  DO i = 1,im
     24     DO j = 1,im
     25        asm( j ) = e( i , im-j+1 )
     26     END DO
     27     DO j = 1,im
     28        e( i,j ) = asm( j )
     29     END DO
     30  END DO
    3131
    32       RETURN
    33       END
     32  RETURN
     33END SUBROUTINE eigen
Note: See TracChangeset for help on using the changeset viewer.