Changeset 5474


Ignore:
Timestamp:
2015-06-24T15:50:59+02:00 (5 years ago)
Author:
clem
Message:

correct limthd_dif to avoid floating point exception. see ticket #1547

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/LIM_SRC_3/limthd_dif.F90

    r5407 r5474  
    756756 
    757757      ! --- diagnose the change in non-solar flux due to surface temperature change --- ! 
    758       IF ( ln_it_qnsice ) hfx_err_dif_1d(:) = hfx_err_dif_1d(:) - ( qns_ice_1d(:)  - zqns_ice_b(:) ) * a_i_1d(:)  
     758      IF ( ln_it_qnsice ) THEN 
     759         DO ji = kideb, kiut 
     760            hfx_err_dif_1d(ji) = hfx_err_dif_1d(ji) - ( qns_ice_1d(ji)  - zqns_ice_b(ji) ) * a_i_1d(ji)  
     761         END DO 
     762      END IF 
    759763 
    760764      ! --- diag conservation imbalance on heat diffusion - PART 2 --- ! 
Note: See TracChangeset for help on using the changeset viewer.