cleanup, bugfix

scripts/apply_obs_op_dart.py

@@ -32,10 +32,7 @@ if __name__ == '__main__':
 
     savedir = cluster.archivedir+'/obs_seq_final_1min/'
 
-    first_obscfg = exp.observations[0]
-    aso.set_DART_nml(cov_loc_radius_km=first_obscfg['cov_loc_radius_km'],
-                        cov_loc_vert_km=first_obscfg.get('cov_loc_vert_km', False), 
-                        just_prior_values=True)
+    aso.set_DART_nml(just_prior_values=True)
     osq.create_obsseqin_alltypes(time, exp.observations, obs_errors=False)
 
     # prepare dummy nature (this Hx is irrelevant)
@@ -49,33 +46,3 @@ if __name__ == '__main__':
     # only the prior state values are of interest in this file
     # observation and truth is wrong in this file (dummy)
     aso.archive_osq_final(time, posterior_1min=True)
-
-    # sys.exit()  # multi stage below
-    # n_stages = len(exp.observations)
-    # for istage, obscfg in enumerate(exp.observations):
-
-    #     n_obs = obscfg['n_obs']
-    #     sat_channel = obscfg.get('sat_channel', False)
-    #     obscfg['folder_obs_coords'] = False
-
-    #     aso.set_DART_nml(cov_loc_radius_km=obscfg['cov_loc_radius_km'],
-    #                      cov_loc_vert_km=obscfg.get('cov_loc_vert_km', False), 
-    #                      just_prior_values=True)
-
-    #     osq.create_obsseq_in(time, obscfg)
-
-
-
-    #     # prepare dummy nature (this Hx is irrelevant)
-    #     os.chdir(cluster.dartrundir)
-    #     os.system('cp ./advance_temp1/wrfout_d01 ./wrfout_d01')
-    #     wrfout_add_geo.run(cluster.dartrundir+'/geo_em.d01.nc',
-    #                        cluster.dartrundir+'/wrfout_d01')
-    #     aso.run_perfect_model_obs()
-    #     aso.assimilate(nproc=96)
-
-    #     # only the prior state values are of interest in this file
-    #     # observation and truth is wrong in this file (dummy)
-    #     archive_stage = savedir+'/assim_stage'+str(istage)
-    #     aso.archive_diagnostics(archive_stage, time)
-