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.
Changeset 14673 – NEMO

Changeset 14673


Ignore:
Timestamp:
2021-04-02T15:06:01+02:00 (3 years ago)
Author:
jchanut
Message:

Update ICE_AGRIF child domain as for VORTEX (position dependent on ghost cell number), #2638

Location:
NEMO/branches/2021/dev_r14608_AGRIF_domcfg/tests/ICE_AGRIF
Files:
2 deleted
3 edited

Legend:

Unmodified
Added
Removed

  • NEMO/branches/2021/dev_r14608_AGRIF_domcfg/tests/ICE_AGRIF/EXPREF/AGRIF_FixedGrids.in

    r13286 r14673  
    111 
    2 33 62 33 62 3 3 3 
     234 63 34 63 3 3 3 
    330 

  • NEMO/branches/2021/dev_r14608_AGRIF_domcfg/tests/ICE_AGRIF/EXPREF/make_INITICE.py

    r10516 r14673  
    2020# Reading coordinates file 
    2121nccoord=netcdf(fcoord,'r') 
    22 nav_lon=nccoord.variables['nav_lon'] 
    23 nav_lat=nccoord.variables['nav_lat'] 
     22nav_lon=nccoord.variables['x'] 
     23nav_lat=nccoord.variables['y'] 
    2424time_counter=1 
    2525LON1= nav_lon.shape[1] 

  • NEMO/branches/2021/dev_r14608_AGRIF_domcfg/tests/ICE_AGRIF/MY_SRC/usrdef_hgr.F90

    r14223 r14673  
    7979      IF( Agrif_Root() ) THEN 
    8080#endif 
    81          zlam0 = -REAL(Ni0glo, wp) * 0.5 * 1.e-3 * rn_dx 
    82          zphi0 = -REAL(Nj0glo, wp) * 0.5 * 1.e-3 * rn_dy 
     81         zlam0 = -REAL(Ni0glo-1, wp) * 0.5 * 1.e-3 * rn_dx 
     82         zphi0 = -REAL(Nj0glo-1, wp) * 0.5 * 1.e-3 * rn_dy 
    8383#if defined key_agrif  
    8484      ELSE 
    … …  
    8686!clem         zlam0  = Agrif_Parent(zlam0) + (Agrif_ix())*Agrif_Parent(rn_dx) * 1.e-5 
    8787!clem         zphi0  = Agrif_Parent(zphi0) + (Agrif_iy())*Agrif_Parent(rn_dy) * 1.e-5 
    88          zlam0 = ( 0.5_wp - REAL(Ni0glo, wp) * 0.5 ) * 1.e-3 * Agrif_irhox() * rn_dx  & 
    89             &  + ( Agrif_Ix() + nbghostcells - 1 ) * Agrif_irhox() * rn_dx * 1.e-3 - ( 0.5_wp + nbghostcells ) * rn_dx * 1.e-3 
    90          zphi0 = ( 0.5_wp - REAL(Nj0glo, wp) * 0.5 ) * 1.e-3 * Agrif_irhoy() * rn_dy  & 
    91             &  + ( Agrif_Iy() + nbghostcells - 1 ) * Agrif_irhoy() * rn_dy * 1.e-3 - ( 0.5_wp + nbghostcells ) * rn_dy * 1.e-3 
     88         zlam0 = ( 0.5_wp - REAL(Agrif_Parent(Ni0glo)-1, wp) * 0.5 ) * 1.e-3 * Agrif_irhox() * rn_dx  & 
     89            &  + ( Agrif_Ix() + Agrif_Parent(nbghostcells_x  ) - 1 ) * Agrif_irhox() * rn_dx * 1.e-3 - ( 0.5_wp + nbghostcells_x  ) * rn_dx * 1.e-3 
     90         zphi0 = ( 0.5_wp - REAL(Agrif_Parent(Nj0glo)-1, wp) * 0.5 ) * 1.e-3 * Agrif_irhoy() * rn_dy  & 
     91            &  + ( Agrif_Iy() + Agrif_Parent(nbghostcells_y_s) - 1 ) * Agrif_irhoy() * rn_dy * 1.e-3 - ( 0.5_wp + nbghostcells_y_s ) * rn_dy * 1.e-3 
    9292      ENDIF 
    9393#endif          
    9494 
    9595      DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 
    96          zti = REAL( mig0(ji), wp ) - 0.5_wp  ! start at i=0.5 in the global grid without halos 
    97          ztj = REAL( mjg0(jj), wp ) - 0.5_wp  ! start at j=0.5 in the global grid without halos 
     96         zti = REAL( mig0(ji)-1, wp ) ! start at i=0. in the global grid without halos 
     97         ztj = REAL( mjg0(jj)-1, wp ) ! start at j=0. in the global grid without halos 
    9898           
    9999         plamt(ji,jj) = zlam0 + rn_dx * 1.e-3 *   zti 
Note: See TracChangeset for help on using the changeset viewer.