Ignore:
Timestamp:
2017-09-01T18:15:45+02:00 (3 years ago)
Author:
clem
Message:

debug icevar.F90

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icecor.F90

    r8486 r8491  
    124124            DO jj = 1, jpj  
    125125               DO ji = 1, jpi 
    126                   zsal = smv_i(ji,jj,jl) 
    127                   smv_i(ji,jj,jl) = MIN(  MAX( rn_simin*v_i(ji,jj,jl) , smv_i(ji,jj,jl) ) , rn_simax*v_i(ji,jj,jl)  ) 
    128                   ! associated salt flux 
    129                   sfx_res(ji,jj) = sfx_res(ji,jj) - ( smv_i(ji,jj,jl) - zsal ) * zzc 
     126                  IF( v_i(ji,jj,jl) > 0._wp ) THEN   ! clem: useless IF ??? 
     127                     zsal = smv_i(ji,jj,jl) 
     128                     smv_i(ji,jj,jl) = MIN(  MAX( rn_simin*v_i(ji,jj,jl) , smv_i(ji,jj,jl) ) , rn_simax*v_i(ji,jj,jl)  ) 
     129                     ! associated salt flux 
     130                     sfx_res(ji,jj) = sfx_res(ji,jj) - ( smv_i(ji,jj,jl) - zsal ) * zzc 
     131                  ENDIF 
    130132               END DO 
    131133            END DO 
Note: See TracChangeset for help on using the changeset viewer.