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

Contents of /trunk/dyn3d/Guide/coordij.f

Parent Directory Parent Directory | Revision Log Revision Log


Revision 115 - (show annotations)
Fri Sep 19 17:36:20 2014 UTC (9 years, 8 months ago) by guez
File size: 936 byte(s)
Extracted code from tau2alpha for first call into new procedure
init_tau2alpha. dxdys, dxdyu, dxdyv are now local variables if guide
computed by init_tau2alpha. This allows us to remove terrible argument
type of tau2alpha: we just give to tau2alpha the right dxdy and
rlat.

In module conf_guide_m, changed default values of tau_min_*, because
0.02 is too small for the default daystep = 240, iperiod = 5. Changed
default values of guide_[uv] to false. Moved variable ok_guide from
conf_gcm_m to conf_guide_m, ok_guide is no longer an input parameter,
it is computed from guide_*. Had then to move test on ok_guide and
day_step from conf_gcm_m to conf_guide_m. Added checks on input
nudging parameters in procedure conf_guide. Upgraded variable factt to
module conf_guide_m in order to check nudging parameters. Bug fix:
variable guide_q was not in namelist conf_guide_nml.

Removed unused variables aire_min, aire_max of MODULE guide_m.

Moved the call to conf_guide from guide to gcm. This was needed to
define ok_guide before getting into guide (since ok_guide is no longer
in conf_gcm_m). Moved test on grossismx and grossismy from tau2alpha
to guide. It is clearer now that only tau_max is used for a regular
grid, and we do not have to repeat this test in each call to
tau2alpha. In guide, we only call writefield when alpha is not a
constant.

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, intent(in):: 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

  ViewVC Help
Powered by ViewVC 1.1.21