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

Diff of /trunk/dyn3d/fyhyp.f

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

revision 149 by guez, Thu Jun 18 12:23:44 2015 UTC revision 212 by guez, Thu Jan 12 12:31:31 2017 UTC
# Line 20  contains Line 20  contains
20      use dynetat0_m, only: clat, grossismy, dzoomy, tauy      use dynetat0_m, only: clat, grossismy, dzoomy, tauy
21      use heavyside_m, only: heavyside      use heavyside_m, only: heavyside
22    
23      REAL, intent(out):: rlatu(jjm + 1)      REAL, intent(out):: rlatu(:) ! (jjm + 1)
24      REAL, intent(out):: rlatv(jjm)      REAL, intent(out):: rlatv(:) ! (jjm)
25      real, intent(out):: rlatu2(jjm), yprimu2(jjm), rlatu1(jjm), yprimu1(jjm)      real, intent(out):: rlatu2(:), yprimu2(:), rlatu1(:), yprimu1(:) ! (jjm)
26    
27      ! Local:      ! Local:
28    
     DOUBLE PRECISION champmin, champmax  
29      INTEGER, PARAMETER:: nmax=30000, nmax2=2*nmax      INTEGER, PARAMETER:: nmax=30000, nmax2=2*nmax
30      REAL dzoom ! distance totale de la zone du zoom (en radians)      REAL dzoom ! distance totale de la zone du zoom (en radians)
31      DOUBLE PRECISION ylat(jjm + 1), yprim(jjm + 1)      DOUBLE PRECISION ylat(jjm + 1), yprim(jjm + 1)
# Line 41  contains Line 40  contains
40      DOUBLE PRECISION pi, pis2, epsilon, pisjm      DOUBLE PRECISION pi, pis2, epsilon, pisjm
41      DOUBLE PRECISION yo1, yi, ylon2, ymoy, yprimin      DOUBLE PRECISION yo1, yi, ylon2, ymoy, yprimin
42      DOUBLE PRECISION yfi, yf1, ffdy      DOUBLE PRECISION yfi, yf1, ffdy
43      DOUBLE PRECISION ypn, deply, y00      DOUBLE PRECISION ypn
44      SAVE y00, deply      DOUBLE PRECISION, save::deply, y00
45    
46      INTEGER i, j, it, ik, iter, jlat      INTEGER i, j, it, ik, iter, jlat, jjpn
47      INTEGER jpn, jjpn      INTEGER, save:: jpn
     SAVE jpn  
48      DOUBLE PRECISION yi2, heavyy0, heavyy0m      DOUBLE PRECISION yi2, heavyy0, heavyy0m
49      DOUBLE PRECISION fa(0:nmax2), fb(0:nmax2)      DOUBLE PRECISION fa(0:nmax2), fb(0:nmax2)
50      REAL y0min, y0max      REAL y0min, y0max
# Line 280  contains Line 278  contains
278      DO j = 1, jjm      DO j = 1, jjm
279         ylat(j) = rlatu(j) - rlatu(j + 1)         ylat(j) = rlatu(j) - rlatu(j + 1)
280      END DO      END DO
     champmin = 1e12  
     champmax = -1e12  
     DO j = 1, jjm  
        champmin = min(champmin, ylat(j))  
        champmax = max(champmax, ylat(j))  
     END DO  
     champmin = champmin*180./pi  
     champmax = champmax*180./pi  
281    
282      DO j = 1, jjm      DO j = 1, jjm
283         IF (rlatu1(j) <= rlatu2(j)) THEN         IF (rlatu1(j) <= rlatu2(j)) THEN
# Line 322  contains Line 312  contains
312      ENDDO      ENDDO
313    
314      print *, 'Latitudes'      print *, 'Latitudes'
315      print 3, champmin, champmax      print 3, minval(ylat(:jjm)) *180d0/pi, maxval(ylat(:jjm))*180d0/pi
316    
317  3   Format(1x, ' Au centre du zoom, la longueur de la maille est', &  3   Format(1x, ' Au centre du zoom, la longueur de la maille est', &
318           ' d environ ', f0.2, ' degres ', /, &           ' d environ ', f0.2, ' degres ', /, &
Legend:
Removed from v.149  
changed lines
  Added in v.212
  ViewVC Help
Powered by ViewVC 1.1.21