Changes in src/FLEXPART_MPI.f90 [20963b1:0c00f1f] in flexpart.git
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
src/FLEXPART_MPI.f90
r20963b1 r0c00f1f 53 53 use conv_mod 54 54 use mpi_mod 55 use netcdf_output_mod, only: writeheader_netcdf 55 56 use random_mod, only: gasdev1 56 57 use class_gribfile 57 58 #ifdef USE_NCF59 use netcdf_output_mod, only: writeheader_netcdf60 #endif61 58 62 59 implicit none … … 67 64 integer :: metdata_format = GRIBFILE_CENTRE_UNKNOWN 68 65 integer :: detectformat 69 integer(selected_int_kind(16)), dimension(maxspec) :: tot_b=0, & 70 & tot_i=0 66 71 67 72 68 … … 207 203 208 204 if (metdata_format.eq.GRIBFILE_CENTRE_ECMWF) then 209 if (lroot)print *,'ECMWF metdata detected'205 print *,'ECMWF metdata detected' 210 206 elseif (metdata_format.eq.GRIBFILE_CENTRE_NCEP) then 211 if (lroot)print *,'NCEP metdata detected'207 print *,'NCEP metdata detected' 212 208 else 213 if (lroot)print *,'Unknown metdata format'209 print *,'Unknown metdata format' 214 210 stop 215 211 endif … … 382 378 383 379 if (mp_measure_time) call mpif_mtime('iotime',0) 384 385 380 if (lroot) then ! MPI: this part root process only 386 #ifdef USE_NCF 387 388 389 390 391 392 393 394 395 396 397 398 399 400 #endif 401 381 382 if (lnetcdfout.eq.1) then 383 call writeheader_netcdf(lnest=.false.) 384 else 385 call writeheader 386 end if 387 388 if (nested_output.eq.1) then 389 if (lnetcdfout.eq.1) then 390 call writeheader_netcdf(lnest=.true.) 391 else 392 call writeheader_nest 393 endif 394 endif 395 396 ! 402 397 if (verbosity.gt.0) then 403 398 print*,'call writeheader' … … 407 402 ! FLEXPART 9.2 ticket ?? write header in ASCII format 408 403 call writeheader_txt 409 404 !if (nested_output.eq.1) call writeheader_nest 410 405 if (nested_output.eq.1.and.surf_only.ne.1) call writeheader_nest 411 406 if (nested_output.eq.1.and.surf_only.eq.1) call writeheader_nest_surf … … 414 409 415 410 if (mp_measure_time) call mpif_mtime('iotime',0) 411 412 !open(unitdates,file=path(2)(1:length(2))//'dates') 416 413 417 414 if (verbosity.gt.0 .and. lroot) then … … 484 481 485 482 ! NIK 16.02.2005 486 if ( mp_partgroup_pid.ge.0) then ! Skip for readwind process487 call MPI_Reduce( tot_blc_count, tot_b, nspec, MPI_INTEGER8, MPI_SUM, id_root, &483 if (lroot) then 484 call MPI_Reduce(MPI_IN_PLACE, tot_blc_count, nspec, MPI_INTEGER8, MPI_SUM, id_root, & 488 485 & mp_comm_used, mp_ierr) 489 call MPI_Reduce( tot_inc_count, tot_i, nspec, MPI_INTEGER8, MPI_SUM, id_root, &486 call MPI_Reduce(MPI_IN_PLACE, tot_inc_count, nspec, MPI_INTEGER8, MPI_SUM, id_root, & 490 487 & mp_comm_used, mp_ierr) 488 else 489 if (mp_partgroup_pid.ge.0) then ! Skip for readwind process 490 call MPI_Reduce(tot_blc_count, 0, nspec, MPI_INTEGER8, MPI_SUM, id_root, & 491 & mp_comm_used, mp_ierr) 492 call MPI_Reduce(tot_inc_count, 0, nspec, MPI_INTEGER8, MPI_SUM, id_root, & 493 & mp_comm_used, mp_ierr) 494 end if 491 495 end if 492 493 496 494 497 if (lroot) then … … 497 500 write(*,*) 'Scavenging statistics for species ', species(i), ':' 498 501 write(*,*) 'Total number of occurences of below-cloud scavenging', & 499 & tot_b(i) 500 ! & tot_blc_count(i) 502 & tot_blc_count(i) 501 503 write(*,*) 'Total number of occurences of in-cloud scavenging', & 502 & tot_i(i) 503 ! & tot_inc_count(i) 504 & tot_inc_count(i) 504 505 write(*,*) '**********************************************' 505 506 end do
Note: See TracChangeset
for help on using the changeset viewer.