[84] | 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 | } |
---|