process_merra2ics_day.sh

#!/bin/bash

#########################################################################
# This file will retrieve archived MERRA2 data from HPSS.
# 
# use ./add_merra2_day.sh YYYYMMDD
#########################################################################
ml hpss

yyyymmdd=${1}
yyyy=$(date -d ${yyyymmdd} +"%Y")


yyyymmddhh=${yyyymmdd}00
# just a print statement
echo ${yyyymmddhh}

# enter directory
pushd ${yyyymmddhh}/gfs/C384/INPUT

# copy python file - change to your copy of merra2_to_ufs_cubesphere.py 
# cp /scratch2/NAGAPE/arl/Barry.Baker/UFS_AERO_CONFIGS/develop/utils/merra2_to_ufs_cubesphere.py .
cp /scratch2/NCEPDEV/stmp3/Barry.Baker/C384_025/MERRA2_UFS_ICS/merra2_to_ufs_cubesphere.py .

# grab MERRA2 from HPSS
htar -xvf /NCEPDEV/emc-naqfc/5year/Barry.Baker/MERRA2_INST_3D_AERO/MERRA2_400.inst_3d_aero_Nv.${yyyy}.nc4 MERRA2_400.inst3_3d_aer_Nv.${yyyymmdd}.nc4

# If you are using your own file for a different day than is archived please replace the path to the file here
MERRA2_PATH='MERRA2_400.inst3_3d_aer_Nv.${yyyymmdd}.nc4'
	    
# Modify the GFS files for UFS-Aersols with MERRA2 ICs
for i in {1..6}; do 
    ./merra2_to_ufs_cubesphere.py -m MERRA2_400.inst3_3d_aer_Nv.${yyyymmdd}.nc4 -c gfs_ctrl.nc -t gfs_data.tile${i}.nc -r C384 
done

if [ -f MERRA2_400.inst3_3d_aer_Nv.${yyyymmdd}.nc4 ]; then
    rm MERRA2_400.inst3_3d_aer_Nv.*.nc4
else
    echo "MERRA2_400.inst3_3d_aer_Nv.${yyyymmdd}.nc4" >> ../../../../getfiles
fi 
    
rm *.old	    
# go back to base directory 
popd