Changeset 1353 for palm/trunk/SOURCE/data_output_spectra.f90
- Timestamp:
- Apr 8, 2014 3:21:23 PM (9 years ago)
- File:
-
- 1 edited
-
palm/trunk/SOURCE/data_output_spectra.f90 (modified) (14 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/data_output_spectra.f90
r1329 r1353 20 20 ! Current revisions: 21 21 ! ------------------ 22 ! 22 ! REAL constants provided with KIND-attribute 23 23 ! 24 24 ! Former revisions: … … 183 183 ! 184 184 !-- Reset spectra values 185 spectrum_x = 0.0 ; spectrum_y = 0.0185 spectrum_x = 0.0_wp; spectrum_y = 0.0_wp 186 186 187 187 ENDIF … … 242 242 243 243 DO i = 1, nx/2 244 frequency = 2.0 * pi * i / ( dx * ( nx + 1 ) )244 frequency = 2.0_wp * pi * i / ( dx * ( nx + 1 ) ) 245 245 netcdf_data_x(i) = frequency * spectrum_x(i,k,nsp) 246 246 ENDDO … … 260 260 261 261 DO i = 1, ny/2 262 frequency = 2.0 * pi * i / ( dy * ( ny + 1 ) )262 frequency = 2.0_wp * pi * i / ( dy * ( ny + 1 ) ) 263 263 netcdf_data_y(i) = frequency * spectrum_y(i,k,nsp) 264 264 ENDDO … … 317 317 LOGICAL :: frame_written !: 318 318 319 REAL(wp) :: frequency = 0.0 !:319 REAL(wp) :: frequency = 0.0_wp !: 320 320 ! 321 321 !-- Variables needed for PROFIL-namelist … … 344 344 LOGICAL :: ylog = .TRUE. !: 345 345 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 !: 352 352 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 !: 356 356 357 357 NAMELIST /RAHMEN/ cranz, datleg, rtext, swap … … 406 406 !-- Write all data and calculate uymi and uyma. They serve to calculate 407 407 !-- the CROSS-parameters uymin and uymax 408 uymi = 999.999 ; uyma = -999.999408 uymi = 999.999_wp; uyma = -999.999_wp 409 409 DO i = 1, nx/2 410 frequency = 2.0 * pi * i / ( dx * ( nx + 1 ) )410 frequency = 2.0_wp * pi * i / ( dx * ( nx + 1 ) ) 411 411 WRITE ( 82, 102 ) frequency, ( frequency * spectrum_x(i,k,m), k = 1, & 412 412 n_sp_x ) … … 431 431 klist(i) = k + klist_x 432 432 ELSE 433 uymi(k) = 999.999 434 uyma(k) = -999.999 433 uymi(k) = 999.999_wp 434 uyma(k) = -999.999_wp 435 435 ENDIF 436 436 k = k + 1 … … 438 438 i = i + 1 439 439 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 442 442 utext = 'x'//utext_char( pr ) 443 IF ( averaging_interval_sp /= 0.0 ) THEN443 IF ( averaging_interval_sp /= 0.0_wp ) THEN 444 444 WRITE ( atext, 104 ) averaging_interval_sp 445 445 utext = TRIM(utext) // ', ' // TRIM( atext ) 446 446 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 ) ) 451 451 ytext = ytext_char( pr ) 452 452 … … 520 520 LOGICAL :: frame_written !: 521 521 522 REAL(wp) :: frequency = 0.0 !:522 REAL(wp) :: frequency = 0.0_wp !: 523 523 524 524 ! … … 548 548 LOGICAL :: ylog = .TRUE. !: 549 549 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 !: 561 561 562 562 NAMELIST /RAHMEN/ cranz, datleg, rtext, swap … … 611 611 !-- Write all data and calculate uymi and uyma. They serve to calculate 612 612 !-- the CROSS-parameters uymin and uymax 613 uymi = 999.999 ; uyma = -999.999613 uymi = 999.999_wp; uyma = -999.999_wp 614 614 DO j = 1, ny/2 615 frequency = 2.0 * pi * j / ( dy * ( ny + 1 ) )615 frequency = 2.0_wp * pi * j / ( dy * ( ny + 1 ) ) 616 616 WRITE ( 84, 102 ) frequency, ( frequency * spectrum_y(j,k,m), & 617 617 k = 1, n_sp_y ) … … 636 636 klist(j) = k + klist_y 637 637 ELSE 638 uymi(k) = 999.999 639 uyma(k) = -999.999 638 uymi(k) = 999.999_wp 639 uyma(k) = -999.999_wp 640 640 ENDIF 641 641 k = k + 1 … … 643 643 j = j + 1 644 644 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 647 647 utext = 'y'//utext_char( pr ) 648 IF ( averaging_interval_sp /= 0.0 ) THEN648 IF ( averaging_interval_sp /= 0.0_wp ) THEN 649 649 WRITE ( atext, 104 ) averaging_interval_sp 650 650 utext = TRIM(utext) // ', ' // TRIM( atext ) 651 651 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 ) ) 656 656 ytext = ytext_char( pr ) 657 657
Note: See TracChangeset
for help on using the changeset viewer.
