Changeset 10349


Ignore:
Timestamp:
2018-11-21T14:59:27+01:00 (17 months ago)
Author:
jchanut
Message:

Fixes volume conservation with isf, #2168

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/OCE/DYN/dynnxt.F90

    r10068 r10349  
    246246                  DO jj = 1, jpj 
    247247                     DO ji = 1, jpi 
    248                         IF( misfkt(ji,jj) <= jk .AND. jk <= misfkb(ji,jj) ) THEN 
    249                            e3t_b(ji,jj,jk) =   e3t_b(ji,jj,jk) - zcoef * ( fwfisf_b(ji,jj) - fwfisf(ji,jj) ) & 
     248                        IF( misfkt(ji,jj) <=jk .and. jk < misfkb(ji,jj) ) THEN 
     249                           e3t_b(ji,jj,jk) = e3t_b(ji,jj,jk) - zcoef * ( fwfisf_b(ji,jj) - fwfisf(ji,jj) ) & 
    250250                                &          * ( e3t_n(ji,jj,jk) * r1_hisf_tbl(ji,jj) ) * tmask(ji,jj,jk) 
     251                        ELSEIF ( jk==misfkb(ji,jj) ) THEN 
     252                           e3t_b(ji,jj,jk) = e3t_b(ji,jj,jk) - zcoef * ( fwfisf_b(ji,jj) - fwfisf(ji,jj) ) & 
     253                                &          * ( e3t_n(ji,jj,jk) * r1_hisf_tbl(ji,jj) ) * ralpha(ji,jj) * tmask(ji,jj,jk) 
    251254                        ENDIF 
    252255                     END DO 
Note: See TracChangeset for help on using the changeset viewer.