dyn3d/Guide/coordij.f


! $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)

  ! =======================================================================

  ! calcul des coordonnees i et j de la maille scalaire dans
  ! laquelle se trouve le point (lon,lat) en radian

  ! =======================================================================

  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)>zlon) THEN
      ilon = i
      GO TO 10
    END IF
  END DO
10 CONTINUE

  j = 0
  DO j = 1, jjm
    IF (rlatv(j)<zlat) THEN
      jlat = j
      GO TO 20
    END IF
  END DO
20 CONTINUE
  IF (j==0) j = jjm + 1

  RETURN
END SUBROUTINE coordij
1
2	! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/coordij.F,v 1.1.1.1 2004/05/19
3	! 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