New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
user/clem/Heat_Flux_air-ice-ocean_interface (diff) – NEMO

Changes between Version 10 and Version 11 of user/clem/Heat_Flux_air-ice-ocean_interface


Ignore:
Timestamp:
2016-05-19T16:52:06+02:00 (8 years ago)
Author:
clem
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • user/clem/Heat_Flux_air-ice-ocean_interface

    v10 v11  
     1== Decision? (Update 2016/05/19) == 
     2If the atm sends to the ocean/ice: {{{Qsnow = Cp(snow) * Tsnow}}}, and {{{Qrain = Cp(rain) * Train}}}, and {{{Qevap(ice) = Cp(ice) * Tice}}}, and {{{Qevap(ocean) = Cp(ocean) * Toce}}}, then heat is conserved and the differences between atm. physical constants and ocean/ice physical constants are unimportant. 
     3If the atm does not provide these fluxes, then we will make a guess and we have to make sure to have the same constants between Atm. models and NEMO. This is not yet the case. 
     4 
     5 
    16== Where is the problem? == 
    27I spent some time sending emails back and forth trying to solve the problem of heat exchanges at the air - ice/ocean interface, especially those associated with mass exchanges (i.e. heat content of precipitation and evaporation). I think there is a misunderstanding between what the atmosphere expects to send/receive (from an ice-ocean perspective) and what is really sent/received.[[BR]] 
     
    611So far, the total non-solar flux (Qns) is imposed by the atmosphere but it excludes (from what I understand) the heat contents of snow precipitation (Qsnow) and rain (Qrain). It also does not take into account the heat content (if any) of evaporated ice and water (Qevap, already in vapor form) that the atmosphere should receive. The total non-solar (conservative) flux should thus be written (including the snow melting when it reaches the ocean):[[BR]] 
    712 
    8 {{{Qns_tot = Qns + Qsnow + Qrain + Qevap + "Snow-melt-in-ocean"}}} 
     13{{{Qns_tot = Qns + Qsnow + Qrain + Qevap + "Snow-melt"}}} 
    914 
    1015 
     
    2126{{{Qevap (ice)   = 0}}} 
    2227 
    23 {{{"Snow-melt-in-ocean" = sprecip * (1 - ice_frac) * ( - Lfus )}}} 
     28{{{"Snow-melt" = sprecip * ( - Lfus )}}} 
    2429 
    2530with sprecip = snow precip, lprecip = liquid precip, ice_frac = ice fraction, Cp = specific heat of sea water, Lfus = Latent heat of fusion of pure ice at 0deg, SST = sea surface temperature.[[BR]]