Changeset 1327 for palm/trunk/SOURCE/lpm_data_output_particles.f90
- Timestamp:
- Mar 21, 2014 11:00:16 AM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/lpm_data_output_particles.f90
r1321 r1327 20 20 ! Current revisions: 21 21 ! ------------------ 22 ! -netcdf output queries 22 23 ! 23 24 ! Former revisions: … … 43 44 44 45 USE control_parameters, & 45 ONLY: netcdf_output,prt_time_count, simulated_time46 ONLY: prt_time_count, simulated_time 46 47 47 48 USE cpulog, & … … 80 81 ! 81 82 !-- Output in netCDF format 82 IF ( netcdf_output ) THEN 83 84 CALL check_open( 108 ) 83 CALL check_open( 108 ) 85 84 86 85 ! 87 !-- 88 86 !-- Update the NetCDF time axis 87 prt_time_count = prt_time_count + 1 89 88 90 91 92 93 89 nc_stat = NF90_PUT_VAR( id_set_prt, id_var_time_prt, & 90 (/ simulated_time /), & 91 start = (/ prt_time_count /), count = (/ 1 /) ) 92 CALL handle_netcdf_error( 'lpm_data_output_particles', 1 ) 94 93 95 94 ! 96 !-- 97 98 99 100 95 !-- Output the real number of particles used 96 nc_stat = NF90_PUT_VAR( id_set_prt, id_var_rnop_prt, & 97 (/ number_of_particles /), & 98 start = (/ prt_time_count /), count = (/ 1 /) ) 99 CALL handle_netcdf_error( 'lpm_data_output_particles', 2 ) 101 100 102 101 ! 103 !-- 104 105 106 107 102 !-- Output all particle attributes 103 nc_stat = NF90_PUT_VAR( id_set_prt, id_var_prt(1), particles%age, & 104 start = (/ 1, prt_time_count /), & 105 count = (/ maximum_number_of_particles /) ) 106 CALL handle_netcdf_error( 'lpm_data_output_particles', 3 ) 108 107 109 110 111 112 108 nc_stat = NF90_PUT_VAR( id_set_prt, id_var_prt(2), particles%dvrp_psize,& 109 start = (/ 1, prt_time_count /), & 110 count = (/ maximum_number_of_particles /) ) 111 CALL handle_netcdf_error( 'lpm_data_output_particles', 4 ) 113 112 114 115 116 117 113 nc_stat = NF90_PUT_VAR( id_set_prt, id_var_prt(3), particles%origin_x, & 114 start = (/ 1, prt_time_count /), & 115 count = (/ maximum_number_of_particles /) ) 116 CALL handle_netcdf_error( 'lpm_data_output_particles', 5 ) 118 117 119 120 121 122 118 nc_stat = NF90_PUT_VAR( id_set_prt, id_var_prt(4), particles%origin_y, & 119 start = (/ 1, prt_time_count /), & 120 count = (/ maximum_number_of_particles /) ) 121 CALL handle_netcdf_error( 'lpm_data_output_particles', 6 ) 123 122 124 125 126 127 123 nc_stat = NF90_PUT_VAR( id_set_prt, id_var_prt(5), particles%origin_z, & 124 start = (/ 1, prt_time_count /), & 125 count = (/ maximum_number_of_particles /) ) 126 CALL handle_netcdf_error( 'lpm_data_output_particles', 7 ) 128 127 129 130 131 132 128 nc_stat = NF90_PUT_VAR( id_set_prt, id_var_prt(6), particles%radius, & 129 start = (/ 1, prt_time_count /), & 130 count = (/ maximum_number_of_particles /) ) 131 CALL handle_netcdf_error( 'lpm_data_output_particles', 8 ) 133 132 134 135 136 137 133 nc_stat = NF90_PUT_VAR( id_set_prt, id_var_prt(7), particles%speed_x, & 134 start = (/ 1, prt_time_count /), & 135 count = (/ maximum_number_of_particles /) ) 136 CALL handle_netcdf_error( 'lpm_data_output_particles', 9 ) 138 137 139 140 141 142 138 nc_stat = NF90_PUT_VAR( id_set_prt, id_var_prt(8), particles%speed_y, & 139 start = (/ 1, prt_time_count /), & 140 count = (/ maximum_number_of_particles /) ) 141 CALL handle_netcdf_error( 'lpm_data_output_particles', 10 ) 143 142 144 145 146 147 143 nc_stat = NF90_PUT_VAR( id_set_prt, id_var_prt(9), particles%speed_z, & 144 start = (/ 1, prt_time_count /), & 145 count = (/ maximum_number_of_particles /) ) 146 CALL handle_netcdf_error( 'lpm_data_output_particles', 11 ) 148 147 149 150 151 152 153 148 nc_stat = NF90_PUT_VAR( id_set_prt,id_var_prt(10), & 149 particles%weight_factor, & 150 start = (/ 1, prt_time_count /), & 151 count = (/ maximum_number_of_particles /) ) 152 CALL handle_netcdf_error( 'lpm_data_output_particles', 12 ) 154 153 155 156 157 158 154 nc_stat = NF90_PUT_VAR( id_set_prt, id_var_prt(11), particles%x, & 155 start = (/ 1, prt_time_count /), & 156 count = (/ maximum_number_of_particles /) ) 157 CALL handle_netcdf_error( 'lpm_data_output_particles', 13 ) 159 158 160 nc_stat = NF90_PUT_VAR( id_set_prt, id_var_prt(12), particles%y, &161 162 163 159 nc_stat = NF90_PUT_VAR( id_set_prt, id_var_prt(12), particles%y, & 160 start = (/ 1, prt_time_count /), & 161 count = (/ maximum_number_of_particles /) ) 162 CALL handle_netcdf_error( 'lpm_data_output_particles', 14 ) 164 163 165 166 167 168 164 nc_stat = NF90_PUT_VAR( id_set_prt, id_var_prt(13), particles%z, & 165 start = (/ 1, prt_time_count /), & 166 count = (/ maximum_number_of_particles /) ) 167 CALL handle_netcdf_error( 'lpm_data_output_particles', 15 ) 169 168 170 171 172 173 169 nc_stat = NF90_PUT_VAR( id_set_prt, id_var_prt(14), particles%class, & 170 start = (/ 1, prt_time_count /), & 171 count = (/ maximum_number_of_particles /) ) 172 CALL handle_netcdf_error( 'lpm_data_output_particles', 16 ) 174 173 175 176 177 178 174 nc_stat = NF90_PUT_VAR( id_set_prt, id_var_prt(15), particles%group, & 175 start = (/ 1, prt_time_count /), & 176 count = (/ maximum_number_of_particles /) ) 177 CALL handle_netcdf_error( 'lpm_data_output_particles', 17 ) 179 178 180 181 182 183 184 179 nc_stat = NF90_PUT_VAR( id_set_prt, id_var_prt(16), & 180 particles%tailpoints, & 181 start = (/ 1, prt_time_count /), & 182 count = (/ maximum_number_of_particles /) ) 183 CALL handle_netcdf_error( 'lpm_data_output_particles', 18 ) 185 184 186 nc_stat = NF90_PUT_VAR( id_set_prt, id_var_prt(17), particles%tail_id, & 187 start = (/ 1, prt_time_count /), & 188 count = (/ maximum_number_of_particles /) ) 189 CALL handle_netcdf_error( 'lpm_data_output_particles', 19 ) 190 191 ENDIF 185 nc_stat = NF90_PUT_VAR( id_set_prt, id_var_prt(17), particles%tail_id, & 186 start = (/ 1, prt_time_count /), & 187 count = (/ maximum_number_of_particles /) ) 188 CALL handle_netcdf_error( 'lpm_data_output_particles', 19 ) 192 189 193 190 #endif
Note: See TracChangeset
for help on using the changeset viewer.