/[lmdze]/trunk/dyn3d/coordij.f90
ViewVC logotype

Diff of /trunk/dyn3d/coordij.f90

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

trunk/dyn3d/coordij.f revision 80 by guez, Fri Nov 15 18:45:49 2013 UTC trunk/dyn3d/coordij.f90 revision 81 by guez, Wed Mar 5 14:38:41 2014 UTC
# Line 1  Line 1 
 !  
 ! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/coordij.F,v 1.1.1.1 2004/05/19 12:53:05 lmdzadmin Exp $  
 !  
       SUBROUTINE coordij(lon,lat,ilon,jlat)  
   
 c=======================================================================  
 c  
 c   calcul des coordonnees i et j de la maille scalaire dans  
 c   laquelle se trouve le point (lon,lat) en radian  
 c  
 c=======================================================================  
   
       use dimens_m  
       use paramet_m  
       use comconst  
       use comgeom  
       use serre  
       use nr_util, only: pi  
       IMPLICIT NONE  
       REAL lon,lat  
       INTEGER ilon,jlat  
       INTEGER i,j  
   
   
       real zlon,zlat  
   
       zlon=lon*pi/180.  
       zlat=lat*pi/180.  
   
       DO i=1,iim+1  
          IF (rlonu(i).GT.zlon) THEN  
             ilon=i  
             GOTO 10  
          ENDIF  
       ENDDO  
 10    CONTINUE  
   
       j=0  
       DO j=1,jjm  
          IF(rlatv(j).LT.zlat) THEN  
             jlat=j  
             GOTO 20  
          ENDIF  
       ENDDO  
 20    CONTINUE  
       IF(j.EQ.0) j=jjm+1  
1    
2        RETURN  ! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/coordij.F,v 1.1.1.1 2004/05/19
3        END  ! 12:53:05 lmdzadmin Exp $
4    
5    SUBROUTINE coordij(lon, lat, ilon, jlat)
6    
7      ! =======================================================================
8    
9      ! calcul des coordonnees i et j de la maille scalaire dans
10      ! laquelle se trouve le point (lon,lat) en radian
11    
12      ! =======================================================================
13    
14      USE dimens_m
15      USE paramet_m
16      USE comconst
17      USE comgeom
18      USE serre
19      USE nr_util, ONLY: pi
20      IMPLICIT NONE
21      REAL lon, lat
22      INTEGER ilon, jlat
23      INTEGER i, j
24    
25    
26      REAL zlon, zlat
27    
28      zlon = lon*pi/180.
29      zlat = lat*pi/180.
30    
31      DO i = 1, iim + 1
32        IF (rlonu(i)>zlon) THEN
33          ilon = i
34          GO TO 10
35        END IF
36      END DO
37    10 CONTINUE
38    
39      j = 0
40      DO j = 1, jjm
41        IF (rlatv(j)<zlat) THEN
42          jlat = j
43          GO TO 20
44        END IF
45      END DO
46    20 CONTINUE
47      IF (j==0) j = jjm + 1
48    
49      RETURN
50    END SUBROUTINE coordij
Legend:
Removed from v.80  
changed lines
  Added in v.81
  ViewVC Help
Powered by ViewVC 1.1.21