Skip to content

Commit

Permalink
Add gbox12 and ice_data_type='smallblock' (CICE-Consortium#51)
Browse files Browse the repository at this point in the history
* add gbox12 and ice_data_type='smallblock'

* update documentation

* update gbox12 default
  • Loading branch information
apcraig authored Jan 25, 2022
1 parent 311302a commit 14f42ab
Show file tree
Hide file tree
Showing 4 changed files with 34 additions and 2 deletions.
16 changes: 15 additions & 1 deletion cicecore/cicedynB/general/ice_init.F90
Original file line number Diff line number Diff line change
Expand Up @@ -2687,7 +2687,7 @@ subroutine set_state_var (nx_block, ny_block, &

if (trim(ice_ic) == 'default') then

if (trim(ice_data_type) == 'box2001') then
if (trim(ice_data_type) == 'box2001' .or. trim(ice_data_type) == 'smallblock') then

hbar = c2 ! initial ice thickness
do n = 1, ncat
Expand Down Expand Up @@ -2764,6 +2764,20 @@ subroutine set_state_var (nx_block, ny_block, &
endif
enddo
enddo

else if (trim(ice_data_type) == 'smallblock') then
! 2x2 ice in center of domain
icells = 0
do j = jlo, jhi
do i = ilo, ihi
if ((iglob(i) == nx_global/2 .or. iglob(i) == nx_global/2+1) .and. &
(jglob(j) == ny_global/2 .or. jglob(j) == ny_global/2+1)) then
icells = icells + 1
indxi(icells) = i
indxj(icells) = j
endif
enddo
enddo

else ! default behavior

Expand Down
11 changes: 11 additions & 0 deletions configuration/scripts/cice_decomp.csh
Original file line number Diff line number Diff line change
Expand Up @@ -66,6 +66,17 @@ else if (${grid} == 'gbox80') then
set blckx = 8; set blcky = 8
endif

else if (${grid} == 'gbox12') then
set nxglob = 12
set nyglob = 12
if (${cicepes} <= 1) then
set blckx = 12; set blcky = 12
else if (${cicepes} <= 8) then
set blckx = 4; set blcky = 4
else
set blckx = 2; set blcky = 2
endif

else if (${grid} == 'gx3') then
set nxglob = 100
set nyglob = 116
Expand Down
5 changes: 5 additions & 0 deletions configuration/scripts/options/set_nml.gbox12
Original file line number Diff line number Diff line change
@@ -0,0 +1,5 @@
ice_ic = 'default'
grid_type = 'rectangular'
atm_data_type = 'box2001'
ocn_data_type = 'calm'
ice_data_type = 'box2001'
4 changes: 3 additions & 1 deletion doc/source/user_guide/ug_case_settings.rst
Original file line number Diff line number Diff line change
Expand Up @@ -577,7 +577,9 @@ forcing_nml
"``highfreq``", "logical", "high-frequency atmo coupling", "``.false.``"
"``ice_data_type``", "``boxslotcyl``", "initialize ice concentration and velocity for :ref:`boxslotcyl` test (:cite:`Zalesak79`)", "``default``"
"", "``box2001``", "initialize ice concentration for :ref:`box2001` test (:cite:`Hunke01`)", ""
"", "``default``", "no special initialization", ""
"", "``default``", "ice dependent on latitude and ocean temperature", ""
"", "``smallblock``", "uniform 2x2 block ice concentration and thickness in center of domain", ""
"", "``uniform``", "uniform ice concentration and thickness across domain", ""
"``iceruf``", "real", "ice surface roughness at atmosphere interface", "0.0005"
"``l_mpond_fresh``", "``.false.``", "release pond water immediately to ocean", "``.false.``"
"", "``true``", "retain (topo) pond water until ponds drain", ""
Expand Down

0 comments on commit 14f42ab

Please sign in to comment.