Ignore:
Timestamp:
2020-05-22T09:05:34+02:00 (8 months ago)
Author:
smasson
Message:

Extra_Halo: additional bugfixes and developments, see #2366

Location:
NEMO/branches/2020/dev_r12558_HPC-08_epico_Extra_Halo/tests
Files:
15 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/dev_r12558_HPC-08_epico_Extra_Halo/tests/BENCH/MY_SRC/usrdef_hgr.F90

    r12939 r12960  
    7676      DO_2D_11_11 
    7777          
    78          zti = REAL( mig0(ji) - 1, wp )   ! start at i=0 in the global grid without halos 
    79          ztj = REAL( mjg0(jj) - 1, wp )   ! start at j=0 in the global grid without halos 
     78         zti = REAL( mig0_oldcmp(ji) - 1, wp )   ! start at i=0 in the global grid without halos 
     79         ztj = REAL( mjg0_oldcmp(jj) - 1, wp )   ! start at j=0 in the global grid without halos 
    8080          
    8181         plamt(ji,jj) = zti 
  • NEMO/branches/2020/dev_r12558_HPC-08_epico_Extra_Halo/tests/BENCH/MY_SRC/usrdef_istate.F90

    r12939 r12960  
    7070      ! Once we decide to forget trunk compatibility, we must simply define z2d as: 
    7171!!$      DO_2D_00_00 
    72 !!$         z2d(ji,jj) = 0.1 * ( 0.5 - REAL( mig0(ji) + (mjg0(jj)-1) * Ni0glo, wp ) / REAL( Ni0glo * Nj0glo, wp ) ) 
     72!!$         z2d(ji,jj) = 0.1 * ( 0.5 - REAL( mig0_oldcmp(ji) + (mjg0_oldcmp(jj)-1) * Ni0glo, wp ) / REAL( Ni0glo * Nj0glo, wp ) ) 
    7373!!$      END_2D 
    7474      igloi = Ni0glo + 2 * COUNT( (/ jperio == 1 .OR. jperio == 4 .OR. jperio == 6 .OR. jperio == 7 /) ) 
    7575      igloj = Nj0glo + 2 * COUNT( (/ jperio == 2 .OR. jperio == 7 /) ) + 1 * COUNT( (/ jperio >= 4 .AND. jperio <= 6 /) ) 
    7676      DO_2D_00_00 
    77          z2d(ji,jj) = 0.1 * ( 0.5 - REAL( mig0(ji) + (mjg0(jj)-1) * igloi, wp ) / REAL( igloi * igloj, wp ) ) 
     77         z2d(ji,jj) = 0.1 * ( 0.5 - REAL( mig0_oldcmp(ji) + (mjg0_oldcmp(jj)-1) * igloi, wp ) / REAL( igloi * igloj, wp ) ) 
    7878      END_2D 
    7979      ! 
  • NEMO/branches/2020/dev_r12558_HPC-08_epico_Extra_Halo/tests/BENCH/MY_SRC/usrdef_sbc.F90

    r12939 r12960  
    110110      ! Once we decide to forget trunk compatibility, we must simply define z2d as: 
    111111!!$      DO_2D_00_00 
    112 !!$         z2d(ji,jj) = 0.1 * ( 0.5 - REAL( mig0(ji) + (mjg0(jj)-1) * Ni0glo, wp ) / REAL( Ni0glo * Nj0glo, wp ) ) 
     112!!$         z2d(ji,jj) = 0.1 * ( 0.5 - REAL( mig0_oldcmp(ji) + (mjg0_oldcmp(jj)-1) * Ni0glo, wp ) / REAL( Ni0glo * Nj0glo, wp ) ) 
    113113!!$      END_2D 
    114114      igloi = Ni0glo + 2 * COUNT( (/ jperio == 1 .OR. jperio == 4 .OR. jperio == 6 .OR. jperio == 7 /) ) 
    115115      igloj = Nj0glo + 2 * COUNT( (/ jperio == 2 .OR. jperio == 7 /) ) + 1 * COUNT( (/ jperio >= 4 .AND. jperio <= 6 /) ) 
    116116      DO_2D_00_00 
    117          z2d(ji,jj) = 0.1 * ( 0.5 - REAL( mig0(ji) + (mjg0(jj)-1) * igloi, wp ) / REAL( igloi * igloj, wp ) ) 
     117         z2d(ji,jj) = 0.1 * ( 0.5 - REAL( mig0_oldcmp(ji) + (mjg0_oldcmp(jj)-1) * igloi, wp ) / REAL( igloi * igloj, wp ) ) 
    118118      END_2D 
    119119      utau_ice(:,:) = 0.1_wp + z2d(:,:) 
  • NEMO/branches/2020/dev_r12558_HPC-08_epico_Extra_Halo/tests/CANAL/MY_SRC/domvvl.F90

    r12939 r12960  
    450450            ELSE 
    451451               ijk_max = MAXLOC( ze3t(:,:,:) ) 
    452                ijk_max(1) = mig0(ijk_max(1)) 
    453                ijk_max(2) = mjg0(ijk_max(2)) 
     452               ijk_max(1) = mig0_oldcmp(ijk_max(1)) 
     453               ijk_max(2) = mjg0_oldcmp(ijk_max(2)) 
    454454               ijk_min = MINLOC( ze3t(:,:,:) ) 
    455                ijk_min(1) = mig0(ijk_min(1)) 
    456                ijk_min(2) = mjg0(ijk_min(2)) 
     455               ijk_min(1) = mig0_oldcmp(ijk_min(1)) 
     456               ijk_min(2) = mjg0_oldcmp(ijk_min(2)) 
    457457            ENDIF 
    458458            IF (lwp) THEN 
  • NEMO/branches/2020/dev_r12558_HPC-08_epico_Extra_Halo/tests/CANAL/MY_SRC/usrdef_hgr.F90

    r12939 r12960  
    9191          
    9292      DO_2D_11_11          
    93          zti = REAL( mig0(ji) - 1, wp )   ! start at i=0 in the global grid without halos 
    94          ztj = REAL( mjg0(jj) - 1, wp )   ! start at j=0 in the global grid without halos 
     93         zti = REAL( mig0_oldcmp(ji) - 1, wp )   ! start at i=0 in the global grid without halos 
     94         ztj = REAL( mjg0_oldcmp(jj) - 1, wp )   ! start at j=0 in the global grid without halos 
    9595          
    9696         plamt(ji,jj) = zlam0 + rn_dx *   zti 
  • NEMO/branches/2020/dev_r12558_HPC-08_epico_Extra_Halo/tests/ICE_ADV1D/MY_SRC/usrdef_hgr.F90

    r12939 r12960  
    7979 
    8080      DO_2D_11_11 
    81          zti = REAL( mig0(ji) - 1, wp )   ! start at i=0 in the global grid without halos 
    82          ztj = REAL( mjg0(jj) - 1, wp )   ! start at j=0 in the global grid without halos 
     81         zti = REAL( mig0_oldcmp(ji) - 1, wp )   ! start at i=0 in the global grid without halos 
     82         ztj = REAL( mjg0_oldcmp(jj) - 1, wp )   ! start at j=0 in the global grid without halos 
    8383          
    8484         plamt(ji,jj) = zlam0 + rn_dx * 1.e-3 *   zti 
  • NEMO/branches/2020/dev_r12558_HPC-08_epico_Extra_Halo/tests/ICE_ADV2D/MY_SRC/usrdef_hgr.F90

    r12939 r12960  
    9191 
    9292      DO_2D_11_11 
    93          zti = REAL( mig0(ji) - 1, wp )   ! start at i=0 in the global grid without halos 
    94          ztj = REAL( mjg0(jj) - 1, wp )   ! start at j=0 in the global grid without halos 
     93         zti = REAL( mig0_oldcmp(ji) - 1, wp )   ! start at i=0 in the global grid without halos 
     94         ztj = REAL( mjg0_oldcmp(jj) - 1, wp )   ! start at j=0 in the global grid without halos 
    9595          
    9696         plamt(ji,jj) = zlam0 + rn_dx * 1.e-3 *   zti 
  • NEMO/branches/2020/dev_r12558_HPC-08_epico_Extra_Halo/tests/ICE_AGRIF/MY_SRC/usrdef_hgr.F90

    r12939 r12960  
    9191 
    9292      DO_2D_11_11 
    93          zti = REAL( mig0(ji) - 1, wp )   ! start at i=0 in the global grid without halos 
    94          ztj = REAL( mjg0(jj) - 1, wp )   ! start at j=0 in the global grid without halos 
     93         zti = REAL( mig0_oldcmp(ji) - 1, wp )   ! start at i=0 in the global grid without halos 
     94         ztj = REAL( mjg0_oldcmp(jj) - 1, wp )   ! start at j=0 in the global grid without halos 
    9595          
    9696         plamt(ji,jj) = zlam0 + rn_dx * 1.e-3 *   zti 
  • NEMO/branches/2020/dev_r12558_HPC-08_epico_Extra_Halo/tests/ISOMIP/MY_SRC/usrdef_hgr.F90

    r12939 r12960  
    7979      DO_2D_11_11 
    8080         !                       ! longitude   (west coast at lon=0°) 
    81          plamt(ji,jj) = rn_e1deg * (  - 0.5 + REAL( mig0(ji)-1 , wp )  )   
    82          plamu(ji,jj) = rn_e1deg * (          REAL( mig0(ji)-1 , wp )  ) 
     81         plamt(ji,jj) = rn_e1deg * (  - 0.5 + REAL( mig0_oldcmp(ji)-1 , wp )  )   
     82         plamu(ji,jj) = rn_e1deg * (          REAL( mig0_oldcmp(ji)-1 , wp )  ) 
    8383         plamv(ji,jj) = plamt(ji,jj) 
    8484         plamf(ji,jj) = plamu(ji,jj) 
    8585         !                       ! latitude   (south coast at lat= 81°) 
    86          pphit(ji,jj) = rn_e2deg * (  - 0.5 + REAL( mjg0(jj)-1 , wp )  ) - 80._wp 
     86         pphit(ji,jj) = rn_e2deg * (  - 0.5 + REAL( mjg0_oldcmp(jj)-1 , wp )  ) - 80._wp 
    8787         pphiu(ji,jj) = pphit(ji,jj) 
    88          pphiv(ji,jj) = rn_e2deg * (          REAL( mjg0(jj)-1 , wp )  ) - 80_wp 
     88         pphiv(ji,jj) = rn_e2deg * (          REAL( mjg0_oldcmp(jj)-1 , wp )  ) - 80_wp 
    8989         pphif(ji,jj) = pphiv(ji,jj) 
    9090      END_2D 
  • NEMO/branches/2020/dev_r12558_HPC-08_epico_Extra_Halo/tests/LOCK_EXCHANGE/MY_SRC/usrdef_hgr.F90

    r12939 r12960  
    7676      DO_2D_11_11 
    7777         !                       ! longitude 
    78          plamt(ji,jj) = zfact * (  - 0.5 + REAL( mig0(ji)-1 , wp )  )   
    79          plamu(ji,jj) = zfact * (          REAL( mig0(ji)-1 , wp )  ) 
     78         plamt(ji,jj) = zfact * (  - 0.5 + REAL( mig0_oldcmp(ji)-1 , wp )  )   
     79         plamu(ji,jj) = zfact * (          REAL( mig0_oldcmp(ji)-1 , wp )  ) 
    8080         plamv(ji,jj) = plamt(ji,jj) 
    8181         plamf(ji,jj) = plamu(ji,jj) 
    8282         !                       ! latitude 
    83          pphit(ji,jj) = zfact * (  - 0.5 + REAL( mjg0(jj)-1 , wp )  ) 
     83         pphit(ji,jj) = zfact * (  - 0.5 + REAL( mjg0_oldcmp(jj)-1 , wp )  ) 
    8484         pphiu(ji,jj) = pphit(ji,jj) 
    85          pphiv(ji,jj) = zfact * (          REAL( mjg0(jj)-1 , wp )  ) 
     85         pphiv(ji,jj) = zfact * (          REAL( mjg0_oldcmp(jj)-1 , wp )  ) 
    8686         pphif(ji,jj) = pphiv(ji,jj) 
    8787      END_2D 
  • NEMO/branches/2020/dev_r12558_HPC-08_epico_Extra_Halo/tests/OVERFLOW/MY_SRC/usrdef_hgr.F90

    r12939 r12960  
    7676      DO_2D_11_11 
    7777         !                       ! longitude 
    78          plamt(ji,jj) = zfact * (  - 0.5 + REAL( mig0(ji)-1 , wp )  )   
    79          plamu(ji,jj) = zfact * (          REAL( mig0(ji)-1 , wp )  ) 
     78         plamt(ji,jj) = zfact * (  - 0.5 + REAL( mig0_oldcmp(ji)-1 , wp )  )   
     79         plamu(ji,jj) = zfact * (          REAL( mig0_oldcmp(ji)-1 , wp )  ) 
    8080         plamv(ji,jj) = plamt(ji,jj) 
    8181         plamf(ji,jj) = plamu(ji,jj) 
    8282         !                       ! latitude 
    83          pphit(ji,jj) = zfact * (  - 0.5 + REAL( mjg0(jj)-1 , wp )  ) 
     83         pphit(ji,jj) = zfact * (  - 0.5 + REAL( mjg0_oldcmp(jj)-1 , wp )  ) 
    8484         pphiu(ji,jj) = pphit(ji,jj) 
    85          pphiv(ji,jj) = zfact * (          REAL( mjg0(jj)-1 , wp )  ) 
     85         pphiv(ji,jj) = zfact * (          REAL( mjg0_oldcmp(jj)-1 , wp )  ) 
    8686         pphif(ji,jj) = pphiv(ji,jj) 
    8787      END_2D 
  • NEMO/branches/2020/dev_r12558_HPC-08_epico_Extra_Halo/tests/STATION_ASF/MY_SRC/nemogcm.F90

    r12939 r12960  
    231231      ! 
    232232      IF( ln_read_cfg ) THEN            ! Read sizes in domain configuration file 
    233          CALL domain_cfg ( cn_cfg, nn_cfg, jpiglo, jpjglo, jpkglo, jperio ) 
     233         CALL domain_cfg ( cn_cfg, nn_cfg, Ni0glo, Nj0glo, jpkglo, jperio ) 
    234234      ELSE                              ! user-defined namelist 
    235          CALL usr_def_nam( cn_cfg, nn_cfg, jpiglo, jpjglo, jpkglo, jperio ) 
     235         CALL usr_def_nam( cn_cfg, nn_cfg, Ni0glo, Nj0glo, jpkglo, jperio ) 
    236236      ENDIF 
    237237      ! 
  • NEMO/branches/2020/dev_r12558_HPC-08_epico_Extra_Halo/tests/VORTEX/MY_SRC/domvvl.F90

    r12939 r12960  
    450450            ELSE 
    451451               ijk_max = MAXLOC( ze3t(:,:,:) ) 
    452                ijk_max(1) = mig0(ijk_max(1)) 
    453                ijk_max(2) = mjg0(ijk_max(2)) 
     452               ijk_max(1) = mig0_oldcmp(ijk_max(1)) 
     453               ijk_max(2) = mjg0_oldcmp(ijk_max(2)) 
    454454               ijk_min = MINLOC( ze3t(:,:,:) ) 
    455                ijk_min(1) = mig0(ijk_min(1)) 
    456                ijk_min(2) = mjg0(ijk_min(2)) 
     455               ijk_min(1) = mig0_oldcmp(ijk_min(1)) 
     456               ijk_min(2) = mjg0_oldcmp(ijk_min(2)) 
    457457            ENDIF 
    458458            IF (lwp) THEN 
  • NEMO/branches/2020/dev_r12558_HPC-08_epico_Extra_Halo/tests/VORTEX/MY_SRC/usrdef_hgr.F90

    r12939 r12960  
    9191          
    9292      DO_2D_11_11 
    93          zti = REAL( mig0(ji) - 1, wp )   ! start at i=0 in the global grid without halos 
    94          ztj = REAL( mjg0(jj) - 1, wp )   ! start at j=0 in the global grid without halos 
     93         zti = REAL( mig0_oldcmp(ji) - 1, wp )   ! start at i=0 in the global grid without halos 
     94         ztj = REAL( mjg0_oldcmp(jj) - 1, wp )   ! start at j=0 in the global grid without halos 
    9595          
    9696         plamt(ji,jj) = zlam0 + rn_dx * 1.e-3 *   zti 
  • NEMO/branches/2020/dev_r12558_HPC-08_epico_Extra_Halo/tests/WAD/MY_SRC/usrdef_hgr.F90

    r12939 r12960  
    7676      DO_2D_11_11 
    7777         !                       ! longitude 
    78          plamt(ji,jj) = zfact * (  - 0.5 + REAL( mig0(ji)-1 , wp )  )   
    79          plamu(ji,jj) = zfact * (          REAL( mig0(ji)-1 , wp )  ) 
     78         plamt(ji,jj) = zfact * (  - 0.5 + REAL( mig0_oldcmp(ji)-1 , wp )  )   
     79         plamu(ji,jj) = zfact * (          REAL( mig0_oldcmp(ji)-1 , wp )  ) 
    8080         plamv(ji,jj) = plamt(ji,jj) 
    8181         plamf(ji,jj) = plamu(ji,jj) 
    8282         !                       ! latitude 
    83          pphit(ji,jj) = zfact * (  - 0.5 + REAL( mjg0(jj)-1 , wp )  ) 
     83         pphit(ji,jj) = zfact * (  - 0.5 + REAL( mjg0_oldcmp(jj)-1 , wp )  ) 
    8484         pphiu(ji,jj) = pphit(ji,jj) 
    85          pphiv(ji,jj) = zfact * (          REAL( mjg0(jj)-1 , wp )  ) 
     85         pphiv(ji,jj) = zfact * (          REAL( mjg0_oldcmp(jj)-1 , wp )  ) 
    8686         pphif(ji,jj) = pphiv(ji,jj) 
    8787      END_2D 
Note: See TracChangeset for help on using the changeset viewer.