CESM ocean forcing - hmix division by zero #585
Comments
|
Good catch. I will work on a fix. It should be a max(c20,work1(i,j,iblk)). |
apcraig
added a commit
to apcraig/CICE
that referenced
this issue
Mar 27, 2021
- update set_nml.gx1prod to match current production system - add apr 1 test case for gx1 - update landice tests, use gx1 - fix hmix default CICE-Consortium#585 - delete old code in ice_forcing.F90 (should have been done in earlier PR)
16 tasks
apcraig
added a commit
that referenced
this issue
Apr 6, 2021
…efault value (#586) * Update gx1 ic - update set_nml.gx1prod to match current production system - add apr 1 test case for gx1 - update landice tests, use gx1 - fix hmix default #585 - delete old code in ice_forcing.F90 (should have been done in earlier PR) * update gx1coreii initial condition * update icepack * add gx1 debug test, expected to fail * Update version number to 6.2.0 * Update gx3 and gx1 input filenames Update CI wget implementation #588 Create variable in ice_forcing.F90 called mixed_layer_depth_default and use it instead of c20 Fix bug in hourly output, created in time manager update #589 Set start year for all runs to 2005 and turn leap years on by default Update documentation in history output section to add information about hist_avg namelist, noted in #566 * update ic filenames, tests, and documentation * update color links in test results wiki * update hmix initialization
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
In icepack_ocean.F90 line 103.
sst = sst + dt * ( &
(fsens_ocn + flat_ocn + flwout_ocn + flw + swabs) * (c1-aice) &
+ fhocn + fswthru) & ! these are aice
/ (cprhohmix)
hmix is zero when using standard ocean forcing (CESM) in standalone simulation.
This is initialized beginning line 4394 in ice_forcing.F90
Below is a hard coded fix (hmix=c20) instead of the value from the forcing file.
do j = 1, ny_block
do i = 1, nx_block
if (n == 2) sss (i,j,:) = c0
if (n == 3) hmix (i,j,:) = c0
if (n == 4) uocn (i,j,:) = c0
if (n == 5) vocn (i,j,:) = c0
if (n == 6) ss_tltx(i,j,:) = c0
if (n == 7) ss_tlty(i,j,:) = c0
if (n == 8) qdp (i,j,:) = c0
do iblk = 1, nblocks
if (hm(i,j,iblk) == c1) then
if (n == 2) sss (i,j,iblk) = work1(i,j,iblk)
if (n == 3) hmix (i,j,iblk) = c20 !!work1(i,j,iblk)
if (n == 4) uocn (i,j,iblk) = work1(i,j,iblk)
if (n == 5) vocn (i,j,iblk) = work1(i,j,iblk)
if (n == 6) ss_tltx(i,j,iblk) = work1(i,j,iblk)
if (n == 7) ss_tlty(i,j,iblk) = work1(i,j,iblk)
if (n == 8) qdp (i,j,iblk) = work1(i,j,iblk)
endif
enddo
enddo
enddo
enddo
The text was updated successfully, but these errors were encountered: