*ID MJRADVCH
*DECLARE ADVSRCE
*D ADVSRCE.391,ADVSRCE.394
          x_div(i,k)=(flux(i,k)-flux(i+1,k))*dxr(i)*recip_cos(j)
          source(i,k)=x_div(i,k)                        
*D ADVSRCE.536,ADVSRCE.542
          y_div(i,k)=(fluxs(i,k)-fluxn(i,k))*dyr(j)  
          source(i,k)=source(i,k)+y_div(i,k)        
*D ADVSRCE.646,ADVSRCE.648
          z_div(i,k)=(flux(i,k+1)-flux(i,k))/dz(k)  
          source(i,k)=source(i,k)+z_div(i,k)       
  
