pblindex(indx) = (j-1)*im + i

       levgather(indx) = levpbl(pblindex(indx),1)

        pigather(indx) =     pz(pblindex(indx),1)

        pkegather(indx,lm+1) = pkht(pblindex(indx),1,lm+1)

        plegather(indx,lm+1) = plze(pblindex(indx),1,lm+1)

         thgather(indx,L) = tz(pblindex(indx),1,L)

         shgather(indx,L) = qz(pblindex(indx),1,L,1)

        pkegather(indx,L) = pkht(pblindex(indx),1,L)

        pkzgather(indx,L) = pkl(pblindex(indx),1,L)

        plegather(indx,L) = plze(pblindex(indx),1,L)

        plzgather(indx,L) = plz(pblindex(indx),1,L)

         dpgather(indx,L) = dpres(pblindex(indx),1,L)

        ugather(indx,L,nt) = qz(pblindex(indx),1,L,nt+ptracer)

        ugather(indx,L,ntracerin-ptracer+1) = uz(pblindex(indx),1,L)

        ugather(indx,L,ntracerin-ptracer+2) = vz(pblindex(indx),1,L)

      call BACK2GRD ( lsp_new,pblindex, lsp_new,im*jm)

      call BACK2GRD (conv_new,pblindex,conv_new,im*jm)

      call BACK2GRD (snow_new,pblindex,snow_new,im*jm)

      call BACK2GRD (raincgath,pblindex,raincgath,im*jm)

      call BACK2GRD (psubcldg ,pblindex,psubcldg ,im*jm)

      call BACK2GRD (psubcldgc,pblindex,psubcldgc,im*jm)

       call BACK2GRD (deltgather(1,L),pblindex, dtmoist(1,1,L)  ,im*jm)

       call BACK2GRD (delqgather(1,L),pblindex, dqmoist(1,1,L,1),im*jm)

       call BACK2GRD (    cpen(1,1,L),pblindex,    cpen(1,1,L)  ,im*jm)

       call BACK2GRD (   cldls(1,1,L),pblindex,   cldls(1,1,L)  ,im*jm)

       call BACK2GRD (cldwater(1,1,L),pblindex,cldwater(1,1,L)  ,im*jm)

       call BACK2GRD ( pkzgather(1,L),pblindex, pkzgather(1,L)  ,im*jm)

       call BACK2GRD(tmpgather(1,L),pblindex,

       call BACK2GRD(pkegather(1,L),pblindex,

       call BACK2GRD(deltrnev(1,L),pblindex,

       call BACK2GRD(delqrnev(1,L),pblindex,

       call BACK2GRD(cldsr(1,1,L),pblindex,

       call BACK2GRD (delugather(1,L,nt),pblindex,

       call BACK2GRD (delugather(1,L,ntracerin-ptracer+1),pblindex,

       call BACK2GRD (delugather(1,L,ntracerin-ptracer+2),pblindex,

      integer  pblindex(im*jm)