Changeset 1320 for palm/trunk/SOURCE/data_output_2d.f90
- Timestamp:
- Mar 20, 2014 8:40:49 AM (10 years ago)
- File:
-
- 1 edited
-
palm/trunk/SOURCE/data_output_2d.f90 (modified) (36 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/data_output_2d.f90
r1319 r1320 20 20 ! Current revisions: 21 21 ! ----------------- 22 ! 22 ! ONLY-attribute added to USE-statements, 23 ! kind-parameters added to all INTEGER and REAL declaration statements, 24 ! kinds are defined in new module kinds, 25 ! old module precision_kind is removed, 26 ! revision history before 2012 removed, 27 ! comment fields (!:) to be used for variable explanations added to 28 ! all variable declaration statements 23 29 ! 24 30 ! Former revisions: … … 65 71 ! +z0h 66 72 ! 67 ! 790 2011-11-29 03:11:20Z raasch68 ! bugfix: calculation of 'pr' must depend on the particle weighting factor69 !70 ! 771 2011-10-27 10:56:21Z heinze71 ! +lpt72 !73 ! 759 2011-09-15 13:58:31Z raasch74 ! Splitting of parallel I/O75 !76 ! 729 2011-05-26 10:33:34Z heinze77 ! Exchange ghost layers for p regardless of used pressure solver (except SOR).78 !79 ! 691 2011-03-04 08:45:30Z maronga80 ! Replaced simulated_time by time_since_reference_point81 !82 ! 673 2011-01-18 16:19:48Z suehring83 ! When using Multigrid or SOR solver an additional CALL exchange_horiz is84 ! is needed for pressure output.85 !86 ! 667 2010-12-23 12:06:00Z suehring/gryschka87 ! nxl-1, nxr+1, nys-1, nyn+1 replaced by nxlg, nxrg, nysg, nyng in loops and88 ! allocation of arrays local_2d and total_2d.89 ! Calls of exchange_horiz are modiefied.90 !91 ! 622 2010-12-10 08:08:13Z raasch92 ! optional barriers included in order to speed up collective operations93 !94 ! 493 2010-03-01 08:30:24Z raasch95 ! netCDF4 support (parallel output)96 !97 ! 367 2009-08-25 08:35:52Z maronga98 ! simulated_time in netCDF output replaced by time_since_reference_point.99 ! Output of netCDF messages with aid of message handling routine.100 ! Bugfix: averaging along z is not allowed for 2d quantities (e.g. u* and z0)101 ! Output of messages replaced by message handling routine.102 ! Output of user defined 2D (XY) arrays at z=nzb+1 is now possible103 ! Bugfix: to_be_resorted => s_av for time-averaged scalars104 ! Calculation of shf* and qsws* added.105 !106 ! 215 2008-11-18 09:54:31Z raasch107 ! Bugfix: no output of particle concentration and radius unless particles108 ! have been started109 !110 ! 96 2007-06-04 08:07:41Z raasch111 ! Output of density and salinity112 !113 ! 75 2007-03-22 09:54:05Z raasch114 ! Output of precipitation amount/rate and roughness length,115 ! 2nd+3rd argument removed from exchange horiz116 !117 ! RCS Log replace by Id keyword, revision history cleaned up118 !119 ! Revision 1.5 2006/08/22 13:50:29 raasch120 ! xz and yz cross sections now up to nzt+1121 !122 ! Revision 1.2 2006/02/23 10:19:22 raasch123 ! Output of time-averaged data, output of averages along x, y, or z,124 ! output of user-defined quantities,125 ! section data are copied from local_pf to local_2d before they are output,126 ! output of particle concentration and mean radius,127 ! Former subroutine plot_2d renamed data_output_2d, pl2d.. renamed do2d..,128 ! anz renamed ngp, ebene renamed section, pl2d_.._anz renamed do2d_.._n129 !130 73 ! Revision 1.1 1997/08/11 06:24:09 raasch 131 74 ! Initial revision … … 140 83 !------------------------------------------------------------------------------! 141 84 142 USE arrays_3d 85 USE arrays_3d, & 86 ONLY: dzw, e, nr, p, pt, q, qc, ql, ql_c, ql_v, ql_vp, qr, qsws, & 87 rho, sa, shf, tend, ts, u, us, v, vpt, w, z0, z0h, zu, zw 88 143 89 USE averaging 144 USE cloud_parameters 145 USE control_parameters 146 USE cpulog 147 USE grid_variables 148 USE indices 90 91 USE cloud_parameters, & 92 ONLY: hyrho, l_d_cp, precipitation_amount, precipitation_rate, prr, & 93 pt_d_t 94 95 USE control_parameters, & 96 ONLY: cloud_physics, data_output_2d_on_each_pe, data_output_xy, & 97 data_output_xz, data_output_yz, do2d, & 98 do2d_xy_last_time, do2d_xy_n, do2d_xy_time_count, & 99 do2d_xz_last_time, do2d_xz_n, do2d_xz_time_count, & 100 do2d_yz_last_time, do2d_yz_n, do2d_yz_time_count, & 101 ibc_uv_b, icloud_scheme, io_blocks, io_group, iso2d_output, & 102 message_string, netcdf_data_format, netcdf_output, & 103 ntdim_2d_xy, ntdim_2d_xz, ntdim_2d_yz, psolver, section, & 104 simulated_time, simulated_time_chr, time_since_reference_point 105 106 USE cpulog, & 107 ONLY: cpu_log, log_point 108 109 USE grid_variables, & 110 ONLY: dx, dy 111 112 USE indices, & 113 ONLY: nbgp, nx, nxl, nxlg, nxr, nxrg, ny, nyn, nyng, nys, nysg, & 114 nz, nzb, nzt 115 116 USE kinds 117 149 118 USE netcdf_control 150 USE particle_attributes 119 120 USE particle_attributes, & 121 ONLY: particle_advection_start, particles, prt_count, & 122 prt_start_index 123 151 124 USE pegrid 152 125 153 126 IMPLICIT NONE 154 127 155 CHARACTER (LEN=2) :: do2d_mode, mode 156 CHARACTER (LEN=4) :: grid 157 CHARACTER (LEN=25) :: section_chr 158 CHARACTER (LEN=50) :: rtext 159 INTEGER :: av, ngp, file_id, i, if, is, iis, j, k, l, layer_xy, n, ns, & 160 psi, s, sender, & 161 ind(4) 162 LOGICAL :: found, resorted, two_d 163 REAL :: mean_r, s_r3, s_r4 164 REAL, DIMENSION(:), ALLOCATABLE :: level_z 165 REAL, DIMENSION(:,:), ALLOCATABLE :: local_2d, local_2d_l 166 REAL, DIMENSION(:,:,:), ALLOCATABLE :: local_pf, local_2d_sections, & 167 local_2d_sections_l 128 CHARACTER (LEN=2) :: do2d_mode !: 129 CHARACTER (LEN=2) :: mode !: 130 CHARACTER (LEN=4) :: grid !: 131 CHARACTER (LEN=25) :: section_chr !: 132 CHARACTER (LEN=50) :: rtext !: 133 134 INTEGER(iwp) :: av !: 135 INTEGER(iwp) :: ngp !: 136 INTEGER(iwp) :: file_id !: 137 INTEGER(iwp) :: i !: 138 INTEGER(iwp) :: if !: 139 INTEGER(iwp) :: is !: 140 INTEGER(iwp) :: iis !: 141 INTEGER(iwp) :: j !: 142 INTEGER(iwp) :: k !: 143 INTEGER(iwp) :: l !: 144 INTEGER(iwp) :: layer_xy !: 145 INTEGER(iwp) :: n !: 146 INTEGER(iwp) :: ns !: 147 INTEGER(iwp) :: psi !: 148 INTEGER(iwp) :: s !: 149 INTEGER(iwp) :: sender !: 150 INTEGER(iwp) :: ind(4) !: 151 152 LOGICAL :: found !: 153 LOGICAL :: resorted !: 154 LOGICAL :: two_d !: 155 156 REAL(wp) :: mean_r !: 157 REAL(wp) :: s_r3 !: 158 REAL(wp) :: s_r4 !: 159 160 REAL(wp), DIMENSION(:), ALLOCATABLE :: level_z !: 161 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: local_2d !: 162 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: local_2d_l !: 163 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: local_pf !: 164 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: local_2d_sections !: 165 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: local_2d_sections_l !: 168 166 #if defined( __parallel ) 169 REAL , DIMENSION(:,:), ALLOCATABLE :: total_2d170 #endif 171 REAL , DIMENSION(:,:,:), POINTER :: to_be_resorted167 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: total_2d !: 168 #endif 169 REAL(wp), DIMENSION(:,:,:), POINTER :: to_be_resorted !: 172 170 173 171 NAMELIST /LOCAL/ rtext … … 184 182 !-- the given end time by the length of the given output interval. 185 183 IF ( netcdf_data_format > 4 ) THEN 186 IF ( mode == 'xy' .AND. do2d_xy_time_count(av) + 1 > &184 IF ( mode == 'xy' .AND. do2d_xy_time_count(av) + 1 > & 187 185 ntdim_2d_xy(av) ) THEN 188 WRITE ( message_string, * ) 'Output of xy cross-sections is not ', &189 'given at t=', simulated_time, '&because the', &186 WRITE ( message_string, * ) 'Output of xy cross-sections is not ', & 187 'given at t=', simulated_time, '&because the', & 190 188 ' maximum number of output time levels is exceeded.' 191 189 CALL message( 'data_output_2d', 'PA0384', 0, 1, 0, 6, 0 ) 192 190 RETURN 193 191 ENDIF 194 IF ( mode == 'xz' .AND. do2d_xz_time_count(av) + 1 > &192 IF ( mode == 'xz' .AND. do2d_xz_time_count(av) + 1 > & 195 193 ntdim_2d_xz(av) ) THEN 196 WRITE ( message_string, * ) 'Output of xz cross-sections is not ', &197 'given at t=', simulated_time, '&because the', &194 WRITE ( message_string, * ) 'Output of xz cross-sections is not ', & 195 'given at t=', simulated_time, '&because the', & 198 196 ' maximum number of output time levels is exceeded.' 199 197 CALL message( 'data_output_2d', 'PA0385', 0, 1, 0, 6, 0 ) 200 198 RETURN 201 199 ENDIF 202 IF ( mode == 'yz' .AND. do2d_yz_time_count(av) + 1 > &200 IF ( mode == 'yz' .AND. do2d_yz_time_count(av) + 1 > & 203 201 ntdim_2d_yz(av) ) THEN 204 WRITE ( message_string, * ) 'Output of yz cross-sections is not ', &205 'given at t=', simulated_time, '&because the', &202 WRITE ( message_string, * ) 'Output of yz cross-sections is not ', & 203 'given at t=', simulated_time, '&because the', & 206 204 ' maximum number of output time levels is exceeded.' 207 205 CALL message( 'data_output_2d', 'PA0386', 0, 1, 0, 6, 0 ) … … 363 361 DO i = nxlg, nxrg 364 362 DO j = nysg, nyng 365 local_pf(i,j,nzb+1) = SUM( ql(nzb:nzt,j,i) * &363 local_pf(i,j,nzb+1) = SUM( ql(nzb:nzt,j,i) * & 366 364 dzw(1:nzt+1) ) 367 365 ENDDO … … 427 425 s_r4 = 0.0 428 426 DO n = psi, psi+prt_count(k,j,i)-1 429 s_r3 = s_r3 + particles(n)%radius**3 * &427 s_r3 = s_r3 + particles(n)%radius**3 * & 430 428 particles(n)%weight_factor 431 s_r4 = s_r4 + particles(n)%radius**4 * &429 s_r4 = s_r4 + particles(n)%radius**4 * & 432 430 particles(n)%weight_factor 433 431 ENDDO … … 499 497 DO i = nxlg, nxrg 500 498 DO j = nysg, nyng 501 local_pf(i,j,nzb+1) = prr_av(nzb+1,j,i) * hyrho(nzb+1) 499 local_pf(i,j,nzb+1) = prr_av(nzb+1,j,i) * & 500 hyrho(nzb+1) 502 501 ENDDO 503 502 ENDDO … … 539 538 DO j = nysg, nyng 540 539 DO k = nzb, nzt+1 541 local_pf(i,j,k) = pt(k,j,i) + l_d_cp * &542 pt_d_t(k) * &540 local_pf(i,j,k) = pt(k,j,i) + l_d_cp * & 541 pt_d_t(k) * & 543 542 ql(k,j,i) 544 543 ENDDO … … 600 599 psi = prt_start_index(k,j,i) 601 600 DO n = psi, psi+prt_count(k,j,i)-1 602 tend(k,j,i) = tend(k,j,i) + &603 particles(n)%weight_factor / &601 tend(k,j,i) = tend(k,j,i) + & 602 particles(n)%weight_factor / & 604 603 prt_count(k,j,i) 605 604 ENDDO … … 824 823 ! 825 824 !-- User defined quantity 826 CALL user_data_output_2d( av, do2d(av,if), found, grid, &825 CALL user_data_output_2d( av, do2d(av,if), found, grid, & 827 826 local_pf, two_d ) 828 827 resorted = .TRUE. … … 837 836 838 837 IF ( .NOT. found ) THEN 839 message_string = 'no output provided for: ' // &838 message_string = 'no output provided for: ' // & 840 839 TRIM( do2d(av,if) ) 841 840 CALL message( 'data_output_2d', 'PA0181', 0, 0, 0, 6, 0 ) … … 881 880 do2d_xy_last_time(av) = simulated_time 882 881 IF ( myid == 0 ) THEN 883 IF ( ( .NOT. data_output_2d_on_each_pe .AND. &884 netcdf_output ) .OR. netcdf_data_format > 4 ) &882 IF ( ( .NOT. data_output_2d_on_each_pe .AND. & 883 netcdf_output ) .OR. netcdf_data_format > 4 ) & 885 884 THEN 886 885 #if defined( __netcdf ) … … 947 946 #if defined( __netcdf ) 948 947 IF ( netcdf_output .AND. myid == 0 ) THEN 949 WRITE ( 21 ) time_since_reference_point, &948 WRITE ( 21 ) time_since_reference_point, & 950 949 do2d_xy_time_count(av), av 951 950 ENDIF … … 981 980 !-- Index limits are received in arbitrary order from 982 981 !-- the PEs. 983 CALL MPI_RECV( ind(1), 4, MPI_INTEGER, &984 MPI_ANY_SOURCE, 0, comm2d, &982 CALL MPI_RECV( ind(1), 4, MPI_INTEGER, & 983 MPI_ANY_SOURCE, 0, comm2d, & 985 984 status, ierr ) 986 985 sender = status(MPI_SOURCE) 987 986 DEALLOCATE( local_2d ) 988 987 ALLOCATE( local_2d(ind(1):ind(2),ind(3):ind(4)) ) 989 CALL MPI_RECV( local_2d(ind(1),ind(3)), ngp, &990 MPI_REAL, sender, 1, comm2d, &988 CALL MPI_RECV( local_2d(ind(1),ind(3)), ngp, & 989 MPI_REAL, sender, 1, comm2d, & 991 990 status, ierr ) 992 991 total_2d(ind(1):ind(2),ind(3):ind(4)) = local_2d … … 1026 1025 ind(1) = nxlg; ind(2) = nxrg 1027 1026 ind(3) = nysg; ind(4) = nyng 1028 CALL MPI_SEND( ind(1), 4, MPI_INTEGER, 0, 0, &1027 CALL MPI_SEND( ind(1), 4, MPI_INTEGER, 0, 0, & 1029 1028 comm2d, ierr ) 1030 1029 ! 1031 1030 !-- Send data to PE0 1032 CALL MPI_SEND( local_2d(nxlg,nysg), ngp, &1031 CALL MPI_SEND( local_2d(nxlg,nysg), ngp, & 1033 1032 MPI_REAL, 0, 1, comm2d, ierr ) 1034 1033 ENDIF … … 1076 1075 ENDIF 1077 1076 IF ( av == 0 ) THEN 1078 rtext = TRIM( do2d(av,if) ) // ' t = ' // &1079 TRIM( simulated_time_chr ) // ' ' // &1077 rtext = TRIM( do2d(av,if) ) // ' t = ' // & 1078 TRIM( simulated_time_chr ) // ' ' // & 1080 1079 TRIM( section_chr ) 1081 1080 ELSE 1082 rtext = TRIM( do2d(av,if) ) // ' averaged t = ' // &1083 TRIM( simulated_time_chr ) // ' ' // &1081 rtext = TRIM( do2d(av,if) ) // ' averaged t = ' // & 1082 TRIM( simulated_time_chr ) // ' ' // & 1084 1083 TRIM( section_chr ) 1085 1084 ENDIF … … 1103 1102 do2d_xz_last_time(av) = simulated_time 1104 1103 IF ( myid == 0 ) THEN 1105 IF ( ( .NOT. data_output_2d_on_each_pe .AND. &1106 netcdf_output ) .OR. netcdf_data_format > 4 ) &1104 IF ( ( .NOT. data_output_2d_on_each_pe .AND. & 1105 netcdf_output ) .OR. netcdf_data_format > 4 ) & 1107 1106 THEN 1108 1107 #if defined( __netcdf ) … … 1130 1129 DO j = nys, nyn 1131 1130 DO i = nxlg, nxrg 1132 local_2d_l(i,k) = local_2d_l(i,k) + &1131 local_2d_l(i,k) = local_2d_l(i,k) + & 1133 1132 local_pf(i,j,k) 1134 1133 ENDDO … … 1139 1138 !-- Now do the averaging over all PEs along y 1140 1139 IF ( collective_wait ) CALL MPI_BARRIER( comm2d, ierr ) 1141 CALL MPI_ALLREDUCE( local_2d_l(nxlg,nzb), &1142 local_2d(nxlg,nzb), ngp, MPI_REAL, &1140 CALL MPI_ALLREDUCE( local_2d_l(nxlg,nzb), & 1141 local_2d(nxlg,nzb), ngp, MPI_REAL, & 1143 1142 MPI_SUM, comm1dy, ierr ) 1144 1143 #else … … 1167 1166 !-- sections reside. Cross sections averaged along y are 1168 1167 !-- output on the respective first PE along y (myidy=0). 1169 IF ( ( section(is,s) >= nys .AND. &1170 section(is,s) <= nyn ) .OR. &1168 IF ( ( section(is,s) >= nys .AND. & 1169 section(is,s) <= nyn ) .OR. & 1171 1170 ( section(is,s) == -1 .AND. myidy == 0 ) ) THEN 1172 1171 #if defined( __netcdf ) … … 1192 1191 #if defined( __netcdf ) 1193 1192 IF ( netcdf_output .AND. myid == 0 ) THEN 1194 WRITE ( 22 ) time_since_reference_point, &1193 WRITE ( 22 ) time_since_reference_point, & 1195 1194 do2d_xz_time_count(av), av 1196 1195 ENDIF … … 1198 1197 DO i = 0, io_blocks-1 1199 1198 IF ( i == io_group ) THEN 1200 IF ( ( section(is,s) >= nys .AND. &1201 section(is,s) <= nyn ) .OR. &1202 ( section(is,s) == -1 .AND. &1203 nys-1 == -1 ) ) &1199 IF ( ( section(is,s) >= nys .AND. & 1200 section(is,s) <= nyn ) .OR. & 1201 ( section(is,s) == -1 .AND. & 1202 nys-1 == -1 ) ) & 1204 1203 THEN 1205 1204 WRITE (22) nxlg, nxrg, nzb, nzt+1 … … 1239 1238 !-- Index limits are received in arbitrary order from 1240 1239 !-- the PEs. 1241 CALL MPI_RECV( ind(1), 4, MPI_INTEGER, &1242 MPI_ANY_SOURCE, 0, comm2d, &1240 CALL MPI_RECV( ind(1), 4, MPI_INTEGER, & 1241 MPI_ANY_SOURCE, 0, comm2d, & 1243 1242 status, ierr ) 1244 1243 ! … … 1247 1246 sender = status(MPI_SOURCE) 1248 1247 DEALLOCATE( local_2d ) 1249 ALLOCATE( local_2d(ind(1):ind(2), &1248 ALLOCATE( local_2d(ind(1):ind(2), & 1250 1249 ind(3):ind(4)) ) 1251 1250 CALL MPI_RECV( local_2d(ind(1),ind(3)), ngp, & 1252 1251 MPI_REAL, sender, 1, comm2d, & 1253 1252 status, ierr ) 1254 total_2d(ind(1):ind(2),ind(3):ind(4)) = &1253 total_2d(ind(1):ind(2),ind(3):ind(4)) = & 1255 1254 local_2d 1256 1255 ENDIF … … 1291 1290 ind(3) = -9999; ind(4) = -9999 1292 1291 ENDIF 1293 CALL MPI_SEND( ind(1), 4, MPI_INTEGER, 0, 0, &1292 CALL MPI_SEND( ind(1), 4, MPI_INTEGER, 0, 0, & 1294 1293 comm2d, ierr ) 1295 1294 ! 1296 1295 !-- If applicable, send data to PE0. 1297 1296 IF ( ind(1) /= -9999 ) THEN 1298 CALL MPI_SEND( local_2d(nxlg,nzb), ngp, &1297 CALL MPI_SEND( local_2d(nxlg,nzb), ngp, & 1299 1298 MPI_REAL, 0, 1, comm2d, ierr ) 1300 1299 ENDIF … … 1335 1334 ENDIF 1336 1335 IF ( av == 0 ) THEN 1337 rtext = TRIM( do2d(av,if) ) // ' t = ' // &1338 TRIM( simulated_time_chr ) // ' ' // &1336 rtext = TRIM( do2d(av,if) ) // ' t = ' // & 1337 TRIM( simulated_time_chr ) // ' ' // & 1339 1338 TRIM( section_chr ) 1340 1339 ELSE 1341 rtext = TRIM( do2d(av,if) ) // ' averaged t = ' // &1342 TRIM( simulated_time_chr ) // ' ' // &1340 rtext = TRIM( do2d(av,if) ) // ' averaged t = ' // & 1341 TRIM( simulated_time_chr ) // ' ' // & 1343 1342 TRIM( section_chr ) 1344 1343 ENDIF … … 1355 1354 do2d_yz_last_time(av) = simulated_time 1356 1355 IF ( myid == 0 ) THEN 1357 IF ( ( .NOT. data_output_2d_on_each_pe .AND. &1358 netcdf_output ) .OR. netcdf_data_format > 4 ) &1356 IF ( ( .NOT. data_output_2d_on_each_pe .AND. & 1357 netcdf_output ) .OR. netcdf_data_format > 4 ) & 1359 1358 THEN 1360 1359 #if defined( __netcdf ) … … 1382 1381 DO j = nysg, nyng 1383 1382 DO i = nxl, nxr 1384 local_2d_l(j,k) = local_2d_l(j,k) + &1383 local_2d_l(j,k) = local_2d_l(j,k) + & 1385 1384 local_pf(i,j,k) 1386 1385 ENDDO … … 1391 1390 !-- Now do the averaging over all PEs along x 1392 1391 IF ( collective_wait ) CALL MPI_BARRIER( comm2d, ierr ) 1393 CALL MPI_ALLREDUCE( local_2d_l(nysg,nzb), &1394 local_2d(nysg,nzb), ngp, MPI_REAL, &1392 CALL MPI_ALLREDUCE( local_2d_l(nysg,nzb), & 1393 local_2d(nysg,nzb), ngp, MPI_REAL, & 1395 1394 MPI_SUM, comm1dx, ierr ) 1396 1395 #else … … 1419 1418 !-- sections reside. Cross sections averaged along x are 1420 1419 !-- output on the respective first PE along x (myidx=0). 1421 IF ( ( section(is,s) >= nxl .AND. &1422 section(is,s) <= nxr ) .OR. &1420 IF ( ( section(is,s) >= nxl .AND. & 1421 section(is,s) <= nxr ) .OR. & 1423 1422 ( section(is,s) == -1 .AND. myidx == 0 ) ) THEN 1424 1423 #if defined( __netcdf ) … … 1444 1443 #if defined( __netcdf ) 1445 1444 IF ( netcdf_output .AND. myid == 0 ) THEN 1446 WRITE ( 23 ) time_since_reference_point, &1445 WRITE ( 23 ) time_since_reference_point, & 1447 1446 do2d_yz_time_count(av), av 1448 1447 ENDIF … … 1450 1449 DO i = 0, io_blocks-1 1451 1450 IF ( i == io_group ) THEN 1452 IF ( ( section(is,s) >= nxl .AND. &1453 section(is,s) <= nxr ) .OR. &1454 ( section(is,s) == -1 .AND. &1455 nxl-1 == -1 ) ) &1451 IF ( ( section(is,s) >= nxl .AND. & 1452 section(is,s) <= nxr ) .OR. & 1453 ( section(is,s) == -1 .AND. & 1454 nxl-1 == -1 ) ) & 1456 1455 THEN 1457 1456 WRITE (23) nysg, nyng, nzb, nzt+1 … … 1491 1490 !-- Index limits are received in arbitrary order from 1492 1491 !-- the PEs. 1493 CALL MPI_RECV( ind(1), 4, MPI_INTEGER, &1494 MPI_ANY_SOURCE, 0, comm2d, &1492 CALL MPI_RECV( ind(1), 4, MPI_INTEGER, & 1493 MPI_ANY_SOURCE, 0, comm2d, & 1495 1494 status, ierr ) 1496 1495 ! … … 1499 1498 sender = status(MPI_SOURCE) 1500 1499 DEALLOCATE( local_2d ) 1501 ALLOCATE( local_2d(ind(1):ind(2), &1500 ALLOCATE( local_2d(ind(1):ind(2), & 1502 1501 ind(3):ind(4)) ) 1503 1502 CALL MPI_RECV( local_2d(ind(1),ind(3)), ngp, & 1504 1503 MPI_REAL, sender, 1, comm2d, & 1505 1504 status, ierr ) 1506 total_2d(ind(1):ind(2),ind(3):ind(4)) = &1505 total_2d(ind(1):ind(2),ind(3):ind(4)) = & 1507 1506 local_2d 1508 1507 ENDIF … … 1543 1542 ind(3) = -9999; ind(4) = -9999 1544 1543 ENDIF 1545 CALL MPI_SEND( ind(1), 4, MPI_INTEGER, 0, 0, &1544 CALL MPI_SEND( ind(1), 4, MPI_INTEGER, 0, 0, & 1546 1545 comm2d, ierr ) 1547 1546 ! 1548 1547 !-- If applicable, send data to PE0. 1549 1548 IF ( ind(1) /= -9999 ) THEN 1550 CALL MPI_SEND( local_2d(nysg,nzb), ngp, &1549 CALL MPI_SEND( local_2d(nysg,nzb), ngp, & 1551 1550 MPI_REAL, 0, 1, comm2d, ierr ) 1552 1551 ENDIF … … 1587 1586 ENDIF 1588 1587 IF ( av == 0 ) THEN 1589 rtext = TRIM( do2d(av,if) ) // ' t = ' // &1590 TRIM( simulated_time_chr ) // ' ' // &1588 rtext = TRIM( do2d(av,if) ) // ' t = ' // & 1589 TRIM( simulated_time_chr ) // ' ' // & 1591 1590 TRIM( section_chr ) 1592 1591 ELSE 1593 rtext = TRIM( do2d(av,if) ) // ' averaged t = ' // &1594 TRIM( simulated_time_chr ) // ' ' // &1592 rtext = TRIM( do2d(av,if) ) // ' averaged t = ' // & 1593 TRIM( simulated_time_chr ) // ' ' // & 1595 1594 TRIM( section_chr ) 1596 1595 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.
