Changeset 12212


Ignore:
Timestamp:
2019-12-12T16:00:01+01:00 (5 months ago)
Author:
dancopsey
Message:

Reinstate fhld_1d into icethd_dh.F90 but set it to zero in icethd.F90. It is used later on in icethd.F90 so it is best to set it to zero early on.

Location:
NEMO/branches/UKMO/NEMO_4.0.1_remove_leads_botm/src/ICE
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/UKMO/NEMO_4.0.1_remove_leads_botm/src/ICE/icethd.F90

    r11715 r12212  
    164164            ! If the grid cell is fully covered by ice (no leads) => transfer energy from the lead budget to the ice bottom budget 
    165165            IF( ( zqld >= 0._wp .AND. at_i(ji,jj) > 0._wp ) .OR. at_i(ji,jj) >= (1._wp - epsi10) ) THEN 
    166                fhld (ji,jj) = rswitch * zqld * r1_rdtice / MAX( at_i(ji,jj), epsi10 ) ! divided by at_i since this is (re)multiplied by a_i in icethd_dh.F90 
     166               fhld (ji,jj) = 0._wp 
    167167               qlead(ji,jj) = 0._wp 
    168168            ELSE 
  • NEMO/branches/UKMO/NEMO_4.0.1_remove_leads_botm/src/ICE/icethd_dh.F90

    r12188 r12212  
    144144      ! 
    145145      DO ji = 1, npti 
    146          zf_tt(ji)         = qcn_ice_bot_1d(ji) + qsb_ice_bot_1d(ji)  
     146         zf_tt(ji)         = qcn_ice_bot_1d(ji) + qsb_ice_bot_1d(ji) + fhld_1d(ji)  
    147147         zq_bot(ji)        = MAX( 0._wp, zf_tt(ji) * rdt_ice ) 
    148148      END DO 
Note: See TracChangeset for help on using the changeset viewer.