anchor.diff:          if (associated(this%soil_totalmoisture%data_2d)) this%soil_totalmoisture%data_2d = 0.2 * 2
           if (associated(this%soil_temperature%data_3d)) this%soil_temperature%data_3d = 280
           if (associated(this%soil_deep_temperature%data_2d)) this%soil_deep_temperature%data_2d = 280
  +        if (associated(this%skin_temperature%data_2d)) this%skin_temperature%data_2d = 280
           if (associated(this%roughness_z0%data_2d)) this%roughness_z0%data_2d = 0.001
  +        if (associated(this%sensible_heat%data_2d)) this%sensible_heat%data_2d=0
  +        if (associated(this%latent_heat%data_2d)) this%latent_heat%data_2d=0
  +        if (associated(this%u_10m%data_2d)) this%u_10m%data_2d=0
  +        if (associated(this%v_10m%data_2d)) this%v_10m%data_2d=0
  +        if (associated(this%temperature_2m%data_2d)) this%temperature_2m%data_2d=280
  +        if (associated(this%humidity_2m%data_2d)) this%humidity_2m%data_2d=0.001
  +        if (associated(this%surface_pressure%data_2d)) this%surface_pressure%data_2d=102000
  +        if (associated(this%longwave_up%data_2d)) this%longwave_up%data_2d=0
  +        if (associated(this%ground_heat_flux%data_2d)) this%ground_heat_flux%data_2d=0
  +        if (associated(this%soil_totalmoisture%data_2d)) this%soil_totalmoisture%data_2d=0.5
   
           call standardize_coordinates(this%geo_u)
           call standardize_coordinates(this%geo_v)
  
positive.code: 
          call setup_geo(this%geo,   this%latitude%data_2d,   this%longitude%data_2d,   this%z%data_3d)
  
      end subroutine initialize_variables
  
      !> -------------------------------
      !! Populare the metadata structure in the domain for later output
  
negative.code:     !! Initialize various domain variables, mostly z, dz, etc.
      !!
      !! -------------------------------
      subroutine initialize_core_variables(this, options)
          implicit none
          class(domain_t), intent(inout)  :: this
          class(options_t),intent(in)     :: options
  