Changeset 20690 for issm/trunk-jpl/src/c/analyses/SealevelriseAnalysis.cpp
- Timestamp:
- 06/03/16 10:05:11 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/analyses/SealevelriseAnalysis.cpp
r20345 r20690 31 31 32 32 /*Create inputs: */ 33 iomodel->FetchDataToInput(elements, MaskIceLevelsetEnum);34 iomodel->FetchDataToInput(elements, MaskOceanLevelsetEnum);35 iomodel->FetchDataToInput(elements, MaskLandLevelsetEnum);36 iomodel->FetchDataToInput(elements, SealevelriseDeltathicknessEnum);37 iomodel->FetchDataToInput(elements, SealevelEnum,0);33 iomodel->FetchDataToInput(elements,"md.mask.ice_levelset",MaskIceLevelsetEnum); 34 iomodel->FetchDataToInput(elements,"md.mask.ocean_levelset",MaskOceanLevelsetEnum); 35 iomodel->FetchDataToInput(elements,"md.mask.land_levelset",MaskLandLevelsetEnum); 36 iomodel->FetchDataToInput(elements,"md.sealevelrise.deltathickness",SealevelriseDeltathicknessEnum); 37 iomodel->FetchDataToInput(elements,"md.slr.sealevel",SealevelEnum,0); 38 38 39 39 }/*}}}*/ … … 60 60 61 61 /*some constant parameters: */ 62 parameters->AddObject(iomodel->CopyConstantObject( SealevelriseReltolEnum));63 parameters->AddObject(iomodel->CopyConstantObject( SealevelriseAbstolEnum));64 parameters->AddObject(iomodel->CopyConstantObject( SealevelriseMaxiterEnum));65 parameters->AddObject(iomodel->CopyConstantObject( SealevelriseRigidEnum));66 parameters->AddObject(iomodel->CopyConstantObject( SealevelriseElasticEnum));67 68 iomodel->FetchData(&elastic, SealevelriseElasticEnum);62 parameters->AddObject(iomodel->CopyConstantObject("md.slr.reltol",SealevelriseReltolEnum)); 63 parameters->AddObject(iomodel->CopyConstantObject("md.slr.abstol",SealevelriseAbstolEnum)); 64 parameters->AddObject(iomodel->CopyConstantObject("md.slr.maxiter",SealevelriseMaxiterEnum)); 65 parameters->AddObject(iomodel->CopyConstantObject("md.slr.rigid",SealevelriseRigidEnum)); 66 parameters->AddObject(iomodel->CopyConstantObject("md.slr.elastic",SealevelriseElasticEnum)); 67 68 iomodel->FetchData(&elastic,"md.slr.elastic"); 69 69 if(elastic){ 70 70 71 71 /*love numbers: */ 72 iomodel->FetchData(&love_h,&nl,NULL, SealevelriseLoveHEnum);73 iomodel->FetchData(&love_k,&nl,NULL, SealevelriseLoveKEnum);72 iomodel->FetchData(&love_h,&nl,NULL,"md.slr.love_h"); 73 iomodel->FetchData(&love_k,&nl,NULL,"md.slr.love_k"); 74 74 75 75 /*compute elastic green function for a range of angles*/ 76 iomodel->FetchData(°acc, SealevelriseDegaccEnum);76 iomodel->FetchData(°acc,"md.slr.degacc"); 77 77 M=reCast<int,IssmDouble>(180./degacc+1.); 78 78 G_elastic=xNew<IssmDouble>(M); … … 133 133 134 134 /*Transitions: */ 135 iomodel->FetchData(&transitions,&transitions_M,&transitions_N,&ntransitions, SealevelriseTransitionsEnum);135 iomodel->FetchData(&transitions,&transitions_M,&transitions_N,&ntransitions,"md.slr.transitions"); 136 136 if(transitions){ 137 137 parameters->AddObject(new DoubleMatArrayParam(SealevelriseTransitionsEnum,transitions,ntransitions,transitions_M,transitions_N)); … … 147 147 148 148 /*Requested outputs*/ 149 iomodel->F etchData(&requestedoutputs,&numoutputs,SealevelriseRequestedOutputsEnum);149 iomodel->FindConstant(&requestedoutputs,&numoutputs,"md.slr.requested_outputs"); 150 150 parameters->AddObject(new IntParam(SealevelriseNumRequestedOutputsEnum,numoutputs)); 151 151 if(numoutputs)parameters->AddObject(new StringArrayParam(SealevelriseRequestedOutputsEnum,requestedoutputs,numoutputs)); 152 iomodel->DeleteData(&requestedoutputs,numoutputs, SealevelriseRequestedOutputsEnum);152 iomodel->DeleteData(&requestedoutputs,numoutputs,"md.slr.requested_outputs"); 153 153 154 154
Note:
See TracChangeset
for help on using the changeset viewer.