Skip to content

Commit

Permalink
Browse files Browse the repository at this point in the history
Defined default perturbation parameters in registry.fire consistent w…
…ith the original version

Brought back original namelist.input for the hill case
Cosmetic modification in module_initialize_fire.F
  • Loading branch information
adamk0 committed Feb 27, 2010
1 parent ef550db commit 993cd88
Show file tree
Hide file tree
Showing 3 changed files with 28 additions and 22 deletions.
10 changes: 5 additions & 5 deletions wrfv2_fire/Registry/registry.fire
Expand Up @@ -147,11 +147,11 @@ rconfig real fire_mountain_start_x namelist,fire max_domains
rconfig real fire_mountain_start_y namelist,fire max_domains 100. - "fire_mountain_start_y" "y coord of start of the mountain" "m"
rconfig real fire_mountain_end_x namelist,fire max_domains 100. - "fire_mountain_end_x" "x coord of end of the mountain" "m"
rconfig real fire_mountain_end_y namelist,fire max_domains 100. - "fire_mountain_end_y" "y coord of end of the mountain" "m"
rconfig real delt_perturbation namelist,fire max_domains 0. - "delt_perturbation" "temperature perturbation for cold (-) /warm (+) bubble" "K"
rconfig real xrad_perturbation namelist,fire max_domains 0. - "xrad_perturbation" "horizontal radius of the perturbation in E-W direction" "m"
rconfig real yrad_perturbation namelist,fire max_domains 0. - "yrad_perturbation" "horizontal radius of the perturbation in N-S direction" "m"
rconfig real zrad_perturbation namelist,fire max_domains 0. - "zrad_perturbation" "vertical radius of the perturbation (bubble) direction" "m"
rconfig real hght_perturbation namelist,fire max_domains 0. - "hght_perturbation" "height at which the perturbation (bubble) will be suspended" "m"
rconfig real delt_perturbation namelist,fire max_domains 3.0 - "delt_perturbation" "temperature perturbation for cold (-) /warm (+) bubble" "K"
rconfig real xrad_perturbation namelist,fire max_domains 10000.0 - "xrad_perturbation" "horizontal radius of the perturbation in E-W direction" "m"
rconfig real yrad_perturbation namelist,fire max_domains 10000.0 - "yrad_perturbation" "horizontal radius of the perturbation in N-S direction" "m"
rconfig real zrad_perturbation namelist,fire max_domains 1500.0 - "zrad_perturbation" "vertical radius of the perturbation (bubble) direction" "m"
rconfig real hght_perturbation namelist,fire max_domains 1500.0 - "hght_perturbation" "height at which the perturbation (bubble) will be suspended" "m"



Expand Down
35 changes: 23 additions & 12 deletions wrfv2_fire/dyn_em/module_initialize_fire.F
Expand Up @@ -92,6 +92,7 @@ SUBROUTINE init_domain_rk ( grid &
INTEGER :: icm,jcm, ii, im1, jj, jm1, loop, error, fid, nxc, nyc
REAL :: u_mean,v_mean, f0, p_surf, p_level, qvf, z_at_v, z_at_u
REAL :: z_scale, xrad, yrad, zrad, rad, delt, cof1, cof2
REAL :: x_rad, y_rad, z_rad, hght_pert
! REAL, EXTERNAL :: interp_0
REAL :: hm
REAL :: pi
Expand Down Expand Up @@ -173,7 +174,17 @@ SUBROUTINE init_domain_rk ( grid &

CALL boundary_condition_check( config_flags, bdyzone, error, grid%id )

delt = config_flags%delt_perturbation
delt = config_flags%delt_perturbation
x_rad = config_flags%xrad_perturbation
y_rad = config_flags%yrad_perturbation
z_rad = config_flags%zrad_perturbation
hght_pert = config_flags%hght_perturbation

! delt = 3.0
! x_rad = 10000.0
! y_rad = 10000.0
! z_rad = 1500.0
! hght_pert = 1500.0

moisture_init = .true.

Expand Down Expand Up @@ -532,23 +543,23 @@ SUBROUTINE init_domain_rk ( grid &
ENDDO

! checking if the perturbation (bubble) is to be applied
IF ((delt/=0.) .and. (config_flags%xrad_perturbation > 0.) &
.and. (config_flags%yrad_perturbation > 0.) &
.and. (config_flags%zrad_perturbation > 0.)) THEN
IF ((delt/=0.) .and. (x_rad > 0.) &
.and. (y_rad > 0.) &
.and. (z_rad > 0.)) THEN
! thermal perturbation to kick off convection

write(6,*) ' nxc, nyc for perturbation ',nxc,nyc
write(6,*) ' delt for perturbation ',delt
write(6,*) ' x radius of the perturbation ' ,config_flags%xrad_perturbation
write(6,*) ' y radius of the perturbation ' ,config_flags%yrad_perturbation
write(6,*) ' z radius of the perturbation ' ,config_flags%zrad_perturbation
write(6,*) ' height of the perturbation ' ,config_flags%hght_perturbation
write(6,('A23,f18.16')) ' delt for perturbation ',delt
write(6,('A30,f18.12')) ' x radius of the perturbation ' ,x_rad
write(6,('A30,f18.12')) ' y radius of the perturbation ' ,y_rad
write(6,('A30,f18.12')) ' z radius of the perturbation ' ,z_rad
write(6,('A30,f18.12')) ' height of the perturbation ' ,hght_pert

DO J = jts, min(jde-1,jte)
yrad = config_flags%dy*float(j-nyc)/config_flags%yrad_perturbation
yrad = config_flags%dy*float(j-nyc)/y_rad
! yrad = 0.
DO I = its, min(ide-1,ite)
xrad = config_flags%dx*float(i-nxc)/config_flags%xrad_perturbation
xrad = config_flags%dx*float(i-nxc)/x_rad
! xrad = 0.
DO K = 1, kte-1

Expand All @@ -558,7 +569,7 @@ SUBROUTINE init_domain_rk ( grid &

zrad = 0.5*(grid%ph_1(i,k,j)+grid%ph_1(i,k+1,j) &
+grid%phb(i,k,j)+grid%phb(i,k+1,j))/g
zrad = (zrad-config_flags%hght_perturbation)/config_flags%zrad_perturbation
zrad = (zrad-hght_pert)/z_rad
RAD=SQRT(xrad*xrad+yrad*yrad+zrad*zrad)
IF(RAD <= 1.) THEN
grid%t_1(i,k,j)=grid%t_1(i,k,j)+delt*COS(.5*PI*RAD)**2
Expand Down
5 changes: 0 additions & 5 deletions wrfv2_fire/test/em_fire/hill/namelist.input
Expand Up @@ -162,11 +162,6 @@
fire_mountain_end_x=1400., ! (m) coord of end of the mountain from lower left corder (just like ignition)
fire_mountain_end_y=1600., ! (m) coord of end of the mountain from lower left corder (just like ignition)
fire_topo_from_atm=0, ! 0 = fire mesh topo set from fine-res data, 1 = populate by interpolating from atmosphere
delt_perturbation = 3.0, ! Temperature perturbation for creating cold (negative) / warm (positive) bubble [K], 0 turns it off
xrad_perturbation = 10000.0, ! Horizontal radius of the bubble in E-W direction [m]
yrad_perturbation = 10000.0, ! Horizontal radius of the bubble in N-S direction [m]
zrad_perturbation = 1500.0, ! Vertical radius of the bubble [m]
hght_perturbation = 1500.0, ! Perturbation height - height at which the warm/cold bubble will be suspended [m]

!
! method switches for developers only, do not change!
Expand Down

0 comments on commit 993cd88

Please sign in to comment.