/[lmdze]/trunk/dyn3d/fxhyp.f
ViewVC logotype

Diff of /trunk/dyn3d/fxhyp.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 151 by guez, Tue Jun 23 15:14:20 2015 UTC revision 156 by guez, Thu Jul 16 17:39:10 2015 UTC
# Line 24  contains Line 24  contains
24      use principal_cshift_m, only: principal_cshift      use principal_cshift_m, only: principal_cshift
25      use tanh_cautious_m, only: tanh_cautious      use tanh_cautious_m, only: tanh_cautious
26    
27      REAL, intent(out):: xprimm025(:), rlonv(:), xprimv(:) ! (iim + 1)      REAL, intent(out):: xprimm025(:) ! (iim + 1)
28      real, intent(out):: rlonu(:), xprimu(:), xprimp025(:) ! (iim + 1)  
29        REAL, intent(out):: rlonv(:) ! (iim + 1)
30        ! longitudes of points of the "scalar" and "v" grid, in rad
31    
32        REAL, intent(out):: xprimv(:) ! (iim + 1)
33        ! 2 pi / iim * (derivative of the longitudinal zoom function)(rlonv)
34    
35        real, intent(out):: rlonu(:) ! (iim + 1)
36        ! longitudes of points of the "u" grid, in rad
37    
38        real, intent(out):: xprimu(:) ! (iim + 1)
39        ! 2 pi / iim * (derivative of the longitudinal zoom function)(rlonu)
40    
41        real, intent(out):: xprimp025(:) ! (iim + 1)
42    
43      ! Local:      ! Local:
44      real rlonm025(iim + 1), rlonp025(iim + 1), d_rlonv(iim)      real rlonm025(iim + 1), rlonp025(iim + 1), d_rlonv(iim)
45      REAL delta, step      REAL delta, h
46      DOUBLE PRECISION, dimension(0:nmax):: xtild, fhyp, G, Xf, ffdx      DOUBLE PRECISION, dimension(0:nmax):: xtild, fhyp, G, Xf, ffdx
47      DOUBLE PRECISION beta      DOUBLE PRECISION beta
48      INTEGER i, is2      INTEGER i, is2
# Line 40  contains Line 53  contains
53      print *, "Call sequence information: fxhyp"      print *, "Call sequence information: fxhyp"
54    
55      if (grossismx == 1.) then      if (grossismx == 1.) then
56         step = twopi / iim         h = twopi / iim
57    
58         xprimm025(:iim) = step         xprimm025(:iim) = h
59         xprimp025(:iim) = step         xprimp025(:iim) = h
60         xprimv(:iim) = step         xprimv(:iim) = h
61         xprimu(:iim) = step         xprimu(:iim) = h
62    
63         rlonv(:iim) = arth(- pi + clon, step, iim)         rlonv(:iim) = arth(- pi + clon, h, iim)
64         rlonm025(:iim) = rlonv(:iim) - 0.25 * step         rlonm025(:iim) = rlonv(:iim) - 0.25 * h
65         rlonp025(:iim) = rlonv(:iim) + 0.25 * step         rlonp025(:iim) = rlonv(:iim) + 0.25 * h
66         rlonu(:iim) = rlonv(:iim) + 0.5 * step         rlonu(:iim) = rlonv(:iim) + 0.5 * h
67      else      else
68         delta = dzoomx * twopi_d         delta = dzoomx * twopi_d
69         xtild = arth(0d0, pi_d / nmax, nmax + 1)         xtild = arth(0d0, pi_d / nmax, nmax + 1)
Legend:
Removed from v.151  
changed lines
  Added in v.156
  ViewVC Help
Powered by ViewVC 1.1.21