Ignore:
Timestamp:
2019-11-13T13:27:19+01:00 (11 months ago)
Author:
cguiavarch
Message:

Change to add basal melt only if the SST is above the freezing point

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/UKMO/NEMO_4.0.1_ICB_melting_temperature/src/OCE/ICB/icbutl.F90

    r11715 r11897  
    7474      ff_e(1:jpi,1:jpj) = ff_f (:,:)  
    7575      tt_e(1:jpi,1:jpj) = sst_m(:,:) 
     76      ss_e(1:jpi,1:jpj) = sss_m(:,:) 
    7677      fr_e(1:jpi,1:jpj) = fr_i (:,:) 
    7778      ua_e(1:jpi,1:jpj) = utau (:,:) * umask(:,:,1) ! maybe mask useless because mask applied in sbcblk 
     
    8586      CALL lbc_lnk_icb( 'icbutl', fr_e, 'T', +1._wp, 1, 1 ) 
    8687      CALL lbc_lnk_icb( 'icbutl', tt_e, 'T', +1._wp, 1, 1 ) 
     88      CALL lbc_lnk_icb( 'icbutl', ss_e, 'T', +1._wp, 1, 1 ) 
    8789#if defined key_si3 
    8890      hi_e(1:jpi, 1:jpj) = hm_i (:,:)   
     
    107109   SUBROUTINE icb_utl_interp( pi, pe1, puo, pui, pua, pssh_i,   & 
    108110      &                       pj, pe2, pvo, pvi, pva, pssh_j,   & 
    109       &                       psst, pcn, phi, pff            ) 
     111      &                       psst, pcn, phi, pff, psss        ) 
    110112      !!---------------------------------------------------------------------- 
    111113      !!                  ***  ROUTINE icb_utl_interp  *** 
     
    128130      REAL(wp), INTENT(  out) ::   puo, pvo, pui, pvi, pua, pva   ! ocean, ice and wind speeds 
    129131      REAL(wp), INTENT(  out) ::   pssh_i, pssh_j                 ! ssh i- & j-gradients 
    130       REAL(wp), INTENT(  out) ::   psst, pcn, phi, pff            ! SST, ice concentration, ice thickness, Coriolis 
     132      REAL(wp), INTENT(  out) ::   psst, pcn, phi, pff, psss      ! SST, ice concentration, ice thickness, Coriolis, SSS 
    131133      ! 
    132134      REAL(wp) ::   zcd, zmod       ! local scalars 
     
    139141      pvo  = icb_utl_bilin_h( vo_e, pi, pj, 'V', .false.  ) 
    140142      psst = icb_utl_bilin_h( tt_e, pi, pj, 'T', .true.   )    ! SST 
     143      psss = icb_utl_bilin_h( ss_e, pi, pj, 'T', .true.   )    ! SSS 
    141144      pcn  = icb_utl_bilin_h( fr_e, pi, pj, 'T', .true.   )    ! ice concentration 
    142145      pff  = icb_utl_bilin_h( ff_e, pi, pj, 'F', .false.  )    ! Coriolis parameter 
Note: See TracChangeset for help on using the changeset viewer.