Changeset 15279 for utils/tools/DOMAINcfg/make_namelist.py

Timestamp:

2021-09-23T12:00:23+02:00 (3 years ago)

Author:

jchanut

Message:

#2222 and #2638: Enable creating agrif meshes with different vertical grids (geopotential only as a start)

File:

• utils/tools/DOMAINcfg/make_namelist.py (modified) (2 diffs)

utils/tools/DOMAINcfg/make_namelist.py

@@ -62 +62 @@
     print(int(Ni0glo_parent), int(Nj0glo_parent))
 
-    nbghostcells_x = nbghostcells
-    nbghostcells_y = nbghostcells
-    nbghostcells_y_n = nbghostcells_y
-    nbghostcells_y_s = nbghostcells_y
+    nbghostcells_x_e = nbghostcells
+    nbghostcells_x_w = nbghostcells
+    nbghostcells_y_n = nbghostcells
+    nbghostcells_y_s = nbghostcells
     if (int(grid[cnt-1][2]) == 1 ):
         nbghostcells_y_s = 1 
     if int(grid[cnt-1][3]) == int(Nj0glo_parent)-1 :
         nbghostcells_y_n = 1 
+    if (int(grid[cnt-1][0]) == 1 ):
+        nbghostcells_x_w = 1 
+    if int(grid[cnt-1][1]) == int(Ni0glo_parent)-1 :
+        nbghostcells_x_e = 1 
+    if int(grid[cnt-1][3]) == int(Nj0glo_parent) :
+        nbghostcells_y_n = 0 
     if int(grid[cnt-1][1]) - int(grid[cnt-1][0]) == int(Ni0glo_parent) :
-        nbghostcells_x = 0
+        nbghostcells_x_w = 0
+        nbghostcells_x_e = 0
 
-    Ni0glo = (int(grid[cnt-1][1])-int(grid[cnt-1][0]))*int(grid[cnt-1][4]) + 2*nbghostcells_x
+    Ni0glo = (int(grid[cnt-1][1])-int(grid[cnt-1][0]))*int(grid[cnt-1][4]) + nbghostcells_x_w  + nbghostcells_x_e
     Nj0glo = (int(grid[cnt-1][3])-int(grid[cnt-1][2]))*int(grid[cnt-1][5]) + nbghostcells_y_n  + nbghostcells_y_s
     #print( "Grid "+str(cnt)+" : jpiglo = "+cnt(jpiglo)+ "  jpjglo = "+str(jpjglo) ) 
     print(int(grid[cnt-1][0]), int(grid[cnt-1][1]), int(grid[cnt-1][2]),int(grid[cnt-1][3]))
-    print(nbghostcells_x, nbghostcells_y_s, nbghostcells_y_n)
+    print(nbghostcells_x_w, nbghostcells_x_e,  nbghostcells_y_s, nbghostcells_y_n)
     print('Grid {:1d} : Ni0glo = {:3d} , Nj0glo = {:3d}'.format(cnt, Ni0glo, Nj0glo))
 
@@ -91 +98 @@
                else:
                    line = "   jperio = 0\n"
+               if nbghostcells_y_n == 0:
+                   line = "   jperio = 4\n"
            if line.strip().startswith('nn_bathy'):
                   line = "   nn_bathy = 2\n"