From 3f78267b38e0f0a36da49a58b8f00d00789ee67a Mon Sep 17 00:00:00 2001 From: Jan Mandel Date: Mon, 17 Jan 2011 18:40:38 -0700 Subject: [PATCH] fixed error in heat flux instertion code that was there from the beginning --- wrfv2_fire/phys/module_fr_sfire_atm.F | 11 +++++++---- 1 file changed, 7 insertions(+), 4 deletions(-) diff --git a/wrfv2_fire/phys/module_fr_sfire_atm.F b/wrfv2_fire/phys/module_fr_sfire_atm.F index 47a6c7a0..8b0e6323 100644 --- a/wrfv2_fire/phys/module_fr_sfire_atm.F +++ b/wrfv2_fire/phys/module_fr_sfire_atm.F @@ -65,6 +65,7 @@ SUBROUTINE fire_tendency( & REAL :: xlv_i REAL :: z_w REAL :: fact_g, fact_c + REAL :: alfg_i, alfc_i REAL, DIMENSION( its:ite,kts:kte,jts:jte ) :: hfx,qfx @@ -85,6 +86,8 @@ SUBROUTINE fire_tendency( & cp_i = 1./cp ! inverse of specific heat xlv_i = 1./xlv ! inverse of latent heat + alfg_i = 1./alfg + alfc_i = 1./alfc !!write(msg,'(8e11.3)')cp,cp_i,xlv,xlv_i,alfg,alfc,z1can !!call message(msg) @@ -113,11 +116,11 @@ SUBROUTINE fire_tendency( & ! --- heat flux - fact_g = cp_i * EXP( - alfg * z_w ) + fact_g = cp_i * EXP( - alfg_i * z_w ) IF ( z_w < z1can ) THEN fact_c = cp_i ELSE - fact_c = cp_i * EXP( - alfc * (z_w - z1can) ) + fact_c = cp_i * EXP( - alfc_i * (z_w - z1can) ) END IF hfx(i,k,j) = fact_g * grnhfx(i,j) + fact_c * canhfx(i,j) @@ -127,11 +130,11 @@ SUBROUTINE fire_tendency( & ! --- vapor flux - fact_g = xlv_i * EXP( - alfg * z_w ) + fact_g = xlv_i * EXP( - alfg_i * z_w ) IF (z_w < z1can) THEN fact_c = xlv_i ELSE - fact_c = xlv_i * EXP( - alfc * (z_w - z1can) ) + fact_c = xlv_i * EXP( - alfc_i * (z_w - z1can) ) END IF qfx(i,k,j) = fact_g * grnqfx(i,j) + fact_c * canqfx(i,j)