Search found 7 matches

by pwallhead
Thu Jun 02, 2022 7:50 am
Forum: ROMS Discussion
Topic: How to calculate turbulent shear rate / energy dissipation
Replies: 3
Views: 4888

Re: How to calculate turbulent shear rate / energy dissipation

Thanks John and John! I checked the ROMS code to compare how the turbulent energy dissipation is calculated under MY25_MIXING (my25_corstep.F) vs. GLS_MIXING (gls_corstep.F). It seems that the MY25_MIXING code is indeed consistent with the GLS_MIXING formulation for epsilon (Warner et al. 2005 Eqn 1...
by pwallhead
Tue May 31, 2022 8:36 am
Forum: ROMS Discussion
Topic: How to calculate turbulent shear rate / energy dissipation
Replies: 3
Views: 4888

How to calculate turbulent shear rate / energy dissipation

Hi I would like to calculate (and output) the turbulent shear rate [s-1], which I think means first calculating the turbulent energy dissipation (then dividing by kinematic viscosity and taking square root). I can see how this is done in the COAWST sediment dynamics code, using Warner et al. (2005) ...
by pwallhead
Thu Mar 07, 2019 9:09 am
Forum: ROMS Bugs
Topic: A small bug in ana_srflux.h
Replies: 12
Views: 9223

Re: A small bug in ana_srflux.h

Hi Strange that our line numbers don't coincide - I refer to the Rutgers code, svn-updated 05/03/2019, and the ana_srflux.h is from ROMS/Functionals/. Anyway here is line 176: zenith=cff1+cff2*COS(Hangle-lonr(i,j)*deg2rad/15.0_r8) and lines 205, 213: & (cff1+cff2*COS(Hangle-lonr(i,j)*deg2rad))) ...
by pwallhead
Wed Mar 06, 2019 10:25 am
Forum: ROMS Bugs
Topic: A small bug in ana_srflux.h
Replies: 12
Views: 9223

Re: A small bug in ana_srflux.h

For me the division by 15.0 on line 176 is a definite bug: 1) It breaks periodicity over 360 degrees change in lonr 2) It is inconsistent with lines 205 and 213 3) As mitya points out, the division is in effect already accounted for in deg2rad Changes to account for local time definitions should occ...
by pwallhead
Wed Mar 06, 2019 9:11 am
Forum: ROMS Discussion
Topic: srflx in bulk_flux.F
Replies: 3
Views: 1876

Re: srflx in bulk_flux.F

Hi Kate I checked ccsm_flux.F in your latest code and it seems that issues (2) and (3) still apply. (2) is the primary concern: srflx is apparently not corrected for ice/snow albedo (variable albedo_ice). The radiative heating of the ice (qswi) is corrected for albedo_ice, and therefore so are qi_o_...
by pwallhead
Thu Dec 20, 2018 8:21 pm
Forum: ROMS Discussion
Topic: srflx in bulk_flux.F
Replies: 0
Views: 1168

srflx in bulk_flux.F

I have a few questions about the Hedstrom code, perhaps someone can help me. Apologies if I missed something obvious: 1) It seems that ice does not affect the skin heating (SRad). Why not? 2) It seems that the ice/snow albedo is used to correct the shortwave entering the ice, but not the shortwave e...
by pwallhead
Thu Dec 20, 2018 8:21 pm
Forum: ROMS Discussion
Topic: srflx in bulk_flux.F
Replies: 3
Views: 1876

srflx in bulk_flux.F

I have a few questions about the Hedstrom code, perhaps someone can help me. Apologies if I missed something obvious: 1) It seems that ice does not affect the skin heating (SRad). Why not? 2) It seems that the ice/snow albedo is used to correct the shortwave entering the ice, but not the shortwave e...