anchor.diff:  c**** THIS ROUTINE WRITES THE MATRIX AND VECTOR FOR PRESSURE RETRIEVAL
   c
         subroutine ajusp(pa,pp,ax,ay,az,x,v,p,cp0,a,b,c,d
  -     &                 ,tetav0,rho0,ip,n1,n2,n3)
  +     &                 ,tetav0,rho0,ip,n1,n2,n3,n3p1)
   c
         real ax(n1,n2,n3),ay(n1,n2,n3),az(n1,n2,n3)
        &    ,pa(n1,n2,n3),pp(n1,n2,n3)
        &    ,x(27,n1,n2,n3),v(n1,n2,n3),p(n1,n2,n3)
        &    ,a(n1,n2,n3),b(n1,n2,n3),c(n1,n2,n3),d(n1,n2,n3)
  -     &    ,tetav0(40),rho0(40)
  +     &    ,tetav0(n3p1),rho0(n3p1)
         integer ip(n1,n2,n3,7)
   c
         common/grid/hx,hy,hz,nx,ny,nz,i2d
  
positive.code:      +           fnu,fnv,fnw,fnz,addu,addv,hfac,itv12,xmv1,ymv1,
       +           xmv2,ymv2,icenter)
  c    NOTE! Modifying the parameter of dimension if the domain is changed.
        parameter (n1=70,n2=70,n3=36,n3p1=n3+1)
  c
        integer   iday(6),ih(6),jh(6),it(n1,n3,14),kk(n3)
        real      u(n1,n2,n3),v(n1,n2,n3),w(n1,n2,n3)
  
negative.code: c
  c  for ced file
  c
        if(icall.ge.0)then
           ibad=id(67)
           nprec=id(96)
           nz=id(106)
  