Changeset 13234
- Timestamp:
- 09/04/12 10:57:36 (13 years ago)
- Location:
- issm/trunk-jpl/test/Par
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/Par/Pig.par
r13233 r13234 2 2 3 3 %Geometry and observation 4 x =transpose(ncread('../Data/Pig.nc','x'));5 y =transpose(ncread('../Data/Pig.nc','y'));6 vx_obs =transpose(ncread('../Data/Pig.nc','vx_obs'));7 vy_obs =transpose(ncread('../Data/Pig.nc','vy_obs'));8 index =transpose(ncread('../Data/Pig.nc','index'));9 surface =transpose(ncread('../Data/Pig.nc','surface'));10 thickness =transpose(ncread('../Data/Pig.nc','thickness'));4 x = transpose(ncread('../Data/Pig.nc','x')); 5 y = transpose(ncread('../Data/Pig.nc','y')); 6 vx_obs = transpose(ncread('../Data/Pig.nc','vx_obs')); 7 vy_obs = transpose(ncread('../Data/Pig.nc','vy_obs')); 8 index = transpose(ncread('../Data/Pig.nc','index')); 9 surface = transpose(ncread('../Data/Pig.nc','surface')); 10 thickness = transpose(ncread('../Data/Pig.nc','thickness')); 11 11 md.inversion.vx_obs =InterpFromMeshToMesh2d(index,x,y,vx_obs,md.mesh.x,md.mesh.y); 12 12 md.inversion.vy_obs =InterpFromMeshToMesh2d(index,x,y,vy_obs,md.mesh.x,md.mesh.y); -
issm/trunk-jpl/test/Par/SquareSheetConstrained.par
r13233 r13234 11 11 12 12 %Initial velocity 13 x =transpose(ncread('../Data/SquareSheetConstrained.nc','x'));14 y =transpose(ncread('../Data/SquareSheetConstrained.nc','y'));15 vx =transpose(ncread('../Data/SquareSheetConstrained.nc','vx'));16 vy =transpose(ncread('../Data/SquareSheetConstrained.nc','vy'));17 index =transpose(ncread('../Data/SquareSheetConstrained.nc','index'));13 x = transpose(ncread('../Data/SquareSheetConstrained.nc','x')); 14 y = transpose(ncread('../Data/SquareSheetConstrained.nc','y')); 15 vx = transpose(ncread('../Data/SquareSheetConstrained.nc','vx')); 16 vy = transpose(ncread('../Data/SquareSheetConstrained.nc','vy')); 17 index = transpose(ncread('../Data/SquareSheetConstrained.nc','index')); 18 18 md.initialization.vx=InterpFromMeshToMesh2d(index,x,y,vx,md.mesh.x,md.mesh.y); 19 19 md.initialization.vy=InterpFromMeshToMesh2d(index,x,y,vy,md.mesh.x,md.mesh.y); -
issm/trunk-jpl/test/Par/SquareSheetShelf.par
r13233 r13234 14 14 15 15 %Initial velocity 16 x =transpose(ncread('../Data/SquareSheetShelf.nc','x'));17 y =transpose(ncread('../Data/SquareSheetShelf.nc','y'));18 vx =transpose(ncread('../Data/SquareSheetShelf.nc','vx'));19 vy =transpose(ncread('../Data/SquareSheetShelf.nc','vy'));20 index =transpose(ncread('../Data/SquareSheetShelf.nc','index'));16 x = transpose(ncread('../Data/SquareSheetShelf.nc','x')); 17 y = transpose(ncread('../Data/SquareSheetShelf.nc','y')); 18 vx = transpose(ncread('../Data/SquareSheetShelf.nc','vx')); 19 vy = transpose(ncread('../Data/SquareSheetShelf.nc','vy')); 20 index = transpose(ncread('../Data/SquareSheetShelf.nc','index')); 21 21 md.initialization.vx=InterpFromMeshToMesh2d(index,x,y,vx,md.mesh.x,md.mesh.y); 22 22 md.initialization.vy=InterpFromMeshToMesh2d(index,x,y,vy,md.mesh.x,md.mesh.y); -
issm/trunk-jpl/test/Par/SquareShelf.par
r13233 r13234 11 11 12 12 %Initial velocity and pressure 13 x =transpose(ncread('../Data/SquareShelf.nc','x'));14 y =transpose(ncread('../Data/SquareShelf.nc','y'));15 vx =transpose(ncread('../Data/SquareShelf.nc','vx'));16 vy =transpose(ncread('../Data/SquareShelf.nc','vy'));17 index =transpose(ncread('../Data/SquareShelf.nc','index'));13 x = transpose(ncread('../Data/SquareShelf.nc','x')); 14 y = transpose(ncread('../Data/SquareShelf.nc','y')); 15 vx = transpose(ncread('../Data/SquareShelf.nc','vx')); 16 vy = transpose(ncread('../Data/SquareShelf.nc','vy')); 17 index = transpose(ncread('../Data/SquareShelf.nc','index')); 18 18 md.initialization.vx=InterpFromMeshToMesh2d(index,x,y,vx,md.mesh.x,md.mesh.y); 19 19 md.initialization.vy=InterpFromMeshToMesh2d(index,x,y,vy,md.mesh.x,md.mesh.y); -
issm/trunk-jpl/test/Par/SquareShelfConstrained.par
r13233 r13234 11 11 12 12 %Initial velocity 13 x =transpose(ncread('../Data/SquareShelfConstrained.nc','x'));14 y =transpose(ncread('../Data/SquareShelfConstrained.nc','y'));15 vx =transpose(ncread('../Data/SquareShelfConstrained.nc','vx'));16 vy =transpose(ncread('../Data/SquareShelfConstrained.nc','vy'));17 index =transpose(ncread('../Data/SquareShelfConstrained.nc','index'));13 x = transpose(ncread('../Data/SquareShelfConstrained.nc','x')); 14 y = transpose(ncread('../Data/SquareShelfConstrained.nc','y')); 15 vx = transpose(ncread('../Data/SquareShelfConstrained.nc','vx')); 16 vy = transpose(ncread('../Data/SquareShelfConstrained.nc','vy')); 17 index = transpose(ncread('../Data/SquareShelfConstrained.nc','index')); 18 18 md.initialization.vx=InterpFromMeshToMesh2d(index,x,y,vx,md.mesh.x,md.mesh.y); 19 19 md.initialization.vy=InterpFromMeshToMesh2d(index,x,y,vy,md.mesh.x,md.mesh.y);
Note:
See TracChangeset
for help on using the changeset viewer.