anchor.diff:                      do i=1,size(this%albedo%data_3d, 2)
                           this%albedo%data_3d(:,i,:) = temporary_data_3d(this%grid%ims:this%grid%ime, this%grid%jms:this%grid%jme,i)
                       enddo
  +
  +                    if (maxval(temporary_data_3d) > 1) then
  +                        if (this_image()==1) print*, "Changing input ALBEDO % to fraction"
  +                        this%albedo%data_3d = this%albedo%data_3d / 100
  +                    endif
                   endif
  -
  -                if (maxval(temporary_data_3d) > 1) then
  -                    if (this_image()==1) print*, "Changing input ALBEDO % to fraction"
  -                    this%albedo%data_3d = this%albedo%data_3d / 100
  -                endif
  -
               else
                   call io_read(options%parameters%init_conditions_file,   &
                                  options%parameters%albedo_var,          &
  
positive.code:         var_to_update%data_3d = var_to_update%data_3d + (var_to_update%dqdt_3d * dt%seconds())
  
          if (associated(this%external_precipitation%data_2d)) then
              this%accumulated_precipitation%data_2d = this%accumulated_precipitation%data_2d + (this%external_precipitation%data_2d * dt%seconds())
          endif
  
  
  
negative.code:                            temporary_data)
              if (allocated(this%land_mask)) then
                  this%land_mask = temporary_data(this%grid%ims:this%grid%ime, this%grid%jms:this%grid%jme)
                  where(this%land_mask==0) this%land_mask = kLC_WATER
              endif
          endif
  
  