--> //initializations

--> h1=7.0 //m (total head at extremities of floor 1)
 h1  = 

   7.


--> h2=0.5 //m (total head at extremities of floor 2)
 h2  = 

   0.5


--> Yw=9.81 //KN/m^3 (unit weight of water)
 Yw  = 

   9.81


--> ds=1 //m (dimension of flow field)
 ds  = 

   1.


--> L=28.5 //m (length)
 L  = 

   28.5


--> t=1 //m (thickness of floor)
 t  = 

   1.


--> nd=15 
 nd  = 

   15.


--> x=7.5 //m (height from datum)
 x  = 

   7.5


--> //calculations

--> dh=x/nd // (equipotential drop in the flow field)
 dh  = 

   0.5


--> A=L*t // (area)
 A  = 

   28.5


--> U=((h1+h2)*Yw*A)/2 // (uplift pressure)
 U  = 

   1048.4438


--> i=dh/ds // (gradient)
 i  = 

   0.5


--> //results

--> printf('uplift pressure on the impervious concrete floor=%f',U) //KN/m
uplift pressure on the impervious concrete floor=1048.443750
--> printf('exit gradient=%f',i) 
exit gradient=0.500000