Changeset 5802 for branches/NERC/dev_r5589_is_oce_cpl/NEMOGCM/NEMO/OPA_SRC/lib_fortran.F90

Timestamp:

2015-10-16T19:19:43+02:00 (9 years ago)

Author:

mathiot

Message:

ice sheet coupling: reproducibility fixed in MISOMIP configuration + contrain bathy to be constant during all the run

File:

• branches/NERC/dev_r5589_is_oce_cpl/NEMOGCM/NEMO/OPA_SRC/lib_fortran.F90 (modified) (4 diffs)

branches/NERC/dev_r5589_is_oce_cpl/NEMOGCM/NEMO/OPA_SRC/lib_fortran.F90

@@ -172 +172 @@
       !!-----------------------------------------------------------------------
       !
-      glob_sum_full_2d = SUM( ptab(:,:) )
+      glob_sum_full_2d = SUM( ptab(:,:) * tmask_h(:,:) )
       IF( lk_mpp )   CALL mpp_sum( glob_sum_full_2d )
       !
@@ -194 +194 @@
       glob_sum_full_3d = 0.e0
       DO jk = 1, ijpk
-         glob_sum_full_3d = glob_sum_full_3d + SUM( ptab(:,:,jk) )
+         glob_sum_full_3d = glob_sum_full_3d + SUM( ptab(:,:,jk) * tmask_h(:,:) )
       END DO
       IF( lk_mpp )   CALL mpp_sum( glob_sum_full_3d )
@@ -376 +376 @@
       DO jj = 1, jpj
          DO ji =1, jpi
-         ztmp =  ptab(ji,jj)
+         ztmp =  ptab(ji,jj) * tmask_h(ji,jj) 
          CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp )
          END DO
@@ -407 +407 @@
          DO jj = 1, jpj
             DO ji =1, jpi
-            ztmp =  ptab(ji,jj,jk)
+            ztmp =  ptab(ji,jj,jk) * tmask_h(ji,jj)
             CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp )
             END DO