# lift forces

Hi Dear All,
I facing one problem about lift coefficient. The drag value is good but even i run the code for Re=100 i npt get the value of lift. Can any one have such experience before and tell about the posibilities that why it’s happened. The forces i calculate like this

void get_force(int wall[NX][NY], int wb[NX][NY], double f_x, double f_y)
{
int i,j; f_x = 0.0; f_y = 0.0;
for(i=1;i<=NX-2;i++){ int ip = i+1; int im = i-1;
for(j=1;j<=NY-2;j++){ int jp = j+1; int jm = j-1;
if(wb[i][j]==1){
f_x = f_x+(f_1[ip][j]+f_1[i][j]) (1-wall[ip][j]) - (f_1[im][j]+f_1[i][j]) (1-wall[im][j])
+ (f_1[ip][jp]+f_1[i][j])
(1-wall[ip][jp]) - (f_1[im][jp]+f_1[i][j])
(1-wall[im][jp])
- (f_1[im][jm]+f_1[i][j])
(1-wall[im][jm]) + (f_1[ip][jm]+f_1[i][j])
(1-wall[ip][jm]);

``````	        *f_y = *f_y+(f_1[i][jp]+f_1[i][j])*(1-wall[i][jp])   	- (f_1[i][jm]+f_1[i][j])*(1-wall[i][jm])
+ (f_1[ip][jp]+f_1[i][j])*(1-wall[ip][jp])+ (f_1[im][jp]+f_1[i][j])*(1-wall[im][jp])
- (f_1[im][jm]+f_1[i][j])*(1-wall[im][jm]) - (f_1[ip][jm]+f_1[i][j])*(1-wall[ip][jm]);

}		}	}
``````

where wb i used for the solid objet like circular cylinder and wall means the computational domain NX and NY.

Then i calculate the drag and lift