Changeset 19682

Timestamp:

10/28/15 12:48:57 (9 years ago)

Author:

Eric.Larour

Message:

CHG: added esmfbinders to return surface output.

File:

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

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

@@ -9 +9 @@
 const int GCMForcingNumTerms = 1;
 const int GCMForcingTerms[GCMForcingNumTerms]= { SMBgcmEnum}; 
+
+const int ISSMOutputNumTerms = 1;
+const int ISSMOutputTerms[ISSMOutputNumTerms]= { SurfaceEnum };
 
 extern "C" {
@@ -96 +99 @@
                 /*}}}*/
 
+                /*Retrieve ISSM outputs and pass them back to the Gcm : {{{*/
+                for (int f=0;f<ISSMOutputNumTerms;f++){
+
+                        int output_type=ISSMOutputTerms[f];
+
+                        for (int i=0;i<femmodel->elements->Size();i++){
+                                Element* element=dynamic_cast<Element*>(femmodel->elements->GetObjectByOffset(i));
+
+                                switch(output_type){
+                                        case SurfaceEnum:
+                                                /*{{{*/
+                                                {
+
+                                                IssmDouble surface;
+                                                
+                                                /*Recover surface from the ISSM element: */
+                                                Input* surface_input = element->GetInput(SurfaceEnum); _assert_(surface_input);
+                                                surface_input->GetInputAverage(&surface);
+                        
+                                                *(issmoutputs+f*numberofelements+i) = surface;
+
+                                                }
+                                                /*}}}*/
+                                                break; 
+                                        default: 
+                                                { _error_("Unknown output type " << output_type << "\n"); }
+                                                break;
+                                }
+                        }
+                }
+
+                /*}}}*/
+
                 /*Before running, setup the time interval: */
                 femmodel->parameters->FindParam(&start_time,TimesteppingStartTimeEnum);