Changeset 13642


Ignore:
Timestamp:
2020-10-20T00:58:34+02:00 (4 months ago)
Author:
clem
Message:

4.0-HEAD: use the transmitted solar for the ice budget. See #2545

Location:
NEMO/releases/r4.0/r4.0-HEAD/src/ICE
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • NEMO/releases/r4.0/r4.0-HEAD/src/ICE/ice1d.F90

    r13284 r13642  
    145145   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) ::   sst_1d 
    146146   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) ::   sss_1d 
     147   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) ::   frq_m_1d 
    147148 
    148149   ! convergence check 
     
    225226      ! 
    226227      ii = ii + 1 
    227       ALLOCATE( sst_1d(jpij) , sss_1d(jpij) , STAT=ierr(ii) ) 
     228      ALLOCATE( sst_1d(jpij) , sss_1d(jpij) , frq_m_1d(jpij) , STAT=ierr(ii) ) 
    228229      ! 
    229230      ii = ii + 1 
  • NEMO/releases/r4.0/r4.0-HEAD/src/ICE/icethd.F90

    r13640 r13642  
    442442         CALL tab_2d_1d( npti, nptidx(1:npti), sst_1d(1:npti), sst_m ) 
    443443         CALL tab_2d_1d( npti, nptidx(1:npti), sss_1d(1:npti), sss_m ) 
     444         CALL tab_2d_1d( npti, nptidx(1:npti), frq_m_1d(1:npti), frq_m ) 
    444445         ! 
    445446         ! to update ice age 
  • NEMO/releases/r4.0/r4.0-HEAD/src/ICE/icethd_dh.F90

    r13589 r13642  
    139139      ! 
    140140      DO ji = 1, npti 
    141          zf_tt(ji)         = qcn_ice_bot_1d(ji) + qsb_ice_bot_1d(ji) + fhld_1d(ji)  
     141         zf_tt(ji)         = qcn_ice_bot_1d(ji) + qsb_ice_bot_1d(ji) + fhld_1d(ji) + qtr_ice_bot_1d(ji) * frq_m_1d(ji)  
    142142         zq_bot(ji)        = MAX( 0._wp, zf_tt(ji) * rdt_ice ) 
    143143      END DO 
  • NEMO/releases/r4.0/r4.0-HEAD/src/ICE/iceupdate.F90

    r13589 r13642  
    138138               !                                        solar flux transmitted thru the 1st level of the ocean (i.e. not used by sea-ice) 
    139139               qsr(ji,jj) = ( 1._wp - at_i_b(ji,jj) ) * qsr_oce(ji,jj) * ( 1._wp - frq_m(ji,jj) ) & 
    140                   !                                   + solar flux transmitted thru ice (also not used by sea-ice) 
    141                   &             + SUM( a_i_b(ji,jj,:) * qtr_ice_bot(ji,jj,:) ) 
     140                  !                                   + solar flux transmitted thru ice and the 1st ocean level (also not used by sea-ice) 
     141                  &             + SUM( a_i_b(ji,jj,:) * qtr_ice_bot(ji,jj,:) ) * ( 1._wp - frq_m(ji,jj) ) 
    142142               ! 
    143143            ELSE                                                       !-- cooling or no ice left 
Note: See TracChangeset for help on using the changeset viewer.