anchor.diff:          character(len=16) :: two_d_t_dimensions(3)              = [character(len=16) :: "lon_x","lat_y","time"]
           character(len=16) :: two_d_u_dimensions(2)              = [character(len=16) :: "lon_u","lat_y"]
           character(len=16) :: two_d_v_dimensions(2)              = [character(len=16) :: "lon_x","lat_v"]
  -        character(len=16) :: three_d_soil_dimensions(3)         = [character(len=16) :: "lon_x","lat_y","nsoil"]
           character(len=16) :: three_d_t_soil_dimensions(4)       = [character(len=16) :: "lon_x","lat_y","nsoil","time"]
  +        character(len=16) :: three_d_t_snow_dimensions(4)       = [character(len=16) :: "lon_x","lat_y","nsnow","time"]
  +        character(len=16) :: three_d_t_snowsoil_dimensions(4)   = [character(len=16) :: "lon_x","lat_y","nsnowsoil","time"]
  +        character(len=16) :: three_d_soilcomp_dimensions(3)     = [character(len=16) :: "lon_x","lat_y","nsoil_composition"]
  +        character(len=16) :: three_d_crop_dimensions(3)         = [character(len=16) :: "lon_x","lat_y","crop"]
  +        character(len=16) :: three_d_t_gecros_dimensions(4)     = [character(len=16) :: "lon_x","lat_y","gecros","time"]
  +        character(len=16) :: two_d_month_dimensions(3)          = [character(len=16) :: "lon_x","lat_y","month"]
   
           if (allocated(var_meta)) deallocate(var_meta)
   
  
positive.code:         !!------------------------------------------------------------
          associate(var=>var_meta(kVARS%snow_layer_liquid_water))
              var%name        = "snow_layer_liquid_water"
              var%dimensions  = three_d_t_dimensions
              var%unlimited_dim=.True.
              var%attributes  = [attribute_t("non_standard_name", "snow_layer_liquid_water_content"), &
                                 attribute_t("units",         "mm"),                                  &
  
negative.code:         !!------------------------------------------------------------
          associate(var=>var_meta(kVARS%soil_sand_and_clay))
              var%name        = "soil_sand_and_clay_composition"
              var%dimensions  = three_d_soilcomp_dimensions
              var%attributes  = [attribute_t("non_standard_name", "soil_sand_and_clay_composition"),    &
                                 attribute_t("units",         "1"),                                     &
                                 attribute_t("coordinates",   "lat lon")]
  