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

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

Parent Directory Parent Directory | Revision Log Revision Log


Revision 115 - (hide annotations)
Fri Sep 19 17:36:20 2014 UTC (9 years, 9 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 guez 3
2 guez 81 ! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/coordij.F,v 1.1.1.1 2004/05/19
3     ! 12:53:05 lmdzadmin Exp $
4 guez 3
5 guez 81 SUBROUTINE coordij(lon, lat, ilon, jlat)
6 guez 3
7 guez 81 ! =======================================================================
8 guez 3
9 guez 81 ! calcul des coordonnees i et j de la maille scalaire dans
10     ! laquelle se trouve le point (lon,lat) en radian
11 guez 3
12 guez 81 ! =======================================================================
13 guez 3
14 guez 81 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 guez 115 REAL, intent(in):: lon, lat
22 guez 81 INTEGER ilon, jlat
23     INTEGER i, j
24 guez 3
25    
26 guez 81 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