Skip to content

Commit

Permalink
Update ptrack4 notes
Browse files Browse the repository at this point in the history
  • Loading branch information
josephzhang8 committed Nov 19, 2024
1 parent 65e2bc8 commit 5b4a850
Show file tree
Hide file tree
Showing 2 changed files with 11 additions and 13 deletions.
8 changes: 4 additions & 4 deletions src/Utility/Particle_Tracking/particle.bp
Original file line number Diff line number Diff line change
Expand Up @@ -2,18 +2,18 @@ Input for ptrack*
1 nscreen
0 mod_part (0: passive; 1: oil spill)
0 ibiofoul (0: no biofouling; 1: with biofouling)
1 ibf
0 istiff
1 ibf !(1: forward; -1: backward)
0 istiff !1: fixed distance from surface
0 ibnd_beh (behavior near bnd or wet/dry; 0: reflect off; 1: slide)
1 -122.6 37.38 ics slam0 sfea0
1 -122.6 37.38 ics slam0 sfea0 (from param.nml)
0.01 140. 120. 30 720 40 !h0,rnday,dtm,nspool,ihfskip,ndeltp
2 !# of particles
1 0.432000000E+05 0.334902000E+06 0.293906000E+06 0. !particle ID, starting time in sec, starting x,y,z (<=0, relative to the instantaneous surface)
2 0.864000000E+05 0.334902000E+06 0.293906000E+06 0.
!Oil spill parameters needed only if mod_part=1
1 3.0 0.2 !ihdf(0: constant diffusivity as hdc; 1: Smagorinsky),hdc,horcon (Smagorinsky coefficient)
1 0 !ibuoy (buoyancy off (0) or on (1)),iwind (wind effect off (0) or on (1))
20.0 !pbeach - % of stranding when particles reach shore
20.0 !set minimum percentage of stranding on shore (particles may be stranded if the random # exceeds this threshold)
!start of biofouling paramaters, needed only if ibiofoul=1
2.5e-3 0. 6.e-6 0.83 1.38 !bio_R0,bio_BT0,bio_BR,bio_den0 (\rho_0),bio_den_biolayer (\rho_D)

16 changes: 7 additions & 9 deletions src/Utility/Particle_Tracking/ptrack4.f90
Original file line number Diff line number Diff line change
Expand Up @@ -55,7 +55,13 @@
! (8) nparticle: # of particles;
! (9) idp(i),st_p(i),xpar(i),ypar(i),zpar0(i): particle id, start time (sec),
! starting x,y, and z relative to the instant f.s. (<=0).
! (10) additional parameters for oil spill
! (10) (additional parameters for oil spill) comment line
! (11) ihdf : turn on Smagorinsky algorithm - off(0), on(1)
! (12) ibuoy,iwind: turn on(1)/off(0) buoyancy of particle; wind effect on(1)/off(0)
! (13) pbeach: set minimum percentage of stranding on shore (particles may be stranded if the random # exceeds this threshold)
! (14) (optional for biofouling) comment line
! (15) bio_R0,bio_BT0,bio_BR,bio_den0,bio_den_biolayer: init plastic radius (m),init thickness of fouling layer (m),
! bio growth (m/day),init density of plastic particle (kg/m^3), density of biofoul layer (kg/m^3)
!
! Output: particle.pth, particle.pth.more (more info), fort.11 (fatal errors).
!
Expand Down Expand Up @@ -246,11 +252,6 @@ program ptrack
enddo !i

!... Additional parameters for oil spill
! ... Description of parameters
! ihdf : turn on Smagorinsky algorithm - off(0), on(1)
! ibuoy : turn buoyancy of particle - off(0), on(1)
! iwind : turn wind effect - off(0), on(1)
! pbeach : set percentage of stranding on shore
! ........................................................
if(mod_part==1) then
read(95,*) !comment line
Expand All @@ -261,9 +262,6 @@ program ptrack

if(ibiofoul/=0) then
read(95,*) !comment line
!init plastic radius (m),init thickness of fouling layer (m),
!bio growth (m/day),init density of plastic particle (kg/m^3),
!density of biofoul layer (kg/m^3)
read(95,*)bio_R0,bio_BT0,bio_BR,bio_den0,bio_den_biolayer
endif !ibiofoul/

Expand Down

0 comments on commit 5b4a850

Please sign in to comment.