Commit a857886f authored by Lars Ackermann's avatar Lars Ackermann

add dynamic downscaling in iterative coupling

parent 9b967620
......@@ -179,9 +179,9 @@ ic_atm2pism_regrid_dyn_downscale() {
ic_atm2pism_regrid_dyn_downscale_generate_elevation_difference
ic_atm2pism_regrid_dyn_downscale_temperature
ic_atm2pism_regrid_dyn_downscale_precipitation
ic_atm2pism_regrid_dyn_downscale_shortwave_downward_radiation
ic_atm2pism_regrid_dyn_downscale_longwave_downward_radiation
ic_atm2pism_regrid_dyn_downscale_cloud_cover
#ic_atm2pism_regrid_dyn_downscale_shortwave_downward_radiation
#ic_atm2pism_regrid_dyn_downscale_longwave_downward_radiation
#ic_atm2pism_regrid_dyn_downscale_cloud_cover
echo; echo -e "\t\t...done."; echo
}
......@@ -201,9 +201,9 @@ ic_atm2pism_regrid_downscale() {
ic_atm2pism_regrid_downscale_generate_elevation_difference
ic_atm2pism_regrid_downscale_temperature
ic_atm2pism_regrid_downscale_precipitation
ic_atm2pism_regrid_downscale_shortwave_downward_radiation
ic_atm2pism_regrid_downscale_longwave_downward_radiation
ic_atm2pism_regrid_downscale_cloud_cover
#ic_atm2pism_regrid_downscale_shortwave_downward_radiation
#ic_atm2pism_regrid_downscale_longwave_downward_radiation
#ic_atm2pism_regrid_downscale_cloud_cover
echo; echo -e "\t\t...done."; echo
}
......@@ -451,26 +451,26 @@ ic_atm2pism_regrid_dyn_downscale_generate_elevation_difference() {
}
ic_atm2pism_regrid_dyn_downscale_temperature() {
echo ""; echo " * downscaling temperature"
DYN_DOWNSCALING_LAPSE_RATE=${DYN_DOWNSCALING_LAPSE_RATE:-"5"}
if (( ${DOWNSCALING_LAPSE_RATE} <= 0 )); then
echo "YOUR DYNAMIC DOWNSCALING LAPSE RATE IS NEGATIVE!"
echo "IT WAS: $DYN_DOWNSCALING_LAPSE_RATE"
exit 6
fi
echo " - Lapse rate applied was Γ=${DYN_DOWNSCALING_LAPSE_RATE} (K/km)"
echo""; echo -e "\t\t* downscaling temperature"
#echo ""; echo " * downscaling temperature"
#DYN_DOWNSCALING_LAPSE_RATE=${DYN_DOWNSCALING_LAPSE_RATE:-"5"}
#if (( ${DOWNSCALING_LAPSE_RATE} <= 0 )); then
# echo "YOUR DYNAMIC DOWNSCALING LAPSE RATE IS NEGATIVE!"
# echo "IT WAS: $DYN_DOWNSCALING_LAPSE_RATE"
# exit 6
#fi
#echo " - Lapse rate applied was Γ=${DYN_DOWNSCALING_LAPSE_RATE} (K/km)"
echo""; echo -e "\t\t* dynamical downscaling temperature"
echo -e "\t\t- falling back to bilinear interpolation for temperature!"
mv ${COUPLE_DIR}/atmosphere_file_for_ice_for_downscaling_${atmosphere_name_temperature}.nc ${COUPLE_DIR}/temperature.nc
}
ic_atm2pism_regrid_dyn_downscale_precipitation() {
ic_atm2pism_regrid_downscale_temperature
ic_atm2pism_regrid_downscale_precipitation
#ic_atm2pism_regrid_downscale_temperature
#ic_atm2pism_regrid_downscale_precipitation
#echo""; echo -e "\t\t* downscaling precipitation"
#echo -e "\t\t- falling back to bilinear interpolation for precipitation!"
#mv ${COUPLE_DIR}/atmosphere_file_for_ice_for_downscaling_${atmosphere_name_precipitation}.nc ${COUPLE_DIR}/precipitation.nc
echo""; echo -e "\t\t* dynamical downscaling precipitation"
echo -e "\t\t- falling back to bilinear interpolation for precipitation!"
mv ${COUPLE_DIR}/atmosphere_file_for_ice_for_downscaling_${atmosphere_name_precipitation}.nc ${COUPLE_DIR}/precipitation.nc
}
ic_atm2pism_regrid_dyn_downscale_shortwave_downward_radiation() {
......@@ -993,8 +993,8 @@ ic_atm2pism_ablation_pdd_set_options_atmo_given() {
${COUPLE_DIR}/pism_coupler_${CHUNK_START_DATE_pism}-${CHUNK_END_DATE_pism}.dat
echo "pism_set_forcing_file___atmosphere_lapse_rate=${COUPLE_DIR}/surface_altitude_coarse_resolution.nc" >> \
${COUPLE_DIR}/pism_forcing_${CHUNK_START_DATE_pism}-${CHUNK_END_DATE_pism}.dat
echo "pism_set_config_value___temp_lapse_rate=${DYN_DOWNSCALING_LAPSE_RATE}" >> \
${COUPLE_DIR}/pism_config_${CHUNK_START_DATE_pism}-${CHUNK_END_DATE_pism}.dat
#echo "pism_set_config_value___temp_lapse_rate=${DYN_DOWNSCALING_LAPSE_RATE}" >> \
#${COUPLE_DIR}/pism_config_${CHUNK_START_DATE_pism}-${CHUNK_END_DATE_pism}.dat
else
echo "pism_set_coupler___atmosphere=given" >> \
${COUPLE_DIR}/pism_coupler_${CHUNK_START_DATE_pism}-${CHUNK_END_DATE_pism}.dat
......@@ -1129,8 +1129,8 @@ ic_debm2pism_ablation_debm_set_options_surface_given() {
${COUPLE_DIR}/pism_coupler_${CHUNK_START_DATE_pism}-${CHUNK_END_DATE_pism}.dat
echo "pism_set_forcing_file___atmosphere_lapse_rate_file=${COUPLE_DIR}/surface_altitude_coarse_resolution.nc" >> \
${COUPLE_DIR}/pism_forcing_${CHUNK_START_DATE_pism}-${CHUNK_END_DATE_pism}.dat
echo "pism_set_config_value___temp_lapse_rate=${DYN_DOWNSCALING_LAPSE_RATE}" >> \
${COUPLE_DIR}/pism_config_${CHUNK_START_DATE_pism}-${CHUNK_END_DATE_pism}.dat
#echo "pism_set_config_value___temp_lapse_rate=${DYN_DOWNSCALING_LAPSE_RATE}" >> \
# ${COUPLE_DIR}/pism_config_${CHUNK_START_DATE_pism}-${CHUNK_END_DATE_pism}.dat
else
echo "pism_set_coupler___atmosphere=given" >> \
${COUPLE_DIR}/pism_coupler_${CHUNK_START_DATE_pism}-${CHUNK_END_DATE_pism}.dat
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment