Changeset 4676

Timestamp:

2014-06-19T09:59:53+02:00 (10 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:

• branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limrst.F90 (modified) (2 diffs)

branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limrst.F90

@@ -142 +142 @@
          CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
       END DO
-
+      
       DO jl = 1, jpl 
          WRITE(zchar,'(I1)') jl
@@ -520 +520 @@
       END DO
       !
+      ! clem: I do not understand why the following IF is needed
+      !       I suspect something inconsistent in the main code with option num_sal=1
+      IF( num_sal == 1 ) THEN
+         DO jl = 1, jpl 
+            sm_i(:,:,jl) = bulk_sal
+            DO jk = 1, nlay_i 
+               s_i(:,:,jk,jl) = bulk_sal
+            END DO
+         END DO
+      ENDIF
+      !
       !CALL iom_close( numrir ) !clem: closed in sbcice_lim.F90
       !