common.py

import os

if os.environ['USER'] == 'ksargsy' or os.environ['USER'] == 'root':
    print('Hello Khachik')
    oscm_dir = os.environ['HOME'] + '/research/OSCM_SciDAC/'
elif os.environ['USER'] == 'csafta':
    print('Hello Cosmin')
    oscm_dir = os.environ['HOME'] + '/Projects/OSCM_SciDAC.dmr/'
else:
    oscm_dir = '../../'


pdict = {'gdd_crit': 500.0, 'crit_dayl': 39300., 'ndays_on': 30, 'ndays_off': 15,
         'nue_tree': 15.0, 'nue_grass': 9.0, 'slatop_everg': 0.01, 'slatop_decid': 0.03,
         'livewdcn': 50, 'leafcn': 25, 'frootcn': 42,
         'fstor2tran': 0.5, 'stem_leaf': 2.7, 'croot_stem': 0.3, 'f_livewd': 0.1,
         'froot_leaf': 1.0,
         'rg_frac': 0.3, 'br_mr': 2.52e-6, 'q10_mr': 1.5, 'cstor_tau': 3.0,
         'r_mort': 0.02, 'lwtop_ann': 0.7, 'leaf_long': 3.0, 'froot_long': 3.0,
         'q10_hr': 1.5, 'k_l1': 1.2039728, 'k_l2': 0.0725707, 'k_l3': 0.0140989,
         'k_s1': 0.0725707, 'k_s2': 0.0140989244, 'k_s3': 0.00140098, 'k_s4': 0.0001,
         'k_frag': 0.0010005, 'rf_l1s1': 0.39, 'rf_l2s2': 0.55, 'rf_l3s3': 0.29,
         'rf_s1s2': 0.28, 'rf_s2s3': 0.46, 'rf_s3s4': 0.55, 'soil4ci': 1000.,
         'cwd_flig': 0.24, 'fr_flig': 0.25, 'lf_flig': 0.25, 'fr_flab': 0.25, 'lf_flab': 0.25,
         'fpi': 0.1, 'fpg': 0.9}


pmin = {}
pmax = {}

for p in pdict:
    if (p == 'crit_dayl'):
        pmin[p] = 36000.
        pmax[p] = 43000.
    elif (p == 'gdd_crit'):
        pmin[p] = 100.0
        pmax[p] = 700.0
    elif (p == 'fpg'):
        pmin[p] = 0.7  # 0.50
        pmax[p] = 0.95  # 1.00
    else:
        pmin[p] = pdict[p] * 0.5
        pmax[p] = pdict[p] * 1.5


dim = 47


outvars = ['gpp', 'npp', 'gr', 'mr', 'hr', 'nee', 'lai',
           'leafc', 'leafc_stor', 'frootc', 'frootc_stor', 'livestemc',
           'deadstemc', 'livecrootc', 'deadcrootc', 'ctcpools', 'totecosysc',
           'totsomc', 'totlitc', 'cstor']