anchor.diff:          if (0<opt%vars_to_allocate( kVARS%terrain) )                    call setup(this%forcing_terrain,          this%grid2d,    forcing_var=opt%parameters%hgtvar, list=this%variables_to_force)
           if (0<opt%vars_to_allocate( kVARS%terrain) )                    call setup(this%terrain_u,                this%u_grid2d) 
           if (0<opt%vars_to_allocate( kVARS%terrain) )                    call setup(this%terrain_v,                this%v_grid2d)
  -        ! if (0<opt%vars_to_allocate( kVARS%terrain) )                    call setup(this%delta_terrain,            this%grid2d)  terrain_u, terrain_v,dz_scl, max_level, H               
           if (0<opt%vars_to_allocate( kVARS%sensible_heat) )              call setup(this%sensible_heat,            this%grid2d)
           if (0<opt%vars_to_allocate( kVARS%latent_heat) )                call setup(this%latent_heat,              this%grid2d)
           if (0<opt%vars_to_allocate( kVARS%u_10m) )                      call setup(this%u_10m,                    this%grid2d)
  
positive.code:         class(domain_t),  intent(inout) :: this
          type(options_t), intent(in)     :: options
  
          real, allocatable :: forcing_terrain_u(:,:), forcing_terrain_v(:,:), delta_terrain(:,:), delta_dzdx_lc(:,:,:)
          real, allocatable :: zf_interface(:,:,:), dzf_interface(:,:,:), zf(:,:,:), dzf_mass(:,:,:), dzfdx(:,:,:), dzfdy(:,:,:), delta_dzdx(:,:,:)
          real :: s
          integer :: i
  
  
negative.code: 
  
                  ! write and compare 
          ! if ((this_image()==1).and.(options%parameters%debug)) then  ! Print some diagnostics. USeful for development.         
            if ((this_image()==1)) then  ! Print some diagnostics. USeful for development.         
            call io_write("forcing_terrain.nc", "forcing_terrain", forcing_terrain(:,:) ) ! check in plot
            call io_write("terrain.nc", "terrain", terrain(:,:) ) ! check in plot        
            call io_write("delta_dzdx_sc.nc", "delta_dzdx_sc", delta_dzdx_sc(:,:,:) )
            ! call io_write("delta_dzdx_lc.nc", "delta_dzdx_lc", delta_dzdx(:,:,:) )
            call io_write("dzdx.nc", "dzdx", dzdx(:,:,:) )
            call io_write("zfr_u.nc", "zfr_u", zfr_u(:,:,:) ) ! check in plot
          endif
  