Changeset 8176


Ignore:
Timestamp:
2017-06-14T18:33:09+02:00 (3 years ago)
Author:
vancop
Message:

Bugfix on ocean-to-ice sensible heat flux found during SIMIP diags

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2016/v3_6_CMIP6_ice_diagnostics/NEMOGCM/NEMO/LIM_SRC_3/limthd.F90

    r8169 r8176  
    151151 
    152152            ! --- Above-freezing sensible heat content (J/m2 grid) 
    153             !-> target zqfr_neg = tmask(ji,jj,1) * rau0 * rcp * fse3t_m(ji,jj) * MIN( ( t_bo(ji,jj) - ( sst_m(ji,jj) + rt0 ) ), 0._wp ) 
     153            zqfr_neg = tmask(ji,jj,1) * rau0 * rcp * fse3t_m(ji,jj) * MIN( ( t_bo(ji,jj) - ( sst_m(ji,jj) + rt0 ) ), 0._wp ) 
    154154 
    155155            ! --- Sensible ocean-to-ice heat flux (W/m2) 
    156156            zfric_u      = MAX( SQRT( ust2s(ji,jj) ), zfric_umin )  
    157             fhtur(ji,jj) = MAX( 0._wp, rswitch * rau0 * rcp * zch  * zfric_u * ( ( sst_m(ji,jj) + rt0 ) - t_bo(ji,jj) ) ) ! W.m-2 
    158             !-> target fhtur(ji,jj) = rswitch * rau0 * rcp * zch  * zfric_u * ( ( sst_m(ji,jj) + rt0 ) - t_bo(ji,jj) ) ! W.m-2 
    159  
    160             fhtur(ji,jj) = rswitch * MIN( fhtur(ji,jj), - zqfr * r1_rdtice / MAX( at_i(ji,jj), epsi10 ) ) 
    161             !-> target fhtur(ji,jj) = rswitch * MIN( fhtur(ji,jj), - zqfr_neg * r1_rdtice / MAX( at_i(ji,jj), epsi10 ) ) 
     157            fhtur(ji,jj) = rswitch * rau0 * rcp * zch  * zfric_u * ( ( sst_m(ji,jj) + rt0 ) - t_bo(ji,jj) ) ! W.m-2 
     158 
     159            fhtur(ji,jj) = rswitch * MIN( fhtur(ji,jj), - zqfr_neg * r1_rdtice / MAX( at_i(ji,jj), epsi10 ) ) 
    162160            ! upper bound for fhtur: the heat retrieved from the ocean must be smaller than the heat necessary to reach  
    163161            !                        the freezing point, so that we do not have SST < T_freeze 
Note: See TracChangeset for help on using the changeset viewer.