Changeset 22376

Timestamp:

01/29/18 10:05:44 (7 years ago)

Author:

seroussi

Message:

NEW: implementing coupling with MITgcm

File:

• issm/trunk-jpl/src/c/main/issm_ocean.cpp (modified) (3 diffs)

issm/trunk-jpl/src/c/main/issm_ocean.cpp

@@ -11 +11 @@
         int    icecommsize;
         int    rankzeros[2];
-        int    my_rank;
+        int    my_rank,my_local_rank,my_size,my_local_size;
         ISSM_MPI_Comm worldcomm;
-
-//      int    modelid;
-//      int    count=0;
         ISSM_MPI_Comm modelcomm;
         ISSM_MPI_Comm frommitgcm;
@@ -28 +25 @@
         /*What is my rank?:*/
         ISSM_MPI_Comm_rank(worldcomm,&my_rank);
+        ISSM_MPI_Comm_size(worldcomm,&my_size);
 
         /*First model is ice, second is ocean*/
@@ -37 +35 @@
         rankzeros[1]=icecommsize; 
 
+        printf("Ice: My global rank: %i My world size: %i \n",my_rank,my_size);
+
         /*Split world into sub-communicators for each and every model:*/
         ISSM_MPI_Comm_split(worldcomm,0, my_rank, &modelcomm);
+
+        ISSM_MPI_Comm_rank(modelcomm,&my_local_rank);
+        ISSM_MPI_Comm_size(modelcomm,&my_local_size);
+
+        printf("Ice: My global rank: %i My local rank: %i My world size: %i My local size: %i\n",my_rank,my_local_rank,my_size,my_local_size);
 
         ISSM_MPI_Intercomm_create( modelcomm, 0, worldcomm, rankzeros[1], 0, &tomitgcmcomm);