Skip to content

Bug fix for box2001: aiu now multiplies the air stress components #525

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

Merged
apcraig merged 1 commit into from
Oct 22, 2020
Merged

Bug fix for box2001: aiu now multiplies the air stress components #525

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
14 changes: 11 additions & 3 deletions cicecore/cicedynB/general/ice_forcing.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -5059,23 +5059,31 @@ end subroutine ocn_data_ispol_init
subroutine box2001_data

! wind and current fields as in Hunke, JCP 2001
! these are defined at the u point
! authors: Elizabeth Hunke, LANL

use ice_domain, only: nblocks
use ice_domain_size, only: max_blocks
use ice_blocks, only: nx_block, ny_block, nghost
use ice_flux, only: uocn, vocn, uatm, vatm, wind, rhoa, strax, stray
use ice_grid, only: uvm
use ice_grid, only: uvm, to_ugrid
use ice_state, only: aice_init

! local parameters

integer (kind=int_kind) :: &
iblk, i,j ! loop indices

real (kind=dbl_kind), dimension (nx_block,ny_block,max_blocks) :: &
aiu ! ice fraction on u-grid

real (kind=dbl_kind) :: &
secday, pi , puny, period, pi2, tau
call icepack_query_parameters(pi_out=pi, pi2_out=pi2, puny_out=puny)
call icepack_query_parameters(secday_out=secday)

call to_ugrid(aice_init, aiu)

period = c4*secday

do iblk = 1, nblocks
Expand Down Expand Up @@ -5106,8 +5114,8 @@ subroutine box2001_data
! wind stress
wind(i,j,iblk) = sqrt(uatm(i,j,iblk)**2 + vatm(i,j,iblk)**2)
tau = rhoa(i,j,iblk) * 0.0012_dbl_kind * wind(i,j,iblk)
strax(i,j,iblk) = tau * uatm(i,j,iblk)
stray(i,j,iblk) = tau * vatm(i,j,iblk)
strax(i,j,iblk) = aiu(i,j,iblk) * tau * uatm(i,j,iblk)
stray(i,j,iblk) = aiu(i,j,iblk) * tau * vatm(i,j,iblk)

! initialization test
! Diagonal wind vectors 1
Expand Down