1 | forward(YREAL x, YREAL y) |
---|
2 | { |
---|
3 | //Yj = longitude |
---|
4 | //Yk = latitude |
---|
5 | int j_ll=(int)x; |
---|
6 | int k_ll=(int)y; |
---|
7 | if (j_ll==Yj && k_ll==Yk) |
---|
8 | { |
---|
9 | //printf("j_ll=%d ; k_ll=%d\n",j_ll,k_ll); |
---|
10 | YS1_uinter(Yi,YTemps)=YS1_u(Yj,Yk);//la 1ere sortie du module uinter pour le floater i est la valeur du module u en Yj Yk |
---|
11 | YS2_uinter(Yi,YTemps)=YS1_u(Yj+1,Yk); |
---|
12 | YS3_uinter(Yi,YTemps)=YS1_u(Yj,Yk+1); |
---|
13 | YS4_uinter(Yi,YTemps)=YS1_u(Yj+1,Yk+1); |
---|
14 | YS5_uinter(Yi,YTemps)=YS1_v(Yj,Yk); |
---|
15 | YS6_uinter(Yi,YTemps)=YS1_v(Yj+1,Yk); |
---|
16 | YS7_uinter(Yi,YTemps)=YS1_v(Yj,Yk+1); |
---|
17 | YS8_uinter(Yi,YTemps)=YS1_v(Yj+1,Yk+1); |
---|
18 | YS1_xinter(Yi,YTemps)=(YREAL)Yj; |
---|
19 | YS2_xinter(Yi,YTemps)=(YREAL)Yj+1; |
---|
20 | YS3_xinter(Yi,YTemps)=(YREAL)Yj; |
---|
21 | YS4_xinter(Yi,YTemps)=(YREAL)Yj+1; |
---|
22 | YS5_xinter(Yi,YTemps)=(YREAL)Yk; |
---|
23 | YS6_xinter(Yi,YTemps)=(YREAL)Yk; |
---|
24 | YS7_xinter(Yi,YTemps)=(YREAL)Yk+1; |
---|
25 | YS8_xinter(Yi,YTemps)=(YREAL)Yk+1; |
---|
26 | |
---|
27 | } |
---|
28 | } |
---|
29 | |
---|
30 | backward(YREAL x, YREAL y) |
---|
31 | { |
---|
32 | int j_ll=(int)x; |
---|
33 | int k_ll=(int)y; |
---|
34 | if (j_ll==Yj && k_ll==Yk) |
---|
35 | { |
---|
36 | YG1_u(Yj,Yk)+=YG1_uinter(Yi,YTemps); |
---|
37 | YG1_u(Yj+1,Yk)+=YG2_uinter(Yi,YTemps); |
---|
38 | YG1_u(Yj,Yk+1)+=YG3_uinter(Yi,YTemps); |
---|
39 | YG1_u(Yj+1,Yk+1)+=YG4_uinter(Yi,YTemps); |
---|
40 | YG1_v(Yj,Yk)+=YG5_uinter(Yi,YTemps); |
---|
41 | YG1_v(Yj+1,Yk)+=YG6_uinter(Yi,YTemps); |
---|
42 | YG1_v(Yj,Yk+1)+=YG7_uinter(Yi,YTemps); |
---|
43 | YG1_v(Yj+1,Yk+1)+=YG8_uinter(Yi,YTemps); |
---|
44 | |
---|
45 | } |
---|
46 | |
---|
47 | |
---|
48 | |
---|
49 | } |
---|