Skip to content
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

Merged
merged 11 commits into from
Nov 17, 2023
Merged

RUC LSM update #116

merged 11 commits into from
Nov 17, 2023

Conversation

tanyasmirnova
Copy link
Collaborator

This PR includes:

  1. Adding leaf area index, wilting point and soil field capacity to the history files.
  2. Adding a placeholder for the fire heat flux provided by smoke-dust code.
  3. A change to the transpiration function to take into account the non-linear distribution of roots.

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.
!water and sea ice
smcref2 (i) = one
smcwlt2 (i) = zero
! reassign smcref2 and smcwlt2 to RUC values at kdt=1
Copy link
Collaborator

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.

Copy link
Collaborator Author

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.

@@ -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
Copy link
Collaborator

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?

Copy link
Collaborator Author

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.

Copy link
Collaborator

@grantfirl grantfirl left a 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.
Copy link
Collaborator

@Qingfu-Liu Qingfu-Liu left a 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.

@tanyasmirnova
Copy link
Collaborator Author

@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.
@dustinswales dustinswales merged commit 5c8a8c1 into ufs-community:ufs/dev Nov 17, 2023
3 checks passed
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

4 participants