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

Diff of /trunk/dyn3d/fxy.f

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

revision 76 by guez, Fri Nov 15 18:45:49 2013 UTC revision 78 by guez, Wed Feb 5 17:51:07 2014 UTC
# Line 87  contains Line 87  contains
87    !******************************************************    !******************************************************
88    
89    REAL FUNCTION fx(ri)    REAL FUNCTION fx(ri)
90    
91      ! stretching in x      ! stretching in x
92      USE nr_util, ONLY : pi      USE nr_util, ONLY : pi
93      USE serre, ONLY : alphax, pxo, transx      USE serre, ONLY : alphax, pxo, transx
94      REAL, INTENT (IN) :: ri      REAL, INTENT (IN) :: ri
95    
96      fx = ripx(ri) + transx + alphax * SIN(ripx(ri) + transx - pxo) - pi      fx = ripx(ri) + transx + alphax * SIN(ripx(ri) + transx - pxo) - pi
97    
98    end function fx    end function fx
99    
100    !******************************************************    !******************************************************
# Line 114  contains Line 116  contains
116    !******************************************************    !******************************************************
117    
118    REAL FUNCTION bigy(rj)    REAL FUNCTION bigy(rj)
119    
120      ! stretching in y      ! stretching in y
121      USE nr_util, ONLY : pi      USE nr_util, ONLY : pi
122      USE dimens_m, ONLY : jjm      USE dimens_m, ONLY : jjm
123      REAL, INTENT (IN) :: rj      REAL, INTENT (IN) :: rj
124    
125      bigy = 2 * (REAL(jjm + 1) - rj) * pi / jjm      bigy = 2 * (REAL(jjm + 1) - rj) * pi / jjm
126    
127    end function bigy    end function bigy
128    
129    !******************************************************    !******************************************************
130    
131    REAL FUNCTION fy(rj)    REAL FUNCTION fy(rj)
132    
133      ! stretching in y      ! stretching in y
134      USE nr_util, ONLY : pi      USE nr_util, ONLY : pi
135      USE serre, ONLY : alphay, pyo, transy      USE serre, ONLY : alphay, pyo, transy
# Line 132  contains Line 137  contains
137    
138      fy = (bigy(rj) + transy + alphay * SIN(bigy(rj) + transy - pyo)) / 2 &      fy = (bigy(rj) + transy + alphay * SIN(bigy(rj) + transy - pyo)) / 2 &
139           - pi / 2           - pi / 2
140    
141    end function fy    end function fy
142    
143    !******************************************************    !******************************************************
Legend:
Removed from v.76  
changed lines
  Added in v.78
  ViewVC Help
Powered by ViewVC 1.1.21