gmao-cda / cycle-oletkf

controlling scripts for doing cycling DA with GEOS
0 stars 0 forks source link

test 0d25 flow #24

Closed gmao-cda closed 1 year ago

gmao-cda commented 1 year ago

python: /discover/nobackup/cda/projects/cylc_0d25 geos_run: /discover/nobackup/cda/projects/0d25_run_2

gmao-cda commented 1 year ago

create a new fcst case with old approach

two baseline for this case

clone

cd /discover/nobackup/projects/gmao/scda_iesa/cda/GEOSgcm_08Nov2022/install/bin ./gcm_setup

prepare init

cd /discover/nobackup/cda/projects/0d25_benchmark python3 ../mytools/prep_restart_files_GEOSgcm.py . /discover/nobackup/cda/projects/0d25_run_2/restarts/restarts.e20200101_00z.tar

run

cd /discover/nobackup/cda/projects/0d25_benchmark sbatch gcm_run.j

Time profiling

gmao-cda commented 1 year ago

Generate clean ensemble output directory

./create_empty_exp_directory.py /discover/nobackup/cda/projects2/odas3_0d25 --members 3 --startDate 2017010100 ./prep_init_ensemble_from_restarts.py --members 3 --startDate 2017010100 --cfgFile oletkf_config_0d25.yaml /discover/nobackup/cda/projects2/odas3_0d25

gmao-cda commented 1 year ago

Need to use prep_geos_fcst to prepare the directory

C4MEM="0001"
CYCLE_DATE="20170101T00"

CYCLE_LETKF_SRCS="/discover/nobackup/cda/develop_space/cycle-oletkf"
SHARED_TMP_DIR="/discover/nobackup/cda/projects2/tmp_0d25"
EXP_DIR="/discover/nobackup/cda/projects2/odas3_0d25"
FCST_HRS=6
FWD_EXEC="/discover/nobackup/projects/gmao/scda_iesa/cda/GEOSgcm_08Nov2022/install/bin/GEOSgcm.x"

 echo "prep_geos_fcst ${C4MEM} $CYCLE_DATE"

 echo "${CYCLE_LETKF_SRCS}/prep_geos_fcst.py ${SHARED_TMP_DIR}/${CYCLE_DATE}_fcst${C4MEM} \
 --expTplDir ${CYCLE_LETKF_SRCS}/exp_template/exp_template_0d25deg  \
 --cycleBkgdDir ${EXP_DIR}/bkgd/${C4MEM}/${CYCLE_DATE} \
 --cycleAnalDir ${EXP_DIR}/anal/${C4MEM}/${CYCLE_DATE} \
 --fwdExec ${FWD_EXEC} \
 --fcstStartDate ${CYCLE_DATE} \
 --fcstHrs ${FCST_HRS}"

reprepare the initial conditions

cd ${SHARED_TMP_DIR}/${CYCLE_DATE}_fcst${C4MEM} rm -rf *_rst rm RESTART python3 ../mytools/prep_restart_files_GEOSgcm.py . /discover/nobackup/cda/projects/0d25_run_2/restarts/restarts.e20200101_00z.tar modify CAP.rc: BEG_DATE: 20200101 000000 END_DATE: 20200101 060000 Modify cap_restart: 20200101 000000

Run the run using gcm_run_lib_v2.py

modify cylc/fcst_0d25_v2 and run (md5sum check passed)

gmao-cda commented 1 year ago
mkdir -p /discover/nobackup/cda/projects2/odas3_0d25/obs/20170101T06
mkdir -p /discover/nobackup/cda/projects2/odas3_0d25/obs/20170101T12
cp /discover/nobackup/cda/projects2/obs/SMAP_L2B_SSS_36950_20220101T005200_R18240_V5.0.h5 /discover/nobackup/cda/projects2/odas3_0d25/obs/20170101T00/20170101T00.sss.h5
cp /discover/nobackup/cda/projects2/obs/SMAP_L2B_SSS_36950_20220101T005200_R18240_V5.0.h5 /discover/nobackup/cda/projects2/odas3_0d25/obs/20170101T06/20170101T06.sss.h5
cp /discover/nobackup/cda/projects2/obs/SMAP_L2B_SSS_36950_20220101T005200_R18240_V5.0.h5 /discover/nobackup/cda/projects2/odas3_0d25/obs/20170101T12/20170101T12.sss.h5
gmao-cda commented 1 year ago

L2SST exp

./create_empty_exp_directory.py /discover/nobackup/cda/projects2/odas3_0d25_viirs --members 3 --startDate 2019090100 ./prep_init_ensemble_from_restarts.py --members 3 --startDate 2019090100 --cfgFile oletkf_config_0d25_viirs.yaml /discover/nobackup/cda/projects2/odas3_0d25_viirs

gmao-cda commented 1 year ago

mkdir -p /discover/nobackup/cda/projects2/odas3_0d25_viirs/obs/20190901T06

cp /discover/nobackup/cda/projects2/viirs_sst/superob_viirs_sst.2019090106.dat /discover/nobackup/cda/projects2/odas3_0d25_viirs/obs/20190901T06/20190901T06.sst.bin