In [1]:

import numpy as np
from matplotlib import pyplot as plt
import netCDF4 as nc
import re
import shutil
%matplotlib inline

In [2]:

src='/data/eolson/results/MEOPAR/tmp/SalishSea_03032640_restart_trc.nc'

In [3]:

f=nc.Dataset(src,'r+')

In [5]:

print(f.variables.keys())
for var in f.variables.keys():
    if 'MYRI' in var:
        print(var)

dict_keys(['nav_lon', 'nav_lat', 'nav_lev', 'time_counter', 'kt', 'ndastp', 'adatrj', 'rnf_pis_NO3_b', 'rnf_pis_NH4_b', 'rnf_pis_Si_b', 'rnf_pis_DIAT_b', 'rnf_pis_PHY_b', 'rnf_pis_MYRI_b', 'rnf_pis_MICZ_b', 'rnf_pis_DON_b', 'rnf_pis_PON_b', 'rnf_pis_bSi_b', 'rnf_pis_TRA_b', 'rnf_pis_DIC_b', 'rnf_pis_TA_b', 'rnf_pis_O2_b', 'sbc_NO3_b', 'sbc_NH4_b', 'sbc_Si_b', 'sbc_DIAT_b', 'sbc_PHY_b', 'sbc_MYRI_b', 'sbc_MICZ_b', 'sbc_DON_b', 'sbc_PON_b', 'sbc_bSi_b', 'sbc_TRA_b', 'sbc_DIC_b', 'sbc_TA_b', 'sbc_O2_b', 'rdttrc1', 'TRNNO3', 'TRNNH4', 'TRNSi', 'TRNDIAT', 'TRNPHY', 'TRNMYRI', 'TRNMICZ', 'TRNDON', 'TRNPON', 'TRNbSi', 'TRNTRA', 'TRNDIC', 'TRNTA', 'TRNO2', 'TRBNO3', 'TRBNH4', 'TRBSi', 'TRBDIAT', 'TRBPHY', 'TRBMYRI', 'TRBMICZ', 'TRBDON', 'TRBPON', 'TRBbSi', 'TRBTRA', 'TRBDIC', 'TRBTA', 'TRBO2'])
rnf_pis_MYRI_b
sbc_MYRI_b
TRNMYRI
TRBMYRI

In [ ]:

for var in ('rnf_pis_MYRI_b','sbc_MYRI_b','TRNMYRI','TRBMYRI'):
    f.variables[var][:]=0.0

In [ ]:

for var in ('rnf_pis_MYRI_b','sbc_MYRI_b','TRNMYRI','TRBMYRI'):
    f.variables[var][:]=0.0

In [ ]:

f.close()

In [ ]:

f=nc.Dataset(dst)

In [ ]:

np.max(f.variables['TRNPHY']), np.min(f.variables['TRNPHY'])

In [ ]:

f.close()

In [ ]: