Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Download subset fix data with python script #3400

Open
wants to merge 29 commits into
base: develop
Choose a base branch
from
Open
Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
29 commits
Select commit Hold shift + click to select a range
f071994
change noaacloud NodeName
weihuang-jedi Feb 25, 2025
11accd8
Merge branch 'develop' of github.com:NOAA-EPIC/global-workflow-cloud …
weihuang-jedi Feb 26, 2025
a4771df
add download fix subset data, also serve as a test PR to trigger CI t…
weihuang-jedi Feb 26, 2025
4f1e679
fix pynorm error
weihuang-jedi Feb 26, 2025
2241227
fix pynorm error 2
weihuang-jedi Feb 26, 2025
9679617
fix pynorm error3
weihuang-jedi Feb 26, 2025
fe88751
fix pynorm error 4
weihuang-jedi Feb 26, 2025
0284fb7
fix pynorm error 5
weihuang-jedi Feb 26, 2025
27ca2f1
fix pynorm error 6
weihuang-jedi Feb 26, 2025
f39bfd9
fix pynorm error 7
weihuang-jedi Feb 27, 2025
8f7fc23
fix pynorm error 8
weihuang-jedi Feb 27, 2025
c248e12
fix pynorm error 9
weihuang-jedi Feb 27, 2025
85120f3
fix pynorm error 10
weihuang-jedi Feb 27, 2025
2a7f1c9
fix a syntax error
weihuang-jedi Feb 27, 2025
9f8ee2d
use just /lutre/jenkins for CI testing work directory
weihuang-jedi Feb 27, 2025
e47a913
Merge branch 'NOAA-EMC:develop' into download-subset-fix-data
weihuang-jedi Feb 28, 2025
0373d6f
Update
kbooker79 Feb 28, 2025
ff79845
Update.
kbooker79 Feb 28, 2025
50bb907
Updating label.
kbooker79 Feb 28, 2025
02c6fa4
Update
kbooker79 Feb 28, 2025
08a2ee5
Update
kbooker79 Feb 28, 2025
2c6bd87
Update
kbooker79 Feb 28, 2025
64769ab
Update
kbooker79 Feb 28, 2025
cfb4132
Updating.
kbooker79 Feb 28, 2025
7feb03f
Updating.
kbooker79 Feb 28, 2025
61b57f1
Update.
kbooker79 Feb 28, 2025
36e00a2
Updating.
kbooker79 Feb 28, 2025
4149ec0
Reverting changes.
kbooker79 Feb 28, 2025
1017414
Updating.
kbooker79 Feb 28, 2025
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
12 changes: 8 additions & 4 deletions ci/Jenkinsfile4AWS
Original file line number Diff line number Diff line change
Expand Up @@ -5,9 +5,9 @@ def HOMEgfs = 'none'
def CI_CASES = ''
def GH = 'none'
// Location of the custom workspaces for each machine in the CI system. They are persistent for each iteration of the PR.
def NodeName = [hera: 'Hera-EMC', orion: 'Orion-EMC', hercules: 'Hercules-EMC', gaea: 'Gaea', noaacloud: 'awsepicglobalworkflow']
def custom_workspace = [hera: '/scratch1/NCEPDEV/global/CI', orion: '/work2/noaa/stmp/CI/ORION', hercules: '/work2/noaa/global/CI/HERCULES', gaea: '/gpfs/f5/epic/proj-shared/global/CI', noaacloud: '/lustre/jenkins']
def repo_url = 'git@github.com:NOAA-EMC/global-workflow.git'
def NodeName = [hera: 'Hera-EMC', orion: 'Orion-EMC', hercules: 'Hercules-EMC', gaea: 'Gaea', awsepicglobalworkflow: 'Awsepicglobalworkflow']
def custom_workspace = [hera: '/scratch1/NCEPDEV/global/CI', orion: '/work2/noaa/stmp/CI/ORION', hercules: '/work2/noaa/global/CI/HERCULES', gaea: '/gpfs/f5/epic/proj-shared/global/CI', awsepicglobalworkflow: '/lustre/jenkins']
def repo_url = 'git@github.com:NOAA-EPIC/global-workflow-cloud.git'
def STATUS = 'Passed'

pipeline {
Expand Down Expand Up @@ -42,11 +42,15 @@ pipeline {
} else {
echo "This is parent job so getting list of nodes matching labels:"
for (label in pullRequest.labels) {
print label
if (label.matches("CI-(.*?)-Ready")) {
def machine_name = label.split('-')[1].toString().toLowerCase()
print machine_name
jenkins.model.Jenkins.get().computers.each { c ->
if (c.node.selfLabel.name == NodeName[machine_name]) {
run_nodes.add(c.node.selfLabel.name)
print "condition matched"
print run_nodes.size()
}
}
}
Expand All @@ -56,7 +60,7 @@ pipeline {
run_nodes.init().each { node ->
def machine_name = node.split('-')[0].toLowerCase()
echo "Spawning job on node: ${node} with machine name: ${machine_name}"
build job: "/global-workflow/EMC-Global-Pipeline/PR-${env.CHANGE_ID}", parameters: [
build job: "/global-workflow/EPIC-AWS-Global-Workflow-Pipeline-Test/PR-${env.CHANGE_ID}", parameters: [
string(name: 'machine', value: machine_name),
string(name: 'Node', value: node) ],
wait: false
Expand Down
351 changes: 351 additions & 0 deletions ush/fetch-fix-data.py
Original file line number Diff line number Diff line change
@@ -0,0 +1,351 @@
#!/usr/bin/env python
# cfetch-fix-data.py
# wei.huang@noaa.gov
# 2025-02-26
# script to download a subset of FIX data to local machines.
import os
import time
import sys
import getopt
import subprocess
from pathlib import Path

# ------------------------------------------------------------------------------


class FetchFIXdata():

def __init__(self, atmgridarray=['C48'], ocngridarray=['500'], localdir=None, verbose=0):

self.aws_fix_bucket = f's3://noaa-nws-global-pds/fix'
self.aws_cp = f'aws --no-sign-request s3 cp'
self.aws_sync = f'aws --no-sign-request s3 sync'

self.atmgridarray = atmgridarray
self.ocngridarray = ocngridarray
self.localdir = localdir
self.verbose = verbose

# if (os.path.isdir(localdir)):
# print('Prepare to download FIX data for %s and %s to %s' %(atmgrid, ocngrid, localdir))
# else:
# print(f'local dir: <{localdir}> does not exist. Stop')
# sys.exit(-1)

self.verdict = {}
self.s3dict = {}
self.s3dict['raworog'] = f'raw/orog'

if (self.localdir.find('fix') < 0):
self.targetdir = f'{self.localdir}/fix.subset'
else:
self.targetdir = self.localdir

# --------------------------------------------------------------------------
def update_s3dict(self):

self.update_s3dick_grid_independent()
self.add_grid_data()

if (self.verbose):
self.printinfo()

# --------------------------------------------------------------------------
def update_s3dick_grid_independent(self):

for key in self.fix_ver_dict.keys():
val = self.fix_ver_dict[key]
if (key == 'aer_ver'):
self.s3dict['aer'] = f'aer/{val}'
elif (key == 'am_ver'):
self.s3dict['am'] = f'am/{val}'
elif (key == 'chem_ver'):
self.s3dict['fimdata_chem'] = f'chem/{val}/fimdata_chem'
self.s3dict['Emission_data'] = f'chem/{val}/Emission_data'
elif (key == 'datm_ver'):
self.s3dict['cfsr'] = f'datm/{val}/cfsr'
self.s3dict['gefs'] = f'datm/{val}/gefs'
self.s3dict['gfs'] = f'datm/{val}/gfs'
self.s3dict['mom6'] = f'datm/{val}/mom6'
elif (key == 'glwu_ver'):
self.s3dict['glwu'] = f'glwu/{val}'
elif (key == 'gsi_ver'):
self.s3dict['gsi'] = f'gsi/{val}'
elif (key == 'lut_ver'):
self.s3dict['lut'] = f'lut/{val}'
elif (key == 'mom6_ver'):
self.s3dict['mom6post'] = f'mom6/{val}/post'
elif (key == 'reg2grb2_ver'):
self.s3dict['reg2grb2'] = f'reg2grb2/{val}'
elif (key == 'sfc_climb_ver'):
self.s3dict['sfc_climo'] = f'sfc_climo/{val}'
elif (key == 'verif_ver'):
self.s3dict['verif'] = f'verif/{val}'
elif (key == 'wave_ver'):
self.s3dict['wave'] = f'wave/{val}'

# --------------------------------------------------------------------------
def add_grid_data(self):

for key in self.fix_ver_dict.keys():
val = self.fix_ver_dict[key]
if (key == 'orog_ver'):
self.add_atmgrid2s3dict('orog', key, val)
elif (key == 'ugwd_ver'):
self.add_atmgrid2s3dict('ugwd', key, val)
elif (key == 'mom6_ver'):
self.add_ocngrid2s3dict('mom6', key, val)
elif (key == 'cice_ver'):
self.add_ocngrid2s3dict('cice', key, val)
elif (key == 'cpl_ver'):
self.add_cpl2s3dict('cpl', key, val)

# --------------------------------------------------------------------------
def add_atmgrid2s3dict(self, varname, key, val):

for atmgrid in self.atmgridarray:
newkey = f'{key}_{atmgrid}'
self.s3dict[newkey] = f'{varname}/{val}/{atmgrid}'

# -------------------------------------------------------------------------
def add_ocngrid2s3dict(self, varname, key, val):

for ocngrid in self.ocngridarray:
newkey = f'{key}_{atmgrid}'
self.s3dict[newkey] = f'{varname}/{val}/{ocngrid}'

# -------------------------------------------------------------------------
def add_cpl2s3dict(self, varname, key, val):

for atmgrid in self.atmgridarray:
for ocngrid in self.ocngridarray:
newkey = f'{key}_a{atmgrid}o{ocngrid}'
self.s3dict[newkey] = f'{varname}/{val}/a{atmgrid}o{ocngrid}'

# -------------------------------------------------------------------------
def printinfo(self):

print(f'Preparing to fetch')
print(f'ATM grid: {self.atmgridarray}')
print(f'ONC grid: {self.ocngridarray}')
print(f'From: {self.aws_fix_bucket}')
print(f'To: {self.targetdir}')
for key in self.s3dict.keys():
val = self.s3dict[key]
print(f'{key}: {val}')

# -------------------------------------------------------------------------
def fetchdata(self):

if (self.verbose):
print(f'Create local fix dir: {self.targetdir}')

path = Path(self.targetdir)
path.mkdir(parents=True, exist_ok=True)

self.fetch_ugwp_limb_tau()

for key in self.s3dict.keys():
self.fetch_dir(self.s3dict[key])

# -------------------------------------------------------------------------
def fetch_dir(self, dir):

remotedir = f'{self.aws_fix_bucket}/{dir}'
localdir = f'{self.targetdir}/{dir}'
cmd = f'{self.aws_sync} {remotedir} {localdir}'
self.download_dir(cmd, localdir)

# --------------------------------------------------------------------------
def download_dir(self, cmd, localdir):

# returned_value = os.system(cmd) # returns the exit code in unix
# print('returned value:', returned_value)

if (os.path.isdir(localdir)):
print(f'{localdir} already exist. skip')
else:
parentdir, dirname = os.path.split(localdir)
if (self.verbose):
print(f'Create local {parentdir} dir:')
path = Path(parentdir)
path.mkdir(parents=True, exist_ok=True)
if (self.verbose):
print(cmd)
print(f'Downloading {localdir}')
returned_value = subprocess.call(cmd, shell=True) # returns the exit code in unix
if (self.verbose):
print('returned value:', returned_value)

# --------------------------------------------------------------------------
def fetch_ugwp_limb_tau(self):

ugwp_limb_tau_remotepath = f'{self.aws_fix_bucket}/ugwd/{self.fix_ver_dict['ugwd_ver']}/ugwp_limb_tau.nc'
ugwp_limb_tau_localdir = f'{self.targetdir}/ugwd/{self.fix_ver_dict['ugwd_ver']}'
filename = f'{ugwp_limb_tau_localdir}/ugwp_limb_tau.nc'
path = Path(ugwp_limb_tau_localdir)
path.mkdir(parents=True, exist_ok=True)
cmd = f'{self.aws_cp} {ugwp_limb_tau_remotepath} {filename}'
self.download_file(cmd, filename)

# -------------------------------------------------------------------------
def download_file(self, cmd, filename):

# returned_value = os.system(cmd) # returns the exit code in unix
# print('returned value:', returned_value)

if (os.path.isfile(filename)):
print(f'{filename} already exist. skip')
else:
if (self.verbose):
print(cmd)
print(f'Downloading {filename}')
returned_value = subprocess.call(cmd, shell=True) # returns the exit code in unix
if (self.verbose):
print('returned value:', returned_value)

# --------------------------------------------------------------------------
def set_fix_ver_from_gwhome(self, gwhome, verdict):

fix_ver_file = f'{gwhome}/versions/fix.ver'
self.fix_ver_dict = verdict
if (os.path.isfile(fix_ver_file)):
with open(fix_ver_file, "r") as file:
for line in file.readlines():
if (line.find('export ') >= 0):
headstr, _, value = line.strip().partition('=')
exphead, _, key = headstr.partition(' ')
self.fix_ver_dict[key] = value
else:
print(f'fix_ver_file: {ix_ver_file}s does not exist.')

# ------------------------------------------------------------------------
def set_default_fix_ver(self, verdict):

self.fix_ver_dict = verdict

# -----------------------------------------------------------------------------


def print_usage(verdict):

print('Usage: python fetch-fix-data.py \\')
print(' --atmgrid=AtmospericGrid (for multiple grids, separate with ",") \\')
print(' --ocngrid=OceanGrid (for multiple grids, separate with ",") \\')
print(' --localdir=Your-local-fix-dir \\')
print(' [options]')
print('options are:')
print('\t--gwhome=xxxx (Global-Workflow directory)')

for key in verdict.keys():
print(f'\t--{key}=yyyymmdd default: {verdict[key]}')

# ------------------------------------------------------------------------------


if __name__ == '__main__':

atmgridlist = ['C48', 'C96', 'C192', 'C384', 'C768', 'C1152']
ocngridlist = ['500', '100', '050', '025']

verbose = 0
atmgrid = f'C48'
ocngrid = f'500'
localdir = f'/contrib/global-workflow-shared-data'

# default fix-version
verdict = {}
verdict['aer_ver'] = f'20220805'
verdict['am_ver'] = f'20220805'
verdict['chem_ver'] = f'20220805'
verdict['cice_ver'] = f'20240416'
verdict['cpl_ver'] = f'20230526'
verdict['datm_ver'] = f'20220805'
verdict['glwu_ver'] = f'20220805'
verdict['gsi_ver'] = f'20240208'
verdict['lut_ver'] = f'20220805'
verdict['mom6_ver'] = f'20240416'
verdict['orog_ver'] = f'20231027'
verdict['reg2grb2_ver'] = f'20220805'
verdict['sfc_climo_ver'] = f'20220805'
verdict['ugwd_ver'] = f'20240624'
verdict['verif_ver'] = f'20220805'
verdict['wave_ver'] = f'20240105'

gwhome = None

opts, args = getopt.getopt(sys.argv[1:], '', ['help', 'atmgrid=', 'ocngrid=',
'verbose=', 'localdir=',
'gwhome=',
'aer_ver=',
'am_ver=',
'chem_ver=',
'cice_ver=',
'cpl_ver=',
'datm_ver=',
'glwu_ver=',
'gsi_ver=',
'lut_ver=',
'mom6_ver=',
'orog_ver=',
'reg2grb2_ver=',
'sfc_climo_ver=',
'ugwd_ver=',
'verif_ver=',
'wave_ver='])
for o, a in opts:
# print(f'o: {o}, a: {a}')
if o in ['--help']:
print_usage(verdict)
sys.exit(0)
elif o in ['--verbose']:
verbose = int(a)
elif o in ['--atmgrid']:
atmgrid = a
elif o in ['--ocngrid']:
ocngrid = a
elif o in ['--localdir']:
localdir = a
elif o in ['--gwhome']:
gwhome = a
else:
_, vername = o.split('--')
print(f'vername: <{vername}>')
verdict[vername] = a

if (atmgrid.find(',') > 0):
atmgridarray = atmgrid.split(',')
else:
atmgridarray = [atmgrid]

for grid in atmgridarray:
if (grid not in atmgridlist):
print('atmgrid: ', grid)
print('is not in supported grids: ', atmgridlist)
print_usage(verdict)
sys.exit(-1)

if (ocngrid.find(',') > 0):
ocngridarray = ocngrid.split(',')
else:
ocngridarray = [ocngrid]

for grid in ocngridarray:
if (grid not in ocngridlist):
print('ocngrid: ', grid)
print('is not in supported grids: ', ocngridlist)
print_usage(verdict)
sys.exit(-1)

# ------------------------------------------------------------------
ffd = FetchFIXdata(atmgridarray=atmgridarray,
ocngridarray=ocngridarray,
localdir=localdir, verbose=verbose)

if (gwhome is None):
ffd.set_default_fix_ver(verdict)
else:
ffd.set_fix_ver_from_gwhome(gwhome, verdict)

ffd.update_s3dict()
ffd.fetchdata()