Changeset 274 for palm/trunk/SOURCE/advec_particles.f90
- Timestamp:
- Mar 26, 2009 3:11:21 PM (15 years ago)
- File:
-
- 1 edited
-
palm/trunk/SOURCE/advec_particles.f90 (modified) (15 diffs)
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 )
Note: See TracChangeset
for help on using the changeset viewer.
