anchor.diff:          end associate
           associate(var=>var_meta(kVARS%water_vapor))
               var%name        = "qv"
  -            var%dimensions  = [character(len=16) :: "lon_x","lat_x","level"]
  +            var%dimensions  = three_d_dimensions
               var%attributes  = [attribute_t("standard_name", "mass_fraction_of_water_vapor_in_air"), &
                                  attribute_t("long_name",     "Water Vapor Mixing Ratio"),            &
                                  attribute_t("units",         "kg kg-1"),                             &
  
positive.code:         end associate
          associate(var=>var_meta(kVARS%pressure))
              var%name        = "pressure"
              var%dimensions  = [character(len=16) :: "lon_x","lat_x","level"]
              var%attributes  = [attribute_t("standard_name", "air_pressure"),                    &
                                 attribute_t("long_name",     "Pressure"),                        &
                                 attribute_t("units",         "Pa"),                              &
  
negative.code:         end associate
          associate(var=>var_meta(kVARS%cloud_water))
              var%name        = "qc"
              var%dimensions  = three_d_dimensions
              var%attributes  = [attribute_t("standard_name", "cloud_liquid_water_mixing_ratio"),     &
                                 attribute_t("units",         "kg kg-1"),                             &
                                 attribute_t("coordinates",   "lat lon")]
          end associate
          associate(var=>var_meta(kVARS%cloud_number_concentration))
              var%name        = "nc"
              var%dimensions  = three_d_dimensions
              var%attributes  = [attribute_t("non_standard_name", "number_concentration_of_cloud_droplets_in_air"), &
                                 attribute_t("units",         "cm-3"),                                              &
                                 attribute_t("coordinates",   "lat lon")]
          end associate
          associate(var=>var_meta(kVARS%cloud_ice))
              var%name        = "qi"
              var%dimensions  = three_d_dimensions
              var%attributes  = [attribute_t("standard_name", "cloud_ice_mixing_ratio"),              &
                                 attribute_t("units",         "kg kg-1"),                             &
                                 attribute_t("coordinates",   "lat lon")]
          end associate
          associate(var=>var_meta(kVARS%ice_number_concentration))
              var%name        = "ni"
              var%dimensions  = three_d_dimensions
              var%attributes  = [attribute_t("non_standard_name", "number_concentration_of_ice_crystals_in_air"), &
                                 attribute_t("units",         "cm-3"),                                            &
                                 attribute_t("coordinates",   "lat lon")]
          end associate
          associate(var=>var_meta(kVARS%rain_in_air))
              var%name        = "qr"
              var%dimensions  = three_d_dimensions
              var%attributes  = [attribute_t("standard_name", "mass_fraction_of_rain_in_air"),        &
                                 attribute_t("units",         "kg kg-1"),                             &
                                 attribute_t("coordinates",   "lat lon")]
          end associate
          associate(var=>var_meta(kVARS%rain_number_concentration))
              var%name        = "nr"
              var%dimensions  = three_d_dimensions
              var%attributes  = [attribute_t("non_standard_name", "number_concentration_of_rain_particles_in_air"), &
                                 attribute_t("units",         "cm-3"),                                              &
                                 attribute_t("coordinates",   "lat lon")]
          end associate
          associate(var=>var_meta(kVARS%snow_in_air))
              var%name        = "qs"
              var%dimensions  = three_d_dimensions
              var%attributes  = [attribute_t("standard_name", "mass_fraction_of_snow_in_air"),        &
                                 attribute_t("units",         "kg kg-1"),                             &
                                 attribute_t("coordinates",   "lat lon")]
          end associate
          associate(var=>var_meta(kVARS%snow_number_concentration))
              var%name        = "ns"
              var%dimensions  = three_d_dimensions
              var%attributes  = [attribute_t("non_standard_name", "number_concentration_of_snow_particles_in_air"), &
                                 attribute_t("units",         "cm-3"),                                              &
                                 attribute_t("coordinates",   "lat lon")]
          end associate
          associate(var=>var_meta(kVARS%graupel_in_air))
              var%name        = "qg"
              var%dimensions  = three_d_dimensions
              var%attributes  = [attribute_t("standard_name", "mass_fraction_of_graupel_in_air"),     &
                                 attribute_t("units",         "kg kg-1"),                             &
                                 attribute_t("coordinates",   "lat lon")]
          end associate
          associate(var=>var_meta(kVARS%graupel_number_concentration))
              var%name        = "ng"
              var%dimensions  = three_d_dimensions
              var%attributes  = [attribute_t("non_standard_name", "number_concentration_of_graupel_particles_in_air"), &
                                 attribute_t("units",         "cm-3"),                                                 &
                                 attribute_t("coordinates",   "lat lon")]
          end associate
          ! associate(var=>var_meta(kVARS%precip_rate))
          !     var%name        = "precip_rate"
          !     var%dimensions  = two_d_dimensions
          !     var%attributes  = [attribute_t("standard_name",   "precipitation_flux"),      &
          !                        attribute_t("units",           "kg m-2 s-1"),              &
          !                        attribute_t("coordinates",     "lat lon")]
          ! end associate
          associate(var=>var_meta(kVARS%precipitation))
              var%name        = "precipitation"
              var%dimensions  = two_d_dimensions
              var%attributes  = [attribute_t("standard_name", "precipitation_amount"),                &
                                 attribute_t("units",         "kg m-2"),                              &
                                 attribute_t("coordinates",   "lat lon")]
          end associate
          associate(var=>var_meta(kVARS%snowfall))
              var%name        = "snowfall"
              var%dimensions  = two_d_dimensions
              var%attributes  = [attribute_t("standard_name", "snowfall_amount"),                     &
                                 attribute_t("units",         "kg m-2"),                              &
                                 attribute_t("coordinates",   "lat lon")]
          end associate
          associate(var=>var_meta(kVARS%graupel))
              var%name        = "graupel"
              var%dimensions  = two_d_dimensions
              var%attributes  = [attribute_t("standard_name", "graupel_amount"),                      &
                                 attribute_t("units",         "kg m-2"),                              &
                                 attribute_t("coordinates",   "lat lon")]
          end associate
          associate(var=>var_meta(kVARS%exner))
              var%name        = "exner"
              var%dimensions  = three_d_dimensions
              var%attributes  = [attribute_t("non_standard_name", "exner_function_result"),           &
                                 attribute_t("units",         "K K-1"),                               &
                                 attribute_t("coordinates",   "lat lon")]
          end associate
          associate(var=>var_meta(kVARS%density))
              var%name        = "density"
              var%dimensions  = three_d_dimensions
              var%attributes  = [attribute_t("standard_name", "air_density"),                         &
                                 attribute_t("units",         "kg m-3"),                              &
                                 attribute_t("coordinates",   "lat lon")]
          end associate
          associate(var=>var_meta(kVARS%z))
              var%name        = "z"
              var%dimensions  = three_d_dimensions
              var%attributes  = [attribute_t("standard_name", "height_above_reference_ellipsoid"),    &
                                 attribute_t("units",         "m"),                                   &
                                 attribute_t("coordinates",   "lat lon")]
          end associate
          associate(var=>var_meta(kVARS%z_interface))
              var%name        = "z_i"
              var%dimensions  = three_d_interface_dimensions
              var%attributes  = [attribute_t("standard_name", "height_above_reference_ellipsoid"),    &
                                 attribute_t("units",         "m"),                                   &
                                 attribute_t("coordinates",   "lat lon")]
          end associate
          associate(var=>var_meta(kVARS%dz))
              var%name        = "dz"
              var%dimensions  = three_d_dimensions
              var%attributes  = [attribute_t("non_standard_name", "layer_thickness"),                 &
                                 attribute_t("units",         "m"),                                   &
                                 attribute_t("coordinates",   "lat lon")]
          end associate
          associate(var=>var_meta(kVARS%dz_interface))
              var%name        = "dz_i"
              var%dimensions  = three_d_interface_dimensions
              var%attributes  = [attribute_t("non_standard_name", "layer_thickness"),                 &
                                 attribute_t("units",         "m"),                                   &
                                 attribute_t("coordinates",   "lat lon")]
          end associate
          associate(var=>var_meta(kVARS%shortwave))
              var%name        = "rsds"
              var%dimensions  = two_d_dimensions
              var%attributes  = [attribute_t("standard_name", "surface_downwelling_shortwave_flux_in_air"), &
                                 attribute_t("units",         "W m-2"),                                     &
                                 attribute_t("coordinates",   "lat lon")]
          end associate
          associate(var=>var_meta(kVARS%longwave))
              var%name        = "rlds"
              var%dimensions  = two_d_dimensions
              var%attributes  = [attribute_t("standard_name", "surface_downwelling_longwave_flux_in_air"), &
                                 attribute_t("units",         "W m-2"),                                    &
                                 attribute_t("coordinates",   "lat lon")]
          end associate
          associate(var=>var_meta(kVARS%longwave_up))
              var%name        = "rlus"
              var%dimensions  = two_d_dimensions
              var%attributes  = [attribute_t("standard_name", "surface_upwelling_longwave_flux_in_air"),   &
                                 attribute_t("units",         "W m-2"),                                    &
                                 attribute_t("coordinates",   "lat lon")]
          end associate
          associate(var=>var_meta(kVARS%ground_heat_flux))
              var%name        = "hfgs"
              var%dimensions  = two_d_dimensions
              var%attributes  = [attribute_t("standard_name", "upward_heat_flux_at_ground_level_in_soil"), &
                                 attribute_t("units",         "W m-2"),                                    &
                                 attribute_t("coordinates",   "lat lon")]
          end associate
          associate(var=>var_meta(kVARS%vegetation_fraction))
              var%name        = "vegetation_fraction"
              var%dimensions  = two_d_dimensions
              var%attributes  = [attribute_t("standard_name", "vegetation_fraction"),                 &
                                 attribute_t("units",         "m2 m-2"),                              &
                                 attribute_t("coordinates",   "lat lon")]
          end associate
          associate(var=>var_meta(kVARS%land_cover))
              var%name        = "land_cover"
              var%dimensions  = two_d_dimensions
              var%attributes  = [attribute_t("non_standard_name", "land_cover_type"),                 &
                                 attribute_t("units",      ""),                                       &
                                 attribute_t("coordinates",   "lat lon")]
          end associate
          associate(var=>var_meta(kVARS%lai))
              var%name        = "lai"
              var%dimensions  = two_d_dimensions
              var%attributes  = [attribute_t("non_standard_name", "leaf_area_index"),                 &
                                 attribute_t("units",         "m2 m-2"),                              &
                                 attribute_t("coordinates",   "lat lon")]
          end associate
          associate(var=>var_meta(kVARS%canopy_water))
              var%name        = "canopy_water"
              var%dimensions  = two_d_dimensions
              var%attributes  = [attribute_t("standard_name", "canopy_water_amount"),                 &
                                 attribute_t("units",         "kg m-2"),                              &
                                 attribute_t("coordinates",   "lat lon")]
          end associate
          associate(var=>var_meta(kVARS%snow_water_equivalent))
              var%name        = "swe"
              var%dimensions  = two_d_dimensions
              var%attributes  = [attribute_t("standard_name", "surface_snow_amount"),                 &
                                 attribute_t("units",         "kg m-2"),                              &
                                 attribute_t("coordinates",   "lat lon")]
          end associate
          associate(var=>var_meta(kVARS%soil_water_content))
              var%name        = "soil_water_content"
              var%dimensions  = three_d_soil_dimensions
              var%attributes  = [attribute_t("standard_name", "soil_moisture_content"),               &
                                 attribute_t("units",         "kg m-2"),                              &
                                 attribute_t("coordinates",   "lat lon")]
          end associate
          associate(var=>var_meta(kVARS%soil_temperature))
              var%name        = "soil_temperature"
              var%dimensions  = three_d_soil_dimensions
              var%attributes  = [attribute_t("standard_name", "soil_temperature"),                    &
                                 attribute_t("units",         "K"),                                   &
                                 attribute_t("coordinates",   "lat lon")]
          end associate
          associate(var=>var_meta(kVARS%air2m_temperature))
              var%name        = "ta2m"
              var%dimensions  = two_d_dimensions
              var%attributes  = [attribute_t("standard_name", "air_temperature"),                     &
                                 attribute_t("long_name",     "Bulk air temperature at 2m"),          &
                                 attribute_t("units",         "K"),                                   &
  