Changeset 25810

Timestamp:

12/02/20 19:58:02 (4 years ago)

Author:

Mathieu Morlighem

Message:

BUG: ESMF not even compiling

File:

• issm/trunk-jpl/src/c/main/esmfbinders.cpp (modified) (6 diffs)

issm/trunk-jpl/src/c/main/esmfbinders.cpp ¶

@@ -16 +16 @@
         FemModel *femmodel;
 
+        /*GEOS 5*/
         void InitializeISSM(int argc, char** argv, int** pelementsonlocalrank, int* pnumberofelements, ISSM_MPI_Comm comm_init){ /*{{{*/
 
@@ -70 +71 @@
                                                 /*{{{*/
                                                 {
-
                                                 /*Recover rho_ice: */
                                                 rho_ice=element->FindParam(MaterialsRhoIceEnum);
@@ -78 +78 @@
 
                                                 /*Recover smb forcing from the gcm forcings: */
-                                                IssmDouble smbforcing=*(gcmforcings+f*numberofelements+i); 
+                                                _error_("not implemented yet, line below may not work");
+                                                IssmDouble smbforcing=gcmforcings[f*numberofelements+element->Sid()];
 
                                                 /*Convert to SI. The smbforcing from GEOS-5 in kg/s, and we transform it into m/s: */
@@ -84 +85 @@
 
                                                 /*Add into the element as new forcing :*/
-                                                element->inputs->AddInput(new DoubleInput(SmbMassBalanceEnum,smbforcing));
-
+                                                element->AddInput(SmbMassBalanceEnum,&smbforcing,P0Enum);
                                                 }
                                                 /*}}}*/
@@ -116 +116 @@
                                                 Input* surface_input = element->GetInput(SurfaceEnum); _assert_(surface_input);
                                                 surface_input->GetInputAverage(&surface);
-
-                                                *(issmoutputs+f*numberofelements+i) = surface;
+                                                issmoutputs[f*numberofelements+element->Sid()] = surface;
 
                                                 }
@@ -154 +153 @@
         } /*}}}*/
 
+        /*FISOC*/
+
 }