Skip to content

Commit

Permalink
how did this happen - droped a whole calc
Browse files Browse the repository at this point in the history
  • Loading branch information
JhanSrbinovsky committed Dec 27, 2024
1 parent b702c69 commit 73bba86
Showing 1 changed file with 4 additions and 0 deletions.
4 changes: 4 additions & 0 deletions src/science/soilsnow/cbl_snowAccum.F90
Original file line number Diff line number Diff line change
Expand Up @@ -47,6 +47,8 @@ SUBROUTINE snow_accum ( dels, canopy, met, ssnow, soil )

ssnow%tgg(i,1) = ssnow%tgg(i,1) + canopy%precis(i) * CHLF &
/ ( REAL( ssnow%gammzz(i,1) ) + Ccswat *canopy%precis(i) )
ssnow%dtmlt(i,1) = ssnow%dtmlt(i,1) + canopy%precis(i) * CHLF &
/ ( REAL( ssnow%gammzz(i,1) ) + Ccswat *canopy%precis(i) )
! change density due to water being added
ssnow%ssdn(i,1) = MIN( max_ssdn, MAX( 120.0, ssnow%ssdn(i,1) &
* ssnow%osnowd(i) / MAX( 0.01, ssnow%snowd(i) ) + Cdensity_liq &
Expand Down Expand Up @@ -87,6 +89,8 @@ SUBROUTINE snow_accum ( dels, canopy, met, ssnow, soil )

ssnow%tggsn(i,1) = ssnow%tggsn(i,1) + canopy%precis(i) * CHLF &
* osm(i) / (sgamm(i) * ssnow%osnowd(i) )
ssnow%dtmlt(i,1) = ssnow%dtmlt(i,1) + canopy%precis(i) * CHLF &
* osm(i) / (sgamm(i) * ssnow%osnowd(i) )
ssnow%smass(i,1) = ssnow%smass(i,1) + canopy%precis(i) &
* osm(i)/ssnow%osnowd(i)

Expand Down

0 comments on commit 73bba86

Please sign in to comment.