Changes between Version 26 and Version 27 of ticket/0665_mass_heat_salt_fluxes


Ignore:
Timestamp:
2010-05-14T13:29:49+02:00 (10 years ago)
Author:
gm
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • ticket/0665_mass_heat_salt_fluxes

    v26 v27  
    1010'''Branch''' : [https://forge.ipsl.jussieu.fr/nemo/browser/branches/DEV_r1837_mass_heat_salt_fluxes  DEV_r1837_mass_heat_salt_fluxes ]  
    1111---- 
     12 
    1213useful commands  
    1314{{{ 
     
    211212==== (2) LIM-3 specific changes  see '''revision XXXX''' ====  
    212213 
    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 ith 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   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.   
    214215 
    215216   Five key elements have been changed : 
     
    257258 
    258259 
    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 !!! :-) 
    264263 
    265264 
     
    270269   Introduce betas in the SBC, no more in LIM-3 alone. 
    271270 
    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 
    275272 
    276273     • '''sbc_oce.F90,sbcmod.F90, ice.F90 ''' : in namelist namsbc, add the following variable: 
     
    312309 
    313310            ! partitioning of snownfall : the snowfall over icy area is modified 
    314             zzc = za_o**betas 
     311            zzc = za_o**rn_psnw 
    315312            zzc_ice = sprecip(ji,jj) * ( 1.e0 - zzc ) / MAX( za_i, 1.e-20 ) *          zinda  
    316313            zzc_oce = sprecip(ji,jj) *          zzc   / MAX( za_o, 1.e-20 ) * ( 1.e0 - zinda ) 
     
    326323!!gm ?         &             + fdtcn(ji,jj)                *   za_i           &   ! turbulent ice-ocean heat 
    327324               &             + 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) 
    328328!!gm end 
    329329}}} 
     
    334334         '''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''' 
    335335 
     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: 
    336360 
    337361