/[lmdze]/trunk/dyn3d/Guide/tau2alpha.f
ViewVC logotype

Diff of /trunk/dyn3d/Guide/tau2alpha.f

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

revision 103 by guez, Fri Aug 29 13:00:05 2014 UTC revision 109 by guez, Wed Sep 17 10:08:00 2014 UTC
# Line 1  Line 1 
1  module tau2alpha_m  module tau2alpha_m
2    
   USE paramet_m, ONLY : iip1, jjp1  
   USE dimens_m, ONLY : jjm  
   
3    IMPLICIT NONE    IMPLICIT NONE
4    
   private iip1, jjp1, jjm  
   
   REAL dxdys(iip1, jjp1), dxdyu(iip1, jjp1), dxdyv(iip1, jjm)  
   
5  contains  contains
6    
7    SUBROUTINE tau2alpha(type, factt, taumin, taumax, alpha)    SUBROUTINE tau2alpha(type, factt, taumin, taumax, alpha)
8    
9      USE comgeom, ONLY : cu_2d, cv_2d, rlatu, rlatv      USE comgeom, ONLY : cu_2d, cv_2d, rlatu, rlatv
10        USE dimens_m, ONLY : jjm
11      use conf_guide_m, only: lat_min_guide, lat_max_guide      use conf_guide_m, only: lat_min_guide, lat_max_guide
12      USE dimens_m, ONLY : iim      USE dimens_m, ONLY : iim
     use dump2d_m, only: dump2d  
13      USE nr_util, ONLY : pi      USE nr_util, ONLY : pi
14        USE paramet_m, ONLY : iip1, jjp1
15      USE serre, ONLY : clat, clon, grossismx, grossismy      USE serre, ONLY : clat, clon, grossismx, grossismy
16        use writefield_m, only: writefield
17    
18      INTEGER, intent(in):: type      INTEGER, intent(in):: type
19      REAL, intent(in):: factt, taumin, taumax      REAL, intent(in):: factt, taumin, taumax
# Line 33  contains Line 28  contains
28      LOGICAL:: first = .TRUE.      LOGICAL:: first = .TRUE.
29      REAL zdx(iip1, jjp1), zdy(iip1, jjp1)      REAL zdx(iip1, jjp1), zdy(iip1, jjp1)
30      REAL zlat      REAL zlat
31        REAL dxdys(iip1, jjp1), dxdyu(iip1, jjp1), dxdyv(iip1, jjm)
32    
33      !------------------------------------------------------------      !------------------------------------------------------------
34    
# Line 54  contains Line 50  contains
50            zdy(i, 1) = zdy(i, 2)            zdy(i, 1) = zdy(i, 2)
51            zdy(i, jjp1) = zdy(i, jjm)            zdy(i, jjp1) = zdy(i, jjm)
52         END DO         END DO
53    
54         DO j = 1, jjp1         DO j = 1, jjp1
55            DO i = 1, iip1            DO i = 1, iip1
56               dxdys(i, j) = sqrt(zdx(i, j)**2 + zdy(i, j)**2)               dxdys(i, j) = sqrt(zdx(i, j)**2 + zdy(i, j)**2)
57            END DO            END DO
58         END DO         END DO
59           CALL writefield("dxdys", dxdys)
60    
61         DO j = 1, jjp1         DO j = 1, jjp1
62            DO i = 1, iim            DO i = 1, iim
63               dxdyu(i, j) = 0.5 * (dxdys(i, j) + dxdys(i + 1, j))               dxdyu(i, j) = 0.5 * (dxdys(i, j) + dxdys(i + 1, j))
# Line 71  contains Line 70  contains
70            END DO            END DO
71         END DO         END DO
72    
        CALL dump2d(iip1, jjp1, dxdys, 'DX2DY2 SCAL ')  
        CALL dump2d(iip1, jjp1, dxdyu, 'DX2DY2 U ')  
        CALL dump2d(iip1, jjp1, dxdyv, 'DX2DY2 v ')  
   
73         ! coordonnees du centre du zoom         ! coordonnees du centre du zoom
74         CALL coordij(clon, clat, ilon, ilat)         CALL coordij(clon, clat, ilon, ilat)
75         ! aire de la maille au centre du zoom         ! aire de la maille au centre du zoom
Legend:
Removed from v.103  
changed lines
  Added in v.109
  ViewVC Help
Powered by ViewVC 1.1.21