anchor.diff:                    dzdx                  => this%dzdx,                           &
                     dzdy                  => this%dzdy,                           &
                     dz_scl                  => this%dz_scl,                           &
  -                  max_level             => this%max_level                       &
  +                  max_level             => this%max_level,                       &
                     z_level_ratio         => this%z_level_ratio,                  &
                     zr_u                  => this%zr_u,                           &
                     zr_v                  => this%zr_v)
  
positive.code:                                     this% kms : this% kme, &
                                      this% jms : this% jme) )
  
              print*,"maxval(forcing_terrain): " maxval(forcing_terrain)     
              do i = this%grid%kms, this%grid%kme
                
                if (i<=max_level) then
  
negative.code:             call io_write("delta_dzdx_sc.nc", "delta_dzdx", delta_dzdx(:,:,:) )
              call io_write("delta_dzdx_lc.nc", "delta_dzdx_lc", delta_dzdx_lc(:,:,:) )
              print*, "sinh term ", sin_h(:) 
              call io_write("dzdx.nc", "dzdx", dzdx(:,:,:) )
              call io_write("dzfdx.nc", "dzfdx", dzfdx(:,:,:) )
              call io_write("zf.nc", "zf", zf(:,:,:) )
  
      
          !_________ 2. Calculate the ratio bewteen z levels from hi-res and forcing data for wind acceleration  _________                                
  