Changeset 12766 for NEMO/branches/UKMO/dev_r12745_HPC-02_Daley_Tiling_trial_structure/src/OCE/do_loop_substitute.h90

Timestamp:

2020-04-17T14:54:46+02:00 (4 years ago)

Author:

hadcv

Message:

tra_ldf_iso trial using structures

File:

• NEMO/branches/UKMO/dev_r12745_HPC-02_Daley_Tiling_trial_structure/src/OCE/do_loop_substitute.h90 (modified) (1 diff)

NEMO/branches/UKMO/dev_r12745_HPC-02_Daley_Tiling_trial_structure/src/OCE/do_loop_substitute.h90

@@ -102 +102 @@
 #define END_2D   END DO   ;   END DO
 #define END_3D   END DO   ;   END DO   ;   END DO
+
+#define __kIs_T     ktile % ntsi
+#define __kJs_T     ktile % ntsj
+#define __kIsm1_T   ktile % ntsim1
+#define __kJsm1_T   ktile % ntsjm1
+
+#define __kIe_T     ktile % ntei
+#define __kJe_T     ktile % ntej
+#define __kIep1_T   ktile % nteip1
+#define __kJep1_T   ktile % ntejp1
+
+#define IND_2D          __kIsm1_T:__kIep1_T,__kJsm1_T:__kJep1_T
+
+#define DO_2D_00_00_T   DO jj = __kJs_T, __kJe_T   ;   DO ji = __kIs_T, __kIe_T
+#define DO_2D_00_01_T   DO jj = __kJs_T, __kJe_T   ;   DO ji = __kIs_T, __kIep1_T
+#define DO_2D_00_10_T   DO jj = __kJs_T, __kJe_T   ;   DO ji = __kIsm1_T, __kIe_T
+#define DO_2D_00_11_T   DO jj = __kJs_T, __kJe_T   ;   DO ji = __kIsm1_T, __kIep1_T
+
+#define DO_2D_01_00_T   DO jj = __kJs_T, __kJep1_T   ;   DO ji = __kIs_T, __kIe_T
+#define DO_2D_01_01_T   DO jj = __kJs_T, __kJep1_T   ;   DO ji = __kIs_T, __kIep1_T
+#define DO_2D_01_10_T   DO jj = __kJs_T, __kJep1_T   ;   DO ji = __kIsm1_T, __kIe_T
+#define DO_2D_01_11_T   DO jj = __kJs_T, __kJep1_T   ;   DO ji = __kIsm1_T, __kIep1_T
+
+#define DO_2D_10_00_T   DO jj = __kJsm1_T, __kJe_T   ;   DO ji = __kIs_T, __kIe_T
+#define DO_2D_10_10_T   DO jj = __kJsm1_T, __kJe_T   ;   DO ji = __kIsm1_T, __kIe_T
+#define DO_2D_10_11_T   DO jj = __kJsm1_T, __kJe_T   ;   DO ji = __kIsm1_T, __kIep1_T
+
+#define DO_2D_11_00_T   DO jj = __kJsm1_T, __kJep1_T   ;   DO ji = __kIs_T, __kIe_T
+#define DO_2D_11_01_T   DO jj = __kJsm1_T, __kJep1_T   ;   DO ji = __kIs_T, __kIep1_T
+#define DO_2D_11_10_T   DO jj = __kJsm1_T, __kJep1_T   ;   DO ji = __kIsm1_T, __kIe_T
+#define DO_2D_11_11_T   DO jj = __kJsm1_T, __kJep1_T   ;   DO ji = __kIsm1_T, __kIep1_T
+
+#define DO_3D_00_00_T(ks,ke)   DO jk = ks, ke   ;   DO_2D_00_00_T
+#define DO_3D_00_10_T(ks,ke)   DO jk = ks, ke   ;   DO_2D_00_10_T
+
+#define DO_3D_01_01_T(ks,ke)   DO jk = ks, ke   ;   DO_2D_01_01_T
+
+#define DO_3D_10_00_T(ks,ke)   DO jk = ks, ke   ;   DO_2D_10_00_T
+#define DO_3D_10_10_T(ks,ke)   DO jk = ks, ke   ;   DO_2D_10_10_T
+#define DO_3D_10_11_T(ks,ke)   DO jk = ks, ke   ;   DO_2D_10_11_T
+
+#define DO_3D_11_11_T(ks,ke)   DO jk = ks, ke   ;   DO_2D_11_11_T
+
+#define DO_3DS_00_00_T(ks,ke,ki)   DO jk = ks, ke, ki   ;   DO_2D_00_00_T
+#define DO_3DS_01_01_T(ks,ke,ki)   DO jk = ks, ke, ki   ;   DO_2D_01_01_T
+#define DO_3DS_10_10_T(ks,ke,ki)   DO jk = ks, ke, ki   ;   DO_2D_10_10_T
+#define DO_3DS_11_11_T(ks,ke,ki)   DO jk = ks, ke, ki   ;   DO_2D_11_11_T