Lines Matching refs:blower
41 type(blowerData), allocatable :: blower(:) variable
63 if(.not. allocated(blower)) then
64 allocate(blower(nBlower))
68 blower(i)%mode = blowerMode(i)
69 blower(i)%surfID = surfID(i)
70 blower(i)%enable = (enable(i) == 1)
71 blower(i)%t_cycle = t_cycle(i)
72 blower(i)%t_fullOn = t_fullOn(i)
73 blower(i)%t_riseTime = t_riseTime(i)
74 blower(i)%t_fallTime = t_fallTime(i)
75 blower(i)%vmax = vmax(i)
76 blower(i)%vmin = vmin(i)
77 blower(i)%T = T(i)
78 blower(i)%nu = nu(i)
79 blower(i)%deltaBL = deltaBL(i)
80 blower(i)%deltaBLsclr = deltaBLsclr(i)
108 if(.not. blower(iBlower)%enable .or.
109 & .not. blower(iBlower)%updateBC ) cycle
114 vmax = blower(iBlower)%vmax
115 vmin = blower(iBlower)%vmin
117 if(blower(iBlower)%mode == 1) then !Trapezoid
127 t1 = t0 + blower(iBlower)%t_riseTime
128 t2 = t1 + blower(iBlower)%t_fullOn
129 t3 = t2 + blower(iBlower)%t_fallTime
130 t4 = blower(iBlower)%t_cycle
134 t = mod(BC_t - blower(iBlower)%t0, t4)
145 elseif(blower(iBlower)%mode == 2) then !sinusoid
154 t = (BC_t - blower(iBlower)%t0)/blower(iBlower)%t_cycle
156 elseif(blower(iBlower)%mode == 0) then !Constant
163 do i = 1,blower(iBlower)%nmap
164 imapped = blower(iBlower)%map(i)
165 vmag = vmagref*blower(iBlower)%vscale(i) !scale the reference
183 deallocate(blower(iBlower)%map)
184 deallocate(blower(iBlower)%vscale)
188 deallocate(blower) !note that blower gets allocated regardless
210 blower(i)%t0 = 0.0d00
226 blower(i)%t0 = -phi*blower(i)%t_cycle
240 call MPI_BCAST(blower(i)%t0, 1, MPI_DOUBLE_PRECISION,
267 if(blower(i)%updateBC) then
268 phi = mod((BC_t - blower(i)%t0)/blower(i)%t_cycle,1.0d0)
273 write(1003, "(i9, F18.14)") blower(i)%surfID, phi
338 if(blower(iBlower)%enable) then
357 allocate(blower(iBlower)%map(nshg))
358 call sfID2np(blower(iBlower)%surfID, !analog to sfID2np, except
359 & blower(iBlower)%nmap, !that map is allocated in
360 & blower(iBlower)%map) !the subroutine call.
364 nNodeMap = blower(iBlower)%nmap
367 allocate(blower(iBlower)%vscale(nNodeMap))
369 if(blower(iBlower)%deltaBL > 0.0) then !use the provided BL
370 deltaBLinv = 1/blower(iBlower)%deltaBL
378 imapped = blower(iBlower)%map(i)
379 blower(iBlower)%vscale(i) =
388 if(blower(iBlower)%T > 0.0) then
390 imapped = blower(iBlower)%map(i)
391 BC(imapped, 2) = blower(iBlower)%T
396 if(blower(iBlower)%nu >= 0.0) then
397 if(blower(iBlower)%deltaBLsclr > 0.0) then !use the provided BL
398 deltaBLinv = 1/blower(iBlower)%deltaBLsclr
404 imapped = blower(iBlower)%map(i)
405 BC(imapped, 7) = blower(iBlower)%nu *
411 blower(iBlower)%updateBC = .true. !update velocity later with a call to BC_iter.
417 if(blower(iBlower)%t_cycle >= 1.0 .or. !Unreasonable cycle
418 & blower(iBlower)%t_cycle <= 0.0 .or. !periods are provided,
419 & blower(iBlower)%mode == 0) then !so assume const.
420 blower(iBlower)%updateBC = .false.
422 blower(iBlower)%updateBC = .true.