/[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 105 by guez, Thu Sep 4 10:40:24 2014 UTC revision 112 by guez, Thu Sep 18 13:36:51 2014 UTC
# Line 25  contains Line 25  contains
25    
26      REAL, intent(in):: dzooma ! distance totale de la zone du zoom      REAL, intent(in):: dzooma ! distance totale de la zone du zoom
27    
28      REAL, intent(in):: tau      REAL, intent(in):: tau
29      ! raideur de la transition de l'intérieur à l'extérieur du zoom      ! raideur de la transition de l'intérieur à l'extérieur du zoom
30    
31      ! arguments de sortie      ! arguments de sortie
# Line 80  contains Line 80  contains
80         STOP 1         STOP 1
81      ELSE      ELSE
82         dzoom = dzooma * pi/180.         dzoom = dzooma * pi/180.
83      ENDIF      END IF
84    
85      print *, ' xzoom(rad), grossism, tau, dzoom (rad):'      print *, ' xzoom(rad), grossism, tau, dzoom (rad):'
86      print *, xzoom, grossism, tau, dzoom      print *, xzoom, grossism, tau, dzoom
87    
88      DO i = 0, nmax2      DO i = 0, nmax2
89         xtild(i) = - pi + FLOAT(i) * depi /nmax2         xtild(i) = - pi + REAL(i) * depi /nmax2
90      ENDDO      ENDDO
91    
92      DO i = nmax, nmax2      DO i = nmax, nmax2
# Line 107  contains Line 107  contains
107            ELSE            ELSE
108               fhyp (i) = TANH (fa/fb)               fhyp (i) = TANH (fa/fb)
109            ENDIF            ENDIF
110         ENDIF         END IF
111    
112         IF (xtild(i) == 0.) fhyp(i) = 1.         IF (xtild(i) == 0.) fhyp(i) = 1.
113         IF (xtild(i) == pi) fhyp(i) = -1.         IF (xtild(i) == pi) fhyp(i) = -1.
114      ENDDO      END DO
115    
116      ! Calcul de beta      ! Calcul de beta
117    
# Line 150  contains Line 150  contains
150         print *, 'Attention ! La valeur beta calculée dans fxhyp est mauvaise.'         print *, 'Attention ! La valeur beta calculée dans fxhyp est mauvaise.'
151         print *, 'Modifier les valeurs de grossismx, tau ou dzoomx et relancer.'         print *, 'Modifier les valeurs de grossismx, tau ou dzoomx et relancer.'
152         STOP 1         STOP 1
153      ENDIF      END IF
154    
155      ! calcul de Xprimt      ! calcul de Xprimt
156    
157      DO i = nmax, nmax2      DO i = nmax, nmax2
158         Xprimt(i) = beta + (grossism - beta) * fhyp(i)         Xprimt(i) = beta + (grossism - beta) * fhyp(i)
159      ENDDO      END DO
160    
161      DO i = nmax + 1, nmax2      DO i = nmax + 1, nmax2
162         Xprimt(nmax2 - i) = Xprimt(i)         Xprimt(nmax2 - i) = Xprimt(i)
163      ENDDO      END DO
164    
165      ! Calcul de Xf      ! Calcul de Xf
166    
# Line 218  contains Line 218  contains
218         ENDIF         ENDIF
219    
220         DO i = ii1, ii2         DO i = ii1, ii2
221            xlon2 = - pi + (FLOAT(i) + xuv - decalx) * depi / FLOAT(iim)            xlon2 = - pi + (REAL(i) + xuv - decalx) * depi / REAL(iim)
222            Xfi = xlon2            Xfi = xlon2
223    
224            it = nmax2            it = nmax2
# Line 264  contains Line 264  contains
264            end if            end if
265    
266    
267            xxprim(i) = depi/ (FLOAT(iim) * Xprimin)            xxprim(i) = depi/ (REAL(iim) * Xprimin)
268            xvrai(i) = xi + xzoom            xvrai(i) = xi + xzoom
269         end DO         end DO
270    
Legend:
Removed from v.105  
changed lines
  Added in v.112
  ViewVC Help
Powered by ViewVC 1.1.21