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

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

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

trunk/libf/dyn3d/coordij.f revision 39 by guez, Tue Jan 25 15:11:05 2011 UTC trunk/dyn3d/Guide/coordij.f revision 127 by guez, Tue Feb 10 17:58:56 2015 UTC
# Line 1  Line 1 
1  !  module coordij_m
 ! $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  
2    
3        RETURN    IMPLICIT NONE
4        END  
5    contains
6    
7      SUBROUTINE coordij(lon, lat, ilon, jlat)
8    
9        ! From LMDZ4/libf/dyn3d/coordij.F, version 1.1.1.1 2004/05/19 12:53:05
10    
11        ! calcul des coordonnees i et j de la maille scalaire dans
12        ! laquelle se trouve le point (lon, lat) en radian
13    
14        USE dimens_m, only: iim, jjm
15        USE comgeom, only: rlonu, rlatv
16        USE nr_util, ONLY: pi
17    
18        REAL, intent(in):: lon, lat
19        INTEGER ilon, jlat
20        INTEGER i, j
21    
22        DO i = 1, iim + 1
23           IF (rlonu(i)>lon) THEN
24              ilon = i
25              exit
26           END IF
27        END DO
28    
29        j = 0
30        DO j = 1, jjm
31           IF (rlatv(j)<lat) THEN
32              jlat = j
33              exit
34           END IF
35        END DO
36    
37        IF (j==0) j = jjm + 1
38    
39      END SUBROUTINE coordij
40    
41    end module coordij_m
Legend:
Removed from v.39  
changed lines
  Added in v.127
  ViewVC Help
Powered by ViewVC 1.1.21