Changeset 4676


Ignore:
Timestamp:
2014-06-19T09:59:53+02:00 (7 years ago)
Author:
clem
Message:

LIM3 in NEMO3.4: redefine s_i and sm_i after restart in case of constant salinity (num_sal=1) to avoid unexpected crashing

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limrst.F90

    r4634 r4676  
    142142         CALL iom_rstput( iter, nitrst, numriw, znam , z2d ) 
    143143      END DO 
    144  
     144       
    145145      DO jl = 1, jpl  
    146146         WRITE(zchar,'(I1)') jl 
     
    520520      END DO 
    521521      ! 
     522      ! clem: I do not understand why the following IF is needed 
     523      !       I suspect something inconsistent in the main code with option num_sal=1 
     524      IF( num_sal == 1 ) THEN 
     525         DO jl = 1, jpl  
     526            sm_i(:,:,jl) = bulk_sal 
     527            DO jk = 1, nlay_i  
     528               s_i(:,:,jk,jl) = bulk_sal 
     529            END DO 
     530         END DO 
     531      ENDIF 
     532      ! 
    522533      !CALL iom_close( numrir ) !clem: closed in sbcice_lim.F90 
    523534      ! 
Note: See TracChangeset for help on using the changeset viewer.