New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
#2534 (Extra halo rows in GYRE) – NEMO

Opened 4 years ago

Closed 3 years ago

Last modified 2 years ago

#2534 closed Defect (fixed)

Extra halo rows in GYRE

Reported by: hadcv Owned by: hadcv
Priority: low Milestone: Documentation
Component: USR Version: trunk
Severity: minor Keywords: GYRE doc
Cc:

Description

GYRE includes hard-coded halo rows, see usrdef_nam.F90:

kpi = 30 * nn_GYRE + 2       !                     
kpj = 20 * nn_GYRE + 2

and usrdef_hgr.F90:

zlam0 = zlam1 + zcos_alpha * ze1deg * REAL( Ni0glo - 2, wp )
zphi0 = zphi1 + zsin_alpha * ze1deg * REAL( Nj0glo - 2, wp )

Although a trivial issue, the code and documentation now differ:

! usrdef_nam.F90
Ni0glo = 30 * nn_GYRE + 2                 
Nj0glo = 20 * nn_GYRE + 2

! mppini.F90
jpiglo = Ni0glo + 2 * nn_hls
jpjglo = Nj0glo + 2 * nn_hls

! Documentation
jpiglo = 30 * nn_GYRE + 2                 
jpjglo = 20 * nn_GYRE + 2

Recommendation

Remove these halo points from usrdef_nam.F90:

- kpi = 30 * nn_GYRE + 2       !                     
- kpj = 20 * nn_GYRE + 2
+ kpi = 30 * nn_GYRE       !                     
+ kpj = 20 * nn_GYRE 

and usrdef_hgr.F90:

- zlam0 = zlam1 + zcos_alpha * ze1deg * REAL( Ni0glo - 2, wp )
- zphi0 = zphi1 + zsin_alpha * ze1deg * REAL( Nj0glo - 2, wp )
+ zlam0 = zlam1 + zcos_alpha * ze1deg * REAL( Ni0glo, wp )
+ zphi0 = zphi1 + zsin_alpha * ze1deg * REAL( Nj0glo, wp )

and update the documentation:

jpiglo = 30 * nn_GYRE + 2 * nn_hls
jpjglo = 20 * nn_GYRE + 2 * nn_hls

Commit History (1)

ChangesetAuthorTimeChangeLog
14196hadcv2020-12-16T15:38:32+01:00

#2534: Minor update to GYRE documentation

Change History (4)

comment:1 Changed 3 years ago by hadcv

  • Status changed from new to assigned

This will not be fixed, as the fix would decrease the size of the internal domain (Ni0glo/Nj0glo) and change the results. The documentation still needs updating, however.

jpiglo = 30 * nn_GYRE + 2 + 2 * nn_hls
jpjglo = 20 * nn_GYRE + 2 + 2 * nn_hls
Last edited 3 years ago by hadcv (previous) (diff)

comment:2 Changed 3 years ago by hadcv

In 14196:

Error: Failed to load processor CommitTicketReference
No macro or processor named 'CommitTicketReference' found

comment:3 Changed 3 years ago by hadcv

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

comment:4 Changed 2 years ago by nemo

  • Keywords doc added
Note: See TracTickets for help on using tickets.