Changeset 7398 for branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/OPA_SRC/IOM/iom.F90

Timestamp:

2016-11-30T15:44:11+01:00 (8 years ago)

Author:

cbricaud

Message:

coarsening branch: first implementation of coarsening in PISCES

File:

• branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/OPA_SRC/IOM/iom.F90 (modified) (3 diffs)

branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/OPA_SRC/IOM/iom.F90

@@ -50 +50 @@
    USE ioipsl, ONLY :  ju2ymds    ! for calendar
    USE crs             ! Grid coarsening
+   USE crslbclnk
 
    IMPLICIT NONE
@@ -823 +824 @@
                         jstartrow = MAX(1,jstartrow)
                      ENDIF
-                     istart(1:2) = (/ mig(1), mjg(1) + jstartrow - 1 /)  ! icnt(1:2) done below
+                     IF( lk_crs .AND. jpi == jpi_crs .AND. jpi .NE. 0  )THEN
+                        istart(1:2) = (/ mig_crs(1), mjg_crs(1) + jstartrow - 1 /)  ! icnt(1:2) done below
+                     ELSE
+                        istart(1:2) = (/ mig(1), mjg(1) + jstartrow - 1 /)  ! icnt(1:2) done below
+                     ENDIF
                   ELSEIF( idom == jpdom_global  ) THEN ; istart(1:2) = (/ nimpp , njmpp  /)  ! icnt(1:2) done below
                   ENDIF
@@ -922 +927 @@
             !--- overlap areas and extra hallows (mpp)
             IF(     PRESENT(pv_r2d) .AND. idom /= jpdom_unknown ) THEN
-               CALL lbc_lnk( pv_r2d,'Z',-999.,'no0' )
+               IF( lk_crs .AND. jpi == jpi_crs .AND. jpi .NE. 0 )THEN
+                  CALL crs_lbc_lnk( pv_r2d,'Z',-999.,'no0')
+               ELSE
+                  CALL lbc_lnk( pv_r2d,'Z',-999.,'no0' )
+               ENDIF
             ELSEIF( PRESENT(pv_r3d) .AND. idom /= jpdom_unknown ) THEN
                ! this if could be simplified with the new lbc_lnk that works with any size of the 3rd dimension
                IF( icnt(3) == jpk ) THEN
-                  CALL lbc_lnk( pv_r3d,'Z',-999.,'no0' )
+                  IF( lk_crs .AND. jpi == jpi_crs .AND. jpi .NE. 0 )THEN
+                     CALL crs_lbc_lnk( pv_r3d,'Z',-999.,'no0')
+                  ELSE
+                     CALL lbc_lnk( pv_r3d,'Z',-999.,'no0' )
+                  ENDIF
                ELSE   ! put some arbitrary value (a call to lbc_lnk will be done later...)
                   DO jj = nlcj+1, jpj   ;   pv_r3d(1:nlci, jj, :) = pv_r3d(1:nlci, nlej, :)   ;   END DO