Opened 9 months ago

Closed 9 months ago

Last modified 9 months ago

#2360 closed Bug (fixed)

masking issue in ldftra/ldf_eiv_trp if isf cavities open

Reported by: mathiot Owned by: mathiot
Priority: low Milestone:
Component: LDF Version: trunk
Severity: minor Keywords:
Cc:

Description

Context

zpsi_uw and zpsi_vw seems to be defined at WU and WV points instead of U and V points in ldftra/ldf_eiv_trp. The use of u-,v-mask is not an issue as long as there is no ice shelf cavities.

Analysis

use wumask and wvmask instead of u/v mask.

Fix

Change mask in ldftra/ldf_eiv_trp. sette results do not change as expected.

  • src/OCE/LDF/ldftra.F90

     
    770770         DO jj = 1, jpjm1 
    771771            DO ji = 1, fs_jpim1   ! vector opt. 
    772772               zpsi_uw(ji,jj,jk) = - r1_4 * e2u(ji,jj) * ( wslpi(ji,jj,jk  ) + wslpi(ji+1,jj,jk) )   & 
    773                   &                                    * ( aeiu (ji,jj,jk-1) + aeiu (ji  ,jj,jk) ) * umask(ji,jj,jk) 
     773                  &                                    * ( aeiu (ji,jj,jk-1) + aeiu (ji  ,jj,jk) ) * wumask(ji,jj,jk) 
    774774               zpsi_vw(ji,jj,jk) = - r1_4 * e1v(ji,jj) * ( wslpj(ji,jj,jk  ) + wslpj(ji,jj+1,jk) )   & 
    775                   &                                    * ( aeiv (ji,jj,jk-1) + aeiv (ji,jj  ,jk) ) * vmask(ji,jj,jk) 
     775                  &                                    * ( aeiv (ji,jj,jk-1) + aeiv (ji,jj  ,jk) ) * wvmask(ji,jj,jk) 
    776776            END DO 
    777777         END DO 
    778778      END DO 

Commit History (3)

ChangesetAuthorTimeChangeLog
12298mathiot2019-12-30T20:39:49+01:00

fix #2360

12297mathiot2019-12-30T20:38:11+01:00

fix #2360

12296mathiot2019-12-30T20:34:01+01:00

fix #2360

Change History (3)

comment:1 Changed 9 months ago by mathiot

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

In 12296:

fix #2360

comment:2 Changed 9 months ago by mathiot

In 12297:

fix #2360

comment:3 Changed 9 months ago by mathiot

In 12298:

fix #2360

Note: See TracTickets for help on using tickets.