anchor.diff:                  call subset_array(temp_offset, this%geo_v%lat, g)
               end associate
           endif
  -
  -        if (this_image()==1) print*, "WARNING Soil variables not properly initialized"
  -        if (associated(this%vegetation_fraction%data_3d)) this%vegetation_fraction%data_3d = 0.6
  -        if (associated(this%soil_water_content%data_3d)) this%soil_water_content%data_3d = 0.2
  -        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:         if (0<opt%vars_to_allocate( kVARS%canopy_water) )               call setup(this%canopy_water,             this%grid2d )
          if (0<opt%vars_to_allocate( kVARS%snow_water_equivalent) )      call setup(this%snow_water_equivalent,    this%grid2d )
          if (0<opt%vars_to_allocate( kVARS%skin_temperature) )           call setup(this%skin_temperature,         this%grid2d,   forcing_var=opt%parameters%sst_var,     list=this%variables_to_force)
          if (0<opt%vars_to_allocate( kVARS%soil_water_content) )         call setup(this%soil_water_content,       this%grid_soil,forcing_var=opt%parameters%soil_t_var,  list=this%variables_to_force)
          if (0<opt%vars_to_allocate( kVARS%soil_temperature) )           call setup(this%soil_temperature,         this%grid_soil,forcing_var=opt%parameters%soil_vwc_var,list=this%variables_to_force)
          if (0<opt%vars_to_allocate( kVARS%latitude) )                   call setup(this%latitude,                 this%grid2d)
          if (0<opt%vars_to_allocate( kVARS%longitude) )                  call setup(this%longitude,                this%grid2d)
          if (0<opt%vars_to_allocate( kVARS%u_latitude) )                 call setup(this%u_latitude,               this%u_grid2d)
  
negative.code:         call create_variables(this, options)
  
          call initialize_core_variables(this, options)
  
          call read_land_variables(this, options)
  
          call setup_meta_data(this, options)
  
  