Changeset 2228 for LMDZ5/trunk/libf/dyn3d_common

Timestamp:

Mar 12, 2015, 2:15:18 PM (9 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.

Location:

LMDZ5/trunk/libf/dyn3d_common

Files:

• coefpoly_m.F90 (modified) (1 diff)
• fxhyp_m.F90 (modified) (3 diffs)
• fyhyp_m.F90 (modified) (3 diffs)
• invert_zoom_x_m.F90 (modified) (5 diffs)

LMDZ5/trunk/libf/dyn3d_common/coefpoly_m.F90

@@ -25 +25 @@
     ! a0, a1, a2, a3.
 
-    DOUBLE PRECISION, intent(in):: xf1, xf2, xprim1, xprim2, xtild1, xtild2
-    DOUBLE PRECISION, intent(out):: a0, a1, a2, a3
+    use nrtype, only: k8
+
+    REAL(K8), intent(in):: xf1, xf2, xprim1, xprim2, xtild1, xtild2
+    REAL(K8), intent(out):: a0, a1, a2, a3
 
     ! Local:
-    DOUBLE PRECISION xtil1car, xtil2car, derr, x1x2car
+    REAL(K8) xtil1car, xtil2car, derr, x1x2car
 
     !------------------------------------------------------------

LMDZ5/trunk/libf/dyn3d_common/fxhyp_m.F90

@@ -20 +20 @@
     use arth_m, only: arth
     use invert_zoom_x_m, only: invert_zoom_x, nmax
-    use nrtype, only: pi, pi_d, twopi, twopi_d
+    use nrtype, only: pi, pi_d, twopi, twopi_d, k8
     use principal_cshift_m, only: principal_cshift
 
@@ -36 +36 @@
     REAL dzoom, step
     real d_rlonv(iim)
-    DOUBLE PRECISION xtild(0:2 * nmax)
-    DOUBLE PRECISION fhyp(nmax:2 * nmax), ffdx, beta, Xprimt(0:2 * nmax)
-    DOUBLE PRECISION Xf(0:2 * nmax), xxpr(2 * nmax)
-    DOUBLE PRECISION fa, fb
+    REAL(K8) xtild(0:2 * nmax)
+    REAL(K8) fhyp(nmax:2 * nmax), ffdx, beta, Xprimt(0:2 * nmax)
+    REAL(K8) Xf(0:2 * nmax), xxpr(2 * nmax)
+    REAL(K8) fa, fb
     INTEGER i, is2
-    DOUBLE PRECISION xmoy, fxm
+    REAL(K8) xmoy, fxm
 
     !----------------------------------------------------------------------
@@ -175 +175 @@
 
           call invert_zoom_x(xf, xtild, Xprimt, rlonm025(:iim), &
-               xprimm025(:iim), xuv = - 0.25d0)
+               xprimm025(:iim), xuv = - 0.25_k8)
           call invert_zoom_x(xf, xtild, Xprimt, rlonv(:iim), xprimv(:iim), &
-               xuv = 0d0)
+               xuv = 0._k8)
           call invert_zoom_x(xf, xtild, Xprimt, rlonu(:iim), xprimu(:iim), &
-               xuv = 0.5d0)
+               xuv = 0.5_k8)
           call invert_zoom_x(xf, xtild, Xprimt, rlonp025(:iim), &
-               xprimp025(:iim), xuv = 0.25d0)
+               xprimp025(:iim), xuv = 0.25_k8)
        end if test_grossismx
 

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
 
     !-------------------------------------------------------------------

LMDZ5/trunk/libf/dyn3d_common/invert_zoom_x_m.F90

@@ -10 +10 @@
 
     use coefpoly_m, only: coefpoly
-    use nrtype, only: pi, pi_d, twopi_d
+    use nrtype, only: pi, pi_d, twopi_d, k8
 
     include "dimensions.h"
@@ -18 +18 @@
     ! for clon
 
-    DOUBLE PRECISION, intent(in):: Xf(0:), xtild(0:), Xprimt(0:) ! (0:2 * nmax)
+    REAL(K8), intent(in):: Xf(0:), xtild(0:), Xprimt(0:) ! (0:2 * nmax)
     real, intent(out):: xlon(:), xprimm(:) ! (iim)
 
-    DOUBLE PRECISION, intent(in):: xuv
+    REAL(K8), intent(in):: xuv
     ! 0. si calcul aux points scalaires 
     ! 0.5 si calcul aux points U 
 
     ! Local:
-    DOUBLE PRECISION xo1, Xfi, a0, a1, a2, a3, Xf1, Xprimin
+    REAL(K8) xo1, Xfi, a0, a1, a2, a3, Xf1, Xprimin
     integer i, it, iter
-    DOUBLE PRECISION, parameter:: my_eps = 1d-6
+    REAL(K8), parameter:: my_eps = 1e-6_k8
 
-    DOUBLE PRECISION xxprim(iim), xvrai(iim) 
+    REAL(K8) xxprim(iim), xvrai(iim) 
     ! intermediary variables because xlon and xprimm are simple precision
 
@@ -36 +36 @@
 
     DO i = 1, iim
-       Xfi = - pi_d + (i + xuv - 0.75d0) * twopi_d / iim
+       Xfi = - pi_d + (i + xuv - 0.75_k8) * twopi_d / iim
 
        it = 2 * nmax
@@ -54 +54 @@
             xtild(it), xtild(it + 1), a0, a1, a2, a3)
        Xf1 = Xf(it)
-       Xprimin = a1 + xvrai(i) * (2d0 * a2 + xvrai(i) * 3d0 * a3)
+       Xprimin = a1 + xvrai(i) * (2._k8 * a2 + xvrai(i) * 3._k8 * a3)
        xo1 = xvrai(i)
        iter = 1
@@ -63 +63 @@
           xo1 = xvrai(i)
           Xf1 = a0 + xvrai(i) * (a1 + xvrai(i) * (a2 + xvrai(i) * a3))
-          Xprimin = a1 + xvrai(i) * (2d0 * a2 + xvrai(i) * 3d0 * a3)
+          Xprimin = a1 + xvrai(i) * (2._k8 * a2 + xvrai(i) * 3._k8 * a3)
        end DO