Skip to content
Snippets Groups Projects
Select Git revision
  • consistent_config default protected
1 result
Blame Permalink
calc_linear_posterior.py 2.71 KiB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72 






import os, sys, shutil, warnings, glob
import datetime as dt

from dartwrf.exp_config import exp
from dartwrf.server_config import cluster
from dartwrf import assimilate as aso
from dartwrf import dart_nml
from dartwrf.obs import create_obsseq_out as osq_out

pattern_osf_linear = cluster.archivedir + "/diagnostics/%Y-%m-%d_%H:%M_obs_seq.final-linear"  # how an obs_seq.final file is archived

def calc_lin_posterior(assim_time):
     """Run filter to get the linear, optimal-interpolation posterior in obs-space

     Args:
          prior_init_time (datetime): initialization of prior forecast
          assim_time (datetime): assimilation time

     Returns:
          None
     """
     # what is the prior for this assim time?
     # read link_to_prior
     link_to_prior = cluster.archivedir + assim_time.strftime('/%Y-%m-%d_%H:%M/')+'link_to_prior.txt'
     with open(link_to_prior) as f:
          lines = [l.strip() for l in f.readlines()]
          prior_exp = lines[0]
          prior_init_time = dt.datetime.strptime(lines[1].strip('/'), '%Y-%m-%d_%H:%M')
     
     aso.prepare_prior_ensemble(assim_time, prior_init_time=prior_init_time, prior_valid_time=assim_time, prior_path_exp=prior_exp)
     #aso.prepare_nature_dart(assim_time) 
     
     dart_nml.write_namelist()

     # does an observation exist at this time?
     #f_oso = assim_time.strftime(aso.pattern_obs_seq_out)
     f_oso = cluster.archivedir+assim_time.strftime("/diagnostics/%Y-%m-%d_%H:%M_obs_seq.out")

     if os.path.exists(f_oso):
          # use the existing file
          shutil.copy(f_oso, cluster.dart_rundir+'/obs_seq.out')
     else:
          raise NotImplementedError(f_oso+' does not exist!')
     
     # link in the linear filter version
     from dartwrf.utils import symlink
     path_linear_filter = '/jetfs/home/lkugler/data/compile/DART/DART-10.8.3-kalmanposterior/models/wrf/work/filter'
     symlink(path_linear_filter,
             os.path.join(cluster.dart_rundir, 'filter'))

     aso.filter(nproc=12)
     aso.archive_filter_diagnostics(assim_time, pattern_osf_linear)


if __name__ == "__main__":
     """Evaluate the ensemble forecast in observation space at a given time, apart from the analysis time.

     Note: Observations are not assimilated. This is only for evaluation purposes.

     Usage: python3 evaluate_obs_space.py init1,valid1 init2,valid2 ...
     """
     args = sys.argv[1:]
     #arg_tuples = [a.split(',') for a in args]

     # we need an existing run_DART folder
     aso.prepare_run_DART_folder()

     for assim_time in args:
          #prior_init_time = dt.datetime.strptime(prior_init_time, "%Y-%m-%d_%H:%M")
          assim_time = dt.datetime.strptime(assim_time, "%Y-%m-%d_%H:%M")
          calc_lin_posterior(assim_time)