anchor.diff:      
       call model_init(options)
       
  +    ! read in the atmospheric predictor data (e.g. GCM or GEFS)
  +    print*, ""
  +    print*, "=========================================="
  +    print*, "Reading predictor"
  +    predictions  = read_atm(options%prediction)
  +    
       ! read in the training atmospheric data (e.g. reanalysis or GEFS)
       print*, ""
       print*, "=========================================="
  
positive.code:     call geo_LUT(training_obs, training_atm)
      call geo_LUT(training_obs, predictions)
      
      call io_write("obs.nc","data",training_obs%variables(1)%data)
      call io_write("training.nc","data",training_atm%variables(1)%data)
      call io_write("predictor.nc","data",predictions%variables(1)%data)
      
      print*, ""
      print*, "=========================================="
  
negative.code:     print*, "=========================================="
      print*, "Reading obs"
      training_obs = read_obs(options%obs)
      
      print*, ""
      print*, "=========================================="
  