Ignore:
Timestamp:
2019-09-18T17:52:11+02:00 (14 months ago)
Author:
mathiot
Message:

ENHANCE-02_ISF: fix compatibility issue between rn_hmin and ice shelf (ticket #2142)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/ENHANCE-02_ISF_domcfg/src/domisf.F90

    r11308 r11568  
    122122      ! 1.1 ground ice shelf if water column less than rn_glhw_min m  
    123123      ! => set the grounding line position 
    124       WHERE( bathy(:,:) - risfdep(:,:) < rn_glhw_min )  
     124      WHERE( bathy(:,:) - risfdep(:,:) < rn_glhw_min .AND. risfdep(:,:) > 0.0_wp )  
    125125         bathy  (:,:) = 0._wp ; risfdep(:,:) = 0._wp 
    126126      END WHERE 
     
    130130      DO jj = 1, jpj 
    131131         DO ji = 1, jpi 
    132             IF ( bathy(ji,jj) - risfdep(ji,jj) < rn_isfhw_min ) THEN 
     132            IF ( bathy(ji,jj) - risfdep(ji,jj) < rn_isfhw_min .AND. risfdep(ji,jj) > 0.0_wp ) THEN 
    133133               risfdep(ji,jj) = bathy(ji,jj) - rn_isfhw_min 
    134134               ! sanity check on risfdep (if < zisfdep_min) 
     
    145145      IF (ln_isfchannel) THEN 
    146146         imask(:,:) = 0 
    147          WHERE ( bathy(:,:) > 0._wp .AND. risfdep(:,:) > 0._wp ) 
     147         WHERE ( bathy(:,:) > 0._wp ) 
    148148            imask(:,:) = 1 
    149149         END WHERE 
    150150         DO jj = 2, jpjm1 
    151151            DO ji = 2, jpim1 
    152                IF( imask(ji,jj) == 1 ) THEN 
     152               IF( imask(ji,jj) == 1 .AND. risfdep(ji,jj) > 0._wp) THEN 
    153153                  ! define connexion 
    154154                  imskip1=imask(ji,jj)*imask(ji+1,jj  )  ! 1 = connexion 
     
    507507               gdept_0(ji,jj,ik-1) = gdept_0(ji,jj,ik) - e3w_0(ji,jj,ik) 
    508508               gdepw_0(ji,jj,ik-1) = gdepw_0(ji,jj,ik) - e3t_0(ji,jj,ik-1) 
     509               IF ( ik > 2 ) THEN 
     510                  e3w_0  (ji,jj,ik-1) = gdept_0(ji,jj,ik-1) - gdept_0(ji,jj,ik-2) 
     511               ELSE 
     512                  e3w_0  (ji,jj,ik-1) = 2.0 * gdept_0(ji,jj,ik-1) 
     513               END IF 
    509514            ENDIF 
    510515         END DO 
Note: See TracChangeset for help on using the changeset viewer.