Changeset 11654

Timestamp:

2019-10-04T15:15:19+02:00 (4 years ago)

Author:

mathiot

Message:

ENHANCE-03_domcfg: forced the tools to write data over the EW overlap region. Issue in some eORCA1 decomposition without (ticket #2143)

Location:

NEMO/branches/2019/ENHANCE-03_domcfg/src

Files:

• iom.F90 (modified) (3 diffs)
• iom_nf90.F90 (modified) (2 diffs)

NEMO/branches/2019/ENHANCE-03_domcfg/src/iom.F90

@@ -613 +613 @@
       ! start halo size for x,y dimensions
       ! end halo size for x,y dimensions
+      !
+      INTEGER ::   nldi_save, nlei_save    !:patch before we remove periodicity and close boundaries in output files
+      INTEGER ::   nldj_save, nlej_save    !:
+      !
       !---------------------------------------------------------------------
       ! Initializations and control
@@ -619 +623 @@
       clinfo = '                    iom_open ~~~  '
       istop = nstop
+
+      ! use patch to force the writing off periodicity and close boundaries
+      ! without this, issue in some model decomposition
+      ! seb: patch before we remove periodicity and close boundaries in output files
+      nldi_save = nldi   ;   nlei_save = nlei
+      nldj_save = nldj   ;   nlej_save = nlej
+      IF( nimpp           ==      1 ) nldi = 1
+      IF( nimpp + jpi - 1 == jpiglo ) nlei = jpi
+      IF( njmpp           ==      1 ) nldj = 1
+      IF( njmpp + jpj - 1 == jpjglo ) nlej = jpj
+
       ! if iom_open is called for the first time: initialize iom_file(:)%nfid to 0
       ! (could be done when defining iom_file in f95 but not in f90)
@@ -743 +758 @@
          CALL iom_nf90_open( clname, kiomid, llwrt, llok, idompar, kdlev = kdlev )
       ENDIF
+
+      nldi = nldi_save   ;   nlei = nlei_save
+      nldj = nldj_save   ;   nlej = nlej_save
       !
    END SUBROUTINE iom_open

NEMO/branches/2019/ENHANCE-03_domcfg/src/iom_nf90.F90

@@ -529 +529 @@
       INTEGER               :: idlv                 ! local variable
       INTEGER               :: idim3                ! id of the third dimension
+      !
+      INTEGER ::   nldi_save, nlei_save    !:patch before we remove periodicity and close boundaries in output files
+      INTEGER ::   nldj_save, nlej_save    !:
       !---------------------------------------------------------------------
       !
       clinfo = '          iom_nf90_rp0123d, file: '//TRIM(iom_file(kiomid)%name)//', var: '//TRIM(cdvar)
       if90id = iom_file(kiomid)%nfid
+      !
+      ! use patch to force the writing off periodicity and close boundaries
+      ! without this, issue in some model decomposition
+      ! seb: patch before we remove periodicity and close boundaries in output files
+      nldi_save = nldi   ;   nlei_save = nlei
+      nldj_save = nldj   ;   nlej_save = nlej
+      IF( nimpp           ==      1 ) nldi = 1
+      IF( nimpp + jpi - 1 == jpiglo ) nlei = jpi
+      IF( njmpp           ==      1 ) nldj = 1
+      IF( njmpp + jpj - 1 == jpjglo ) nlej = jpj
       !
       ! define dimension variables if it is not already done
@@ -705 +718 @@
          IF(lwp) WRITE(numout,*) TRIM(clinfo)//' written ok'
       ENDIF
+      !
+      nldi = nldi_save   ;   nlei = nlei_save
+      nldj = nldj_save   ;   nlej = nlej_save
       !     
    END SUBROUTINE iom_nf90_rp0123d