Changeset 18786 for issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateOutputDefinitions.cpp

Timestamp:

11/14/14 16:42:09 (10 years ago)

Author:

Eric.Larour

Message:

NEW: introduced new Masscon output definition, relying on a levelset definition of the masscon shape.

File:

• issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateOutputDefinitions.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateOutputDefinitions.cpp

@@ -114 +114 @@
                                 /*}}}*/
                         }
+                        else if (output_definition_enums[i]==MassconEnum){
+                                /*Deal with masscons: {{{*/
+                                
+                                /*masscon variables: */
+                                int          nummasscons;
+                                char**       masscon_name_s             = NULL;    
+                                IssmDouble** masscon_levelset_s           = NULL;
+                                int*         masscon_levelset_M_s    = NULL;
+                                int*         masscon_levelset_N_s    = NULL;
 
+                                /*Fetch name and levelset, etc ... (see src/m/classes/masscon.m): */
+                                iomodel->FetchMultipleData(&masscon_name_s,&nummasscons,MassconNameEnum);
+                                iomodel->FetchMultipleData(&masscon_levelset_s,&masscon_levelset_M_s,&masscon_levelset_N_s,&nummasscons,MassconLevelsetEnum);
+                                for(j=0;j<nummasscons;j++){
+
+                                        /*Create a masscon object: */
+                                        output_definitions->AddObject(new Masscon(masscon_name_s[j],masscon_levelset_s[j],masscon_levelset_M_s[j]));
+
+                                }
+
+                                /*Free ressources:*/
+                                for(j=0;j<nummasscons;j++){
+                                        char* string=NULL;
+                                        IssmDouble* matrix = NULL;
+
+                                        string = masscon_name_s[j];    xDelete<char>(string);
+                                        matrix = masscon_levelset_s[j]; xDelete<IssmDouble>(matrix);
+                                }
+                                xDelete<char*>(masscon_name_s);
+                                xDelete<IssmDouble*>(masscon_levelset_s);
+                                xDelete<int>(masscon_levelset_M_s);
+                                xDelete<int>(masscon_levelset_N_s);
+                                /*}}}*/
+                        }
                         else _error_("output definition enum " << output_definition_enums[i] << "not supported yet!");
                 }