Changeset 2228 for LMDZ5/trunk/libf/dyn3d_common/fyhyp_m.F90

Timestamp:

Mar 12, 2015, 2:15:18 PM (10 years ago)

Author:

lguez

Message:

Correcting a problem from revision 2218. The type double precision
with option "-fdefault-real-8" of gfortran is promoted to 16-byte
precision and there is no specific procedure in arth with this
precision. Could not add a specific procedure in arth with double
precision because, with ifort, the option "-real-size 64" does not
promote the double precision, so that would make two identical
specific procedures in arth.

In module nrtype, replaced double precision by a parameterized real
kind so that the effective precision does not depend on a compiler
option.

In coefpoly, fxhyp, fyhyp and invert_zoom_x, use the parameterized
real kind defined in nrtype, instead of double precision.

Also, in module nrtype, removed unused derived types sprs2_sp and
sprs2_dp.

File:

• LMDZ5/trunk/libf/dyn3d_common/fyhyp_m.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/dyn3d_common/fyhyp_m.F90

@@ -17 +17 @@
 
     use coefpoly_m, only: coefpoly
+    use nrtype, only: k8
 
     include "dimensions.h"
@@ -30 +31 @@
     ! Local: 
 
-    DOUBLE PRECISION champmin, champmax
+    REAL(K8) champmin, champmax
     INTEGER, PARAMETER:: nmax=30000, nmax2=2*nmax
     REAL dzoom ! distance totale de la zone du zoom (en radians)
-    DOUBLE PRECISION ylat(jjm + 1), yprim(jjm + 1)
-    DOUBLE PRECISION yuv
-    DOUBLE PRECISION, save:: yt(0:nmax2)
-    DOUBLE PRECISION fhyp(0:nmax2), beta
-    DOUBLE PRECISION, save:: ytprim(0:nmax2)
-    DOUBLE PRECISION fxm(0:nmax2)
-    DOUBLE PRECISION, save:: yf(0:nmax2)
-    DOUBLE PRECISION yypr(0:nmax2)
-    DOUBLE PRECISION yvrai(jjm + 1), yprimm(jjm + 1), ylatt(jjm + 1)
-    DOUBLE PRECISION pi, pis2, epsilon, y0, pisjm
-    DOUBLE PRECISION yo1, yi, ylon2, ymoy, yprimin
-    DOUBLE PRECISION yfi, yf1, ffdy
-    DOUBLE PRECISION ypn, deply, y00
+    REAL(K8) ylat(jjm + 1), yprim(jjm + 1)
+    REAL(K8) yuv
+    REAL(K8), save:: yt(0:nmax2)
+    REAL(K8) fhyp(0:nmax2), beta
+    REAL(K8), save:: ytprim(0:nmax2)
+    REAL(K8) fxm(0:nmax2)
+    REAL(K8), save:: yf(0:nmax2)
+    REAL(K8) yypr(0:nmax2)
+    REAL(K8) yvrai(jjm + 1), yprimm(jjm + 1), ylatt(jjm + 1)
+    REAL(K8) pi, pis2, epsilon, y0, pisjm
+    REAL(K8) yo1, yi, ylon2, ymoy, yprimin
+    REAL(K8) yfi, yf1, ffdy
+    REAL(K8) ypn, deply, y00
     SAVE y00, deply
 
@@ -51 +52 @@
     INTEGER jpn, jjpn
     SAVE jpn
-    DOUBLE PRECISION a0, a1, a2, a3, yi2, heavyy0, heavyy0m
-    DOUBLE PRECISION fa(0:nmax2), fb(0:nmax2)
+    REAL(K8) a0, a1, a2, a3, yi2, heavyy0, heavyy0m
+    REAL(K8) fa(0:nmax2), fb(0:nmax2)
     REAL y0min, y0max
 
-    DOUBLE PRECISION heavyside
+    REAL(K8) heavyside
 
     !-------------------------------------------------------------------