Changeset 21829
- Timestamp:
- 07/20/17 12:51:28 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/classes/FemModel.cpp
r21822 r21829 2419 2419 element->GetVerticesCoordinates(&xyz_list); 2420 2420 for(int i=0;i<numvertices;i++){ 2421 x = *(xyz_list+3*i+0);2422 y = *(xyz_list+3*i+1);2421 x = xyz_list[3*i+0]; 2422 y = xyz_list[3*i+1]; 2423 2423 2424 2424 bx=-150.-728.8*std::pow(x/300000.,2)+343.91*std::pow(x/300000.,4)-50.57*std::pow(x/300000.,6); 2425 2425 by=500./(1.+std::exp((-2./4000.)*(y-80000./2.-24000.)))+500./(1.+std::exp((2./4000.)*(y-80000./2.+24000.))); 2426 2426 2427 r[i]=std::max(bx+by,-720.); 2427 r[i] = -720.; 2428 if(bx+by>-720.) r[i] = bx+by; 2428 2429 } 2429 2430 … … 2437 2438 /*Delete*/ 2438 2439 xDelete<IssmDouble>(r); 2439 2440 return;2441 2440 } 2442 2441 /*}}}*/
Note:
See TracChangeset
for help on using the changeset viewer.