- Timestamp:
- 2021-04-06T14:36:08+02:00 (3 years ago)
- Location:
- utils/tools/DOMAINcfg
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
utils/tools/DOMAINcfg/AGRIF_FixedGrids.in_cmems
r14633 r14674 2 2 -10 15 57 75 3 3 3 3 3 30 70 126 170 3 3 3 4 1 1 791 35 3 3 34 1 181 1 35 3 3 3 5 5 0 6 6 0 -
utils/tools/DOMAINcfg/make_namelist.py
r14626 r14674 64 64 if int(grid[cnt-1][2]) == 1: 65 65 nbghostcells_y_s = 0 66 if int(grid[cnt-1][1]) + int(grid[cnt-1][0]) == int(Ni0glo_parent) :66 if int(grid[cnt-1][1]) + int(grid[cnt-1][0]) == int(Ni0glo_parent) + 2 : 67 67 nbghostcells_x = 0 68 69 Ni0glo = (int(grid[cnt-1][1])-int(grid[cnt-1][0]))*int(grid[cnt-1][4]) + 2 + 2*nbghostcells_x 70 Nj0glo = (int(grid[cnt-1][3])-int(grid[cnt-1][2]))*int(grid[cnt-1][5]) + 2 + nbghostcells_y_n + nbghostcells_y_s 68 Ni0glo = (int(grid[cnt-1][1])-int(grid[cnt-1][0]))*int(grid[cnt-1][4]) + 2*nbghostcells_x 69 Nj0glo = (int(grid[cnt-1][3])-int(grid[cnt-1][2]))*int(grid[cnt-1][5]) + nbghostcells_y_n + nbghostcells_y_s 71 70 #print( "Grid "+str(cnt)+" : jpiglo = "+cnt(jpiglo)+ " jpjglo = "+str(jpjglo) ) 72 71 print('Grid {:1d} : Ni0glo = {:3d} , Nj0glo = {:3d}'.format(cnt, Ni0glo, Nj0glo)) … … 79 78 line = fp.readline() 80 79 if line.strip().startswith('jperio'): 81 if int(grid[cnt-1][1]) + int(grid[cnt-1][0]) == int(Ni0glo_parent) :80 if int(grid[cnt-1][1]) + int(grid[cnt-1][0]) == int(Ni0glo_parent) + 2: 82 81 line = " jperio = 1\n" 83 82 else: -
utils/tools/DOMAINcfg/src/agrif_user.F90
r14640 r14674 215 215 !--------------------------------------------------------------------- 216 216 217 ind2 = nn_hls + 2+ nbghostcells_x218 ind3 = nn_hls + 2+ nbghostcells_y_s217 ind2 = nn_hls + 1 + nbghostcells_x 218 ind3 = nn_hls + 1 + nbghostcells_y_s 219 219 220 220 nbghostcellsfine_tot_x=nbghostcells_x+1 … … 1115 1115 IF (.not.agrif_root()) THEN 1116 1116 lk_west = .NOT. ( Agrif_Ix() == 1 ) 1117 lk_east = .NOT. ( Agrif_Ix() + nbcellsx/AGRIF_Irhox() == Agrif_Parent(Ni0glo) -1 )1117 lk_east = .NOT. ( Agrif_Ix() + nbcellsx/AGRIF_Irhox() == Agrif_Parent(Ni0glo) + 1 ) 1118 1118 lk_south = .NOT. ( Agrif_Iy() == 1 ) 1119 lk_north = .NOT. ( Agrif_Iy() + nbcellsy/AGRIF_Irhoy() == Agrif_Parent(Nj0glo) -1)1119 lk_north = .NOT. ( Agrif_Iy() + nbcellsy/AGRIF_Irhoy() == Agrif_Parent(Nj0glo) + 1) 1120 1120 IF (.NOT.lk_south) THEN 1121 1121 nbghostcells_y_s = 0 -
utils/tools/DOMAINcfg/src/mppini.F90
r14623 r14674 109 109 print *,'nbcellsx = ',nbcellsx,nbghostcells_x 110 110 print *,'nbcellsy = ',nbcellsy,nbghostcells_y_s,nbghostcells_y_n 111 IF( Ni0glo /= nbcellsx + 2 + 2*nbghostcells_x ) THEN111 IF( Ni0glo /= nbcellsx + 2*nbghostcells_x ) THEN 112 112 IF(lwp) THEN 113 113 WRITE(numout,*) 114 WRITE(numout,*) 'Ni0glo should be: ', nbcellsx + 2 + 2*nbghostcells_x114 WRITE(numout,*) 'Ni0glo should be: ', nbcellsx + 2*nbghostcells_x 115 115 ENDIF 116 CALL ctl_stop( 'STOP', 'mpp_init: Agrif children requires Ni0glo == nbcellsx + 2 + 2*nbghostcells_x' )116 CALL ctl_stop( 'STOP', 'mpp_init: Agrif children requires Ni0glo == nbcellsx + 2*nbghostcells_x' ) 117 117 ENDIF 118 IF( Nj0glo /= nbcellsy + 2 +nbghostcells_y_s + nbghostcells_y_n ) THEN118 IF( Nj0glo /= nbcellsy + nbghostcells_y_s + nbghostcells_y_n ) THEN 119 119 IF(lwp) THEN 120 120 WRITE(numout,*) 121 WRITE(numout,*) 'Nj0glo shoud be: ', nbcellsy + 2 +nbghostcells_y_s + nbghostcells_y_n121 WRITE(numout,*) 'Nj0glo shoud be: ', nbcellsy + nbghostcells_y_s + nbghostcells_y_n 122 122 ENDIF 123 123 CALL ctl_stop( 'STOP', & 124 'mpp_init: Agrif children requires Nj0glo == nbcellsy + 2 +nbghostcells_y_s + nbghostcells_y_n' )124 'mpp_init: Agrif children requires Nj0glo == nbcellsy + nbghostcells_y_s + nbghostcells_y_n' ) 125 125 ENDIF 126 126 IF( ln_use_jattr ) CALL ctl_stop( 'STOP', 'mpp_init:Agrif children requires ln_use_jattr = .false. ' ) … … 361 361 CALL agrif_nemo_init() 362 362 IF( .NOT. Agrif_Root() ) THEN ! AGRIF children: specific setting (cf. agrif_user.F90) 363 IF( Ni0glo /= nbcellsx + 2 + 2*nbghostcells_x ) THEN363 IF( Ni0glo /= nbcellsx + 2*nbghostcells_x ) THEN 364 364 IF(lwp) THEN 365 365 WRITE(numout,*) 366 WRITE(numout,*) 'Ni0glo should be: ', nbcellsx + 2 + 2*nbghostcells_x366 WRITE(numout,*) 'Ni0glo should be: ', nbcellsx + 2*nbghostcells_x 367 367 ENDIF 368 CALL ctl_stop( 'STOP', 'mpp_init: Agrif children requires Ni0glo == nbcellsx + 2 + 2*nbghostcells_x' )368 CALL ctl_stop( 'STOP', 'mpp_init: Agrif children requires Ni0glo == nbcellsx + 2*nbghostcells_x' ) 369 369 ENDIF 370 IF( Nj0glo /= nbcellsy + 2+ nbghostcells_y_s + nbghostcells_y_n ) THEN370 IF( Nj0glo /= nbcellsy + nbghostcells_y_s + nbghostcells_y_n ) THEN 371 371 IF(lwp) THEN 372 372 WRITE(numout,*) 373 WRITE(numout,*) 'Nj0glo shoud be: ', nbcellsy + 2+ nbghostcells_y_s + nbghostcells_y_n373 WRITE(numout,*) 'Nj0glo shoud be: ', nbcellsy + nbghostcells_y_s + nbghostcells_y_n 374 374 ENDIF 375 375 CALL ctl_stop( 'STOP', & 376 'mpp_init: Agrif children requires Nj0glo == nbcellsy + 2+ nbghostcells_y_s + nbghostcells_y_n' )376 'mpp_init: Agrif children requires Nj0glo == nbcellsy + nbghostcells_y_s + nbghostcells_y_n' ) 377 377 ENDIF 378 378 IF( ln_use_jattr ) CALL ctl_stop( 'STOP', 'mpp_init:Agrif children requires ln_use_jattr = .false. ' )
Note: See TracChangeset
for help on using the changeset viewer.