#2168 closed Bug (fixed)

Volume leak with ice cavities

Reported by: jchanut Owned by: systeam
Priority: low Milestone:
Component: TRA Version: trunk
Severity: minor Keywords: conservation, ISF
Cc:

Description (last modified by nicolasmartin)

Context

Volume is not perfectly conserved with ice shelves cavities. This occurs if the top bbl thickness is greater than the first cell thickness (i.e. rn_hisf_tbl /= 0).

Analysis

MLF corrections on before scale factor are only partly implemented in dynnxt.F90

Fix

Add:

            IF ( ln_isf ) THEN   ! if ice shelf melting
               DO jk = 1, jpkm1 ! Deal with isf separetely, as can be through depth too
                  DO jj = 1, jpj
                     DO ji = 1, jpi
                        IF( misfkt(ji,jj) <=jk .and. jk < misfkb(ji,jj)  ) THEN
                           e3t_b(ji,jj,jk) = e3t_b(ji,jj,jk) - zcoef * ( fwfisf_b(ji,jj) - fwfisf(ji,jj) ) &
                                &          * ( e3t_n(ji,jj,jk) * r1_hisf_tbl(ji,jj) ) * tmask(ji,jj,jk)
                        ELSEIF ( jk==misfkb(ji,jj) ) THEN
                           e3t_b(ji,jj,jk) = e3t_b(ji,jj,jk) - zcoef * ( fwfisf_b(ji,jj) - fwfisf(ji,jj) ) &
                                &          * ( e3t_n(ji,jj,jk) * r1_hisf_tbl(ji,jj) ) * ralpha(ji,jj) * tmask(ji,jj,jk)
                        ENDIF
                     END DO
                  END DO
               END DO
            END IF

instead of

            IF ( ln_isf ) THEN   ! if ice shelf melting
               DO jk = 1, jpkm1 ! Deal with isf separetely, as can be through depth too
                  DO jj = 1, jpj
                     DO ji = 1, jpi
                        IF( misfkt(ji,jj) <=jk .and. jk <= misfkb(ji,jj)  ) THEN
                           e3t_b(ji,jj,jk) = e3t_b(ji,jj,jk) - zcoef * ( fwfisf_b(ji,jj) - fwfisf(ji,jj) ) &
                                &          * ( e3t_n(ji,jj,jk) * r1_hisf_tbl(ji,jj) ) * tmask(ji,jj,jk)
                        ENDIF
                     END DO
                  END DO
               END DO
            END IF


Commit History (1)

ChangesetAuthorTimeChangeLog
10349jchanut2018-11-21T14:59:27+01:00

Fixes volume conservation with isf, #2168

Change History (3)

comment:1 Changed 20 months ago by jchanut

In 10349:

Fixes volume conservation with isf, #2168

comment:2 Changed 20 months ago by nicolasmartin

  • Description modified (diff)

comment:3 Changed 18 months ago by mathiot

  • Resolution set to fixed
  • Status changed from new to closed

I agree with Jerome fix. Has been implemented by Jerome in r10349.

Note: See TracTickets for help on using tickets.