Changeset 2874 for branches/2011/dev_r2802_UKMO8_cice/NEMOGCM/NEMO/OPA_SRC/LBC/mppini.F90

Timestamp:

2011-09-28T12:19:59+02:00 (13 years ago)

Author:

charris

Message:

Code for running NEMO with CICE (for fully coupled mode this should be used in combination with dev_r2802_UKMO8_sbccpl). Changes are described briefly below.

physct: Constants modified to be consistent with CICE

nemogcm / prtctl / mppini: Changes to NEMO decomposition (activated using key_nemocice_decomp) to produce 'square' options in CICE. Can run without this key / code but this requires a global gather / scatter in the NEMO-CICE coupling which gets very slow on large processors numbers.

sbc_ice: CICE options and arrays added

sbcmod: CICE option added, including calls for initialising and finalising CICE.

sbcblk_core: Make sure necessary forcing field are available for CICE

sbcice_cice: Main CICE coupling code.

File:

• branches/2011/dev_r2802_UKMO8_cice/NEMOGCM/NEMO/OPA_SRC/LBC/mppini.F90 (modified) (3 diffs)

branches/2011/dev_r2802_UKMO8_cice/NEMOGCM/NEMO/OPA_SRC/LBC/mppini.F90

@@ -152 +152 @@
 
       IF(  iresti == 0 )   iresti = jpni
+
+#if defined key_nemocice_decomp
+      ! In order to match CICE the size of domains in NEMO has to be changed
+      ! The last line of blocks (west) will have fewer points
+
+      DO jj = 1, jpnj
+         DO ji=1, jpni-1
+            ilcit(ji,jj) = jpi
+         END DO
+         ilcit(jpni,jj) = jpiglo - (jpni - 1) * (jpi - nreci)
+      END DO
+
+#else
+
       DO jj = 1, jpnj
          DO ji = 1, iresti
@@ -161 +175 @@
       END DO
       
+#endif
       IF( irestj == 0 )   irestj = jpnj
+
+#if defined key_nemocice_decomp
+      ! Same change to domains in North-South direction as in East-West. 
+      DO ji=1,jpni
+         DO jj=1,jpnj-1
+            ilcjt(ji,jj) = jpj
+         END DO
+         ilcjt(ji,jpnj) = jpjglo - (jpnj - 1) * (jpj - nrecj)
+      END DO
+
+#else
+
       DO ji = 1, jpni
          DO jj = 1, irestj
@@ -171 +198 @@
       END DO
       
+#endif
       IF(lwp) THEN
          WRITE(numout,*)