/[lmdze]/trunk/Sources/dyn3d/Guide/init_tau2alpha.f
ViewVC logotype

Diff of /trunk/Sources/dyn3d/Guide/init_tau2alpha.f

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

revision 134 by guez, Wed Apr 29 15:47:56 2015 UTC revision 173 by guez, Tue Oct 6 15:57:02 2015 UTC
# Line 8  contains Line 8  contains
8    
9    SUBROUTINE init_tau2alpha(dxdys, dxdyu, dxdyv)    SUBROUTINE init_tau2alpha(dxdys, dxdyu, dxdyv)
10    
11      USE comgeom, ONLY: cu_2d, cv_2d, rlatu      USE comgeom, ONLY: cu_2d, cv_2d
12      use conf_guide_m, only: guide_u, guide_v      use conf_guide_m, only: guide_u, guide_v
13      use coordij_m, only: coordij      use coordij_m, only: coordij
14      USE dimens_m, ONLY: iim, jjm      USE dimens_m, ONLY: iim, jjm
15        USE dynetat0_m, ONLY: clat, clon, grossismx, grossismy, rlatu
16      USE paramet_m, ONLY: iip1, jjp1      USE paramet_m, ONLY: iip1, jjp1
     USE serre, ONLY: clat, clon, grossismx, grossismy  
17      use writefield_m, only: writefield      use writefield_m, only: writefield
18    
19      REAL, intent(out):: dxdys(iip1, jjp1), dxdyu(iip1, jjp1), dxdyv(iip1, jjm)      REAL, intent(out):: dxdys(iip1, jjp1), dxdyu(iip1, jjp1), dxdyv(iip1, jjm)
# Line 72  contains Line 72  contains
72      CALL coordij(clon, clat, ilon, ilat)      CALL coordij(clon, clat, ilon, ilat)
73      ! aire de la maille au centre du zoom      ! aire de la maille au centre du zoom
74      dxdy_min = dxdys(ilon, ilat)      dxdy_min = dxdys(ilon, ilat)
75        print *, "dxdy_min = ", dxdy_min
76    
77      ! dxdy maximal de la maille :      ! dxdy maximal de la maille :
78      dxdy_max = 0.      dxdy_max = maxval(dxdys)
79      DO j = 1, jjp1      print *, "dxdy_max = ", dxdy_max
        DO i = 1, iip1  
           dxdy_max = max(dxdy_max, dxdys(i, j))  
        END DO  
     END DO  
80    
81      IF (abs(grossismx - 1.) < 0.1 .OR. abs(grossismy - 1.) < 0.1) THEN      IF (abs(grossismx - 1.) < 0.1 .OR. abs(grossismy - 1.) < 0.1) THEN
82         PRINT *, 'Attention : modèle peu zoomé.'         PRINT *, 'Attention : modèle peu zoomé.'
Legend:
Removed from v.134  
changed lines
  Added in v.173
  ViewVC Help
Powered by ViewVC 1.1.21