anchor.diff:  
   
               ! write and compare 
  -            call io_write("delta_dzdx.nc", delta_dzdx, delta_dzdx(:,:,:) )
  -            call io_write("delta_dzdx_2.nc", delta_dzdx_2, delta_dzdx_2(:,:,:) )
  +            call io_write("delta_dzdx.nc", "delta_dzdx", delta_dzdx(:,:,:) )
  +            call io_write("delta_dzdx_2.nc", "delta_dzdx_2", delta_dzdx_2(:,:,:) )
   
               ! ! if the latter works, put this call at the end of current module
               ! if (use_delta_terrain) then
  
positive.code: 
              ! #----------------------- calc z levels from forcing terrain -------------------
              print* , " hgt variable forcing_terrain: ", shape(this%forcing_terrain%data_2d)
              call io_write("forcing_terrain.nc", forcing_terrain, forcing_terrain(:,:) ) ! check in plot
  
              ! Allodcate al of the 'f' variables below? required, or just faster?
              ! allocate( zf_interface(this% ims : this% ime, &
                                      ! this% kms : this% kme, &
                                      ! this% jms : this% jme) )
  
              do i = this%grid%kms, this%grid%kme
                
  
negative.code:                                           this% kds : this% kde,   &
                                            this% jds : this% jde)   )
  
  
  
      end subroutine allocate_z_arrays
  