-
Notifications
You must be signed in to change notification settings - Fork 34
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
RUC LSM update #116
RUC LSM update #116
Conversation
1. adding a heat source from fires. It is icontrolled by a flag and not turned on yet. Needs more testing. 2. Adding output of surface paramteres: leaf area index (LAI), wilting point and soil field capacity.
physics/lsm_ruc.F90
Outdated
!water and sea ice | ||
smcref2 (i) = one | ||
smcwlt2 (i) = zero | ||
! reassign smcref2 and smcwlt2 to RUC values at kdt=1 |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Would this be clearer to put the if (kdt==1)
in the else statement for if (.not. land(i))
? It looks like the same thing is done for the non-land points whether or not it is the first time step.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
@grantfirl Agree, the change is made.
physics/module_sf_ruclsm.F90
Outdated
@@ -2049,6 +2068,9 @@ SUBROUTINE SFCTMP (debug_print, delt,ktau,conflx,i,j, & !--- input varia | |||
eeta = eetas*(one-snowfrac) + eeta*snowfrac | |||
qfx = qfxs*(one-snowfrac) + qfx*snowfrac | |||
hfx = hfxs*(one-snowfrac) + hfx*snowfrac | |||
!IF ( add_fire_heat_flux ) then ! JLS |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Any reason to keep this commented-out section here and in line 2119?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
@grantfirl The commented out section is removed.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Looks OK to me with a couple minor comments.
1. Removed commented lines in fire feedback to RUC LSM. 2. The coordinates of a test point are passed to the RUC LSM subroutines to make debug printing more manageable. 3. Reverted change in the TRANSF subroutine to original version. The non-linear root distribution needs scientific evaluation in the retro runs.
…hermal conductivity is made in module_sf_ruclsm.F90.
…hard snow slab thermal conductivity.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
There are several lines of the code (physics/module_sf_ruclsm.F90) are commented out (lines 4083, 4086, 4540, 4543, 5217, 5220, 5793 and 5796, !fact = 5._kind_phys and !fact = 2._kind_phys), those lines can be removed.
@Qingfu-Liu The commented lines are removed. |
Also, added the fraction of grid cell that is burned out.
One more variable is added: fraction of grid cell burned by the fire. This fraction is used to take into account fire's effect on surface albedo. Also, added some debug prints. These changes will change the results only when fire feedback is turned on. The default is .false.
This PR includes: