Parent Directory | Revision Log
"pi" comes from "nr_util". Removed subroutine "initialize" in module "comconst". Copied the content of "fxy_sin.h" into "fxysinus", instead of getting it from an "include" line. Removed file "fxy_sin.h". "ps" has rank 2 in "gcm" and "dynetat0". Assumed-shape for argument "q" of "integrd".
1 | guez | 3 | ! |
2 | ! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/coordij.F,v 1.1.1.1 2004/05/19 12:53:05 lmdzadmin Exp $ | ||
3 | ! | ||
4 | SUBROUTINE coordij(lon,lat,ilon,jlat) | ||
5 | |||
6 | c======================================================================= | ||
7 | c | ||
8 | c calcul des coordonnees i et j de la maille scalaire dans | ||
9 | c laquelle se trouve le point (lon,lat) en radian | ||
10 | c | ||
11 | c======================================================================= | ||
12 | |||
13 | use dimens_m | ||
14 | use paramet_m | ||
15 | use comconst | ||
16 | use comgeom | ||
17 | use serre | ||
18 | guez | 39 | use nr_util, only: pi |
19 | guez | 3 | IMPLICIT NONE |
20 | REAL lon,lat | ||
21 | INTEGER ilon,jlat | ||
22 | INTEGER i,j | ||
23 | |||
24 | |||
25 | real zlon,zlat | ||
26 | |||
27 | zlon=lon*pi/180. | ||
28 | zlat=lat*pi/180. | ||
29 | |||
30 | DO i=1,iim+1 | ||
31 | IF (rlonu(i).GT.zlon) THEN | ||
32 | ilon=i | ||
33 | GOTO 10 | ||
34 | ENDIF | ||
35 | ENDDO | ||
36 | 10 CONTINUE | ||
37 | |||
38 | j=0 | ||
39 | DO j=1,jjm | ||
40 | IF(rlatv(j).LT.zlat) THEN | ||
41 | jlat=j | ||
42 | GOTO 20 | ||
43 | ENDIF | ||
44 | ENDDO | ||
45 | 20 CONTINUE | ||
46 | IF(j.EQ.0) j=jjm+1 | ||
47 | |||
48 | RETURN | ||
49 | END |
ViewVC Help | |
Powered by ViewVC 1.1.21 |