Changeset 274 for palm/trunk/SOURCE
- Timestamp:
- Mar 26, 2009 3:11:21 PM (16 years ago)
- Location:
- palm/trunk/SOURCE
- Files:
-
- 27 edited
-
advec_particles.f90 (modified) (15 diffs)
-
calc_spectra.f90 (modified) (4 diffs)
-
check_for_restart.f90 (modified) (2 diffs)
-
check_open.f90 (modified) (13 diffs)
-
cpu_log.f90 (modified) (3 diffs)
-
data_output_2d.f90 (modified) (1 diff)
-
data_output_3d.f90 (modified) (1 diff)
-
data_output_dvrp.f90 (modified) (2 diffs)
-
data_output_profiles.f90 (modified) (2 diffs)
-
fft_xy.f90 (modified) (4 diffs)
-
flow_statistics.f90 (modified) (1 diff)
-
init_3d_model.f90 (modified) (1 diff)
-
init_dvrp.f90 (modified) (2 diffs)
-
init_grid.f90 (modified) (4 diffs)
-
init_particles.f90 (modified) (3 diffs)
-
init_pegrid.f90 (modified) (9 diffs)
-
netcdf.f90 (modified) (13 diffs)
-
parin.f90 (modified) (2 diffs)
-
poisfft_hybrid.f90 (modified) (1 diff)
-
read_3d_binary.f90 (modified) (7 diffs)
-
read_var_list.f90 (modified) (6 diffs)
-
surface_coupler.f90 (modified) (1 diff)
-
timestep.f90 (modified) (1 diff)
-
user_init_plant_canopy.f90 (modified) (2 diffs)
-
user_parin.f90 (modified) (1 diff)
-
user_read_restart_data.f90 (modified) (1 diff)
-
user_spectra.f90 (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/advec_particles.f90
r264 r274 319 319 ( rho_surface * dx * dy * dz ) 320 320 IF ( ql_c(k,j,i) > 100.0 ) THEN 321 WRITE( message_string, * ) 'k=',k,' j=',j,' i=',i, &321 WRITE( message_string, * ) 'k=',k,' j=',j,' i=',i, & 322 322 ' ql_c=',ql_c(k,j,i), ' &part(',n,')%wf=', & 323 323 particles(n)%weight_factor,' delta_r=',delta_r … … 329 329 IF ( ( new_r - particles(n)%radius ) < 0.0 .AND. new_r < 0.0 ) & 330 330 THEN 331 WRITE( message_string, * ) '#1 k=',k,' j=',j,' i=',i, &332 ' e_s=',e_s, ' e_a=',e_a,' t_int=',t_int, &331 WRITE( message_string, * ) '#1 k=',k,' j=',j,' i=',i, & 332 ' e_s=',e_s, ' e_a=',e_a,' t_int=',t_int, & 333 333 ' &d_radius=',d_radius,' delta_r=',delta_r,& 334 334 ' particle_radius=',particles(n)%radius … … 424 424 effective_coll_efficiency < 0.0 ) & 425 425 THEN 426 WRITE( message_string, * ) 'collision_efficiency ' , & 427 'out of range:' ,effective_coll_efficiency 428 CALL message( 'advec_particles', 'PA0145', 2, 2, 0, 6, 0 ) 426 WRITE( message_string, * ) 'collision_efficiency ' , & 427 'out of range:', effective_coll_efficiency 428 CALL message( 'advec_particles', 'PA0145', & 429 2, 2, 0, 6, 0 ) 429 430 END IF 430 431 … … 1920 1921 maximum_number_of_particles ) THEN 1921 1922 IF ( netcdf_output ) THEN 1922 message_string = 'maximum_number_of_particles ' // &1923 'needs to be increased ' // &1923 message_string = 'maximum_number_of_particles ' // & 1924 'needs to be increased ' // & 1924 1925 '&but this is not allowed with ' // & 1925 1926 'NetCDF output switched on' … … 1940 1941 maximum_number_of_tails ) THEN 1941 1942 IF ( netcdf_output ) THEN 1942 message_string = 'maximum_number_of_tails ' // &1943 'needs to be increased ' // &1943 message_string = 'maximum_number_of_tails ' // & 1944 'needs to be increased ' // & 1944 1945 '&but this is not allowed wi' // & 1945 1946 'th NetCDF output switched on' … … 2325 2326 THEN 2326 2327 IF ( netcdf_output ) THEN 2327 message_string = 'maximum_number_of_particles ' // &2328 'needs to be increased ' // &2328 message_string = 'maximum_number_of_particles ' // & 2329 'needs to be increased ' // & 2329 2330 '&but this is not allowed with ' // & 2330 2331 'NetCDF output switched on' … … 2353 2354 THEN 2354 2355 IF ( netcdf_output ) THEN 2355 message_string = 'maximum_number_of_tails ' // &2356 'needs to be increased ' // &2356 message_string = 'maximum_number_of_tails ' // & 2357 'needs to be increased ' // & 2357 2358 '&but this is not allowed wi'// & 2358 2359 'th NetCDF output switched on' … … 2402 2403 THEN 2403 2404 IF ( netcdf_output ) THEN 2404 message_string = 'maximum_number_of_particles ' // &2405 'needs to be increased ' // &2405 message_string = 'maximum_number_of_particles ' // & 2406 'needs to be increased ' // & 2406 2407 '&but this is not allowed with '// & 2407 2408 'NetCDF output switched on' … … 2430 2431 THEN 2431 2432 IF ( netcdf_output ) THEN 2432 message_string = 'maximum_number_of_tails ' // &2433 'needs to be increased ' // &2433 message_string = 'maximum_number_of_tails ' // & 2434 'needs to be increased ' // & 2434 2435 '&but this is not allowed wi'// & 2435 2436 'th NetCDF output switched on' … … 2766 2767 THEN 2767 2768 IF ( netcdf_output ) THEN 2768 message_string = 'maximum_number_of_particles ' // &2769 'needs to be increased ' // &2769 message_string = 'maximum_number_of_particles ' // & 2770 'needs to be increased ' // & 2770 2771 '&but this is not allowed with '// & 2771 2772 'NetCDF output switched on' … … 2794 2795 THEN 2795 2796 IF ( netcdf_output ) THEN 2796 message_string = 'maximum_number_of_tails ' // &2797 'needs to be increased ' // &2797 message_string = 'maximum_number_of_tails ' // & 2798 'needs to be increased ' // & 2798 2799 '&but this is not allowed wi' // & 2799 2800 'th NetCDF output switched on' … … 2843 2844 THEN 2844 2845 IF ( netcdf_output ) THEN 2845 message_string = 'maximum_number_of_particles ' // &2846 'needs to be increased ' // &2846 message_string = 'maximum_number_of_particles ' // & 2847 'needs to be increased ' // & 2847 2848 '&but this is not allowed with ' // & 2848 2849 'NetCDF output switched on' … … 2871 2872 THEN 2872 2873 IF ( netcdf_output ) THEN 2873 message_string = 'maximum_number_of_tails ' // &2874 'needs to be increased ' // &2874 message_string = 'maximum_number_of_tails ' // & 2875 'needs to be increased ' // & 2875 2876 '&but this is not allowed wi'// & 2876 2877 'th NetCDF output switched on' … … 3661 3662 CALL handle_netcdf_error( 'output_particles_netcdf', 12 ) 3662 3663 3663 nc_stat = NF90_PUT_VAR( id_set_prt, id_var_prt(11), particles%x, 3664 nc_stat = NF90_PUT_VAR( id_set_prt, id_var_prt(11), particles%x, & 3664 3665 start = (/ 1, prt_time_count /), & 3665 3666 count = (/ maximum_number_of_particles /) ) … … 3896 3897 k > nzt ) THEN 3897 3898 WRITE( message_string, * ) ' particle out of range: i=', i, ' j=', & 3898 j, ' k=', k, &3899 ' &nxl=', nxl, ' nxr=', nxr, &3900 ' nys=', nys, ' nyn=', nyn, &3899 j, ' k=', k, & 3900 ' &nxl=', nxl, ' nxr=', nxr, & 3901 ' nys=', nys, ' nyn=', nyn, & 3901 3902 ' nzb=', nzb, ' nzt=', nzt 3902 3903 CALL message( 'sort_particles', 'PA0149', 1, 2, 0, 6, 0 ) -
palm/trunk/SOURCE/calc_spectra.f90
r247 r274 85 85 IF ( bc_lr /= 'cyclic' ) THEN 86 86 87 message_string = 'non-cyclic lateral boundaries along x do not '// &87 message_string = 'non-cyclic lateral boundaries along x do not'// & 88 88 '& allow calculation of spectra along x' 89 89 CALL message( 'calc_spectra', 'PA0160', 1, 2, 0, 6, 0 ) … … 100 100 CALL calc_spectra_x( d, pr, m ) 101 101 #else 102 message_string = 'sorry, calculation of spectra in non parallel mode' // &103 ' & is still not realized'102 message_string = 'sorry, calculation of spectra in non parallel ' // & 103 'mode& is still not realized' 104 104 CALL message( 'calc_spectra', 'PA0161', 1, 2, 0, 6, 0 ) 105 105 #endif … … 116 116 IF ( bc_ns /= 'cyclic' ) THEN 117 117 IF ( myid == 0 ) THEN 118 message_string = 'non-cyclic lateral boundaries along y do not' // &119 ' & allow calculation of spectra along y'118 message_string = 'non-cyclic lateral boundaries along y do' // & 119 ' not & allow calculation of spectra along y' 120 120 CALL message( 'calc_spectra', 'PA0162', 1, 2, 0, 6, 0 ) 121 121 ENDIF … … 129 129 CALL calc_spectra_y( d, pr, m ) 130 130 #else 131 message_string = 'sorry, calculation of spectra in non parallel mode' // &132 ' & is still not realized'131 message_string = 'sorry, calculation of spectra in non parallel' // & 132 'mode& is still not realized' 133 133 CALL message( 'calc_spectra', 'PA0161', 1, 2, 0, 6, 0 ) 134 134 #endif -
palm/trunk/SOURCE/check_for_restart.f90
r249 r274 69 69 !-- Output that job will be terminated 70 70 IF ( terminate_run .AND. myid == 0 ) THEN 71 WRITE( message_string, * ) 'run will be terminated because it is running', &72 ' out of job cpu limit & ',&73 'remaining time: ', remaining_time, ' s',&74 71 WRITE( message_string, * ) 'run will be terminated because it is ', & 72 'running out of job cpu limit & ', & 73 'remaining time: ', remaining_time, ' s', & 74 'termination time needed:', termination_time_needed, ' s' 75 75 CALL message( 'check_for_restart', 'PA0163', 0, 1, 0, 6, 0 ) 76 76 ENDIF … … 113 113 114 114 WRITE( message_string, * ) 'run will be terminated due to user ', & 115 'settings of',&116 '&restart_time / dt_restart',&117 118 CALL message( 'check_for_restart', 'PA0164', 0, 0, 0, 6, 0 ) 115 'settings of', & 116 '&restart_time / dt_restart', & 117 '&new restart time is: ', time_restart, ' s' 118 CALL message( 'check_for_restart', 'PA0164', 0, 0, 0, 6, 0 ) 119 119 120 120 ! -
palm/trunk/SOURCE/check_open.f90
r263 r274 102 102 IF ( file_id == 14 .AND. openfile(file_id)%opened_before ) THEN 103 103 message_string = 're-open of unit ' // & 104 ' 104 '14 is not verified. Please check results!' 105 105 CALL message( 'check_open', 'PA0165', 0, 1, 0, 6, 0 ) 106 106 ENDIF … … 108 108 CASE DEFAULT 109 109 WRITE( message_string, * ) 're-opening of file-id ', file_id, & 110 110 ' is not allowed' 111 111 CALL message( 'check_open', 'PA0166', 0, 1, 0, 6, 0 ) 112 112 … … 144 144 !-- File-ids that are used temporarily in other routines 145 145 WRITE( message_string, * ) 'opening file-id ',file_id, & 146 146 ' is not allowed since it is used otherwise' 147 147 CALL message( 'check_open', 'PA0168', 0, 1, 0, 6, 0 ) 148 148 … … 405 405 IF ( .NOT. avs_coor_file_found .OR. & 406 406 .NOT. avs_data_file_found ) THEN 407 message_string= 'no filename for AVS-data-file ' // &408 'found in MRUN-config-file' // &407 message_string= 'no filename for AVS-data-file ' // & 408 'found in MRUN-config-file' // & 409 409 ' &filename in FLD-file set to "unknown"' 410 410 CALL message( 'check_open', 'PA0169', 0, 1, 0, 6, 0 ) … … 725 725 id_set_xy(av) ) 726 726 #else 727 message_string = 'NetCDF: no 64-bit offset allowed on this machine' 727 message_string = 'NetCDF: no 64-bit offset allowed ' // & 728 'on this machine' 728 729 CALL message( 'check_open', 'PA0171', 0, 1, 0, 6, 0 ) 729 730 … … 786 787 id_set_xz(av) ) 787 788 #else 788 message_string = 'NetCDF: no 64-bit offset allowed on this machine' 789 message_string = 'NetCDF: no 64-bit offset allowed ' // & 790 'on this machine' 789 791 CALL message( 'check_open', 'PA0171', 0, 1, 0, 6, 0 ) 790 792 … … 847 849 id_set_yz(av) ) 848 850 #else 849 message_string = 'NetCDF: no 64-bit offset allowed on this machine' 851 message_string = 'NetCDF: no 64-bit offset allowed ' // & 852 'on this machine' 850 853 CALL message( 'check_open', 'PA0171', 0, 1, 0, 6, 0 ) 851 854 … … 902 905 id_set_pr ) 903 906 #else 904 message_string = 'NetCDF: no 64-bit offset allowed on this machine' 907 message_string = 'NetCDF: no 64-bit offset allowed ' // & 908 'on this machine' 905 909 CALL message( 'check_open', 'PA0171', 0, 1, 0, 6, 0 ) 906 910 … … 957 961 id_set_ts ) 958 962 #else 959 message_string = 'NetCDF: no 64-bit offset allowed on this machine' 963 message_string = 'NetCDF: no 64-bit offset allowed ' // & 964 'on this machine' 960 965 CALL message( 'check_open', 'PA0171', 0, 1, 0, 6, 0 ) 961 966 … … 1019 1024 id_set_3d(av) ) 1020 1025 #else 1021 message_string = 'NetCDF: no 64-bit offset allowed on this machine' 1026 message_string = 'NetCDF: no 64-bit offset allowed ' // & 1027 'on this machine' 1022 1028 CALL message( 'check_open', 'PA0171', 0, 1, 0, 6, 0 ) 1023 1029 … … 1076 1082 id_set_sp ) 1077 1083 #else 1078 message_string = 'NetCDF: no 64-bit offset allowed on this machine' 1084 message_string = 'NetCDF: no 64-bit offset allowed ' // & 1085 'on this machine' 1079 1086 CALL message( 'check_open', 'PA0171', 0, 1, 0, 6, 0 ) 1080 1087 … … 1152 1159 id_set_prt ) 1153 1160 #else 1154 message_string = 'NetCDF: no 64-bit offset allowed on this machine' 1161 message_string = 'NetCDF: no 64-bit offset allowed ' // & 1162 'on this machine' 1155 1163 CALL message( 'check_open', 'PA0171', 0, 1, 0, 6, 0 ) 1156 1164 … … 1208 1216 id_set_pts ) 1209 1217 #else 1210 message_string = 'NetCDF: no 64-bit offset allowed on this machine' 1218 message_string = 'NetCDF: no 64-bit offset allowed ' // & 1219 'on this machine' 1211 1220 CALL message( 'check_open', 'PA0171', 0, 1, 0, 6, 0 ) 1212 1221 -
palm/trunk/SOURCE/cpu_log.f90
r254 r274 82 82 ELSEIF ( modus == 'pause' ) THEN 83 83 IF ( ( mtime - log_event%mtime ) < 0.0 .AND. first ) THEN 84 WRITE( message_string, * ) 'negative time interval occured', &85 86 84 WRITE( message_string, * ) 'negative time interval occured', & 85 ' &PE',myid,' L=PAUSE "',TRIM(log_event%place),'" new=', & 86 mtime,' last=',log_event%mtime 87 87 CALL message( 'cpu_log', 'PA0176', 0, 1, -1, 6, 0 ) 88 88 first = .FALSE. … … 93 93 IF ( ( mtime - log_event%mtime + log_event%isum ) < 0.0 .AND. & 94 94 first ) THEN 95 WRITE( message_string, * ) 'negative time interval occured', &96 97 95 WRITE( message_string, * ) 'negative time interval occured', & 96 ' &PE',myid,' L=STOP "',TRIM(log_event%place),'" new=', & 97 mtime,' last=',log_event%mtime,' isum=',log_event%isum 98 98 CALL message( 'cpu_log', 'PA0177', 0, 1, -1, 6, 0 ) 99 99 first = .FALSE. … … 103 103 log_event%sum = log_event%sum + log_event%mtime 104 104 IF ( log_event%sum < 0.0 .AND. first ) THEN 105 WRITE( message_string, * ) 'negative time interval occured', &106 105 WRITE( message_string, * ) 'negative time interval occured', & 106 ' &PE',myid,' L=STOP "',TRIM(log_event%place),'" sum=', & 107 107 log_event%sum,' mtime=',log_event%mtime 108 108 CALL message( 'cpu_log', 'PA0178', 0, 1, -1, 6, 0 ) -
palm/trunk/SOURCE/data_output_2d.f90
r263 r274 520 520 521 521 IF ( .NOT. found ) THEN 522 message_string = 'no output provided for: ' // TRIM( do2d(av,if) ) 522 message_string = 'no output provided for: ' // & 523 TRIM( do2d(av,if) ) 523 524 CALL message( 'data_output_2d', 'PA0181', 0, 0, 0, 6, 0 ) 524 525 ENDIF -
palm/trunk/SOURCE/data_output_3d.f90
r263 r274 313 313 314 314 IF ( .NOT. found ) THEN 315 message_string = 'no output available for: ' // TRIM( do3d(av,if) ) 315 message_string = 'no output available for: ' // & 316 TRIM( do3d(av,if) ) 316 317 CALL message( 'data_output_3d', 'PA0182', 0, 0, 0, 6, 0 ) 317 318 ENDIF -
palm/trunk/SOURCE/data_output_dvrp.f90
r264 r274 401 401 ENDDO 402 402 ELSE 403 message_string = 'if humidity/passive_scalar = ' // & 404 'FALSE output of ' // TRIM( output_variable ) // & 405 'is not provided' 406 CALL message( 'data_output_dvrp', 'PA0183', 0, 0, 0, 6, 0 ) 403 message_string = 'if humidity/passive_scalar = ' // & 404 'FALSE output of ' // TRIM( output_variable ) // & 405 'is not provided' 406 CALL message( 'data_output_dvrp', 'PA0183',& 407 0, 0, 0, 6, 0 ) 407 408 ENDIF 408 409 … … 417 418 ENDDO 418 419 ELSE 419 message_string = 'if cloud_physics = FALSE ' // & 420 'output of ' // TRIM( output_variable) // & 421 'is not provided' 422 CALL message( 'data_output_dvrp', 'PA0184', 0, 0, 0, 6, 0 ) 420 message_string = 'if cloud_physics = FALSE ' // & 421 'output of ' // TRIM( output_variable) // & 422 'is not provided' 423 CALL message( 'data_output_dvrp', 'PA0184',& 424 0, 0, 0, 6, 0 ) 423 425 ENDIF 424 426 -
palm/trunk/SOURCE/data_output_profiles.f90
r263 r274 289 289 IF ( cross_normx_factor(k,j) == 0.0 .OR. & 290 290 cross_normy_factor(k,j) == 0.0 ) THEN 291 WRITE( message_string, * ) 'data_output_profiles: normalizi', & 292 'ng cross ',j, ' is not possible since one o', & 293 'f the & normalizing factors is zero! & ', & 294 'cross_normx_factor(',k,',',j,') = ', & 295 cross_normx_factor(k,j), & 296 ' & cross_normy_factor(',k,',',j,') = ', & 297 cross_normy_factor(k,j) 298 CALL message( 'data_output_profiles', 'PA0185', 0, 1, 0, 6, 0 ) 291 WRITE( message_string, * ) 'data_output_profiles: ', & 292 'normalizing cross ',j, ' is not possible ', & 293 'since one of the & normalizing factors ', & 294 'is zero! & cross_normx_factor(',k,',',j,') = ', & 295 cross_normx_factor(k,j), & 296 ' & cross_normy_factor(',k,',',j,') = ', & 297 cross_normy_factor(k,j) 298 CALL message( 'data_output_profiles', 'PA0185',& 299 0, 1, 0, 6, 0 ) 299 300 cross_normx_factor(k,j) = 1.0 300 301 cross_normy_factor(k,j) = 1.0 … … 534 535 IF ( cross_normx_factor(k,j) == 0.0 .OR. & 535 536 cross_normy_factor(k,j) == 0.0 ) THEN 536 WRITE( message_string, * ) 'data_output_profiles: normalizi', & 537 'ng cross ',j, ' is not possible since one o', & 538 'f the & normalizing factors is zero! & ', & 539 'cross_normx_factor(',k,',',j,') = ', & 540 cross_normx_factor(k,j), & 541 ' & cross_normy_factor(',k,',',j,') = ', & 542 cross_normy_factor(k,j) 543 CALL message( 'data_output_profiles', 'PA0185', 0, 1, 0, 6, 0 ) 537 WRITE( message_string, * ) 'data_output_profiles: ', & 538 'normalizing cross ',j, ' is not possible ', & 539 'since one of the & normalizing factors ', & 540 'is zero! & cross_normx_factor(',k,',',j,') = ', & 541 cross_normx_factor(k,j), & 542 ' & cross_normy_factor(',k,',',j,') = ', & 543 cross_normy_factor(k,j) 544 CALL message( 'data_output_profiles', 'PA0185',& 545 0, 1, 0, 6, 0 ) 544 546 cross_normx_factor(k,j) = 1.0 545 547 cross_normy_factor(k,j) = 1.0 -
palm/trunk/SOURCE/fft_xy.f90
r254 r274 361 361 #endif 362 362 ELSE 363 message_string = 'fft method "' // TRIM( fft_method) // '" not available' 363 message_string = 'fft method "' // TRIM( fft_method) // & 364 '" not available' 364 365 CALL message( 'fft_x', 'PA0189', 1, 2, 0, 6, 0 ) 365 366 … … 543 544 ELSE 544 545 545 message_string = 'fft method "' // TRIM( fft_method) // '" not available' 546 message_string = 'fft method "' // TRIM( fft_method) // & 547 '" not available' 546 548 CALL message( 'fft_y', 'PA0189', 1, 2, 0, 6, 0 ) 547 549 … … 678 680 ELSE 679 681 680 message_string = 'fft method "' // TRIM( fft_method) // '" not available' 682 message_string = 'fft method "' // TRIM( fft_method) // & 683 '" not available' 681 684 CALL message( 'fft_x_m', 'PA0189', 1, 2, 0, 6, 0 ) 682 685 … … 813 816 ELSE 814 817 815 message_string = 'fft method "' // TRIM( fft_method) // '" not available' 818 message_string = 'fft method "' // TRIM( fft_method) // & 819 '" not available' 816 820 CALL message( 'fft_x_m', 'PA0189', 1, 2, 0, 6, 0 ) 817 821 -
palm/trunk/SOURCE/flow_statistics.f90
r254 r274 95 95 IF ( flow_statistics_called ) THEN 96 96 97 message_string = 'flow_statistics is called two times within one timestep' 97 message_string = 'flow_statistics is called two times within one ' // & 98 'timestep' 98 99 CALL message( 'flow_statistics', 'PA0190', 1, 2, 0, 6, 0 ) 99 100 -
palm/trunk/SOURCE/init_3d_model.f90
r254 r274 1379 1379 IF ( dots_num > dots_max ) THEN 1380 1380 WRITE( message_string, * ) 'number of time series quantities exceeds', & 1381 ' its maximum of dots_max = ', dots_max, &1381 ' its maximum of dots_max = ', dots_max, & 1382 1382 ' &Please increase dots_max in modules.f90.' 1383 1383 CALL message( 'init_3d_model', 'PA0194', 1, 2, 0, 6, 0 ) -
palm/trunk/SOURCE/init_dvrp.f90
r273 r274 145 145 IF ( dvrp_output /= 'ftp' .AND. dvrp_output /= 'rtsp' .AND. & 146 146 dvrp_output /= 'local' ) THEN 147 message_string = 'dvrp_output="' // TRIM( dvrp_output ) // '" not allowed' 147 message_string = 'dvrp_output="' // TRIM( dvrp_output ) // & 148 '" not allowed' 148 149 CALL message( 'init_dvrp', 'PA0196', 1, 2, 0, 6, 0 ) 149 150 ENDIF … … 189 190 mode_dvrp(m)(1:9) /= 'pathlines' ) THEN 190 191 191 message_string = 'mode_dvrp="' // TRIM( mode_dvrp(m) ) // '" not allowed' 192 message_string = 'mode_dvrp="' // TRIM( mode_dvrp(m) ) // & 193 '" not allowed' 192 194 CALL message( 'init_dvrp', 'PA0198', 1, 2, 0, 6, 0 ) 193 195 CALL local_stop -
palm/trunk/SOURCE/init_grid.f90
r254 r274 367 367 IF ( ( bxl < 1 ) .OR. ( bxr > nx-1 ) .OR. ( bxr < bxl+3 ) .OR. & 368 368 ( bys < 1 ) .OR. ( byn > ny-1 ) .OR. ( byn < bys+3 ) ) THEN 369 WRITE( message_string, * ) 'inconsistent building parameters:', &370 371 369 WRITE( message_string, * ) 'inconsistent building parameters:', & 370 '& bxl=', bxl, 'bxr=', bxr, 'bys=', bys, & 371 'byn=', byn, 'nx=', nx, 'ny=', ny 372 372 CALL message( 'init_grid', 'PA0203', 1, 2, 0, 6, 0 ) 373 373 ENDIF … … 417 417 ( ch < 3 ) ) THEN 418 418 WRITE( message_string, * ) 'inconsistent canyon parameters:', & 419 '&cxl=', cxl, 'cxr=', cxr, &420 'cwx=', cwx, &419 '&cxl=', cxl, 'cxr=', cxr, & 420 'cwx=', cwx, & 421 421 'ch=', ch, 'nx=', nx, 'ny=', ny 422 422 CALL message( 'init_grid', 'PA0205', 1, 2, 0, 6, 0 ) … … 426 426 ( ch < 3 ) ) THEN 427 427 WRITE( message_string, * ) 'inconsistent canyon parameters:', & 428 '&cys=', cys, 'cyn=', cyn, &429 'cwy=', cwy, &428 '&cys=', cys, 'cyn=', cyn, & 429 'cwy=', cwy, & 430 430 'ch=', ch, 'nx=', nx, 'ny=', ny 431 431 CALL message( 'init_grid', 'PA0206', 1, 2, 0, 6, 0 ) 432 432 ENDIF 433 433 ENDIF 434 IF ( canyon_width_x /= 9999999.9 .AND. canyon_width_y /= 9999999.9 ) 434 IF ( canyon_width_x /= 9999999.9 .AND. canyon_width_y /= 9999999.9 ) & 435 435 THEN 436 message_string = 'inconsistent canyon parameters:' // &436 message_string = 'inconsistent canyon parameters:' // & 437 437 '&street canyon can only be oriented' // & 438 438 '&either in x- or in y-direction' … … 507 507 !-- Consistency checks 508 508 IF ( MINVAL( nzb_local ) < 0 .OR. MAXVAL( nzb_local ) > nz + 1 ) THEN 509 WRITE( message_string, * ) 'nzb_local values are outside the', &510 'model domain',&511 512 509 WRITE( message_string, * ) 'nzb_local values are outside the', & 510 'model domain', & 511 '&MINVAL( nzb_local ) = ', MINVAL(nzb_local), & 512 '&MAXVAL( nzb_local ) = ', MAXVAL(nzb_local) 513 513 CALL message( 'init_grid', 'PA0210', 1, 2, 0, 6, 0 ) 514 514 ENDIF -
palm/trunk/SOURCE/init_particles.f90
r254 r274 98 98 !-- Check the number of particle groups. 99 99 IF ( number_of_particle_groups > max_number_of_particle_groups ) THEN 100 WRITE( message_string, * ) 'max_number_of_particle_groups =', &101 max_number_of_particle_groups , &100 WRITE( message_string, * ) 'max_number_of_particle_groups =', & 101 max_number_of_particle_groups , & 102 102 '&number_of_particle_groups reset to ', & 103 103 max_number_of_particle_groups … … 153 153 particle_binary_version = '3.0' 154 154 IF ( TRIM( version_on_file ) /= TRIM( particle_binary_version ) ) THEN 155 message_string = 'version mismatch concerning data from prior run' // & 156 '&version on file = "' // TRIM( version_on_file ) // & 157 '&version in program = "' // & 158 TRIM( particle_binary_version ) // '"' 155 message_string = 'version mismatch concerning data from prior ' // & 156 'run &version on file = "' // & 157 TRIM( version_on_file ) // & 158 '&version in program = "' // & 159 TRIM( particle_binary_version ) // '"' 159 160 CALL message( 'init_particles', 'PA0214', 1, 2, 0, 6, 0 ) 160 161 ENDIF … … 269 270 IF ( n > maximum_number_of_particles ) THEN 270 271 WRITE( message_string, * ) 'number of initial', & 271 'particles (', n, ') exceeds',&272 '&maximum_number_of_particles (',&273 272 'particles (', n, ') exceeds', & 273 '&maximum_number_of_particles (', & 274 maximum_number_of_particles, ') on PE ', & 274 275 myid 275 CALL message( 'init_particles', 'PA0216', 2, 2, 0, 6, 0 ) 276 CALL message( 'init_particles', 'PA0216', & 277 2, 2, 0, 6, 0 ) 276 278 ENDIF 277 279 particles(n)%x = pos_x -
palm/trunk/SOURCE/init_pegrid.f90
r254 r274 122 122 !-- must be equal to the number of PEs available to the job 123 123 IF ( ( npex * npey ) /= numprocs ) THEN 124 WRITE( message_string, * ) 'number of PEs of the prescribed topology (', &125 npex*npey,') does not match & the number of ',&126 124 WRITE( message_string, * ) 'number of PEs of the prescribed ', & 125 'topology (', npex*npey,') does not match & the number of ', & 126 'PEs available to the job (', numprocs, ')' 127 127 CALL message( 'init_pegrid', 'PA0221', 1, 2, 0, 6, 0 ) 128 128 ENDIF … … 134 134 !-- If the processor topology is prescribed by the user, the number of 135 135 !-- PEs must be given in both directions 136 message_string = 'if the processor topology is prescribed by the user, ' //&137 '& both values of "npex" and "npey" must be given in the ' //&138 '&NAMELIST-parameter file'136 message_string = 'if the processor topology is prescribed by the, ' // & 137 ' user& both values of "npex" and "npey" must be given ' // & 138 'in the &NAMELIST-parameter file' 139 139 CALL message( 'init_pegrid', 'PA0222', 1, 2, 0, 6, 0 ) 140 140 … … 235 235 236 236 IF ( MOD( nxa+1 , pdims(1) ) /= 0 ) THEN 237 WRITE( message_string, * ) 'x-direction: gridpoint number (',nx+1,') is not an',&238 '& integral divisor of the number of proces',&239 'sors (', pdims(1),')'237 WRITE( message_string, * ) 'x-direction: gridpoint number (',nx+1,') ',& 238 'is not an& integral divisor of the number ', & 239 'processors (', pdims(1),')' 240 240 CALL message( 'init_pegrid', 'PA0225', 1, 2, 0, 6, 0 ) 241 241 ELSE 242 242 nnx = ( nxa + 1 ) / pdims(1) 243 243 IF ( nnx*pdims(1) - ( nx + 1) > nnx ) THEN 244 WRITE( message_string, * ) 'x-direction: nx does not match the', & 245 'requirements given by the number of PEs', & 246 '& used',& 247 '& please use nx = ', nx - ( pdims(1) - ( nnx*pdims(1) & 248 - ( nx + 1 ) ) ), ' instead of nx =', nx 244 WRITE( message_string, * ) 'x-direction: nx does not match the', & 245 'requirements given by the number of PEs &used', & 246 '& please use nx = ', nx - ( pdims(1) - ( nnx*pdims(1) & 247 - ( nx + 1 ) ) ), ' instead of nx =', nx 249 248 CALL message( 'init_pegrid', 'PA0226', 1, 2, 0, 6, 0 ) 250 249 ENDIF … … 262 261 !-- Calculate array bounds in y-direction for every PE. 263 262 IF ( MOD( nya+1 , pdims(2) ) /= 0 ) THEN 264 WRITE( message_string, * ) 'y-direction: gridpoint number (',ny+1,') is not an', &265 '& integral divisor of the number of proces',&266 'sors (', pdims(2),')'263 WRITE( message_string, * ) 'y-direction: gridpoint number (',ny+1,') ', & 264 'is not an& integral divisor of the number of', & 265 'processors (', pdims(2),')' 267 266 CALL message( 'init_pegrid', 'PA0227', 1, 2, 0, 6, 0 ) 268 267 ELSE 269 268 nny = ( nya + 1 ) / pdims(2) 270 269 IF ( nny*pdims(2) - ( ny + 1) > nny ) THEN 271 WRITE( message_string, * ) 'y-direction: ny does not match the',& 272 'requirements given by the number of PEs', & 273 '& used', & 274 '& please use ny = ', ny - ( pdims(2) - ( nnx*pdims(2) & 270 WRITE( message_string, * ) 'y-direction: ny does not match the', & 271 'requirements given by the number of PEs &used ', & 272 '& please use ny = ', ny - ( pdims(2) - ( nnx*pdims(2) & 275 273 - ( ny + 1 ) ) ), ' instead of ny =', ny 276 274 CALL message( 'init_pegrid', 'PA0228', 1, 2, 0, 6, 0 ) … … 324 322 nny_x = nny 325 323 IF ( MOD( nza , pdims(1) ) /= 0 ) THEN 326 WRITE( message_string, * ) 'transposition z --> x:', &327 328 pdims(1)324 WRITE( message_string, * ) 'transposition z --> x:', & 325 '&nz=',nz,' is not an integral divisior of pdims(1)=', & 326 pdims(1) 329 327 CALL message( 'init_pegrid', 'PA0230', 1, 2, 0, 6, 0 ) 330 328 ENDIF … … 356 354 nzt_y = nzt_x 357 355 IF ( MOD( nxa+1 , pdims(2) ) /= 0 ) THEN 358 WRITE( message_string, * ) 'transposition x --> y:', &359 360 356 WRITE( message_string, * ) 'transposition x --> y:', & 357 '&nx+1=',nx+1,' is not an integral divisor of ',& 358 'pdims(2)=',pdims(2) 361 359 CALL message( 'init_pegrid', 'PA0231', 1, 2, 0, 6, 0 ) 362 360 ENDIF … … 382 380 nxr_z = nxr_y 383 381 IF ( MOD( nya+1 , pdims(1) ) /= 0 ) THEN 384 WRITE( message_string, * ) 'transposition y --> z:', &385 386 382 WRITE( message_string, * ) 'transposition y --> z:', & 383 '& ny+1=',ny+1,' is not an integral divisor of ',& 384 'pdims(1)=',pdims(1) 387 385 CALL message( 'init_pegrid', 'PA0232', 1, 2, 0, 6, 0 ) 388 386 ENDIF … … 398 396 !-- x --> y. This condition must be fulfilled for a 1D-decomposition along x 399 397 IF ( MOD( nya+1 , pdims(1) ) /= 0 ) THEN 400 WRITE( message_string, * ) 'transposition x --> y:', &401 402 398 WRITE( message_string, * ) 'transposition x --> y:', & 399 '& ny+1=',ny+1,' is not an integral divisor of ',& 400 'pdims(1)=',pdims(1) 403 401 CALL message( 'init_pegrid', 'PA0233', 1, 2, 0, 6, 0 ) 404 402 ENDIF … … 410 408 IF ( dt_dosp /= 9999999.9 ) THEN 411 409 IF ( MOD( nza, pdims(2) ) /= 0 ) THEN 412 WRITE( message_string, * ) 'direct transposition z --> y (needed for spectra):', &413 '& nz=',nz,' is not an integral divisor of ',&414 410 WRITE( message_string, * ) 'direct transposition z --> y (needed ', & 411 'for spectra):& nz=',nz,' is not an integral divisor of ',& 412 'pdims(2)=',pdims(2) 415 413 CALL message( 'init_pegrid', 'PA0234', 1, 2, 0, 6, 0 ) 416 414 ELSE -
palm/trunk/SOURCE/netcdf.f90
r263 r274 110 110 j = INDEX( netcdf_precision(i), '_' ) 111 111 IF ( j == 0 ) THEN 112 WRITE ( message_string, * ) 'netcdf_precision must contain a "_"', &113 ' netcdf_precision(', i, ')="',&112 WRITE ( message_string, * ) 'netcdf_precision must contain a ', & 113 '"_"netcdf_precision(', i, ')="', & 114 114 TRIM( netcdf_precision(i) ),'"' 115 115 CALL message( 'define_netcdf_header', 'PA0241', 1, 2, 0, 6, 0 ) … … 153 153 154 154 CASE DEFAULT 155 WRITE ( message_string, * ) 'unknown variable in inipar assignment:',&156 'netcdf_precision(', i, ')="',&155 WRITE ( message_string, * ) 'unknown variable in inipar ', & 156 'assignment: netcdf_precision(', i, ')="', & 157 157 TRIM( netcdf_precision(i) ),'"' 158 CALL message( 'define_netcdf_header', 'PA0243', 1, 2, 0, 6, 0 ) 158 CALL message( 'define_netcdf_header', 'PA0243', 1, 2, 0, 6, 0 ) 159 159 160 160 END SELECT … … 418 418 WRITE ( message_string, * ) 'no grid defined for', & 419 419 ' variable ', TRIM( do3d(av,i) ) 420 CALL message( 'define_netcdf_header', 'PA0244', 0, 1, 0, 6, 0 ) 420 CALL message( 'define_netcdf_header', 'PA0244', 0, 1, 0, 6, 0 ) 421 421 ENDIF 422 422 … … 950 950 WRITE ( message_string, * ) 'no grid defined for', & 951 951 ' variable ', TRIM( do2d(av,i) ) 952 CALL message( 'define_netcdf_header', 'PA0244', 0, 1, 0, 6, 0 ) 952 CALL message( 'define_netcdf_header', 'PA0244',& 953 0, 1, 0, 6, 0 ) 953 954 ENDIF 954 955 … … 1189 1190 IF ( section(i,1) /= -1 ) THEN 1190 1191 IF ( zu(section(i,1)) /= netcdf_data(i) ) THEN 1191 message_string = 'NetCDF file for cross-sections ' // & 1192 TRIM( var ) // ' from previuos run found,' // & 1193 '&but this file cannot be extended' // & 1194 ' due to mismatch in cross' // & 1195 '§ion levels.' // & 1196 '&New file is created instead.' 1197 CALL message( 'define_netcdf_header', 'PA0251', 0, 1, 0, 6, 0 ) 1192 message_string = 'NetCDF file for cross-sections ' // & 1193 TRIM( var ) // ' from previuos run found,' // & 1194 '&but this file cannot be extended' // & 1195 ' due to mismatch in cross' // & 1196 '§ion levels.' // & 1197 '&New file is created instead.' 1198 CALL message( 'define_netcdf_header', 'PA0251', & 1199 0, 1, 0, 6, 0 ) 1198 1200 extend = .FALSE. 1199 1201 RETURN … … 1201 1203 ELSE 1202 1204 IF ( -1.0 /= netcdf_data(i) ) THEN 1203 message_string = 'NetCDF file for cross-sections ' // & 1204 TRIM( var ) // ' from previuos run found,' // & 1205 '&but this file cannot be extended' // & 1206 ' due to mismatch in cross' // & 1207 '§ion levels.' // & 1208 '&New file is created instead.' 1209 CALL message( 'define_netcdf_header', 'PA0251', 0, 1, 0, 6, 0 ) 1205 message_string = 'NetCDF file for cross-sections ' // & 1206 TRIM( var ) // ' from previuos run found,' // & 1207 '&but this file cannot be extended' // & 1208 ' due to mismatch in cross' // & 1209 '§ion levels.' // & 1210 '&New file is created instead.' 1211 CALL message( 'define_netcdf_header', 'PA0251',& 1212 0, 1, 0, 6, 0 ) 1210 1213 extend = .FALSE. 1211 1214 RETURN … … 1501 1504 WRITE ( message_string, * ) 'no grid defined for', & 1502 1505 ' variable ', TRIM( do2d(av,i) ) 1503 CALL message( 'define_netcdf_header', 'PA0244', 0, 1, 0, 6, 0 ) 1506 CALL message( 'define_netcdf_header', 'PA0244',& 1507 0, 1, 0, 6, 0 ) 1504 1508 ENDIF 1505 1509 … … 1708 1712 IF ( section(i,2) /= -1 ) THEN 1709 1713 IF ( ( section(i,2) * dy ) /= netcdf_data(i) ) THEN 1710 message_string = 'NetCDF file for cross-sections ' // & 1711 TRIM( var ) // ' from previuos run found,' // & 1712 '&but this file cannot be extended' // & 1713 ' due to mismatch in cross' // & 1714 '§ion levels.' // & 1715 '&New file is created instead.' 1716 CALL message( 'define_netcdf_header', 'PA0251', 0, 1, 0, 6, 0 ) 1714 message_string = 'NetCDF file for cross-sections ' // & 1715 TRIM( var ) // ' from previuos run found,' // & 1716 '&but this file cannot be extended' // & 1717 ' due to mismatch in cross' // & 1718 '§ion levels.' // & 1719 '&New file is created instead.' 1720 CALL message( 'define_netcdf_header', 'PA0251', & 1721 0, 1, 0, 6, 0 ) 1717 1722 extend = .FALSE. 1718 1723 RETURN … … 1720 1725 ELSE 1721 1726 IF ( -1.0 /= netcdf_data(i) ) THEN 1722 message_string = 'NetCDF file for cross-sections ' // & 1723 TRIM( var ) // ' from previuos run found,' // & 1724 '&but this file cannot be extended' // & 1725 ' due to mismatch in cross' // & 1726 '§ion levels.' // & 1727 '&New file is created instead.' 1728 CALL message( 'define_netcdf_header', 'PA0251', 0, 1, 0, 6, 0 ) 1727 message_string = 'NetCDF file for cross-sections ' // & 1728 TRIM( var ) // ' from previuos run found,' // & 1729 '&but this file cannot be extended' // & 1730 ' due to mismatch in cross' // & 1731 '§ion levels.' // & 1732 '&New file is created instead.' 1733 CALL message( 'define_netcdf_header', 'PA0251', & 1734 0, 1, 0, 6, 0 ) 1729 1735 extend = .FALSE. 1730 1736 RETURN … … 2020 2026 WRITE ( message_string, * ) 'no grid defined for', & 2021 2027 ' variable ', TRIM( do2d(av,i) ) 2022 CALL message( 'define_netcdf_header', 'PA0244', 0, 1, 0, 6, 0 ) 2028 CALL message( 'define_netcdf_header', 'PA0244',& 2029 0, 1, 0, 6, 0 ) 2023 2030 ENDIF 2024 2031 … … 2225 2232 IF ( section(i,3) /= -1 ) THEN 2226 2233 IF ( ( section(i,3) * dx ) /= netcdf_data(i) ) THEN 2227 message_string = 'NetCDF file for cross-sections ' // & 2228 TRIM( var ) // ' from previuos run found,' // & 2229 '&but this file cannot be extended' // & 2230 ' due to mismatch in cross' // & 2231 '§ion levels.' // & 2232 '&New file is created instead.' 2233 CALL message( 'define_netcdf_header', 'PA0251', 0, 1, 0, 6, 0 ) 2234 message_string = 'NetCDF file for cross-sections ' // & 2235 TRIM( var ) // ' from previuos run found,' // & 2236 '&but this file cannot be extended' // & 2237 ' due to mismatch in cross' // & 2238 '§ion levels.' // & 2239 '&New file is created instead.' 2240 CALL message( 'define_netcdf_header', 'PA0251', & 2241 0, 1, 0, 6, 0 ) 2234 2242 extend = .FALSE. 2235 2243 RETURN … … 2237 2245 ELSE 2238 2246 IF ( -1.0 /= netcdf_data(i) ) THEN 2239 message_string = 'NetCDF file for cross-sections ' // & 2240 TRIM( var ) // ' from previuos run found,' // & 2241 '&but this file cannot be extended' // & 2242 ' due to mismatch in cross' // & 2243 '§ion levels.' // & 2244 '&New file is created instead.' 2245 CALL message( 'define_netcdf_header', 'PA0251', 0, 1, 0, 6, 0 ) 2247 message_string = 'NetCDF file for cross-sections ' // & 2248 TRIM( var ) // ' from previuos run found,' // & 2249 '&but this file cannot be extended' // & 2250 ' due to mismatch in cross' // & 2251 '§ion levels.' // & 2252 '&New file is created instead.' 2253 CALL message( 'define_netcdf_header', 'PA0251', & 2254 0, 1, 0, 6, 0 ) 2246 2255 extend = .FALSE. 2247 2256 RETURN … … 3216 3225 IF ( averaging_interval_sp /= 0.0 ) THEN 3217 3226 WRITE (time_average_text,'('', '',F7.1,'' s average'')') & 3218 3227 averaging_interval_sp 3219 3228 nc_stat = NF90_PUT_ATT( id_set_sp, NF90_GLOBAL, 'title', & 3220 3229 TRIM( run_description_header ) // & -
palm/trunk/SOURCE/parin.f90
r260 r274 179 179 READ ( 11, inipar, ERR=10, END=11 ) 180 180 GOTO 12 181 10 message_string = 'errors in \$inipar &or no \$inipar-namelist found (CRAY-machines only)' 181 10 message_string = 'errors in \$inipar &or no \$inipar-namelist ' // & 182 'found (CRAY-machines only)' 182 183 CALL message( 'parin', 'PA0271', 1, 2, 0, 6, 0 ) 183 184 … … 229 230 230 231 IF ( nx <= 0 ) THEN 231 WRITE( message_string, * ) 'no value or wrong value given for nx: nx=', nx 232 WRITE( message_string, * ) 'no value or wrong value given for nx: ',& 233 'nx=', nx 232 234 CALL message( 'parin', 'PA0273', 1, 2, 0, 6, 0 ) 233 235 ENDIF 234 236 IF ( ny <= 0 ) THEN 235 WRITE( message_string, * ) 'no value or wrong value given for ny: ny=', ny 237 WRITE( message_string, * ) 'no value or wrong value given for ny: ',& 238 'ny=', ny 236 239 CALL message( 'parin', 'PA0274', 1, 2, 0, 6, 0 ) 237 240 ENDIF 238 241 IF ( nz <= 0 ) THEN 239 WRITE( message_string, * ) 'no value or wrong value given for nz: nz=', nz 242 WRITE( message_string, * ) 'no value or wrong value given for nz: ',& 243 'nz=', nz 240 244 CALL message( 'parin', 'PA0275', 1, 2, 0, 6, 0 ) 241 245 ENDIF -
palm/trunk/SOURCE/poisfft_hybrid.f90
r257 r274 187 187 'tasks_per_logical_node = ', & 188 188 tasks_per_logical_node 189 CALL message( 'poisfft_hybrid_ini', 'PA0281', 0, 0, 0, 6, 0 ) 189 CALL message( 'poisfft_hybrid_ini', 'PA0281', 0, 0, 0, 6, 0 ) 190 190 191 191 ENDIF -
palm/trunk/SOURCE/read_3d_binary.f90
r257 r274 204 204 IF ( files_to_be_opened /= 1 .OR. numprocs /= numprocs_previous_run ) & 205 205 THEN 206 WRITE( message_string, * ) 'number of PEs or virtual PE-grid changed in ',&207 'restart run& PE', myid, ' will read from files ', &208 206 WRITE( message_string, * ) 'number of PEs or virtual PE-grid changed ', & 207 'in restart run& PE', myid, ' will read from files ', & 208 file_list(1:files_to_be_opened) 209 209 CALL message( 'read_3d_binary', 'PA0285', 0, 0, 0, 6, 0 ) 210 210 ENDIF … … 237 237 binary_version = '3.1' 238 238 IF ( TRIM( version_on_file ) /= TRIM( binary_version ) ) THEN 239 WRITE( message_string, * ) 'version mismatch concerning data ', 240 'from prior run',&241 '&version on file = "', TRIM( version_on_file ), '"',&242 239 WRITE( message_string, * ) 'version mismatch concerning data ', & 240 'from prior run', & 241 '&version on file = "', TRIM( version_on_file ), '"', & 242 '&version in program = "', TRIM( binary_version ), '"' 243 243 CALL message( 'read_3d_binary', 'PA0286', 1, 2, 0, 6, 0 ) 244 244 ENDIF … … 251 251 252 252 IF ( nxl_on_file /= hor_index_bounds_previous_run(1,j) ) THEN 253 WRITE( message_string, * ) 'problem with index bound nxl on ', 254 'restart file "', myid_char, '"',&255 '&nxl = ', nxl_on_file, ' but it should be',&256 '&= ', hor_index_bounds_previous_run(1,j),&257 253 WRITE( message_string, * ) 'problem with index bound nxl on ', & 254 'restart file "', myid_char, '"', & 255 '&nxl = ', nxl_on_file, ' but it should be', & 256 '&= ', hor_index_bounds_previous_run(1,j), & 257 '&from the index bound information array' 258 258 CALL message( 'read_3d_binary', 'PA0287', 2, 2, 0, 6, 0 ) 259 259 ENDIF 260 260 261 261 IF ( nxr_on_file /= hor_index_bounds_previous_run(2,j) ) THEN 262 WRITE( message_string, * ) 'problem with index bound nxr on ', 263 264 265 266 262 WRITE( message_string, * ) 'problem with index bound nxr on ', & 263 'restart file "', myid_char, '"' , & 264 '&nxr = ', nxr_on_file, ' but it should be', & 265 '&= ', hor_index_bounds_previous_run(2,j), & 266 '&from the index bound information array' 267 267 CALL message( 'read_3d_binary', 'PA0288', 2, 2, 0, 6, 0 ) 268 268 … … 270 270 271 271 IF ( nys_on_file /= hor_index_bounds_previous_run(3,j) ) THEN 272 WRITE( message_string, * ) 'problem with index bound nys on ', 273 274 275 272 WRITE( message_string, * ) 'problem with index bound nys on ', & 273 'restart file "', myid_char, '"', & 274 '&nys = ', nys_on_file, ' but it should be', & 275 '&= ', hor_index_bounds_previous_run(3,j), & 276 276 '&from the index bound information array' 277 277 CALL message( 'read_3d_binary', 'PA0289', 2, 2, 0, 6, 0 ) … … 279 279 280 280 IF ( nyn_on_file /= hor_index_bounds_previous_run(4,j) ) THEN 281 WRITE( message_string, * ) 'problem with index bound nyn on ', 282 283 284 285 281 WRITE( message_string, * ) 'problem with index bound nyn on ', & 282 'restart file "', myid_char, '"', & 283 '&nyn = ', nyn_on_file, ' but it should be', & 284 '&= ', hor_index_bounds_previous_run(4,j), & 285 '&from the index bound information array' 286 286 CALL message( 'read_3d_binary', 'PA0290', 2, 2, 0, 6, 0 ) 287 287 ENDIF … … 614 614 IF ( k == 1 ) THEN 615 615 IF ( nx_on_file /= nx ) THEN 616 message_string = 'read_3d_binary: spectrum_x ' // & 617 'on restart file ignored because' // & 618 '&total numbers of grid points (nx) ' // & 619 'do not match' 620 CALL message( 'read_3d_binary', 'PA0293', 0, 1, 0, 6, 0 ) 616 message_string = 'read_3d_binary: spectrum_x ' // & 617 'on restart file ignored because' // & 618 '&total numbers of grid points (nx) ' // & 619 'do not match' 620 CALL message( 'read_3d_binary', 'PA0293',& 621 0, 1, 0, 6, 0 ) 621 622 READ ( 13 ) rdummy 622 623 ELSE … … 628 629 IF ( k == 1 ) THEN 629 630 IF ( ny_on_file /= ny ) THEN 630 message_string = 'read_3d_binary: spectrum_y ' // & 631 'on restart file ignored because' // & 632 '&total numbers of grid points (ny) '// & 633 'do not match' 634 CALL message( 'read_3d_binary', 'PA0294', 0, 1, 0, 6, 0 ) 635 READ ( 13 ) rdummy 631 message_string = 'read_3d_binary: spectrum_y ' // & 632 'on restart file ignored because' // & 633 '&total numbers of grid points (ny) '// & 634 'do not match' 635 CALL message( 'read_3d_binary', 'PA0294', & 636 0, 1, 0, 6, 0 ) 637 READ ( 13 ) rdummy 636 638 ELSE 637 639 READ ( 13 ) spectrum_y -
palm/trunk/SOURCE/read_var_list.f90
r257 r274 99 99 READ ( 13 ) variable_chr 100 100 IF ( TRIM( variable_chr ) /= 'numprocs' ) THEN 101 WRITE( message_string, * ) 'numprocs not found in data from prior run', &102 ' on PE ', myid101 WRITE( message_string, * ) 'numprocs not found in data from prior ', & 102 'run on PE ', myid 103 103 CALL message( 'read_var_list', 'PA0297', 1, 2, 0, 6, 0 ) 104 104 ENDIF … … 633 633 READ ( 13 ) variable_chr 634 634 IF ( TRIM( variable_chr ) /= 'numprocs' ) THEN 635 WRITE( message_string, * ) 'numprocs not found in data from prior run', &636 ' on PE ', myid635 WRITE( message_string, * ) 'numprocs not found in data from prior ', & 636 'run on PE ', myid 637 637 CALL message( 'read_parts_of_var_list', 'PA0297', 1, 2, 0, 6, 0 ) 638 638 ENDIF … … 663 663 IF ( nz_on_file /= nz ) THEN 664 664 WRITE( message_string, * ) 'mismatch concerning number of ', & 665 'gridpoints along z', &666 '&nz on file = "', nz_on_file, '"', &667 '&nz from run = "', nz, '"'665 'gridpoints along z', & 666 '&nz on file = "', nz_on_file, '"', & 667 '&nz from run = "', nz, '"' 668 668 CALL message( 'read_parts_of_var_list', 'PA0304', 1, 2, 0, 6, 0 ) 669 669 ENDIF … … 676 676 READ ( 13 ) max_pr_user_on_file 677 677 IF ( max_pr_user_on_file /= max_pr_user ) THEN 678 WRITE( message_string, * ) 'number of user profiles on res', &679 'tart data file differs from the current run', &680 ' &max_pr_user on file = "',&681 max_pr_user_on_file, '"', &682 '&max_pr_user from run = "', &678 WRITE( message_string, * ) 'number of user profiles on res', & 679 'tart data file differs from the current ', & 680 'run&max_pr_user on file = "', & 681 max_pr_user_on_file, '"', & 682 '&max_pr_user from run = "', & 683 683 max_pr_user, '"' 684 684 CALL message( 'read_parts_of_var_list', 'PA0306', 0, 0, 0, 6, 0 ) … … 695 695 READ ( 13 ) statistic_regions_on_file 696 696 IF ( statistic_regions_on_file /= statistic_regions ) THEN 697 WRITE( message_string, * ) 'statistic regions on restart data file differ',&698 ' from the current run',&699 '&statistic regions on file = "', &700 statistic_regions_on_file, '"', &701 '&statistic regions from run = "', &702 statistic_regions, '"', &697 WRITE( message_string, * ) 'statistic regions on restart data file ',& 698 'differ from the current run', & 699 '&statistic regions on file = "', & 700 statistic_regions_on_file, '"', & 701 '&statistic regions from run = "', & 702 statistic_regions, '"', & 703 703 '&statistic data may be lost!' 704 704 CALL message( 'read_parts_of_var_list', 'PA0308', 0, 1, 0, 6, 0 ) … … 725 725 READ ( 13 ) average_count_pr 726 726 IF ( average_count_pr /= 0 ) THEN 727 WRITE( message_string, * ) 'inflow profiles not temporally averaged.',& 728 '&Averaging will be done now using ', & 729 average_count_pr, ' samples.' 730 CALL message( 'read_parts_of_var_list', 'PA0309', 0, 1, 0, 6, 0 ) 727 WRITE( message_string, * ) 'inflow profiles not temporally ', & 728 'averaged. &Averaging will be done now using ', & 729 average_count_pr, ' samples.' 730 CALL message( 'read_parts_of_var_list', 'PA0309', & 731 0, 1, 0, 6, 0 ) 731 732 ENDIF 732 733 -
palm/trunk/SOURCE/surface_coupler.f90
r258 r274 53 53 0, comm_inter, status, ierr ) 54 54 IF ( terminate_coupled_remote > 0 ) THEN 55 WRITE( message_string, * ) 'remote model "', TRIM( coupling_mode_remote ), & 56 '" terminated', & 57 '&with terminate_coupled_remote = ', & 58 terminate_coupled_remote, & 59 '&local model "', TRIM( coupling_mode ), & 60 '" has', & 61 '&terminate_coupled = ', & 55 WRITE( message_string, * ) 'remote model "', & 56 TRIM( coupling_mode_remote ), & 57 '" terminated', & 58 '&with terminate_coupled_remote = ', & 59 terminate_coupled_remote, & 60 '&local model "', TRIM( coupling_mode ), & 61 '" has', & 62 '&terminate_coupled = ', & 62 63 terminate_coupled 63 64 CALL message( 'surface_coupler', 'PA0310', 1, 2, 0, 6, 0 ) -
palm/trunk/SOURCE/timestep.f90
r258 r274 202 202 stop_dt = .TRUE. 203 203 204 WRITE( message_string, * ) 'Time step has reached minimum limit.', 205 206 207 208 209 210 211 212 213 214 215 216 204 WRITE( message_string, * ) 'Time step has reached minimum limit.', & 205 '&dt = ', dt_3d, ' s Simulation is terminated.', & 206 '&old_dt = ', old_dt, ' s', & 207 '&dt_u = ', dt_u, ' s', & 208 '&dt_v = ', dt_v, ' s', & 209 '&dt_w = ', dt_w, ' s', & 210 '&dt_diff = ', dt_diff, ' s', & 211 '&u_max = ', u_max, ' m/s k=', u_max_ijk(1), & 212 ' j=', u_max_ijk(2), ' i=', u_max_ijk(3), & 213 '&v_max = ', v_max, ' m/s k=', v_max_ijk(1), & 214 ' j=', v_max_ijk(2), ' i=', v_max_ijk(3), & 215 '&w_max = ', w_max, ' m/s k=', w_max_ijk(1), & 216 ' j=', w_max_ijk(2), ' i=', w_max_ijk(3) 217 217 CALL message( 'timestep', 'PA0312', 0, 1, 0, 6, 0 ) 218 218 ! -
palm/trunk/SOURCE/user_init_plant_canopy.f90
r258 r274 60 60 !-- After definition, please 61 61 !-- remove the following three lines! 62 message_string = 'canopy_mode "' // canopy_mode // '" not available yet' 62 message_string = 'canopy_mode "' // canopy_mode // & 63 '" not available yet' 63 64 CALL message( 'user_init_plant_canopy', 'UI0007', 0, 1, 0, 6, 0 ) 64 65 … … 66 67 ! 67 68 !-- The DEFAULT case is reached if the parameter canopy_mode contains a 68 !-- wrong character string that is neither recognized in init_3d_model nor 69 !-- wrong character string that is neither recognized in init_3d_model nor 69 70 !-- here in user_init_plant_canopy. 70 71 message_string = 'unknown canopy_mode "' // canopy_mode // '"' -
palm/trunk/SOURCE/user_parin.f90
r258 r274 75 75 IF ( TRIM( initializing_actions ) == 'read_restart_data' ) THEN 76 76 IF ( max_pr_user /= max_pr_user_tmp ) THEN 77 WRITE( message_string, * ) 'the number of user-defined profiles given in', &78 '&data_output_pr (', max_pr_user_tmp, ') doe',&79 'snot match the one ',&80 '&found in the restart file (', max_pr_user,&77 WRITE( message_string, * ) 'the number of user-defined profiles ', & 78 'given in &data_output_pr (', max_pr_user_tmp, ') doe', & 79 'snot match the one ', & 80 '&found in the restart file (', max_pr_user, & 81 81 ')' 82 82 CALL message( 'user_parin', 'UI0009', 1, 2, 0, 6, 0 ) -
palm/trunk/SOURCE/user_read_restart_data.f90
r258 r274 82 82 ! 83 83 ! CASE DEFAULT 84 ! WRITE( message_string, * ) 'unknown variable named "', 85 ! 86 ! 84 ! WRITE( message_string, * ) 'unknown variable named "', & 85 ! TRIM( field_char ), '" found in', & 86 ! '&data from prior run on PE ', myid 87 87 ! CALL message( 'user_spectra', 'UI0012', 1, 2, 0, 6, 0 ) 88 88 ! -
palm/trunk/SOURCE/user_spectra.f90
r258 r274 54 54 CASE DEFAULT 55 55 message_string = 'Spectra of ' // & 56 56 TRIM( data_output_sp(m) ) // ' can not be calculated' 57 57 CALL message( 'user_spectra', 'UI0010', 0, 0, 0, 6, 0 ) 58 58
Note: See TracChangeset
for help on using the changeset viewer.
