Ignore:
Timestamp:
Apr 8, 2014 3:21:23 PM (10 years ago)
Author:
heinze
Message:

REAL constants provided with KIND-attribute

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/cpulog.f90

    r1323 r1353  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! REAL constants provided with KIND-attribute
    2323!
    2424! Former revisions:
     
    100100    END TYPE logpoint
    101101
    102     TYPE(logpoint), DIMENSION(100) ::  log_point = logpoint( 0.0, 0.0, 0.0,   &
    103                                        0.0, 0.0, 0.0, 0.0, 0, ' ' ),          &
    104                                        log_point_s = logpoint( 0.0, 0.0, 0.0, &
    105                                        0.0, 0.0, 0.0, 0.0, 0, ' ' )
     102    TYPE(logpoint), DIMENSION(100) ::  log_point = logpoint( 0.0_wp, 0.0_wp,   &
     103                                       0.0_wp, 0.0_wp, 0.0_wp, 0.0_wp, 0.0_wp, &
     104                                       0, ' ' ),                               &
     105                                       log_point_s = logpoint( 0.0_wp, 0.0_wp, &
     106                                       0.0_wp, 0.0_wp, 0.0_wp, 0.0_wp, 0.0_wp, &
     107                                       0, ' ' )
    106108
    107109    SAVE
     
    120122       LOGICAL, SAVE     ::  first = .TRUE.     !:
    121123       
    122        REAL(wp)          ::  mtime = 0.0        !:
    123        REAL(wp)          ::   mtimevec = 0.0    !:
     124       REAL(wp)          ::  mtime = 0.0_wp     !:
     125       REAL(wp)          ::  mtimevec = 0.0_wp  !:
    124126       TYPE(logpoint)    ::  log_event          !:
    125127
     
    212214          log_event%vector   = log_event%vector + log_event%mtimevec
    213215          log_event%counts   = log_event%counts + 1
    214           log_event%isum     = 0.0
    215           log_event%ivect    = 0.0
     216          log_event%isum     = 0.0_wp
     217          log_event%ivect    = 0.0_wp
    216218       ELSE
    217219          message_string = 'unknown modus of time measurement: ' // TRIM( modus )
     
    239241       INTEGER(iwp)    ::  sender          !:
    240242       REAL(wp)       ::  average_cputime  !:
    241        REAL(wp), SAVE ::  norm = 1.0       !:
     243       REAL(wp), SAVE ::  norm = 1.0_wp    !:
    242244       REAL(wp), DIMENSION(:),   ALLOCATABLE ::  pe_max        !:
    243245       REAL(wp), DIMENSION(:),   ALLOCATABLE ::  pe_min        !:
     
    267269          pe_min = log_point%sum
    268270          pe_max = log_point%sum    ! need to be set in case of 1 PE
    269           pe_rms = 0.0
     271          pe_rms = 0.0_wp
    270272
    271273#if defined( __parallel )
     
    319321             sum = log_point%sum
    320322          ELSEWHERE
    321              sum = -1.0
     323             sum = -1.0_wp
    322324          ENDWHERE
    323325
     
    328330                               REAL( nr_timesteps_this_run, KIND=wp ) * 1E6_wp  ! in micro-sec
    329331          ELSE
    330              average_cputime = -1.0
     332             average_cputime = -1.0_wp
    331333          ENDIF
    332334
     
    354356             ii = MAXLOC( sum )
    355357             i = ii(1)
    356              IF ( sum(i) /= -1.0 )  THEN
     358             IF ( sum(i) /= -1.0_wp )  THEN
    357359                WRITE ( 18, 102 ) &
    358               log_point(i)%place, log_point(i)%sum,                &
    359                    log_point(i)%sum / log_point(1)%sum * 100.0,         &
     360              log_point(i)%place, log_point(i)%sum,                            &
     361                   log_point(i)%sum / log_point(1)%sum * 100.0_wp,             &
    360362                   log_point(i)%counts, pe_min(i), pe_max(i), pe_rms(i)
    361                 sum(i) = -1.0
     363                sum(i) = -1.0_wp
    362364             ELSE
    363365                EXIT
     
    391393          pe_min = log_point_s%sum
    392394          pe_max = log_point_s%sum    ! need to be set in case of 1 PE
    393           pe_rms = 0.0
     395          pe_rms = 0.0_wp
    394396
    395397#if defined( __parallel )
     
    416418!--          Calculate rms
    417419             DO  i = 0, numprocs-1
    418                 pe_rms(iii) = pe_rms(iii) + ( &
     420                pe_rms(iii) = pe_rms(iii) + (                                  &
    419421                                    pe_log_points(iii,i) - log_point_s(iii)%sum &
    420422                                            )**2
     
    441443             sum = log_point_s%sum
    442444          ELSEWHERE
    443              sum = -1.0
     445             sum = -1.0_wp
    444446          ENDWHERE
    445447
     
    450452             ii = MAXLOC( sum )
    451453             i = ii(1)
    452              IF ( sum(i) /= -1.0 )  THEN
    453                 WRITE ( 18, 102 ) &
    454                    log_point_s(i)%place, log_point_s(i)%sum, &
    455                    log_point_s(i)%sum / log_point(1)%sum * 100.0, &
     454             IF ( sum(i) /= -1.0_wp )  THEN
     455                WRITE ( 18, 102 )                                              &
     456                   log_point_s(i)%place, log_point_s(i)%sum,                   &
     457                   log_point_s(i)%sum / log_point(1)%sum * 100.0_wp,          &
    456458                   log_point_s(i)%counts, pe_min(i), pe_max(i), pe_rms(i)
    457                 sum(i) = -1.0
     459                sum(i) = -1.0_wp
    458460             ELSE
    459461                EXIT
Note: See TracChangeset for help on using the changeset viewer.