|
| 1 | +#!/usr/bin/env python3 |
| 2 | +from yaml import load,dump |
| 3 | +import argparse |
| 4 | +import os |
| 5 | +import yaml |
| 6 | +import glob |
| 7 | + |
| 8 | +dims_dict = {"2":"xy", "3":"xyz"} |
| 9 | + |
| 10 | +def get_dims(component_map, name): |
| 11 | + for comp in component_map: |
| 12 | + if name in component_map[comp]: |
| 13 | + dims = component_map[comp][name]["dims"] |
| 14 | + return dims |
| 15 | + |
| 16 | +def get_vars_needed(input_file): |
| 17 | + |
| 18 | + output_list = [] |
| 19 | + f = open(input_file,"r") |
| 20 | + lines = f.readlines() |
| 21 | + f.close() |
| 22 | + for line in lines: |
| 23 | + temp = line.split() |
| 24 | + output_list.append(temp[4]) |
| 25 | + |
| 26 | + return output_list |
| 27 | + |
| 28 | +def get_extdata_map(input_dir): |
| 29 | + |
| 30 | + input_files = glob.glob(input_dir+"/*.yaml") |
| 31 | + export_list = {} |
| 32 | + for input_file in input_files: |
| 33 | + f = open(input_file,'r') |
| 34 | + extdata_def = yaml.safe_load(f) |
| 35 | + f.close() |
| 36 | + for short_name in extdata_def["Exports"]: |
| 37 | + export_list.update({short_name:extdata_def["Exports"][short_name]}) |
| 38 | + return export_list |
| 39 | + |
| 40 | +def get_block(cvs_file): |
| 41 | + |
| 42 | + temp = cvs_file[0].split(' ') |
| 43 | + state_type = temp[1][1:].strip() |
| 44 | + component = temp[4].strip() |
| 45 | + i=2 |
| 46 | + for line in cvs_file[2:]: |
| 47 | + if "spec for" in line: |
| 48 | + break |
| 49 | + i=i+1 |
| 50 | + return component,state_type,i |
| 51 | + |
| 52 | + |
| 53 | +def get_component_map(cvs_file): |
| 54 | + |
| 55 | + i_start = 0 |
| 56 | + i_end = 0 |
| 57 | + n_lines = len(cvs_file) |
| 58 | + components = {} |
| 59 | + while i_start < n_lines-1: |
| 60 | + |
| 61 | + comp_name,state_type,i_end = get_block(cvs_file[i_start:]) |
| 62 | + comp_map = {} |
| 63 | + for i in range(i_end-2): |
| 64 | + line = cvs_file[i_start+2+i] |
| 65 | + values = line.split(',') |
| 66 | + short_name = values[1].strip() |
| 67 | + long_name = values[2].strip() |
| 68 | + units = values[3].strip() |
| 69 | + dims = values[4].strip() |
| 70 | + item_type = values[5].strip() |
| 71 | + comp_map.update({short_name:{"long_name":long_name,"units":units,"item_type":item_type,"dims":dims}}) |
| 72 | + components.update({comp_name+"_"+state_type:comp_map}) |
| 73 | + i_start = i_start + i_end |
| 74 | + |
| 75 | + return components |
| 76 | + |
| 77 | +def parse_args(): |
| 78 | + p = argparse.ArgumentParser(description='Generarte input files for ExtDataDriver to simulate GEOS') |
| 79 | + p.add_argument('extdata_provided',type=str,help='a list of items ExtData should fill',default=None) |
| 80 | + p.add_argument('spec_def',type=str,help='the GEOS gcm import state from the printspec',default=None) |
| 81 | + p.add_argument('extdata_dir',type=str,help='diretory with all the yaml imputs for extdata',default=None) |
| 82 | + p.add_argument('-e','--export',action='store_true',help='also include exports for corresponding imports') |
| 83 | + |
| 84 | + return vars(p.parse_args()) |
| 85 | + |
| 86 | +if __name__ == '__main__': |
| 87 | + |
| 88 | + args = parse_args() |
| 89 | + |
| 90 | + extdata_list = args['extdata_provided'] |
| 91 | + do_exports = args['export'] |
| 92 | + input_file = args['spec_def'] |
| 93 | + f = open(input_file,"r") |
| 94 | + input_rc = f.readlines() |
| 95 | + f.close() |
| 96 | + |
| 97 | + extdata_directory = args['extdata_dir'] |
| 98 | + extdata_def = get_extdata_map(extdata_directory) |
| 99 | + |
| 100 | + f_agcm = open("AGCM.rc",'w') |
| 101 | +# component |
| 102 | + component_map = {} |
| 103 | + component_map = get_component_map(input_rc) |
| 104 | + |
| 105 | + vars_needed = get_vars_needed(extdata_list) |
| 106 | + |
| 107 | + nl = "\n" |
| 108 | + cm = " , " |
| 109 | + |
| 110 | +# Import state |
| 111 | + written = [] |
| 112 | + f_agcm.write("IMPORT_STATE::"+nl) |
| 113 | + |
| 114 | + for item in vars_needed: |
| 115 | + if item in extdata_def: |
| 116 | + long_name = "NA" |
| 117 | + units = "NA" |
| 118 | + dims = get_dims(component_map, item) |
| 119 | + cdims = dims_dict[dims] |
| 120 | + |
| 121 | + if item not in written: |
| 122 | + f_agcm.write(item+cm+long_name+cm+units+cm+cdims+cm+"c"+nl) |
| 123 | + written.append(item) |
| 124 | + |
| 125 | + f_agcm.write("::"+nl) |
| 126 | + |
| 127 | +# Export state |
| 128 | + if do_exports: |
| 129 | + written = [] |
| 130 | + f_agcm.write("EXPORT_STATE::"+nl) |
| 131 | + for item in vars_needed: |
| 132 | + if item in extdata_def: |
| 133 | + long_name = "NA" |
| 134 | + units = "NA" |
| 135 | + dims = get_dims(component_map, item) |
| 136 | + cdims = dims_dict[dims] |
| 137 | + |
| 138 | + if item not in written: |
| 139 | + f_agcm.write(item+cm+long_name+cm+units+cm+cdims+cm+"c"+nl) |
| 140 | + written.append(item) |
| 141 | + |
| 142 | + f_agcm.write("::"+nl) |
| 143 | + |
| 144 | + f_hist = open("HISTORY.rc",'w') |
| 145 | + f_hist.write("GRID_LABELS:"+nl) |
| 146 | + f_hist.write("::"+nl) |
| 147 | + f_hist.write("COLLECTIONS: my_collection"+nl) |
| 148 | + f_hist.write("::"+nl) |
| 149 | + f_hist.write("my_collection.template: 'nc4'"+nl) |
| 150 | + f_hist.write("my_collection.format: 'CFIO'"+nl) |
| 151 | + f_hist.write("my_collection.frequency: '240000'"+nl) |
| 152 | + first = True |
| 153 | + written = [] |
| 154 | + for item in vars_needed: |
| 155 | + if item in extdata_def: |
| 156 | + if item not in written: |
| 157 | + if first: |
| 158 | + first = False |
| 159 | + f_hist.write("my_collection.fields:'"+item+"' , 'Root',"+"\n") |
| 160 | + else: |
| 161 | + f_hist.write("'"+item+"' , 'Root',"+"\n") |
| 162 | + written.append(item) |
| 163 | + f_hist.write("::") |
| 164 | + |
| 165 | + f_hist.close() |
| 166 | + f_agcm.close() |
| 167 | + |
0 commit comments