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.

Sign up for GitHub

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

Add probabilistic grounding scheme for landfast ice #565

Merged
merged 27 commits into from
Feb 22, 2021
Merged

Add probabilistic grounding scheme for landfast ice #565

Changes from 1 commit
Commits
Show all changes
27 commits
Select commit Hold shift + click to select a range
8f49435
In the process of adding seabed2 parameterization
JFLemieux73 Jan 26, 2021
73cdca4
Added calls foe icepack_query_parameters
JFLemieux73 Jan 26, 2021
d7a9eec
Modified modules used in subroutine to match latest code (it now comp…
JFLemieux73 Jan 26, 2021
3a4c2b4
Changed expression basal stress to seabed stress
JFLemieux73 Jan 26, 2021
e0aaf74
Replaced basal by better expression seabed at many places
JFLemieux73 Jan 27, 2021
c9222d0
seabed everywhere...basalstress logical replaced by seabedstress
JFLemieux73 Jan 27, 2021
1f6f09a
Added kseabed in namelist for choosing seabed stress method
JFLemieux73 Jan 27, 2021
41ef64e
Changed gacc to gravit (from icepack) and use of alphab from ice_in
JFLemieux73 Jan 28, 2021
c99a3df
Cleaned up the seabed2 code, added comments
JFLemieux73 Feb 3, 2021
001914f
Changing the doc for new grounding scheme
JFLemieux73 Feb 3, 2021
257e7ac
Almost done with the doc
JFLemieux73 Feb 4, 2021
774d2b1
Minor changes to the doc
JFLemieux73 Feb 8, 2021
8f320e9
Modifs to the doc
JFLemieux73 Feb 11, 2021
a8bcbed
Value of sigmab was not the right one with respect to Dupont et el in…
JFLemieux73 Feb 11, 2021
c6426b3
For some rebase did not work completely...replced basalstress logical…
JFLemieux73 Feb 11, 2021
4bfe465
Other issue with the rebase for ice_init...fixed
JFLemieux73 Feb 11, 2021
16e9b1d
Minor mistake corrected...it compiles
JFLemieux73 Feb 11, 2021
cdc82d1
changed basalstress to seabedstress in options for tests
JFLemieux73 Feb 11, 2021
821ab05
Addressed some of the changes required for the PR
JFLemieux73 Feb 17, 2021
eae232d
Replace kseabed by character string for choice of method
JFLemieux73 Feb 18, 2021
1ddff2e
Updated the doc
JFLemieux73 Feb 18, 2021
f0d5e28
Cosmetic change ice_init
JFLemieux73 Feb 18, 2021
9f689ee
Update doc/source/user_guide/ug_testing.rst
JFLemieux73 Feb 18, 2021
852658e
Cosmetic changes to ice_init to align prints
JFLemieux73 Feb 18, 2021
b6bb0a4
Some modifs requested for PR: same taub after 4 months (gx1)
JFLemieux73 Feb 19, 2021
0a930db
Done with requested changes to new grounding method: same taub after …
JFLemieux73 Feb 19, 2021
173d99e
Modifs to doc and new test of prob method
JFLemieux73 Feb 22, 2021
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
Prev Previous commit
Next Next commit
Changed gacc to gravit (from icepack) and use of alphab from ice_in
  • Loading branch information
@JFLemieux73
JFLemieux73 committed Feb 15, 2021
commit 41ef64e8bbe78fadbf1f584af1cc6433462a48ce
19 changes: 11 additions & 8 deletions cicecore/cicedynB/dynamics/ice_dyn_shared.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -857,7 +857,8 @@ subroutine dyn_finish (nx_block, ny_block, &
end subroutine dyn_finish

!=======================================================================
! Computes seabed (basal) stress factor Tbu (landfast ice)
! Computes seabed (basal) stress factor Tbu (landfast ice) based on mean
! thickness and bathymetry data.
!
! Lemieux, J. F., B. Tremblay, F. Dupont, M. Plante, G.C. Smith, D. Dumont (2015).
! A basal stress parameterization form modeling landfast ice, J. Geophys. Res.
Expand Down Expand Up @@ -931,9 +932,12 @@ subroutine seabed1_stress_factor (nx_block, ny_block, &
end subroutine seabed1_stress_factor

!=======================================================================
! Computes seabed stress due to grounded ridges
! Computes seabed (basal) stress factor Tbu (landfast ice) based on
! probability of contact between the ITD and the seabed.
!
! authors: E. Dumas-Lefebvre, D. Dumont, F. Dupont, JF Lemieux (June 2018)
! Dupont, F. Dumont, D., Lemieux, J. F., Dumas-Lefebvre, E., in prep.
!
! authors: D. Dumont, JF Lemieux, E. Dumas-Lefebvre, F. Dupont
!
subroutine seabed2_stress_factor (nx_block, ny_block, &
icellt, indxti, indxtj, &
Expand Down Expand Up @@ -977,9 +981,7 @@ subroutine seabed2_stress_factor (nx_block, ny_block, &

real (kind=dbl_kind), parameter :: &
max_depth = 50.0_dbl_kind, & ! JFL use threshold_hw!!!!
mu_s = 0.1_dbl_kind, &
gacc = 9.81_dbl_kind, & ! JFL use gravit!!!!
alphab = 20.0_dbl_kind ! alphab=Cb factor in Lemieux et al 2015
mu_s = 0.1_dbl_kind

real (kind=dbl_kind), dimension(ncat_i) :: &
x_k, g_k, P_x ! Center of ice thickness categories
Expand All @@ -1000,11 +1002,12 @@ subroutine seabed2_stress_factor (nx_block, ny_block, &
real (kind=dbl_kind), dimension(ncat_i):: tb_tmp
real (kind=dbl_kind), dimension (nx_block,ny_block):: Tbt
real (kind=dbl_kind) :: atot, x_kmax, x95, x99, x995, x996, x997, x998, x999, xinf
real (kind=dbl_kind) :: cut, rhoi, rhow, pi, puny
real (kind=dbl_kind) :: cut, rhoi, rhow, gravit, pi, puny

character(len=*), parameter :: subname = '(seabed2_stress_coeff)'

call icepack_query_parameters(rhow_out=rhow, rhoi_out=rhoi)
call icepack_query_parameters(gravit_out=gravit)
call icepack_query_parameters(pi_out=pi)
call icepack_query_parameters(puny_out=puny)

Expand Down Expand Up @@ -1081,7 +1084,7 @@ subroutine seabed2_stress_factor (nx_block, ny_block, &
if (ii .eq. 0) then
tb_tmp(n) = c0
else
tb_tmp(n) = max(mu_s*gacc*P_x(n)*sum(P_y(1:ii)*(rhoi*x_k(n) - rhow*y_n(1:ii))), c0)
tb_tmp(n) = max(mu_s*gravit*P_x(n)*sum(P_y(1:ii)*(rhoi*x_k(n) - rhow*y_n(1:ii))),c0)
endif
enddo
Tbt(i,j) = sum(tb_tmp)*exp(-alphab * (c1 - atot))
Expand Down