Changeset 867 for codes/icosagcm/devel/src/time

Timestamp:

05/16/19 23:45:56 (5 years ago)

Author:

dubos

Message:

devel : towards Fortran driver for unstructured/LAM meshes

File:

• codes/icosagcm/devel/src/time/euler_scheme.f90 (modified) (1 diff)

codes/icosagcm/devel/src/time/euler_scheme.f90

@@ -167 +167 @@
 
        u=f_u(ind)
-       DO l=ll_begin,ll_end
-          !DIR$ SIMD
-          DO ij=ij_begin,ij_end
-             u(ij+u_right,l)=u(ij+u_right,l)*de(ij+u_right)
-             u(ij+u_lup,l)=u(ij+u_lup,l)*de(ij+u_lup)
-             u(ij+u_ldown,l)=u(ij+u_ldown,l)*de(ij+u_ldown)
-          ENDDO
-       ENDDO
+       SELECT CASE(grid_type)
+       CASE(grid_ico)
+          DO l=ll_begin,ll_end
+             !DIR$ SIMD
+             DO ij=ij_begin,ij_end
+                u(ij+u_right,l)=u(ij+u_right,l)*de(ij+u_right)
+                u(ij+u_lup,l)=u(ij+u_lup,l)*de(ij+u_lup)
+                u(ij+u_ldown,l)=u(ij+u_ldown,l)*de(ij+u_ldown)
+             ENDDO
+          ENDDO
+       CASE(grid_unst)
+          PRINT *, 'legacy_to_DEC', llm, edge_num, SHAPE(u), SHAPE(de) ! FIXME
+          DO ij=1, edge_num
+             DO l=1,llm
+                u(l,ij) = u(l,ij)*de(ij)
+             END DO
+          END DO
+       CASE DEFAULT
+          STOP 'Unsupported grid_type encountered in legacy_to_DEC'
+       END SELECT
     ENDDO