Ignore:
Timestamp:
2019-08-08T16:02:49+02:00 (15 months ago)
Author:
mathiot
Message:

ENHANCE-02_ISF_nemo : add UKESM ice sheet coupling method (ticket #2142)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/ISF/isfcavmlt.F90

    r11403 r11423  
    198198            ! solve the 2nd order equation to find zsfrz 
    199199            zaqe  = risf_lamb1 * (zeps1 + zeps3) 
    200             zaqer = 0.5_wp/MIN(zaqe,-risf_eps) 
     200            zaqer = 0.5_wp/MIN(zaqe,-zeps) 
    201201            zbqe  = zeps1*zeps6+zeps3*zeps7-zeps2 
    202202            zcqe  = zeps2*pstbl(ji,jj) 
     
    248248      REAL(wp), DIMENSION(jpi,jpj), INTENT(in   ) :: pstbl               ! salinity in tbl 
    249249      !!-------------------------------------------------------------------- 
    250       REAL(wp)                     :: zfwf_fld, zfwf_cpl                 ! total fwf in the forcing fields (pattern) and from the cpl interface (amount) 
     250      REAL(wp)                     :: zfwf_fld, zfwf_oasis               ! total fwf in the forcing fields (pattern) and from the oasis interface (amount) 
    251251      REAL(wp), DIMENSION(jpi,jpj) :: ztfrz                              ! tbl freezing temperature 
    252252      REAL(wp), DIMENSION(jpi,jpj) :: zfwf                               ! 2d fwf map after scaling 
     
    266266      ! 
    267267      ! compute glob sum from atm->oce ice shelf fwf 
    268       zfwf_cpl = glob_sum('isfcav_mlt', fwfisf_cpl(:,:)) 
     268      zfwf_oasis = glob_sum('isfcav_mlt', fwfisf_oasis(:,:)) 
    269269      ! 
    270270      ! scale fwf 
    271       zfwf(:,:) = zfwf(:,:) * zfwf_cpl / zfwf_fld 
     271      zfwf(:,:) = zfwf(:,:) * zfwf_oasis / zfwf_fld 
    272272      !  
    273273      ! define fwf and qoce 
     
    281281   !SUBROUTINE isfmlt_3eq_frz_ktm1 
    282282   ! compute tfrz based on sfrz value at kt-1 (need to be SAVED local array) 
     283   ! => should reduce error due to linarisation 
    283284   ! compute qfwf (eq 24) 
    284285   ! compute zqoce, zqlat, zqcon, zqhc 
Note: See TracChangeset for help on using the changeset viewer.