lg_ind = 4 + 3*(ngrid - 1)
j = 3 + 3*(ngrid + 1)
mnc_f_info(indf,3) = ngrid + 1
integer i,j,k,l, n,n1, ngrid, ds,de
ngrid = mnc_f_info(indf,3)
DO i = 1,ngrid
integer i,j,k,ii,jj,kk, n,nf, indf,indg,indd, fid, ngrid
ng_ind = 4 + 3*ngrid
IF (ngrid .EQ. 0) THEN
INTEGER i,j,k, n, nf, indf,indv, fid, nd, ngrid, is,ie, err
IF (ngrid .LT. 1) THEN