Skip to content

add incr-update with IAU #129

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.

Sign up for GitHub

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

Merged
guillaumevernieres merged 5 commits into from
Sep 20, 2022
Merged

add incr-update with IAU #129

Show file tree
Hide file tree
Changes from 2 commits
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
01f5aea
add incr-update with IAU
hyunchul386 Sep 13, 2022
b840cff
update PR#129
hyunchul386 Sep 14, 2022
9b0c307
pointing total increments
hyunchul386 Sep 14, 2022
b7ffd54
pointing full increments with cleaning
hyunchul386 Sep 14, 2022
0d9891a
cleaning comments
hyunchul386 Sep 14, 2022
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
31 changes: 31 additions & 0 deletions scripts/exgdas_global_marine_analysis_run.sh
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -30,6 +30,30 @@ pwd=$(pwd)
# Utilities
export NLN=${NLN:-"/bin/ln -sf"}

# function that converts a soca incr to something MOM6 is happy with
function socaincr2mom6 {
incr=$1
bkg=$2
grid=$3
incr_out=$4

scratch=scratch_socaincr2mom6
mkdir -p $scratch
cd $scratch

echo "at socaincr2mom6" $bkg $grid

cp $incr inc.nc # TODO: use accumulated incremnet, not outerloop intermediates

@guillaumevernieres guillaumevernieres Sep 14, 2022

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

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

remove comment once you point to the total increment.

@hyunchul386 hyunchul386 Sep 14, 2022

Copy link
Copy Markdown
Collaborator Author

Choose a reason for hiding this comment

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

updated

@guillaumevernieres guillaumevernieres Sep 14, 2022

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

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

doesn't look like it, I still see # TODO: use accumulated incremnet, not outerloop intermediates ... which presumably was done

ncks -A -C -v h $bkg inc.nc # Replace h incrememnt (all 0's) by h background (expected by MOM)
ncrename -d zaxis_1,Layer inc.nc # Rename zaxis_1 to Layer
ncks -A -C -v Layer $bkg inc.nc # Replace dimension-less Layer with dimensional Layer
mv inc.nc inc_tmp.nc # ... dummy copy
ncwa -O -a Time inc_tmp.nc inc.nc # Remove degenerate Time dimension
ncks -A -C -v lon $grid inc.nc # Add longitude
ncks -A -C -v lat $grid inc.nc # Add latitude
mv inc.nc $incr_out
}

function bump_vars()
{
tmp=$(ls -d ${1}_* )
Expand Down Expand Up @@ -108,6 +132,13 @@ concatenate_bump 'bump3d'
clean_yaml var.yaml
$APRUN_SOCAANAL $JEDI_BIN/soca_var.x var.yaml 2>var.err


# increment update for MOM6 iau
# Make sure the last (most recent) increment file put into socaincr2mom6 as $1

@guillaumevernieres guillaumevernieres Sep 14, 2022

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

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

I don't know what this comment means ...

@hyunchul386 hyunchul386 Sep 14, 2022

Copy link
Copy Markdown
Collaborator Author

Choose a reason for hiding this comment

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

updated

( socaincr2mom6 `ls -t ${DATA}/Data/ocn.*iter*.incr* | head -1` ${DATA}/INPUT/MOM.res.nc ${DATA}/soca_gridspec.nc ${DATA}/Data/inc.nc )
Comment thread
hyunchul386 marked this conversation as resolved.
Outdated



################################################################################
set +x
if [ $VERBOSE = "YES" ]; then
Expand Down
1 change: 1 addition & 0 deletions test/soca/runtime_vars.sh
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -47,3 +47,4 @@ export R2D2_OBS_SRC=gdasapp
export R2D2_OBS_WINDOW=24 # R2D2 sampling DB window

export APRUN_SOCAANAL="$MPIEXEC_EXEC $MPIEXEC_NPROC 6"

1 change: 1 addition & 0 deletions test/soca/test_run.sh
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -25,3 +25,4 @@ test_file ${project_binary_dir}/test/soca/3dvar/analysis/Data/ocn.3dvarfgat_pse
test_file ${project_binary_dir}/test/soca/3dvar/analysis/Data/ocn.3dvarfgat_pseudo.an.2018-04-15T09:00:00Z.nc
test_file ${project_binary_dir}/test/soca/3dvar/analysis/Data/ocn.3dvarfgat_pseudo.an.2018-04-15T12:00:00Z.nc
test_file ${project_binary_dir}/test/soca/3dvar/analysis/Data/ocn.3dvarfgat_pseudo.an.2018-04-15T15:00:00Z.nc
test_file ${project_binary_dir}/test/soca/3dvar/analysis/Data/inc.nc