Changeset 688 for codes/icosagcm/devel/src/unstructured/transport_unstructured.F90

Timestamp:

04/09/18 15:22:10 (6 years ago)

Author:

dubos

Message:

devel/unstructured : select double or single precision for physical quantities

File:

• codes/icosagcm/devel/src/unstructured/transport_unstructured.F90 (modified) (2 diffs)

codes/icosagcm/devel/src/unstructured/transport_unstructured.F90

@@ -7 +7 @@
   SAVE
 
-#define BINDC_(thename) BIND(C, name=#thename)
-#define BINDC(thename) BINDC_(dynamico_ ## thename)
-
-#define DBL REAL(C_DOUBLE)
-#define DOUBLE1(m) DBL, DIMENSION(m)
-#define DOUBLE2(m,n) DBL, DIMENSION(m,n)
-#define DOUBLE3(m,n,p) DBL, DIMENSION(m,n,p)
-#define DOUBLE4(m,n,p,q) DBL, DIMENSION(m,n,p,q)
-#define INDEX INTEGER(C_INT)
-
 CONTAINS
 
-#define FIELD_PS      DOUBLE1(primal_num)
-#define FIELD_MASS    DOUBLE2(llm, primal_num)
-#define FIELD_Z       DOUBLE2(llm, dual_num)
-#define FIELD_U       DOUBLE2(llm, edge_num)
-#define FIELD_UL      DOUBLE2(llm+1, edge_num)
-#define FIELD_THETA   DOUBLE3(llm, primal_num, nqdyn)
-#define FIELD_GEOPOT  DOUBLE2(llm+1, primal_num)
-
+#include "unstructured.h90"
 #define HASNAN(field) (ANY(.NOT.ABS(field)<1e20))
 
@@ -40 +23 @@
     INTEGER, DIMENSION(primal_num) :: cur_lev
 
-    DBL          :: rhodz_cum_target, rhodz_cum_lev, rhodz_cum_levp1, X, time1,time2
+    NUM          :: rhodz_cum_target, rhodz_cum_lev, rhodz_cum_levp1, X
+    TIME         :: time1,time2
     INTEGER :: ij, l, level, edge, ij_left, ij_right