Changeset 869 for trunk/NEMO/LIM_SRC_3/limthd.F90

Timestamp:

2008-03-26T10:21:54+01:00 (16 years ago)

Author:

rblod

Message:

Parallelisation of LIM3. This commit seems to ensure the reproducibility mono/mpp. See ticket #77.

File:

• trunk/NEMO/LIM_SRC_3/limthd.F90 (modified) (5 diffs)

trunk/NEMO/LIM_SRC_3/limthd.F90

@@ -30 +30 @@
    USE limvar
    USE prtctl          ! Print control
+   USE lib_mpp
 
    IMPLICIT NONE
@@ -116 +117 @@
       WRITE(numout,*) '~~~~~~'
 
+      IF( numit == nstart  )   CALL lim_thd_sal_init  ! Initialization (first time-step only)
 !------------------------------------------------------------------------------!
 ! 1) Initialization of diagnostic variables                                    !
@@ -180 +182 @@
       ! 1.4) Compute global heat content
       !-----------------------------------
-      qt_i_in(:,:)  = 0.0
-      qt_s_in(:,:)  = 0.0
-      qt_i_fin(:,:)  = 0.0
-      qt_s_fin(:,:)  = 0.0
-      sum_fluxq(:,:) = 0.0
-      fatm(:,:) = 0.0
+      qt_i_in(:,:)  = 0.e0
+      qt_s_in(:,:)  = 0.e0
+      qt_i_fin(:,:)  = 0.e0
+      qt_s_fin(:,:)  = 0.e0
+      sum_fluxq(:,:) = 0.e0
+      fatm(:,:) = 0.e0
 
 ! 2) Partial computation of forcing for the thermodynamic sea ice model.      !
@@ -279 +281 @@
 !------------------------------------------------------------------------------!
 
+         IF( lk_mpp ) CALL mpp_ini_ice(nbpb)
+
          IF (nbpb > 0) THEN  ! If there is no ice, do nothing.
 
@@ -414 +418 @@
             !+++++
 
+         IF( lk_mpp ) CALL mpp_comm_free(ncomm_ice) !RB necessary ??
          ENDIF ! nbpb