Changes between Version 26 and Version 27 of ticket/0665_mass_heat_salt_fluxes
- Timestamp:
- 2010-05-14T13:29:49+02:00 (14 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
ticket/0665_mass_heat_salt_fluxes
v26 v27 10 10 '''Branch''' : [https://forge.ipsl.jussieu.fr/nemo/browser/branches/DEV_r1837_mass_heat_salt_fluxes DEV_r1837_mass_heat_salt_fluxes ] 11 11 ---- 12 12 13 useful commands 13 14 {{{ … … 211 212 ==== (2) LIM-3 specific changes see '''revision XXXX''' ==== 212 213 213 In LIM3 case, the problem is somewhat different. Indeed, LIM3 is written in an energy conserving way (Bitz & Lipscomb 1999). Therefore the heat exchanged i thboth atmosphere and ocean is the internal energy (i.e. the heat flux through the interface + the heat content of the mass flux referenced to rt0). There is no need to introduce rdq_ice and rdq_snw fields.214 In LIM3 case, the problem is somewhat different. Indeed, LIM3 is written in an energy conserving way (Bitz & Lipscomb 1999). Therefore the heat exchanged in both atmosphere and ocean is the internal energy (i.e. the heat flux through the interface + the heat content of the mass flux referenced to rt0). There is no need to introduce rdq_ice and rdq_snw fields. 214 215 215 216 Five key elements have been changed : … … 257 258 258 259 259 QUESTION ==> Martin Sublimation in coupled mode !!!! ou comment??? 260 mettre qen_snw à la valeur du bilan sublim - snow ???? 261 est-ce que ça marche ou limitation ?? 262 263 Coupled mode check the value of cpic used in the atmosphere as well as rcp !!! :-) 260 '''OPEN QUESTION ==>''' Sublimation in coupled mode !!!! ===>>> Seb and Martin help required 261 262 '''OPEN QUESTION ==>''' Coupled mode check the value of cpic used in the atmosphere as well as rcp !!! :-) 264 263 265 264 … … 270 269 Introduce betas in the SBC, no more in LIM-3 alone. 271 270 272 reason: (1) make it also avalable in LIM-2 and ice-IF cases 273 274 (2) make it availablein coupled mode (sbccpl.F90 will have to use it) 271 '''reason (1) :''' make it available in coupled mode (sbccpl.F90 will have to use it) ==>> Following change in LIM3 & sbccpl 275 272 276 273 • '''sbc_oce.F90,sbcmod.F90, ice.F90 ''' : in namelist namsbc, add the following variable: … … 312 309 313 310 ! partitioning of snownfall : the snowfall over icy area is modified 314 zzc = za_o** betas311 zzc = za_o**rn_psnw 315 312 zzc_ice = sprecip(ji,jj) * ( 1.e0 - zzc ) / MAX( za_i, 1.e-20 ) * zinda 316 313 zzc_oce = sprecip(ji,jj) * zzc / MAX( za_o, 1.e-20 ) * ( 1.e0 - zinda ) … … 326 323 !!gm ? & + fdtcn(ji,jj) * za_i & ! turbulent ice-ocean heat 327 324 & + fsbbq(ji,jj) *( 1.- zindb ) * za_o ) ! residual heat from previous step 325 ! 326 sprecip(ji,jj) = zzc_ice * sprecip(ji,jj) ! update snowfall and its internal energy over icy area 327 qen_snw(ji,jj) = zzc_ice * qen_snw(ji,jj) ! (used only after over icy area) 328 328 !!gm end 329 329 }}} … … 334 334 '''POTENTIAL BUG''' It seems to me that the among of heat exchange through turbulent fluxes between the ice and the ocean should be fdtcn(ji,jj) * za_i, and not fdtcn(ji,jj) * za_o. Indeed the ice-ocean surface of contact is proportional to at_i (=za_i), not za_o ! '''===>> TO BE CHECKED with Martin''' 335 335 336 • '''limthd_dh.F90 ''' : remove the snow partitioning: lines 337 {{{ 338 zcoeff = ( 1.0 - ( 1.0 - at_i_b(ji) )**betas ) / at_i_b(ji) 339 zdh_s_pre(ji) = zcoeff * sprecip_1d(ji) * rdt_ice / rhosn 340 }}} 341 becomes : 342 {{{ 343 zdh_s_pre(ji) = sprecip_1d(ji) * rdt_ice / rhosn 344 }}} 345 346 • '''limsbc.F90 ''' : remove the lecture from namicethd namelist. Change the budget over lead from: 347 {{{ 348 & + sprecip(ji,jj) * ( 1. - pfrld(ji,jj)**betas ) & ! remove snow is intercepted by sea-ice (previous frld) 349 }}} 350 becomes : 351 {{{ 352 & + sprecip(ji,jj) * ( 1. - pfrld(ji,jj) ) ! remove snow intercepted by sea-ice (before thermo) 353 }}} 354 355 • '''sbccpl.F90 ''' : add the partitioning of snow : 356 357 '''reason (2) :''' make it also avalable in LIM-2. ==>> the Following change is LIM-2 358 359 • '''limthd_2.F90 ''' : remove the lecture from namicethd namelist. Change the budget over lead from: 336 360 337 361