Changeset c8fc724 in flexpart.git
- Timestamp:
- Feb 2, 2017, 2:23:04 PM (7 years ago)
- Branches:
- master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug, univie
- Children:
- 93786a1
- Parents:
- 4c64400
- Location:
- src
- Files:
-
- 6 edited
Legend:
- Unmodified
- Added
- Removed
-
src/FLEXPART.f90
r4c64400 rc8fc724 414 414 415 415 ! NIK 16.02.2005 416 write(*,*) '**********************************************' 417 write(*,*) 'Total number of occurences of below-cloud scavenging', tot_blc_count 418 write(*,*) 'Total number of occurences of in-cloud scavenging', tot_inc_count 419 write(*,*) '**********************************************' 420 416 do i=1,nspec 417 write(*,*) '**********************************************' 418 write(*,*) 'Scavenging statistics for species ', species(i), ':' 419 write(*,*) 'Total number of occurences of below-cloud scavenging', & 420 & tot_blc_count(i) 421 write(*,*) 'Total number of occurences of in-cloud scavenging', & 422 & tot_inc_count(i) 423 write(*,*) '**********************************************' 424 end do 425 421 426 write(*,*) 'CONGRATULATIONS: YOU HAVE SUCCESSFULLY COMPLETED A FLE& 422 427 &XPART MODEL RUN!' -
src/FLEXPART_MPI.f90
r4c64400 rc8fc724 462 462 ! NIK 16.02.2005 463 463 if (lroot) then 464 call MPI_Reduce(MPI_IN_PLACE, tot_blc_count, 1, MPI_INTEGER8, MPI_SUM, id_root, &464 call MPI_Reduce(MPI_IN_PLACE, tot_blc_count, nspec, MPI_INTEGER8, MPI_SUM, id_root, & 465 465 & mp_comm_used, mp_ierr) 466 call MPI_Reduce(MPI_IN_PLACE, tot_inc_count, 1, MPI_INTEGER8, MPI_SUM, id_root, &466 call MPI_Reduce(MPI_IN_PLACE, tot_inc_count, nspec, MPI_INTEGER8, MPI_SUM, id_root, & 467 467 & mp_comm_used, mp_ierr) 468 468 else 469 469 if (mp_partgroup_pid.ge.0) then ! Skip for readwind process 470 call MPI_Reduce(tot_blc_count, 0, 1, MPI_INTEGER8, MPI_SUM, id_root, &470 call MPI_Reduce(tot_blc_count, 0, nspec, MPI_INTEGER8, MPI_SUM, id_root, & 471 471 & mp_comm_used, mp_ierr) 472 call MPI_Reduce(tot_inc_count, 0, 1, MPI_INTEGER8, MPI_SUM, id_root, &472 call MPI_Reduce(tot_inc_count, 0, nspec, MPI_INTEGER8, MPI_SUM, id_root, & 473 473 & mp_comm_used, mp_ierr) 474 474 end if … … 476 476 477 477 if (lroot) then 478 write(*,*) '**********************************************' 479 write(*,*) 'Total number of occurences of below-cloud scavenging', & 480 & tot_blc_count 481 write(*,*) 'Total number of occurences of in-cloud scavenging', & 482 & tot_inc_count 483 write(*,*) '**********************************************' 478 do i=1,nspec 479 write(*,*) '**********************************************' 480 write(*,*) 'Scavenging statistics for species ', species(i), ':' 481 write(*,*) 'Total number of occurences of below-cloud scavenging', & 482 & tot_blc_count(i) 483 write(*,*) 'Total number of occurences of in-cloud scavenging', & 484 & tot_inc_count(i) 485 write(*,*) '**********************************************' 486 end do 484 487 485 488 write(*,*) 'CONGRATULATIONS: YOU HAVE SUCCESSFULLY COMPLETED A FLE& -
src/com_mod.f90
r4c64400 rc8fc724 140 140 141 141 !NIK 16.02.2015 142 integer(selected_int_kind(16)) :: tot_blc_count=0, tot_inc_count=0 142 integer(selected_int_kind(16)), dimension(maxspec) :: tot_blc_count=0, & 143 &tot_inc_count=0 143 144 144 145 … … 576 577 real :: dxoutn,dyoutn,outlon0n,outlat0n,xoutshiftn,youtshiftn 577 578 !real outheight(maxzgrid),outheighthalf(maxzgrid) 578 logical :: DEP,DRYDEP,DRYDEPSPEC(maxspec),WETDEP,OHREA,ASSSPEC 579 logical :: DEP,DRYDEP,DRYDEPSPEC(maxspec),WETDEP,WETDEPSPEC(maxspec),& 580 & OHREA,ASSSPEC 579 581 580 582 ! numxgrid,numygrid number of grid points in x,y-direction … … 593 595 ! DRYDEPSPEC .true., if dry deposition is switched on for that species 594 596 ! WETDEP .true., if wet deposition is switched on 597 ! WETDEPSPEC .true., if wet deposition is switched on for that species 595 598 ! OHREA .true., if OH reaction is switched on 596 599 ! ASSSPEC .true., if there are two species asscoiated -
src/readreleases.f90
r4c64400 rc8fc724 276 276 do i=1,maxspec 277 277 DRYDEPSPEC(i)=.false. 278 WETDEPSPEC(i)=.false. 278 279 end do 279 280 … … 369 370 &(dquer(i).gt.0. .and. (crain_aero(i) .gt. 0. .or. csnow_aero(i).gt.0.))) then 370 371 WETDEP=.true. 372 WETDEPSPEC(i)=.true. 371 373 if (lroot) then 372 374 write (*,*) ' Below-cloud scavenging: ON' … … 380 382 if (dquer(i).gt.0..and.(ccn_aero(i).gt.0. .or. in_aero(i).gt.0.)) then 381 383 WETDEP=.true. 384 WETDEPSPEC(i)=.true. 382 385 if (lroot) then 383 386 write (*,*) ' In-cloud scavenging: ON' … … 399 402 endif 400 403 401 end do 404 end do ! end loop over species 402 405 403 406 if (WETDEP.or.DRYDEP) DEP=.true. -
src/readspecies.f90
r62e65c7 rc8fc724 197 197 weta_gas(pos_spec)=pweta_gas 198 198 wetb_gas(pos_spec)=pwetb_gas 199 crain_aero =pcrain_aero200 csnow_aero =pcsnow_aero199 crain_aero(pos_spec)=pcrain_aero 200 csnow_aero(pos_spec)=pcsnow_aero 201 201 ccn_aero(pos_spec)=pccn_aero 202 202 in_aero(pos_spec)=pin_aero -
src/wetdepo.f90
r05cf28d rc8fc724 90 90 real :: frac_act, liq_frac, dquer_m 91 91 92 integer :: blc_count, inc_count92 integer(selected_int_kind(16)), dimension(nspec) :: blc_count, inc_count 93 93 real :: Si_dummy, wetscav_dummy 94 94 logical :: readclouds_this_nest … … 107 107 !************************ 108 108 109 blc_count =0110 inc_count =0109 blc_count(:)=0 110 inc_count(:)=0 111 111 112 112 do jpart=1,numpart … … 256 256 do ks=1,nspec ! loop over species 257 257 wetdeposit(ks)=0. 258 wetscav=0. 258 wetscav=0. 259 260 ! Cycle loop if wet deposition for the species is off 261 if (WETDEPSPEC(ks).eqv..false.) cycle 259 262 260 263 if (ngrid.gt.0) then … … 274 277 if ((dquer(ks).le.0.).and.(weta_gas(ks).gt.0..or.wetb_gas(ks).gt.0.)) then 275 278 ! if (weta(ks).gt.0. .or. wetb(ks).gt.0.) then 276 blc_count =blc_count+1279 blc_count(ks)=blc_count(ks)+1 277 280 wetscav=weta_gas(ks)*prec(1)**wetb_gas(ks) 278 281 … … 280 283 !********************************************************************************* 281 284 else if ((dquer(ks).gt.0.).and.(crain_aero(ks).gt.0..or.csnow_aero(ks).gt.0.)) then 282 blc_count =blc_count+1285 blc_count(ks)=blc_count(ks)+1 283 286 284 287 !NIK 17.02.2015 … … 320 323 ! given in species file, or if gas and positive Henry's constant 321 324 if ((ccn_aero(ks).gt.0. .or. in_aero(ks).gt.0.).or.(henry(ks).gt.0.and.dquer(ks).le.0)) then 322 inc_count =inc_count+1325 inc_count(ks)=inc_count(ks)+1 323 326 ! if negative coefficients (turned off) set to zero for use in equation 324 327 if (ccn_aero(ks).lt.0.) ccn_aero(ks)=0. … … 432 435 433 436 ! count the total number of below-cloud and in-cloud occurences: 434 tot_blc_count =tot_blc_count+blc_count435 tot_inc_count =tot_inc_count+inc_count437 tot_blc_count(1:nspec)=tot_blc_count(1:nspec)+blc_count(1:nspec) 438 tot_inc_count(1:nspec)=tot_inc_count(1:nspec)+inc_count(1:nspec) 436 439 437 440 end subroutine wetdepo
Note: See TracChangeset
for help on using the changeset viewer.