Skip to content

sensible+latent heatfluxes using linear bulk formula #633

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
eclare108213 merged 13 commits into from
Sep 27, 2021
Merged

sensible+latent heatfluxes using linear bulk formula #633

Show file tree
Hide file tree
Changes from 4 commits
Commits
Show all changes
13 commits
Select commit Hold shift + click to select a range
db07a38
'heatflux_linear' flag: sensible+latent heatfluxes using traditional …
mhrib Sep 14, 2021
ce7b41d
Add option atmbndy='mixed' boundary layer condition
mhrib Sep 15, 2021
416bd74
New options for 'atmbndy'
mhrib Sep 15, 2021
316538c
For backward compability, rename "atmbndy" = default to similarity
mhrib Sep 16, 2021
bc292e5
Change wording for text related to atmbndy
mhrib Sep 23, 2021
779509a
Change wording for text related to atmbndy
mhrib Sep 23, 2021
b759d49
Change wording for text related to atmbndy
mhrib Sep 23, 2021
d3dd99e
Spelling error etc.
mhrib Sep 23, 2021
3954d4b
Update Icepack with required atmbndy changes
mhrib Sep 24, 2021
fa92c7f
Added test for atmbndy={constant,mixed}
mhrib Sep 27, 2021
dd18f9b
Add test for atmbndy={constant,mixed}
mhrib Sep 27, 2021
9d2aacb
Test for atmbndy={constant,mixed}
mhrib Sep 27, 2021
290744b
clean up. Renamed without underscores
mhrib Sep 27, 2021
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
23 changes: 18 additions & 5 deletions cicecore/cicedynB/general/ice_init.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -448,7 +448,7 @@ subroutine input_data
albsnowv = 0.98_dbl_kind ! cold snow albedo, visible
albsnowi = 0.70_dbl_kind ! cold snow albedo, near IR
ahmax = 0.3_dbl_kind ! thickness above which ice albedo is constant (m)
atmbndy = 'default' ! or 'constant'
atmbndy = 'similarity' ! Atm boundary layer: 'similarity', 'constant' or 'mixed'
default_season = 'winter' ! default forcing data, if data is not read in
fyear_init = 1900 ! first year of forcing cycle
ycycle = 1 ! number of years in forcing cycle
Expand Down Expand Up @@ -1214,6 +1214,14 @@ subroutine input_data
endif
endif

if (trim(atmbndy) == 'default') then
if (my_task == master_task) then
write(nu_diag,*) subname//' WARNING: atmbndy = default'
write(nu_diag,*) subname//' WARNING: For backward compability, set atmbndy = similarity'

@phil-blain phil-blain Sep 16, 2021

Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

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

Small wording suggestion:

Suggested change
write(nu_diag,*) subname//' WARNING: atmbndy = default'
write(nu_diag,*) subname//' WARNING: For backward compability, set atmbndy = similarity'
write(nu_diag,*) subname//' WARNING: "atmbndy = default" is deprecated, use "atmbndy = similarity"'
write(nu_diag,*) subname//' WARNING: setting atmbndy = similarity'

@mhrib mhrib Sep 23, 2021

Copy link
Copy Markdown
Contributor Author

Choose a reason for hiding this comment

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

rewritten

endif
atmbndy = 'similarity'
endif

if (formdrag) then
if (trim(atmbndy) == 'constant') then
if (my_task == master_task) write(nu_diag,*) subname//' ERROR: formdrag=T and atmbndy=constant'
Expand Down Expand Up @@ -1641,13 +1649,18 @@ subroutine input_data
write(nu_diag,1010) ' rotate_wind = ', rotate_wind,' : rotate wind/stress to computational grid'
write(nu_diag,1010) ' formdrag = ', formdrag,' : use form drag parameterization'
write(nu_diag,1000) ' iceruf = ', iceruf, ' : ice surface roughness at atmosphere interface (m)'
if (trim(atmbndy) == 'default') then
tmpstr2 = ' : stability-based boundary layer'
if (trim(atmbndy) == 'constant') then
tmpstr2 = ' : boundary layer uses bulk transfer coefficients'

@eclare108213 eclare108213 Sep 22, 2021

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 would change 'bulk' to 'constant' in this line. If I remember correctly, the overall parameterization we use is considered a 'bulk' formulation, and we're just choosing whether the coefficients are constant or not. Anyone please correct me if that's wrong. (And yes I used 'bulk' in the original comment!)

@mhrib mhrib Sep 23, 2021

Copy link
Copy Markdown
Contributor Author

Choose a reason for hiding this comment

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

Changed / rewritten

elseif (trim(atmbndy) == 'similarity' .or. &
trim(atmbndy) == 'mixed') then
write(nu_diag,1010) ' highfreq = ', highfreq,' : high-frequency atmospheric coupling'
write(nu_diag,1020) ' natmiter = ', natmiter,' : number of atmo boundary layer iterations'
write(nu_diag,1002) ' atmiter_conv = ', atmiter_conv,' : convergence criterion for ustar'
elseif (trim(atmbndy) == 'constant') then
tmpstr2 = ' : boundary layer uses bulk transfer coefficients'
if (trim(atmbndy) == 'similarity') then
tmpstr2 = ' : stability-based boundary layer'
else
tmpstr2 = ' : stability-based boundary layer, but constant for sensible+latent heatfluxes'

@eclare108213 eclare108213 Sep 22, 2021

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.

How about
' : stability-based wind stress, constant-coefficient sensible+latent heat fluxes'

@mhrib mhrib Sep 23, 2021

Copy link
Copy Markdown
Contributor Author

Choose a reason for hiding this comment

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

Changed / rewritten

endif
else
tmpstr2 = ' : unknown value'
endif
Expand Down
2 changes: 1 addition & 1 deletion configuration/scripts/ice_in
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -223,7 +223,7 @@

&forcing_nml
formdrag = .false.
atmbndy = 'default'
atmbndy = 'similarity'
rotate_wind = .true.
calc_strair = .true.
calc_Tsfc = .true.
Expand Down
2 changes: 1 addition & 1 deletion doc/source/cice_index.rst
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -65,7 +65,7 @@ either Celsius or Kelvin units).
"atm_data_dir", "directory for atmospheric forcing data", ""
"atm_data_format", "format of atmospheric forcing files", ""
"atm_data_type", "type of atmospheric forcing", ""
"atmbndy", "atmo boundary layer parameterization (‘default’ or ‘constant’)", ""
"atmbndy", "atmo boundary layer parameterization ('similarity', ‘constant’, or 'mixed')", ""
"avail_hist_fields", "type for history field data", ""
"awtidf", "weighting factor for near-ir, diffuse albedo", "0.36218"
"awtidr", "weighting factor for near-ir, direct albedo", "0.00182"
Expand Down
6 changes: 4 additions & 2 deletions doc/source/user_guide/ug_case_settings.rst
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -528,8 +528,10 @@ forcing_nml
:widths: 15, 15, 30, 15

"", "", "", ""
"``atmbndy``", "``constant``", "bulk transfer coefficients", "``default``"
"", "``default``", "stability-based boundary layer", ""
"``atmbndy``", "string", "bulk transfer coefficients", "``similarity``"
"", "``similarity``", "stability-based boundary layer", ""
"", "``constant``", "constant-based boundary layer", ""
"", "``mixed``", "stability-based, but constant for sensible+latent heatfluxes", ""

@eclare108213 eclare108213 Sep 22, 2021

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.

Likewise,
'stability-based wind stress, constant-coefficient sensible+latent heat fluxes'

@mhrib mhrib Sep 23, 2021

Copy link
Copy Markdown
Contributor Author

Choose a reason for hiding this comment

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

Changed / rewritten

"``atmiter_conv``", "real", "convergence criteria for ustar", "0.0"
"``atm_data_dir``", "string", "path to atmospheric forcing data directory", ""
"``atm_data_format``", "``bin``", "read direct access binary atmo forcing file format", "``bin``"
Expand Down