| 61 | You can also read and write layouts to chombo files using |
| 62 | |
| 63 | {{{ |
| 64 | CALL CreateLayout(GmGlobal, layout) |
| 65 | layout%level = 0 |
| 66 | mB=layout%mB(MPI_ID,:,:) |
| 67 | ALLOCATE(data(mB(1,1):mB(1,2), mB(2,1):mB(2,2),mB(3,1):mB(3,2),4)) |
| 68 | CALL LoadFieldIntoLayout(layout, data, FieldID) |
| 69 | CALL StoreLayoutAsChombo(layout, data, FieldID, name) |
| 70 | }}} |
| 71 | |
| 72 | and |
| 73 | |
| 74 | {{{ |
| 75 | CALL CreateLayout(GmGlobal, layout) |
| 76 | layout%level = 0 |
| 77 | mB=layout%mB(MPI_ID,:,:) |
| 78 | ALLOCATE(data(mB(1,1):mB(1,2), mB(2,1):mB(2,2),mB(3,1):mB(3,2),4)) |
| 79 | CALL LoadLayoutFromChombo(layout, data, name) |
| 80 | CALL UnloadFieldFromLayout(layout, data, FieldID, lHydroPeriodic, levels(n)%gmbc) |
| 81 | }}} |
| 82 | |
| 83 | |
| 84 | |
| 85 | |