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 | ! |
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 | use nr_util, only: pi |
19 | 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 |