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

REAL constants provided with KIND-attribute

File:
1 edited

Legend:

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

    r1329 r1353  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! REAL constants provided with KIND-attribute
    2323!
    2424! Former revisions:
     
    183183!
    184184!--    Reset spectra values
    185        spectrum_x = 0.0; spectrum_y = 0.0
     185       spectrum_x = 0.0_wp; spectrum_y = 0.0_wp
    186186
    187187    ENDIF
     
    242242
    243243          DO  i = 1, nx/2
    244              frequency = 2.0 * pi * i / ( dx * ( nx + 1 ) )
     244             frequency = 2.0_wp * pi * i / ( dx * ( nx + 1 ) )
    245245             netcdf_data_x(i) = frequency * spectrum_x(i,k,nsp)
    246246          ENDDO
     
    260260
    261261          DO  i = 1, ny/2
    262              frequency = 2.0 * pi * i / ( dy * ( ny + 1 ) )
     262             frequency = 2.0_wp * pi * i / ( dy * ( ny + 1 ) )
    263263             netcdf_data_y(i) = frequency * spectrum_y(i,k,nsp)
    264264          ENDDO
     
    317317    LOGICAL            ::  frame_written   !:
    318318   
    319     REAL(wp)           ::  frequency = 0.0 !:
     319    REAL(wp)           ::  frequency = 0.0_wp !:
    320320!
    321321!-- Variables needed for PROFIL-namelist
     
    344344    LOGICAL ::  ylog = .TRUE.   !:
    345345   
    346     REAL(wp) ::  gwid = 0.1    !:
    347     REAL(wp) ::  rlegfak = 0.7 !:
    348     REAL(wp) ::  uxmin         !:
    349     REAL(wp) ::  uxmax         !:
    350     REAL(wp) ::  uymin         !:
    351     REAL(wp) ::  uymax         !:
     346    REAL(wp) ::  gwid = 0.1_wp    !:
     347    REAL(wp) ::  rlegfak = 0.7_wp !:
     348    REAL(wp) ::  uxmin            !:
     349    REAL(wp) ::  uxmax            !:
     350    REAL(wp) ::  uymin            !:
     351    REAL(wp) ::  uymax            !:
    352352     
    353     REAL(wp), DIMENSION(1:100) ::  lwid = 0.6 !:
    354     REAL(wp), DIMENSION(100)   ::  uyma       !:
    355     REAL(wp), DIMENSION(100)   ::  uymi       !:
     353    REAL(wp), DIMENSION(1:100) ::  lwid = 0.6_wp !:
     354    REAL(wp), DIMENSION(100)   ::  uyma          !:
     355    REAL(wp), DIMENSION(100)   ::  uymi          !:
    356356
    357357    NAMELIST /RAHMEN/  cranz, datleg, rtext, swap
     
    406406!-- Write all data and calculate uymi and uyma. They serve to calculate
    407407!-- the CROSS-parameters uymin and uymax
    408     uymi = 999.999; uyma = -999.999
     408    uymi = 999.999_wp; uyma = -999.999_wp
    409409    DO  i = 1, nx/2
    410        frequency = 2.0 * pi * i / ( dx * ( nx + 1 ) )
     410       frequency = 2.0_wp * pi * i / ( dx * ( nx + 1 ) )
    411411       WRITE ( 82, 102 )  frequency, ( frequency * spectrum_x(i,k,m), k = 1, &
    412412                          n_sp_x )
     
    431431             klist(i) = k + klist_x
    432432          ELSE
    433              uymi(k) =  999.999
    434              uyma(k) = -999.999
     433             uymi(k) =  999.999_wp
     434             uyma(k) = -999.999_wp
    435435          ENDIF
    436436          k = k + 1
     
    438438       i = i + 1
    439439    ENDDO
    440     uymi(k:n_sp_x) =  999.999
    441     uyma(k:n_sp_x) = -999.999
     440    uymi(k:n_sp_x) =  999.999_wp
     441    uyma(k:n_sp_x) = -999.999_wp
    442442    utext = 'x'//utext_char( pr )
    443     IF ( averaging_interval_sp /= 0.0 ) THEN
     443    IF ( averaging_interval_sp /= 0.0_wp ) THEN
    444444       WRITE ( atext, 104 )  averaging_interval_sp
    445445       utext = TRIM(utext) // ',  ' // TRIM( atext )
    446446    ENDIF
    447     uxmin = 0.8 * 2.0 * pi        / ( dx * ( nx + 1 ) )
    448     uxmax = 1.2 * 2.0 * pi * nx/2 / ( dx * ( nx + 1 ) )
    449     uymin = 0.8 * MIN (  999.999, MINVAL ( uymi ) )
    450     uymax = 1.2 * MAX ( -999.999, MAXVAL ( uyma ) )
     447    uxmin = 0.8_wp * 2.0_wp * pi        / ( dx * ( nx + 1 ) )
     448    uxmax = 1.2_wp * 2.0_wp * pi * nx/2 / ( dx * ( nx + 1 ) )
     449    uymin = 0.8_wp * MIN (  999.999_wp, MINVAL ( uymi ) )
     450    uymax = 1.2_wp * MAX ( -999.999_wp, MAXVAL ( uyma ) )
    451451    ytext = ytext_char( pr )
    452452
     
    520520    LOGICAL            :: frame_written   !:
    521521   
    522     REAL(wp)           :: frequency = 0.0 !:
     522    REAL(wp)           :: frequency = 0.0_wp !:
    523523
    524524!
     
    548548    LOGICAL ::  ylog = .TRUE.   !:
    549549   
    550     REAL(wp) ::  gwid = 0.1     !:
    551     REAL(wp) ::  rlegfak = 0.7  !:
    552     REAL(wp) ::  uxmin          !:
    553     REAL(wp) ::  uxmax          !:
    554     REAL(wp) ::  uymin          !:
    555     REAL(wp) ::  uymax          !:
    556    
    557     REAL(wp), DIMENSION(1:100) ::  lwid = 0.6 !:
    558    
    559     REAL(wp), DIMENSION(100)   ::  uyma       !:
    560     REAL(wp), DIMENSION(100)   ::  uymi       !:
     550    REAL(wp) ::  gwid = 0.1_wp     !:
     551    REAL(wp) ::  rlegfak = 0.7_wp  !:
     552    REAL(wp) ::  uxmin             !:
     553    REAL(wp) ::  uxmax             !:
     554    REAL(wp) ::  uymin             !:
     555    REAL(wp) ::  uymax             !:
     556   
     557    REAL(wp), DIMENSION(1:100) ::  lwid = 0.6_wp !:
     558   
     559    REAL(wp), DIMENSION(100)   ::  uyma          !:
     560    REAL(wp), DIMENSION(100)   ::  uymi          !:
    561561
    562562    NAMELIST /RAHMEN/  cranz, datleg, rtext, swap
     
    611611!-- Write all data and calculate uymi and uyma. They serve to calculate
    612612!-- the CROSS-parameters uymin and uymax
    613     uymi = 999.999; uyma = -999.999
     613    uymi = 999.999_wp; uyma = -999.999_wp
    614614    DO  j = 1, ny/2
    615        frequency = 2.0 * pi * j / ( dy * ( ny + 1 ) )
     615       frequency = 2.0_wp * pi * j / ( dy * ( ny + 1 ) )
    616616       WRITE ( 84, 102 ) frequency, ( frequency * spectrum_y(j,k,m), &
    617617                                      k = 1, n_sp_y )
     
    636636             klist(j) = k + klist_y
    637637          ELSE
    638              uymi(k) =  999.999
    639              uyma(k) = -999.999
     638             uymi(k) =  999.999_wp
     639             uyma(k) = -999.999_wp
    640640          ENDIF
    641641          k = k + 1
     
    643643       j = j + 1
    644644    ENDDO
    645     uymi(k:n_sp_y) =  999.999
    646     uyma(k:n_sp_y) = -999.999
     645    uymi(k:n_sp_y) =  999.999_wp
     646    uyma(k:n_sp_y) = -999.999_wp
    647647    utext = 'y'//utext_char( pr )
    648     IF ( averaging_interval_sp /= 0.0 )  THEN
     648    IF ( averaging_interval_sp /= 0.0_wp )  THEN
    649649       WRITE ( atext, 104 )  averaging_interval_sp
    650650       utext = TRIM(utext) // ',  ' // TRIM( atext )
    651651    ENDIF
    652     uxmin = 0.8 * 2.0 * pi        / ( dy * ( ny + 1 ) )
    653     uxmax = 1.2 * 2.0 * pi * ny/2 / ( dy * ( ny + 1 ) )
    654     uymin = 0.8 * MIN (  999.999, MINVAL ( uymi ) )
    655     uymax = 1.2 * MAX ( -999.999, MAXVAL ( uyma ) )
     652    uxmin = 0.8_wp * 2.0_wp * pi        / ( dy * ( ny + 1 ) )
     653    uxmax = 1.2_wp * 2.0_wp * pi * ny/2 / ( dy * ( ny + 1 ) )
     654    uymin = 0.8_wp * MIN (  999.999_wp, MINVAL ( uymi ) )
     655    uymax = 1.2_wp * MAX ( -999.999_wp, MAXVAL ( uyma ) )
    656656    ytext = ytext_char( pr )
    657657
Note: See TracChangeset for help on using the changeset viewer.