Changeset 22638

Timestamp:

03/27/18 20:48:00 (7 years ago)

Author:

seroussi

Message:

CHG: fixed coupling bug

File:

• issm/trunk-jpl/src/c/cores/transient_core.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/cores/transient_core.cpp

@@ -105 +105 @@
                         _printf_("ny size : "<< oceannysize<<"\n");
                         oceangridx = xNew<IssmDouble>(oceannxsize*oceannysize);
-                        ISSM_MPI_Recv(&oceangridx,oceannxsize*oceannysize,ISSM_MPI_DOUBLE,0,10001005,tomitgcmcomm,&status);
+                        ISSM_MPI_Recv(oceangridx,oceannxsize*oceannysize,ISSM_MPI_DOUBLE,0,10001005,tomitgcmcomm,&status);
                         _printf_("gridx received \n");
                         oceangridy = xNew<IssmDouble>(oceannxsize*oceannysize);
-                        ISSM_MPI_Recv(&oceangridy,oceannxsize*oceannysize,ISSM_MPI_DOUBLE,0,10001006,tomitgcmcomm,&status);
+                        ISSM_MPI_Recv(oceangridy,oceannxsize*oceannysize,ISSM_MPI_DOUBLE,0,10001006,tomitgcmcomm,&status);
                         _printf_("gridy received \n");
-                        //icebase = xNew<IssmDouble>(oceannxsize*oceannysize);
-                        //ISSM_MPI_Send(&icebase,oceannxsize*oceannysize,ISSM_MPI_DOUBLE,0,10001008,tomitgcmcomm);
+                        icebase = xNew<IssmDouble>(oceannxsize*oceannysize);
+                        for(int i=0;i<oceannxsize;i++){
+                                for(int j=0;j<oceannysize;j++){
+                                        icebase[i*oceannysize+j]=2*oceangridx[i*oceannysize+j];
+                                }
+                        }
+                        ISSM_MPI_Send(icebase,oceannxsize*oceannysize,ISSM_MPI_DOUBLE,0,10001008,tomitgcmcomm);
+                        _printf_("base sent\n");
                         //oceanmelt = xNew<IssmDouble>(oceannxsize*oceannysize);
                         //ISSM_MPI_Recv(&oceanmelt,oceannxsize*oceannysize,ISSM_MPI_DOUBLE,0,10001007,tomitgcmcomm,&status);