Skip to content

Commit

Permalink
trunk:
Browse files Browse the repository at this point in the history
Bug fix to shortwave module, to ensure fractions of snow, ponds and bare ice add to 1.
Not BFB.

Updated project code for LANL institutional computing machines.
  • Loading branch information
eclare108213 committed Apr 8, 2016
1 parent faa7189 commit 8052b8a
Show file tree
Hide file tree
Showing 2 changed files with 6 additions and 3 deletions.
2 changes: 1 addition & 1 deletion cice/input_templates/run_ice.Linux.LANL.conejo
Original file line number Diff line number Diff line change
@@ -1,5 +1,5 @@
#!/bin/tcsh
#MSUB -A s11_climate
#MSUB -A s11_climateacme
#MSUB -l walltime=01:00:00
#MSUB -l nodes=1:ppn=4
#MSUB -N cice
Expand Down
7 changes: 5 additions & 2 deletions cice/source/ice_shortwave.F90
Original file line number Diff line number Diff line change
Expand Up @@ -1345,10 +1345,13 @@ subroutine run_dEdd(ilo,ihi,jlo,jhi, &
hpn(i,j) = hp * tmp
fpn(i,j) = fpn(i,j) * tmp
endif
fsn(i,j) = min(fsn(i,j), c1-fpn(i,j))

endif ! hp > puny

! Zero out fraction of thin ponds for radiation only
if (hpn(i,j) < hpmin) fpn(i,j) = c0
fsn(i,j) = min(fsn(i,j), c1-fpn(i,j))

! endif ! masking by lid ice
apeffn(i,j,n) = fpn(i,j) ! for history
enddo ! ij
Expand Down Expand Up @@ -1762,7 +1765,7 @@ subroutine shortwave_dEdd (nx_block, ny_block, &
coszen(i,j) = max(puny,coszen(i,j))
hi(i,j) = vice(i,j) / aice(i,j)
! if non-zero pond fraction and sufficient pond depth
if( fp(i,j) > puny .and. hp(i,j) > hpmin ) then
if (fp(i,j) > puny) then
icells_DE = icells_DE + 1
indxi_DE(icells_DE) = i
indxj_DE(icells_DE) = j
Expand Down

0 comments on commit 8052b8a

Please sign in to comment.