- Timestamp:
- Aug 13, 2008 5:09:26 PM (16 years ago)
- Location:
- palm/trunk
- Files:
-
- 9 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/DOC/app/chapter_3.7.html
r108 r189 1 1 <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN"> 2 <html> 3 <head> 2 <html><head> 4 3 5 4 <meta http-equiv="CONTENT-TYPE" content="text/html; charset=windows-1252"> 6 <title>PALM chapter 3.7</title>5 7 6 8 < meta name="GENERATOR" content="StarOffice 7 (Win32)">7 <title>PALM chapter 3.7</title><meta name="GENERATOR" content="StarOffice 7 (Win32)"> 9 8 10 9 <meta name="AUTHOR" content="Siegfried Raasch"> … … 20 19 @page { size: 21cm 29.7cm } 21 20 --> 22 </style> 23 </head> 21 </style></head> 24 22 25 23 <body style="direction: ltr;" lang="en-US"> … … 121 119 </td> 122 120 <td> 123 <p>see section 4.2</p>121 <p>see section <a href="chapter_4.2.html#particle_parameters">4.2</a></p> 124 122 125 123 </td> … … 143 141 144 142 <td> 145 <p>see section 4.2</p>143 <p>see section <a href="chapter_4.2.html#dvrp_graphics">4.2</a></p> 146 144 </td> 147 145 </tr> … … 162 160 <td> 163 161 <p>see 164 section <a href="chapter_4.2.html#spectra _package">4.2</a></p>162 section <a href="chapter_4.2.html#spectra">4.2</a></p> 165 163 166 164 </td> … … 176 174 change: </i> $Id$</p> 177 175 178 </body> 179 </html> 176 </body></html> -
palm/trunk/DOC/app/chapter_4.2.html
r144 r189 11180 11180 11181 11181 11182 <td style="vertical-align: top;">I(10 )</td>11182 <td style="vertical-align: top;">I(100)</td> 11183 11183 11184 11184 … … 11202 11202 11203 11203 11204 Spectra can be calculated for up to tenlevels.</td>11204 Spectra can be calculated for up to 100 levels.</td> 11205 11205 11206 11206 … … 11417 11417 11418 11418 11419 <td style="vertical-align: top;">I(10 )</td>11419 <td style="vertical-align: top;">I(100)</td> 11420 11420 11421 11421 -
palm/trunk/DOC/app/chapter_4.6.html
r166 r189 65 65 --> 66 66 </style></head> 67 68 67 <body style="direction: ltr;" lang="en-US"> 69 68 … … 2021 2020 2022 2021 <p>I 2023 (10 )</p>2022 (100)</p> 2024 2023 2025 2024 … … 9494 9493 9495 9494 <p>I 9496 (10 )</p>9495 (100)</p> 9497 9496 9498 9497 -
palm/trunk/SOURCE/CURRENT_MODIFICATIONS
r187 r189 42 42 Changed: 43 43 ------- 44 Allow 100 spectra levels instead of 10 for consistency with 45 define_netcdf_header (calc_spectra, data_output_spectra, header, modules, 46 netcdf) 47 44 48 Modification of the integrated version of the profile function for momentum 45 49 for unstable stratification; more consistent flux definitions. (wall_fluxes, … … 83 87 (init_1d_model) 84 88 85 advec_s_ups, advec_u_ups, advec_v_ups, advec_w_ups, calc_spectra, check_open, check_parameters, cpu_statistics, init_1d_model, init_3d_model, modules, palm, parin, poisfft, production_e, read_var_list, read_3d_binary, transpose, wall_fluxes, write_var_list, write_3d_binary89 advec_s_ups, advec_u_ups, advec_v_ups, advec_w_ups, calc_spectra, check_open, check_parameters, cpu_statistics, data_output_spectra, header, init_1d_model, init_3d_model, modules, netcdf, palm, parin, poisfft, production_e, read_var_list, read_3d_binary, transpose, wall_fluxes, write_var_list, write_3d_binary 86 90 87 91 -
palm/trunk/SOURCE/calc_spectra.f90
r164 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, arguments removed from transpose routines 7 9 ! … … 223 225 REAL, DIMENSION(0:nx) :: work 224 226 REAL, DIMENSION(0:nx/2) :: sums_spectra_l 225 REAL, DIMENSION(0:nx/2,10 ):: sums_spectra227 REAL, DIMENSION(0:nx/2,100):: sums_spectra 226 228 227 229 REAL, DIMENSION(0:nxa,nys_x:nyn_xa,nzb_x:nzt_xa) :: ddd … … 234 236 !-- Loop over all levels defined by the user 235 237 n = 1 236 DO WHILE ( comp_spectra_level(n) /= 999999 .AND. n <= 10 )238 DO WHILE ( comp_spectra_level(n) /= 999999 .AND. n <= 100 ) 237 239 238 240 k = comp_spectra_level(n) … … 333 335 REAL, DIMENSION(0:ny) :: work 334 336 REAL, DIMENSION(0:ny/2) :: sums_spectra_l 335 REAL, DIMENSION(0:ny/2,10 ):: sums_spectra337 REAL, DIMENSION(0:ny/2,100):: sums_spectra 336 338 337 339 REAL, DIMENSION(0:nya,nxl_yd:nxr_yda,nzb_yd:nzt_yda) :: ddd … … 345 347 !-- Loop over all levels defined by the user 346 348 n = 1 347 DO WHILE ( comp_spectra_level(n) /= 999999 .AND. n <= 10 )349 DO WHILE ( comp_spectra_level(n) /= 999999 .AND. n <= 100 ) 348 350 349 351 k = comp_spectra_level(n) -
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') -
palm/trunk/SOURCE/header.f90
r167 r189 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! allow 100 spectra levels instead of 10 for consistency with 7 ! define_netcdf_header 6 8 ! Bugfix in the output of the characteristic levels of potential temperature, 7 9 ! geostrophic wind, scalar concentration, humidity and leaf area density … … 964 966 WRITE ( io, 372 ) ( data_output_sp(i), i = 1,10 ), & 965 967 ( spectra_direction(i), i = 1,10 ), & 966 ( comp_spectra_level(i), i = 1,10 ), &967 ( plot_spectra_level(i), i = 1,10 ), &968 ( comp_spectra_level(i), i = 1,100 ), & 969 ( plot_spectra_level(i), i = 1,100 ), & 968 970 averaging_interval_sp, dt_averaging_input_pr 969 971 ENDIF … … 1492 1494 372 FORMAT (' Arrays: ', 10(A5,',')/ & 1493 1495 ' Directions: ', 10(A5,',')/ & 1494 ' height levels k = ', 9(I3,','),I3,'.'/ & 1496 ' height levels k = ', 20(I3,',')/ & 1497 ' ', 20(I3,',')/ & 1498 ' ', 20(I3,',')/ & 1499 ' ', 20(I3,',')/ & 1500 ' ', 19(I3,','),I3,'.'/ & 1495 1501 ' height levels selected for standard plot:'/ & 1496 ' k = ', 9(I3,','),I3,'.'/ & 1502 ' k = ', 20(I3,',')/ & 1503 ' ', 20(I3,',')/ & 1504 ' ', 20(I3,',')/ & 1505 ' ', 20(I3,',')/ & 1506 ' ', 19(I3,','),I3,'.'/ & 1497 1507 ' Time averaged over ', F7.1, ' s,' / & 1498 1508 ' Profiles for the time averaging are taken every ', & -
palm/trunk/SOURCE/modules.f90
r163 r189 5 5 ! Actual revisions: 6 6 ! ----------------- 7 ! allow 100 spectra levels instead of 10 for consistency with 8 ! define_netcdf_header 7 9 ! +canopy_heat_flux, cthf, lai, 8 10 ! +leaf_surface_concentration, scalar_exchange_coefficient, sec, sls … … 1115 1117 INTEGER :: klist_x = 0, klist_y = 0, n_sp_x = 0, n_sp_y = 0 1116 1118 1117 INTEGER :: comp_spectra_level(10) = 999999, & 1118 lstyles(10) = (/ 0, 7, 3, 10, 1, 4, 9, 2, 6, 8 /), & 1119 plot_spectra_level(10) = 999999 1119 INTEGER :: comp_spectra_level(100) = 999999, & 1120 lstyles(100) = (/ 0, 7, 3, 10, 1, 4, 9, 2, 6, 8, & 1121 0, 7, 3, 10, 1, 4, 9, 2, 6, 8, & 1122 0, 7, 3, 10, 1, 4, 9, 2, 6, 8, & 1123 0, 7, 3, 10, 1, 4, 9, 2, 6, 8, & 1124 0, 7, 3, 10, 1, 4, 9, 2, 6, 8, & 1125 0, 7, 3, 10, 1, 4, 9, 2, 6, 8, & 1126 0, 7, 3, 10, 1, 4, 9, 2, 6, 8, & 1127 0, 7, 3, 10, 1, 4, 9, 2, 6, 8, & 1128 0, 7, 3, 10, 1, 4, 9, 2, 6, 8, & 1129 0, 7, 3, 10, 1, 4, 9, 2, 6, 8 /), & 1130 plot_spectra_level(100) = 999999 1120 1131 1121 1132 REAL :: time_to_start_sp = 0.0 -
palm/trunk/SOURCE/netcdf.f90
r175 r189 7 7 ! Current revisions: 8 8 ! ------------------ 9 ! consistently allow 100 spectra levels instead of 10 9 10 ! bug fix in the determination of the number of output heights for spectra 10 11 ! … … 3059 3060 !-- shall be output 3060 3061 ns = 1 3061 DO WHILE ( comp_spectra_level(ns) /= 999999 .AND. ns <= 10 )3062 DO WHILE ( comp_spectra_level(ns) /= 999999 .AND. ns <= 100 ) 3062 3063 ns = ns + 1 3063 3064 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.