|
578 | 578 | "plt.colorbar(axp, location='right')\n", |
579 | 579 | "\n", |
580 | 580 | "# Define the corners of the subregion of interest\n", |
581 | | - "x1, x2 = 4500, 6000\n", |
582 | | - "y1, y2 = 800, 1400\n", |
| 581 | + "\n", |
| 582 | + "# For Agulhas Current, we can use the following coordinates\n", |
| 583 | + "# These coordinates are approximate and can be adjusted based on the specific area of interest\n", |
| 584 | + "# x1, x2 are longitude coordinates (in pixels)\n", |
| 585 | + "# y1, y2 are latitude coordinates (in pixels)\n", |
| 586 | + "x1, x2 = 800, 2400\n", |
| 587 | + "y1, y2 = 2500, 3500\n", |
| 588 | + "\n", |
| 589 | + "# For Strait of Gibraltar, we can use the following coordinates:\n", |
| 590 | + "x1, x2 = 400, 1600\n", |
| 591 | + "y1, y2 = 4500, 5500\n", |
| 592 | + "\n", |
583 | 593 | "\n", |
584 | 594 | "# Draw a red rectangle around the selected region\n", |
585 | 595 | "regionRect = plt.Rectangle(\n", |
|
686 | 696 | "\n", |
687 | 697 | "- Load the `eastwest_ocean_velocity_u` dataset \n", |
688 | 698 | "- Select a specific **time step**, **resolution**, and **region** covering the Agulhas leakage area \n", |
689 | | - "- Visualize the **east-west velocity component** of ocean currents \n", |
| 699 | + "- Visualize the **east-west velocity component** of ocean currents with proper colormap\n", |
690 | 700 | "- Compute and report:\n", |
691 | 701 | " - The **maximum and minimum velocity values**\n", |
692 | 702 | " - The **structure of eddies** and current boundaries in the region\n", |
|
0 commit comments