Changeset 189 for palm/trunk/SOURCE/data_output_spectra.f90
- Timestamp:
- Aug 13, 2008 5:09:26 PM (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/data_output_spectra.f90
r144 r189 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! allow 100 spectra levels instead of 10 for consistency with 7 ! define_netcdf_header 6 8 ! user-defined spectra 7 9 ! … … 270 272 INTEGER :: cranz, labforx = 3, labfory = 3, legpos = 3, & 271 273 timodex = 1 272 INTEGER, DIMENSION(1:10 ):: cucol = 1, klist = 999999, lstyle = 0274 INTEGER, DIMENSION(1:100):: cucol = 1, klist = 999999, lstyle = 0 273 275 LOGICAL :: datleg = .TRUE., grid = .TRUE., & 274 276 lclose = .TRUE., rand = .TRUE., & … … 278 280 REAL :: gwid = 0.1, rlegfak = 0.7, uxmin, uxmax, & 279 281 uymin, uymax 280 REAL, DIMENSION(1:10 ):: lwid = 0.6281 REAL, DIMENSION(10 ):: uyma, uymi282 REAL, DIMENSION(1:100) :: lwid = 0.6 283 REAL, DIMENSION(100) :: uyma, uymi 282 284 283 285 NAMELIST /RAHMEN/ cranz, datleg, rtext, swap … … 301 303 !-- different quantities u, v, w, pt and q) 302 304 DO k = 1, n_sp_x 303 IF ( k < 10 ) THEN305 IF ( k < 100 ) THEN 304 306 IF ( pr == 3 ) THEN 305 307 WRITE ( 82, 100 ) '#', k, header_char( pr ), & … … 351 353 !-- Calculate klist-values from the available comp_spectra_level values 352 354 i = 1; k = 1 353 DO WHILE ( i <= 10 .AND. plot_spectra_level(i) /= 999999 )355 DO WHILE ( i <= 100 .AND. plot_spectra_level(i) /= 999999 ) 354 356 DO WHILE ( k <= n_sp_x .AND. & 355 357 plot_spectra_level(i) >= comp_spectra_level(k) ) … … 398 400 100 FORMAT (A,I1,1X,A,1X,I4,'m ',A) 399 401 101 FORMAT (A,I2,1X,A,1X,I4,'m ',A) 400 102 FORMAT (E15.7,10 (1X,E15.7))402 102 FORMAT (E15.7,100(1X,E15.7)) 401 403 103 FORMAT ('NEXT') 402 404 104 FORMAT ('time averaged over',F7.1,' s') … … 429 431 INTEGER :: cranz, labforx = 3, labfory = 3, legpos = 3, & 430 432 timodex = 1 431 INTEGER, DIMENSION(1:10 ):: cucol = 1, klist = 999999, lstyle = 0433 INTEGER, DIMENSION(1:100):: cucol = 1, klist = 999999, lstyle = 0 432 434 LOGICAL :: datleg = .TRUE., grid = .TRUE., & 433 435 lclose = .TRUE., rand = .TRUE., & … … 437 439 REAL :: gwid = 0.1, rlegfak = 0.7, uxmin, uxmax, & 438 440 uymin, uymax 439 REAL, DIMENSION(1:10 ):: lwid = 0.6440 REAL, DIMENSION(10 ):: uyma, uymi441 REAL, DIMENSION(1:100) :: lwid = 0.6 442 REAL, DIMENSION(100) :: uyma, uymi 441 443 442 444 NAMELIST /RAHMEN/ cranz, datleg, rtext, swap … … 460 462 !-- different quantities u, v, w, pt and q) 461 463 DO k = 1, n_sp_y 462 IF ( k < 10 ) THEN464 IF ( k < 100 ) THEN 463 465 IF ( pr == 3 ) THEN 464 466 WRITE ( 84, 100 ) '#', k, header_char( pr ), & … … 510 512 !-- Calculate klist-values from the available comp_spectra_level values 511 513 j = 1; k = 1 512 DO WHILE ( j <= 10 .AND. plot_spectra_level(j) /= 999999 )514 DO WHILE ( j <= 100 .AND. plot_spectra_level(j) /= 999999 ) 513 515 DO WHILE ( k <= n_sp_y .AND. & 514 516 plot_spectra_level(j) >= comp_spectra_level(k) ) … … 557 559 100 FORMAT (A,I1,1X,A,1X,I4,'m ',A) 558 560 101 FORMAT (A,I2,1X,A,1X,I4,'m ',A) 559 102 FORMAT (E15.7,10 (1X,E15.7))561 102 FORMAT (E15.7,100(1X,E15.7)) 560 562 103 FORMAT ('NEXT') 561 563 104 FORMAT ('time averaged over',F7.1,' s')
Note: See TracChangeset
for help on using the changeset viewer.