source: altifloat/src/locate.h @ 84

Last change on this file since 84 was 84, checked in by leila_ocean, 11 years ago

initial import
File size: 1.3 KB
Line 
1forward(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
30backward(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}
Note: See TracBrowser for help on using the repository browser.