source: trunk/LMDZ.MARS/libf/dyn3d/coordij.F @ 937

Last change on this file since 937 was 38, checked in by emillour, 14 years ago

Ajout du modè Martien (mon LMDZ.MARS.BETA, du 28/01/2011) dans le rértoire mars, pour pouvoir suivre plus facilement les modifs.
EM
File size: 974 bytes
Line 
1      SUBROUTINE coordij(lon,lat,ilon,jlat)
2
3c=======================================================================
4c
5c   calcul des coordonnees i et j de la maille scalaire dans
6c   laquelle se trouve le point (lon,lat) en radian
7c
8c=======================================================================
9
10      IMPLICIT NONE
11      REAL lon,lat
12      INTEGER ilon,jlat
13      INTEGER i,j
14
15#include "dimensions.h"
16#include "paramet.h"
17#include "comconst.h"
18#include "comgeom.h"
19#include "serre.h"
20
21      real zlon,zlat
22
23      zlon=lon
24      zlat=lat
25
26      DO i=1,iim+1
27        print*,i,' ',rlonu(i),' ',zlon
28         IF (rlonu(i).GT.zlon) THEN
29            ilon=i
30            GOTO 10
31         ENDIF
32      ENDDO
3310    CONTINUE
34
35      j=0
36      DO j=1,jjm
37        print*,j,' ',rlatv(j),' ' ,zlat
38         IF(rlatv(j).LT.zlat) THEN
39           
40            jlat=j
41            GOTO 20
42         ENDIF
43      ENDDO
4420    CONTINUE
45      IF(j.EQ.0) j=jjm+1
46
47      RETURN
48      END
Note: See TracBrowser for help on using the repository browser.