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

update icepack to support snowtable better #1

Merged
merged 1 commit into from
Jun 11, 2021
Merged

update icepack to support snowtable better #1

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
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
187 changes: 157 additions & 30 deletions columnphysics/icepack_parameters.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -7,7 +7,8 @@
module icepack_parameters

use icepack_kinds
use icepack_warnings, only: icepack_warnings_aborted
use icepack_warnings, only: icepack_warnings_aborted, &
icepack_warnings_add, icepack_warnings_setabort

implicit none
private
Expand Down Expand Up @@ -317,7 +318,7 @@ module icepack_parameters

character (len=char_len), public :: &
snwredist = 'none', & ! type of snow redistribution
snw_aging_table = 'test' ! lookup table: 'snicar' or 'test'
snw_aging_table = 'test' ! lookup table: 'snicar' or 'test' or 'file'

logical (kind=log_kind), public :: &
use_smliq_pnd = .false. , & ! use liquid in snow for ponds
Expand All @@ -340,10 +341,14 @@ module icepack_parameters
isnw_rhos ! maximum snow density index

! dry snow aging parameters
real (kind=dbl_kind), dimension(:), allocatable, public :: &
snowage_rhos, & ! snowage table dimension data for rhos (kg/m^3)
snowage_Tgrd, & ! snowage table dimension data for temp gradient (deg K/m)
snowage_T ! snowage table dimension data for temperature (deg K)
real (kind=dbl_kind), dimension(:,:,:), allocatable, public :: &
snowage_tau, & ! (10^-6 m)
snowage_kappa, & !
snowage_drdt0 ! (10^-6 m/hr)
snowage_tau, & ! snowage table 3D data for tau (10^-6 m)
snowage_kappa, & ! snowage table 3D data for kappa (10^-6 m)
snowage_drdt0 ! snowage table 3D data for drdt0 (10^-6 m/hr)

!-----------------------------------------------------------------------
! Parameters for biogeochemistry
Expand Down Expand Up @@ -454,6 +459,7 @@ subroutine icepack_init_parameters( &
snwredist_in, use_smliq_pnd_in, rsnw_fall_in, rsnw_tmax_in, &
rhosnew_in, rhosmin_in, rhosmax_in, windmin_in, drhosdwind_in, &
snwlvlfac_in, isnw_T_in, isnw_Tgrd_in, isnw_rhos_in, &
snowage_rhos_in, snowage_Tgrd_in, snowage_T_in, &
snowage_tau_in, snowage_kappa_in, snowage_drdt0_in, &
snw_aging_table_in)

Expand Down Expand Up @@ -513,7 +519,7 @@ subroutine icepack_init_parameters( &
! 1 = Bitz and Lipscomb 1999
! 2 = mushy layer theory

character (char_len), intent(in), optional :: &
character (len=*), intent(in), optional :: &
conduct_in, & ! 'MU71' or 'bubbly'
fbot_xfer_type_in ! transfer coefficient type for ice-ocean heat flux

Expand All @@ -538,7 +544,7 @@ subroutine icepack_init_parameters( &
phi_c_slow_mode_in , & ! liquid fraction porosity cutoff for slow mode
phi_i_mushy_in ! liquid fraction of congelation ice

character(len=char_len), intent(in), optional :: &
character(len=*), intent(in), optional :: &
tfrz_option_in ! form of ocean freezing temperature
! 'minus1p8' = -1.8 C
! 'linear_salt' = -depressT * sss
Expand All @@ -561,7 +567,7 @@ subroutine icepack_init_parameters( &
awtvdf_in, & ! visible, diffuse
awtidf_in ! near IR, diffuse

character (len=char_len), intent(in), optional :: &
character (len=*), intent(in), optional :: &
shortwave_in, & ! shortwave method, 'ccsm3' or 'dEdd'
albedo_type_in ! albedo parameterization, 'ccsm3' or 'constant'
! shortwave='dEdd' overrides this parameter
Expand Down Expand Up @@ -629,7 +635,7 @@ subroutine icepack_init_parameters( &
qqqocn_in, & ! for qsat over ocn
TTTocn_in ! for qsat over ocn

character (len=char_len), intent(in), optional :: &
character (len=*), intent(in), optional :: &
atmbndy_in ! atmo boundary method, 'default' ('ccsm3') or 'constant'

logical (kind=log_kind), intent(in), optional :: &
Expand Down Expand Up @@ -669,14 +675,14 @@ subroutine icepack_init_parameters( &
logical (kind=log_kind), intent(in), optional :: &
wave_spec_in ! if true, use wave forcing

character (len=char_len), intent(in), optional :: &
character (len=*), intent(in), optional :: &
wave_spec_type_in ! type of wave spectrum forcing

!-----------------------------------------------------------------------
! Parameters for biogeochemistry
!-----------------------------------------------------------------------

character(char_len), intent(in), optional :: &
character (len=*), intent(in), optional :: &
bgc_flux_type_in ! type of ocean-ice piston velocity
! 'constant', 'Jin2006'

Expand Down Expand Up @@ -737,7 +743,7 @@ subroutine icepack_init_parameters( &
hs0_in ! snow depth for transition to bare sea ice (m)

! level-ice ponds
character (len=char_len), intent(in), optional :: &
character (len=*), intent(in), optional :: &
frzpnd_in ! pond refreezing parameterization

real (kind=dbl_kind), intent(in), optional :: &
Expand All @@ -755,7 +761,7 @@ subroutine icepack_init_parameters( &
! Parameters for snow redistribution, metamorphosis
!-----------------------------------------------------------------------

character (len=char_len), intent(in), optional :: &
character (len=*), intent(in), optional :: &
snwredist_in, & ! type of snow redistribution
snw_aging_table_in ! snow aging lookup table

Expand All @@ -771,11 +777,18 @@ subroutine icepack_init_parameters( &
rhosmax_in, & ! maximum snow density (kg/m^3)
windmin_in, & ! minimum wind speed to compact snow (m/s)
drhosdwind_in, & ! wind compaction factor (kg s/m^4)
snwlvlfac_in, & ! fractional increase in snow depth
snwlvlfac_in ! fractional increase in snow depth

integer (kind=int_kind), intent(in), optional :: &
eclare108213 marked this conversation as resolved.
Show resolved Hide resolved
isnw_T_in, & ! maxiumum temperature index
isnw_Tgrd_in, & ! maxiumum temperature gradient index
isnw_rhos_in ! maxiumum snow density index

real (kind=dbl_kind), dimension(:), intent(in), optional :: &
snowage_rhos_in, & ! snowage dimension data
snowage_Tgrd_in, & !
snowage_T_in !

real (kind=dbl_kind), dimension(:,:,:), intent(in), optional :: &
snowage_tau_in, & ! (10^-6 m)
snowage_kappa_in, &!
Expand Down Expand Up @@ -911,9 +924,109 @@ subroutine icepack_init_parameters( &
if (present(isnw_T_in) ) isnw_T = isnw_T_in
if (present(isnw_Tgrd_in) ) isnw_Tgrd = isnw_Tgrd_in
if (present(isnw_rhos_in) ) isnw_rhos = isnw_rhos_in
if (present(snowage_tau_in) ) snowage_tau = snowage_tau_in
if (present(snowage_kappa_in) ) snowage_kappa = snowage_kappa_in
if (present(snowage_drdt0_in) ) snowage_drdt0 = snowage_drdt0_in

! check array sizes and re/allocate if necessary
if (present(snowage_rhos_in) ) then
if (size(snowage_rhos_in) /= isnw_rhos) then
call icepack_warnings_add(subname//' incorrect size of snowage_rhos_in')
call icepack_warnings_setabort(.true.,__FILE__,__LINE__)
elseif (.not.allocated(snowage_rhos)) then
allocate(snowage_rhos(isnw_rhos))
snowage_rhos = snowage_rhos_in
elseif (size(snowage_rhos) /= isnw_rhos) then
deallocate(snowage_rhos)
allocate(snowage_rhos(isnw_rhos))
snowage_rhos = snowage_rhos_in
else
snowage_rhos = snowage_rhos_in
endif
endif

! check array sizes and re/allocate if necessary
if (present(snowage_Tgrd_in) ) then
if (size(snowage_Tgrd_in) /= isnw_Tgrd) then
call icepack_warnings_add(subname//' incorrect size of snowage_Tgrd_in')
call icepack_warnings_setabort(.true.,__FILE__,__LINE__)
elseif (.not.allocated(snowage_Tgrd)) then
allocate(snowage_Tgrd(isnw_Tgrd))
snowage_Tgrd = snowage_Tgrd_in
elseif (size(snowage_Tgrd) /= isnw_Tgrd) then
deallocate(snowage_Tgrd)
allocate(snowage_Tgrd(isnw_Tgrd))
snowage_Tgrd = snowage_Tgrd_in
else
snowage_Tgrd = snowage_Tgrd_in
endif
endif

! check array sizes and re/allocate if necessary
if (present(snowage_T_in) ) then
if (size(snowage_T_in) /= isnw_T) then
call icepack_warnings_add(subname//' incorrect size of snowage_T_in')
call icepack_warnings_setabort(.true.,__FILE__,__LINE__)
elseif (.not.allocated(snowage_T)) then
allocate(snowage_T(isnw_T))
snowage_T = snowage_T_in
elseif (size(snowage_T) /= isnw_T) then
deallocate(snowage_T)
allocate(snowage_T(isnw_T))
snowage_T = snowage_T_in
else
snowage_T = snowage_T_in
endif
endif

! check array sizes and re/allocate if necessary
if (present(snowage_tau_in) ) then
if (size(snowage_tau_in) /= isnw_T*isnw_Tgrd*isnw_rhos) then
call icepack_warnings_add(subname//' incorrect size of snowage_tau_in')
call icepack_warnings_setabort(.true.,__FILE__,__LINE__)
elseif (.not.allocated(snowage_tau)) then
allocate(snowage_tau(isnw_T,isnw_Tgrd,isnw_rhos))
snowage_tau = snowage_tau_in
elseif (size(snowage_tau) /= isnw_T*isnw_Tgrd*isnw_rhos) then
deallocate(snowage_tau)
allocate(snowage_tau(isnw_T,isnw_Tgrd,isnw_rhos))
snowage_tau = snowage_tau_in
else
snowage_tau = snowage_tau_in
endif
endif

! check array sizes and re/allocate if necessary
if (present(snowage_kappa_in) ) then
if (size(snowage_kappa_in) /= isnw_T*isnw_Tgrd*isnw_rhos) then
call icepack_warnings_add(subname//' incorrect size of snowage_kappa_in')
call icepack_warnings_setabort(.true.,__FILE__,__LINE__)
elseif (.not.allocated(snowage_kappa)) then
allocate(snowage_kappa(isnw_T,isnw_Tgrd,isnw_rhos))
snowage_kappa = snowage_kappa_in
elseif (size(snowage_kappa) /= isnw_T*isnw_Tgrd*isnw_rhos) then
deallocate(snowage_kappa)
allocate(snowage_kappa(isnw_T,isnw_Tgrd,isnw_rhos))
snowage_kappa = snowage_kappa_in
else
snowage_kappa = snowage_kappa_in
endif
endif

! check array sizes and re/allocate if necessary
if (present(snowage_drdt0_in) ) then
if (size(snowage_drdt0_in) /= isnw_T*isnw_Tgrd*isnw_rhos) then
call icepack_warnings_add(subname//' incorrect size of snowage_drdt0_in')
call icepack_warnings_setabort(.true.,__FILE__,__LINE__)
elseif (.not.allocated(snowage_drdt0)) then
allocate(snowage_drdt0(isnw_T,isnw_Tgrd,isnw_rhos))
snowage_drdt0 = snowage_drdt0_in
elseif (size(snowage_drdt0) /= isnw_T*isnw_Tgrd*isnw_rhos) then
deallocate(snowage_drdt0)
allocate(snowage_drdt0(isnw_T,isnw_Tgrd,isnw_rhos))
snowage_drdt0 = snowage_drdt0_in
else
snowage_drdt0 = snowage_drdt0_in
endif
endif

if (present(bgc_flux_type_in) ) bgc_flux_type = bgc_flux_type_in
if (present(z_tracers_in) ) z_tracers = z_tracers_in
if (present(scale_bgc_in) ) scale_bgc = scale_bgc_in
Expand Down Expand Up @@ -1009,6 +1122,7 @@ subroutine icepack_query_parameters( &
snwredist_out, use_smliq_pnd_out, rsnw_fall_out, rsnw_tmax_out, &
rhosnew_out, rhosmin_out, rhosmax_out, windmin_out, drhosdwind_out, &
snwlvlfac_out, isnw_T_out, isnw_Tgrd_out, isnw_rhos_out, &
snowage_rhos_out, snowage_Tgrd_out, snowage_T_out, &
snowage_tau_out, snowage_kappa_out, snowage_drdt0_out, &
snw_aging_table_out)

Expand Down Expand Up @@ -1077,7 +1191,7 @@ subroutine icepack_query_parameters( &
! 1 = Bitz and Lipscomb 1999
! 2 = mushy layer theory

character (char_len), intent(out), optional :: &
character (len=*), intent(out), optional :: &
conduct_out, & ! 'MU71' or 'bubbly'
fbot_xfer_type_out ! transfer coefficient type for ice-ocean heat flux

Expand All @@ -1102,7 +1216,7 @@ subroutine icepack_query_parameters( &
phi_c_slow_mode_out , & ! liquid fraction porosity cutoff for slow mode
phi_i_mushy_out ! liquid fraction of congelation ice

character(len=char_len), intent(out), optional :: &
character(len=*), intent(out), optional :: &
tfrz_option_out ! form of ocean freezing temperature
! 'minus1p8' = -1.8 C
! 'linear_salt' = -depressT * sss
Expand All @@ -1125,7 +1239,7 @@ subroutine icepack_query_parameters( &
awtvdf_out, & ! visible, diffuse
awtidf_out ! near IR, diffuse

character (len=char_len), intent(out), optional :: &
character (len=*), intent(out), optional :: &
shortwave_out, & ! shortwave method, 'ccsm3' or 'dEdd'
albedo_type_out ! albedo parameterization, 'ccsm3' or 'constant'
! shortwave='dEdd' overrides this parameter
Expand Down Expand Up @@ -1193,7 +1307,7 @@ subroutine icepack_query_parameters( &
qqqocn_out, & ! for qsat over ocn
TTTocn_out ! for qsat over ocn

character (len=char_len), intent(out), optional :: &
character (len=*), intent(out), optional :: &
atmbndy_out ! atmo boundary method, 'default' ('ccsm3') or 'constant'

logical (kind=log_kind), intent(out), optional :: &
Expand Down Expand Up @@ -1233,15 +1347,15 @@ subroutine icepack_query_parameters( &
logical (kind=log_kind), intent(out), optional :: &
wave_spec_out ! if true, use wave forcing

character (len=char_len), intent(out), optional :: &
character (len=*), intent(out), optional :: &
wave_spec_type_out ! type of wave spectrum forcing

!-----------------------------------------------------------------------
! Parameters for biogeochemistry
!-----------------------------------------------------------------------

character(char_len), intent(out), optional :: &
bgc_flux_type_out ! type of ocean-ice piston velocity
character (len=*), intent(out), optional :: &
bgc_flux_type_out ! type of ocean-ice piston velocity
! 'constant', 'Jin2006'

logical (kind=log_kind), intent(out), optional :: &
Expand Down Expand Up @@ -1301,7 +1415,7 @@ subroutine icepack_query_parameters( &
hs0_out ! snow depth for transition to bare sea ice (m)

! level-ice ponds
character (len=char_len), intent(out), optional :: &
character (len=*), intent(out), optional :: &
frzpnd_out ! pond refreezing parameterization

real (kind=dbl_kind), intent(out), optional :: &
Expand All @@ -1319,7 +1433,7 @@ subroutine icepack_query_parameters( &
! Parameters for snow redistribution, metamorphosis
!-----------------------------------------------------------------------

character (len=char_len), intent(out), optional :: &
character (len=*), intent(out), optional :: &
snwredist_out, & ! type of snow redistribution
snw_aging_table_out ! snow aging lookup table

Expand All @@ -1335,11 +1449,18 @@ subroutine icepack_query_parameters( &
rhosmax_out, & ! maximum snow density (kg/m^3)
windmin_out, & ! minimum wind speed to compact snow (m/s)
drhosdwind_out, & ! wind compaction factor (kg s/m^4)
snwlvlfac_out, & ! fractional increase in snow depth
snwlvlfac_out ! fractional increase in snow depth

integer (kind=int_kind), intent(out), optional :: &
isnw_T_out, & ! maxiumum temperature index
isnw_Tgrd_out, & ! maxiumum temperature gradient index
isnw_rhos_out ! maxiumum snow density index

real (kind=dbl_kind), dimension(:), intent(out), optional :: &
snowage_rhos_out, & ! snowage dimension data
snowage_Tgrd_out, & !
snowage_T_out !

real (kind=dbl_kind), dimension(:,:,:), intent(out), optional :: &
snowage_tau_out, & ! (10^-6 m)
snowage_kappa_out, &!
Expand Down Expand Up @@ -1515,6 +1636,9 @@ subroutine icepack_query_parameters( &
if (present(isnw_T_out) ) isnw_T_out = isnw_T
if (present(isnw_Tgrd_out) ) isnw_Tgrd_out = isnw_Tgrd
if (present(isnw_rhos_out) ) isnw_rhos_out = isnw_rhos
if (present(snowage_rhos_out) ) snowage_rhos_out = snowage_rhos
if (present(snowage_Tgrd_out) ) snowage_Tgrd_out = snowage_Tgrd
if (present(snowage_T_out) ) snowage_T_out = snowage_T
if (present(snowage_tau_out) ) snowage_tau_out = snowage_tau
if (present(snowage_kappa_out) ) snowage_kappa_out= snowage_kappa
if (present(snowage_drdt0_out) ) snowage_drdt0_out= snowage_drdt0
Expand Down Expand Up @@ -1713,9 +1837,12 @@ subroutine icepack_write_parameters(iounit)
write(iounit,*) " isnw_T = ", isnw_T
write(iounit,*) " isnw_Tgrd = ", isnw_Tgrd
write(iounit,*) " isnw_rhos = ", isnw_rhos
write(iounit,*) " snowage_tau = ", snowage_tau
write(iounit,*) " snowage_kappa = ", snowage_kappa
write(iounit,*) " snowage_drdt0 = ", snowage_drdt0
write(iounit,*) " snowage_rhos = ", snowage_rhos(1)
write(iounit,*) " snowage_Tgrd = ", snowage_Tgrd(1)
write(iounit,*) " snowage_T = ", snowage_T(1)
write(iounit,*) " snowage_tau = ", snowage_tau(1,1,1)
write(iounit,*) " snowage_kappa = ", snowage_kappa(1,1,1)
write(iounit,*) " snowage_drdt0 = ", snowage_drdt0(1,1,1)
write(iounit,*) " bgc_flux_type = ", bgc_flux_type
write(iounit,*) " z_tracers = ", z_tracers
write(iounit,*) " scale_bgc = ", scale_bgc
Expand Down
1 change: 1 addition & 0 deletions columnphysics/icepack_shortwave.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1050,6 +1050,7 @@ subroutine run_dEdd(dt, ncat, &
hs0, hsnlvl, &
rhosnwn(:), rsnwn(:), &
rsnow(:,n))
if (icepack_warnings_aborted(subname)) return
endif ! snwredist

fpn = c0 ! fraction of ice covered in pond
Expand Down
Loading