Skip to content

Commit 223675e

Browse files
Merge pull request #461 from NCAR/CAM_bias_correction
CAM reanalysis bias correction for CLM-DART
2 parents c389957 + 2e0c46c commit 223675e

11 files changed

+2026
-1
lines changed

CHANGELOG.rst

+7
Original file line numberDiff line numberDiff line change
@@ -22,6 +22,13 @@ individual files.
2222

2323
The changes are now listed with the most recent at the top.
2424

25+
**May 8 2023 :: CLM-DART: CAM reanalysis site-level bias correction tool. Tag v10.7.1**
26+
27+
- Initial version of bias correction for CAM reanalysis forcing for
28+
site-level assimilation with CLM-DART.
29+
- mkmf change: make clean removes .mod files.
30+
- bug-fix: readthedocs yaml file for online documentation build.
31+
2532
**April 21 2023 :: MOM6. Tag v10.7.0**
2633

2734
- CESM-MOM6 interface added to DART.

conf.py

+1-1
Original file line numberDiff line numberDiff line change
@@ -21,7 +21,7 @@
2121
author = 'Data Assimilation Research Section'
2222

2323
# The full version, including alpha/beta/rc tags
24-
release = '10.7.0'
24+
release = '10.7.1'
2525
master_doc = 'README'
2626

2727
# -- General configuration ---------------------------------------------------

models/clm/CAM_bias_correction/CAM4_bias_correct.m

+511
Large diffs are not rendered by default.
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,52 @@
1+
#!/bin/bash
2+
3+
# Extracts site level grid cell from CAM4 reanalysis
4+
5+
# CAM4 reanalysis
6+
CAM4_grid_path="/glade/collections/rda/data/ds199.1/"
7+
# Site Level grid cell from CAM4 renanalysis
8+
CAM4_site_path="/glade/work/bmraczka/CAM4_NR1/"
9+
10+
# The US-NR1 flux tower site is used as an example
11+
# NR1 location (40.03, -105.55) or (40.03, 254.45)
12+
# CAM4 grid is 1.95x2.5 or 288 longitude grid points and 96 latitude grid points
13+
# This lat/lon corresponds with doma_lon= 102 (255 degrees east) ; doma_lat= 69 ; (40.7368 degrees)
14+
15+
16+
for YEAR in {1997..2010..1}
17+
18+
do
19+
20+
echo " "
21+
echo "entering YEAR loop where value of YEAR is:"
22+
echo ${YEAR}
23+
echo " "
24+
for NINST in {0001..0080..1} # Will create 0001, 0002, ... 0080
25+
do
26+
echo " "
27+
echo "entering NINST loop where value of YEAR is:"
28+
echo ${NINST}
29+
echo " "
30+
31+
#Create output directory if necessary
32+
OUTDIR=${CAM4_site_path}${NINST}
33+
if [[ ! -d ${OUTDIR} ]] ; then
34+
mkdir ${OUTDIR}
35+
fi
36+
37+
38+
ncks -d doma_nx,102,102 -d doma_ny,69,69 -d a2x6h_nx,102,102 -d a2x6h_ny,69,69 \
39+
${CAM4_grid_path}/CAM_DATM.cpl_${NINST}.ha2x1dx6h.${YEAR}.nc \
40+
${CAM4_site_path}${NINST}/CAM4_NR1.cpl_${NINST}.ha2x1dx6h.${YEAR}.nc
41+
42+
echo " "
43+
echo "Completed CAM4 NR1 extract"
44+
echo " "
45+
46+
done
47+
done
48+
49+
50+
exit 0
51+
52+

0 commit comments

Comments
 (0)