Skip to content

Add VIIRS obs processing to dev/gcafs.v1#4473

Merged
CoryMartin-NOAA merged 101 commits intoNOAA-EMC:dev/gcafs.v1from
CoryMartin-NOAA:feature/more-ecflow
Feb 5, 2026
Merged

Add VIIRS obs processing to dev/gcafs.v1#4473
CoryMartin-NOAA merged 101 commits intoNOAA-EMC:dev/gcafs.v1from
CoryMartin-NOAA:feature/more-ecflow

Conversation

@CoryMartin-NOAA
Copy link
Contributor

Description

For GCAFSv1 instead of using ObsForge for upstream observation processing, we will include the VIIRS AOD preprocessing from DCOM directly in the GCAFS workflow. This PR adds the scripts borrowed from ObsForge into the dev branch.

Type of change

  • Bug fix (fixes something broken)
  • New feature (adds functionality)
  • Maintenance (code refactor, clean-up, new CI test, etc.)

Change characteristics

  • Is this change expected to change outputs (e.g. value changes to existing outputs, new files stored in COM, files removed from COM, filename changes, additions/subtractions to archives)? YES/NO (If YES, please indicate to which system(s))
    • GFS
    • GEFS
    • SFS
    • GCAFS
  • Is this a breaking change (a change in existing functionality)? YES/NO
  • Does this change require a documentation update? YES/NO
  • Does this change require an update to any of the following submodules? YES/NO (If YES, please add a link to any PRs that are pending.)
    • EMC verif-global
    • GDAS
    • GFS-utils
    • GSI
    • GSI-monitor
    • GSI-utils
    • UFS-utils
    • UFS-weather-model
    • wxflow

How has this been tested?

Checklist

  • Any dependent changes have been merged and published
  • My code follows the style guidelines of this project
  • I have performed a self-review of my own code
  • I have commented my code, particularly in hard-to-understand areas
  • I have documented my code, including function, input, and output descriptions
  • My changes generate no new warnings
  • New and existing tests pass with my changes
  • This change is covered by an existing CI test or a new one has been added
  • Any new scripts have been added to the .github/CODEOWNERS file with owners
  • I have made corresponding changes to the system documentation if necessary

# COMIN_OBSPROC:COM_OBSPROC_TMPL \
# COMIN_ATMOS_ANALYSIS:COM_ATMOS_ANALYSIS_TMPL
export COMIN_OBSPROC=/lfs/h1/ops/prod/com/obsproc/v1.2/gdas.${PDY}/${cyc}/atmos/
export COMIN_ATMOS_ANALYSIS=/lfs/h1/ops/prod/com/gfs/v16.3/gdas.${PDY}/${cyc}/atmos/
Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@aerorahul @DavidHuber-NOAA while I think the OBSPROC one can be rectified by run_envir=nco this COMIN_ATMOS_ANALYSIS is a good example I think of why this declare_from_tmpl may not be ideal

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Or is the convention here that this should be COMINgfs?

Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

if [[ "${RUN_ENVIR:-emc}" == "nco" ]]; then
COM_OBSPROC_TMPL=$(compath.py "${envir}/obsproc/${obsproc_ver}")'/${RUN}.${YMD}/${HH}/atmos'
COM_OBSFORGE_TMPL=$(compath.py "${envir}/obsforge/${obsforge_ver}")'/${RUN}.${YMD}/${HH}/'
COM_RTOFS_TMPL=$(compath.py "${envir}/rtofs/${rtofs_ver}")
COM_TCVITAL_TMPL=$(compath.py "${envir}/gfs/${gfs_ver}")'/${RUN}.${YMD}/${HH}/atmos'
else
COM_OBSPROC_TMPL='${DMPDIR}/${RUN}${DUMP_SUFFIX}.${YMD}/${HH}/atmos'
COM_OBSFORGE_TMPL='${IODADIR}/${RUN}${DUMP_SUFFIX}.${YMD}/${HH}'
COM_RTOFS_TMPL='${DMPDIR}'
COM_TCVITAL_TMPL='${DMPDIR}/${RUN}.${YMD}/${HH}/atmos'
fi

Comment on lines 29 to 32
result = subprocess.run([nc2ioda_exe, nc2ioda_yaml],
cwd=task_config['DATA'],
capture_output=True,
text=True)
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Use wxflow.executable

logger.info("Copying ioda files to destination COMROOT directory")
yyyymmdd = self.task_config['PDY'].strftime('%Y%m%d')

comout = os.path.join(self.task_config['COMROOT'],
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

COMIN/COMOUT should be declared in the j-jobs

"""
"""
for platform in self.task_config.platforms:
print(f"========= platform: {platform}")
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

use logging here and elsewhere

# COMIN_OBSPROC:COM_OBSPROC_TMPL \
# COMIN_ATMOS_ANALYSIS:COM_ATMOS_ANALYSIS_TMPL
export COMIN_OBSPROC=/lfs/h1/ops/prod/com/obsproc/v1.2/gdas.${PDY}/${cyc}/atmos/
export COMIN_ATMOS_ANALYSIS=/lfs/h1/ops/prod/com/gfs/v16.3/gdas.${PDY}/${cyc}/atmos/
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

if [[ "${RUN_ENVIR:-emc}" == "nco" ]]; then
COM_OBSPROC_TMPL=$(compath.py "${envir}/obsproc/${obsproc_ver}")'/${RUN}.${YMD}/${HH}/atmos'
COM_OBSFORGE_TMPL=$(compath.py "${envir}/obsforge/${obsforge_ver}")'/${RUN}.${YMD}/${HH}/'
COM_RTOFS_TMPL=$(compath.py "${envir}/rtofs/${rtofs_ver}")
COM_TCVITAL_TMPL=$(compath.py "${envir}/gfs/${gfs_ver}")'/${RUN}.${YMD}/${HH}/atmos'
else
COM_OBSPROC_TMPL='${DMPDIR}/${RUN}${DUMP_SUFFIX}.${YMD}/${HH}/atmos'
COM_OBSFORGE_TMPL='${IODADIR}/${RUN}${DUMP_SUFFIX}.${YMD}/${HH}'
COM_RTOFS_TMPL='${DMPDIR}'
COM_TCVITAL_TMPL='${DMPDIR}/${RUN}.${YMD}/${HH}/atmos'
fi

@CoryMartin-NOAA
Copy link
Contributor Author

Thanks @aerorahul I'll make those suggested changes later today. What do you suggest for the COMIN_ATMOS_ANALYSIS? Rename it to COMINgfs? Just define it manually and do not use the tmpl function? Something else?

@aerorahul
Copy link
Contributor

Thanks @aerorahul I'll make those suggested changes later today. What do you suggest for the COMIN_ATMOS_ANALYSIS? Rename it to COMINgfs? Just define it manually and do not use the tmpl function? Something else?

We need a COMINgfs_ATMOS_ANALYSIS variable that makes it clear that the atmos analysis comes from GFS.
In the parm/config/gcafs/config.com, I think we should add a variable COM_GFS_TMPL=$(compath.py "${envir}/gfs/${gfs_ver}")'/${RUN}.${YMD}/${HH}' that should be used to create COMINgfs.
Once that is created, COMINgfs_ATMOS_ANALYSIS=$COMINgfs/analysis/atmos.

@CoryMartin-NOAA CoryMartin-NOAA merged commit 8f1c321 into NOAA-EMC:dev/gcafs.v1 Feb 5, 2026
5 checks passed
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

7 participants