Skip to content
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.

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

Bug fixes and changes necessary to run the Ens. var #46

Merged
merged 13 commits into from
Dec 2, 2024
17 changes: 15 additions & 2 deletions algorithm/marine/soca_2cice_global.yaml.j2
Original file line number Diff line number Diff line change
Expand Up @@ -29,7 +29,13 @@ variable change:
sno_lev: 1
cice output:
restart: {{ ice_rst }}
output variables: [sea_water_potential_temperature, sea_water_salinity, sea_water_cell_thickness, sea_ice_area_fraction, sea_ice_thickness, sea_ice_snow_thickness]
output variables:
- sea_water_potential_temperature
- sea_water_salinity
- sea_water_cell_thickness
- sea_ice_area_fraction
- sea_ice_thickness
- sea_ice_snow_thickness

states:
- input:
Expand All @@ -38,7 +44,14 @@ states:
ocn_filename: {{ ocn_ana }}
ice_filename: {{ ice_ana }}
date: '{{ fcst_begin }}'
state variables: [sea_water_potential_temperature, sea_water_salinity, sea_water_cell_thickness, sea_ice_area_fraction, sea_ice_thickness, sea_ice_snow_thickness]
state variables:
- sea_water_potential_temperature
- sea_water_salinity
- sea_water_cell_thickness
- sea_ice_area_fraction
- sea_ice_thickness
- sea_ice_snow_thickness

output:
datadir: Data
exp: soca2cice
Expand Down
15 changes: 13 additions & 2 deletions algorithm/marine/soca_diagb.yaml.j2
Original file line number Diff line number Diff line change
Expand Up @@ -18,9 +18,20 @@ background error:
type: incr

variables:
name: [sea_water_potential_temperature, sea_water_salinity, eastward_sea_water_velocity, northward_sea_water_velocity, sea_water_cell_thickness, sea_surface_height_above_geoid, sea_ice_area_fraction, sea_ice_thickness, sea_ice_snow_thickness, mom6_mld]
name:
- sea_water_potential_temperature
- sea_water_salinity
- eastward_sea_water_velocity
- northward_sea_water_velocity
- sea_water_cell_thickness
- sea_surface_height_above_geoid
- sea_ice_area_fraction
- sea_ice_thickness
- sea_ice_snow_thickness
- mom6_mld

rescale: 2.0 # rescales the filtered std. dev. by "rescale"
vertical e-folding scale: 500.0
rescale: 1.0 # rescales the filtered std. dev. by "rescale"
min sst: 0.0 # Added to sst bkg. err.
max ssh: 0.0 # Limits the amplitude of the unbalanced bkg err
min depth: 500.0 # zero out the bkg. error. at less than min depth
Expand Down
31 changes: 25 additions & 6 deletions algorithm/marine/soca_ensb.yaml.j2
Original file line number Diff line number Diff line change
Expand Up @@ -5,11 +5,20 @@ geometry:

date: '{{ marine_window_begin_iso }}'

layers variable: [hocn]
layers variable: [sea_water_cell_thickness]

increment variables: [tocn, socn, uocn, vocn, ssh, hocn, cicen, hicen, hsnon]
increment variables:
- sea_water_potential_temperature
- sea_water_salinity
- eastward_sea_water_velocity
- northward_sea_water_velocity
- sea_surface_height_above_geoid
- sea_water_cell_thickness
- sea_ice_area_fraction
- sea_ice_thickness
- sea_ice_snow_thickness

set increment variables to zero: [ssh]
set increment variables to zero: [sea_surface_height_above_geoid]

vertical geometry:
date: '{{ marine_window_begin_iso }}'
Expand Down Expand Up @@ -81,16 +90,26 @@ linear variable change:
- linear variable change name: BalanceSOCA

trajectory:
state variables: [tocn, socn, uocn, vocn, ssh, hocn, layer_depth, mld, cicen, hicen, hsnon]
state variables:
- sea_water_potential_temperature
- sea_water_salinity
- eastward_sea_water_velocity
- northward_sea_water_velocity
- sea_surface_height_above_geoid
- sea_water_cell_thickness
- sea_water_depth
- ocean_mixed_layer_thickness
- sea_ice_area_fraction
- sea_ice_thickness
- sea_ice_snow_thickness
date: '{{ marine_window_begin_iso }}'
basename: ./INPUT/
ocn_filename: MOM.res.nc
ice_filename: cice.res.nc
read_from_file: 1

output increment:
# TODO: Revert this when fms can take more than 128 charactres file names
datadir: '{{ enspert_relpath }}/enspert/'
datadir: '{{ enspert_relpath }}/'
date: '{{ marine_window_begin_iso }}'
exp: trash
type: incr
Expand Down
13 changes: 11 additions & 2 deletions algorithm/marine/soca_ensweights.yaml.j2
Original file line number Diff line number Diff line change
Expand Up @@ -5,8 +5,17 @@ geometry:
date: '{{ marine_window_middle_iso }}'

variables:
ice: [cicen, hicen, hsnon]
ocean: [tocn, socn, uocn, vocn, ssh]
ice:
- sea_ice_area_fraction
- sea_ice_thickness
- sea_ice_snow_thickness

ocean:
- sea_water_potential_temperature
- sea_water_salinity
- eastward_sea_water_velocity
- northward_sea_water_velocity
- sea_surface_height_above_geoid

background:
date: '{{ marine_window_begin_iso }}'
Expand Down
13 changes: 8 additions & 5 deletions algorithm/marine/soca_parameters_diffusion_hz.yaml.j2
Original file line number Diff line number Diff line change
Expand Up @@ -30,8 +30,11 @@ background error:
write:
filepath: ./staticb/hz_ocean

- horizontal:
as gaussian: true
fixed value: 50000.0
write:
filepath: ./staticb/hz_ice
# TODO(G): OK to not use what's below but it will need
# to be fixed when we add ice thickness and
# snow depth
# - horizontal:
# as gaussian: true
# fixed value: 50000.0
# write:
# filepath: ./staticb/hz_ice
35 changes: 29 additions & 6 deletions model/marine/marine_background_error_static_diffusion.yaml.j2
Original file line number Diff line number Diff line change
Expand Up @@ -3,15 +3,22 @@ saber central block:
saber block name: diffusion
read:
groups:
- variables: [sea_water_potential_temperature, sea_water_salinity, sea_surface_height_above_geoid]
- variables:
- sea_water_potential_temperature
- sea_water_salinity
- sea_surface_height_above_geoid
- sea_ice_area_fraction
horizontal:
filepath: ./staticb/hz_ocean
vertical:
levels: {{marine_vt_levels}}
filepath: ./staticb/vt_ocean
- variables: [sea_ice_area_fraction]
horizontal:
filepath: ./staticb/hz_ice
# TODO(G): OK to not use what's below but it will need
# to be fixed when we add ice thickness and
# snow depth
# - variables: [sea_ice_area_fraction]
# horizontal:
# filepath: ./staticb/hz_ice

saber outer blocks:
- saber block name: StdDev
Expand All @@ -24,7 +31,23 @@ saber outer blocks:
read_from_file: 3

linear variable change:
input variables: [sea_ice_area_fraction, sea_ice_thickness, sea_ice_snow_thickness, sea_water_salinity, sea_water_potential_temperature, eastward_sea_water_velocity, northward_sea_water_velocity, sea_surface_height_above_geoid]
output variables: [sea_ice_area_fraction, sea_ice_thickness, sea_ice_snow_thickness, sea_water_salinity, sea_water_potential_temperature, eastward_sea_water_velocity, northward_sea_water_velocity, sea_surface_height_above_geoid]
input variables:
- sea_ice_area_fraction
- sea_ice_thickness
- sea_ice_snow_thickness
- sea_water_salinity
- sea_water_potential_temperature
- eastward_sea_water_velocity
- northward_sea_water_velocity
- sea_surface_height_above_geoid
output variables:
- sea_ice_area_fraction
- sea_ice_thickness
- sea_ice_snow_thickness
- sea_water_salinity
- sea_water_potential_temperature
- eastward_sea_water_velocity
- northward_sea_water_velocity
- sea_surface_height_above_geoid
linear variable changes:
- linear variable change name: BalanceSOCA
Loading