{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": true }, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "%matplotlib inline" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "#### Increasing Wind Stress" ] }, { "cell_type": "raw", "metadata": {}, "source": [ "IF( kt <= nn_tau000 ) THEN ! Increase the stress to its nominal value \n", " ! ! during the first nn_tau000 time-steps\n", " zfact = 0.5 * ( 1. - COS( rpi * REAL( kt, wp ) / REAL( nn_tau000, wp ) ) )\n", " zcoef = 1. / ( zrhoa * zcdrag ) \n", " ztx = zfact * rn_utau0\n", " zty = zfact * rn_vtau0\n", " zmod = SQRT( ztx * ztx + zty * zty )\n", " utau(:,:) = ztx\n", " vtau(:,:) = zty\n", " taum(:,:) = zmod\n", " zmod = SQRT( zmod * zcoef )\n", " wndm(:,:) = zmod\n", " ENDIF" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "#### Decreasing Wind Stress" ] }, { "cell_type": "raw", "metadata": {}, "source": [ "nn_tau002 = start of decrease period\n", "nn_tau001 = length of decrease period\n", "\n", "IF(( nn_tau002 + nn_tau001)>= kt >= nn_tau002) THEN ! Decrease the stress to its nominal value \n", " ! ! during the first nn_tau000 time-steps\n", " zfact = -0.5 * ( 1. - COS( rpi * (REAL( kt, wp ) - nn_tau001) / REAL( nn_tau001, wp ) ) ) + 1\n", " zcoef = 1. / ( zrhoa * zcdrag ) \n", " ztx = zfact * rn_utau0\n", " zty = zfact * rn_vtau0\n", " zmod = SQRT( ztx * ztx + zty * zty )\n", " utau(:,:) = ztx\n", " vtau(:,:) = zty\n", " taum(:,:) = zmod\n", " zmod = SQRT( zmod * zcoef )\n", " wndm(:,:) = zmod\n", " ENDIF" ] }, { "cell_type": "code", "execution_count": 54, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD8CAYAAACMwORRAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl4VdW9//H3FxAHQBDRVAWFClZRRI2iKDUJCAJ6xXng\n59Rqeahya2/b6/D8br0dvW3tYK3Wua1Dr1yvIlJEASGJWqsFRBmcQJzAKiqVNk5M398f68RfjCQ5\nOTnnrH32+byeJ084Jzvh4zb5ZLPO2muZuyMiIunSKXYAERHJP5W7iEgKqdxFRFJI5S4ikkIqdxGR\nFFK5i4ikkMpdRCSFVO4iIimkchcRSaEusf7iPn36eP/+/XP63A8++IBu3brlN1AeJDUXJDebcrWP\ncrVPGnMtWrToXXffpc0D3T3KW2VlpeeqtrY2588tpKTmck9uNuVqH+VqnzTmAhZ6Fh2rYRkRkRRS\nuYuIpJDKXUQkhVTuIiIppHIXEUmhNsvdzH5nZmvNbFkLHzczu9bMVprZEjM7JP8xRUSkPbK5cv8D\nMLaVj48DBmXeJgE3dDyWiIh0RJs3Mbn7o2bWv5VDJgB3ZOZfPmlmvcxsN3f/W54ySg7c4Zln4Ikn\n4K23YP/9t4sdSdJo0ya48kr6vv9+eDxiBHSJdm+kNGGexR6qmXKf6e4HbOVjM4GfuPvjmcfzgMvc\nfeFWjp1EuLqnoqKicurUqTmFbmhooHv37jl9biElJddTT/Xmllu+yMsvhyxmzg9/+CRHHfVJ5GSf\nl5Rz1pxyZafL+vUcdfLJ2JYtAHy8yy68dvbZ/O3446FT/Jf0kna+GnUkV01NzSJ3P7TNA7O50wno\nDyxr4WMzgRFNHs8DDm3ra+oO1fz7+GP3SZPcwX3QIPebbnJ/4w33zZvjZ2uJcrVPInNt2eKPT5/u\nfu+97iNGhG/Amhr3tWtjJ0vm+fLSuUN1DdCvyeO+meekiD74AI47Dm6+GS67DJYtg0mToG/fRFxA\nSZqZsbFnTzjlFHj0Ubj1VvjLX+DII+H112OnK1v5+LGfAZybmTVzBLDeNd5eVBs3hp+r2lq4/Xb4\nyU+ga9fYqaQsmcEFF8D8+fDOOzBqFKxdGztVWcpmKuTdwF+AL5nZajO7wMwmm9nkzCGzgFXASuAW\n4KKCpZWtuuQSmD0bbrkFzj03dhoRYPhweOghWLMGTj4ZNmyInajsZDNb5qw2Pu7AxXlLJO3yxz/C\nDTfApZfCV78aO41IE8OHwx/+AGecAZdfDr/8ZexEZUWjsSXstdfg618Ps89+/OPYaUS24vTT4eKL\n4Ve/gkceiZ2mrKjcS5Q7TJ4c3t9xh6YWS4JdfTV86Uvwta+FV/6lKFTuJeree+Hhh+Gqq2DAgNhp\nRFqx/fZhBs2rr8IPfhA7TdlQuZegDz+Eb38bDj4YLtLL11IKRoyA888PwzMrVsROUxZU7iXommvg\njTfC+86dY6cRydJ//VeYo3v55bGTlAWVe4lZtw5++lM44QQ4+ujYaUTa4QtfCNO6pk2DBQtip0k9\nlXuJ+cUv4J//hB/9KHYSkRz8279Bnz5w5ZWxk6Seyr2ErFsHv/kNnHoqDBkSO41IDnr0gH//9zAb\n4K9/jZ0m1VTuJeT668NV+3/8R+wkIh1w0UXQu3cYg5eCUbmXiI8+gmuvDYuDHXhg7DQiHdC9O0yZ\nAtOnwwsvxE6TWir3EnHnnfDuu+FftCIlb8oU2G67MOVLCkLlXgLc4de/hkMO0QwZSYlddoGzzw7L\nmL73Xuw0qaRyLwG1tfDcc2H1R7PYaUTy5BvfgI8/httui50klVTuJeD662HnncMaTCKpMWQIVFWF\nZU0z2/RJ/qjcE+7NN+GBB8L+B9tpj2tJm69/Paw5M2dO7CSpo3JPuN//HjZvDgvqiaTOSSeF8feb\nb46dJHVU7gm2ZUsYjqypgYEDY6cRKYCuXeG88+BPf4K3346dJlVU7glWXw+vvBKGZERS64ILYNOm\nMN9X8kblnmC33w477hi2oBRJrX33hcMPD9/w7rHTpIbKPaE++CBsyHH66WGvA5FUO/98WLYMnn46\ndpLUULkn1P33h4I/99zYSUSK4PTTw/j7H/8YO0lqqNwT6q67YK+94KijYicRKYLevcPCSXffHaaH\nSYep3BNo7dqwUfzEidBJ/4ekXEycCG+9BXV1sZOkgqojge69N1y8nHVW7CQiRXTccWG99//+79hJ\nUkHlnkBTp8LgwdqQQ8rM9tvDhAlhG74NG2KnKXkq94RZswYefxzOOCN2EpEIzjgD3n8f5s6NnaTk\nqdwTZtq0MNX3tNNiJxGJYMwY6NkzjE1Kh6jcE+bee2H//WG//WInEYmga9cwNDN9uoZmOkjlniBv\nvw2PPRY2wBYpW6eeGoZmamtjJylpKvcEeeCBMCRzyimxk4hENHo0dOsW7uSTnGVV7mY21sxeNLOV\nZnb5Vj7e08z+ZGbPmtlyM/tK/qOm37RpsPfecMABsZOIRLTddmFa5PTpuqGpA9osdzPrDFwPjAMG\nA2eZ2eBmh10MPOfuQ4Fq4Bdm1jXPWVNt/XqYPz8sb62t9KTsnXhiGKd86qnYSUpWNlfuw4CV7r7K\n3TcAU4EJzY5xoIeZGdAdWAdsymvSlHvoIdi4MXxPi5S98eNhm23C1bvkJJty3wN4o8nj1ZnnmroO\n2A94E1gKXOLu2hSxHWbMCBvSHHFE7CQiCdCzJ1RXhxeiJCfmbayfbGanAmPd/cLM43OAw919SrNj\njgK+BewNzAWGuvs/mn2tScAkgIqKisqpU6fmFLqhoYHu3bvn9LmFlGuuTZuME088iqOPfodLL32x\nAMnSd84KTbnapxC5dr//fva59lqeuv12Ptpzz8TkyoeO5KqpqVnk7oe2eaC7t/oGDAdmN3l8BXBF\ns2MeBL7c5PF8YFhrX7eystJzVVtbm/PnFlKuuebNcwf3++/Pb56m0nbOCk252qcguV57Lfxg/PSn\nOX+JNJ4vYKG30dvuntWwzAJgkJkNyLxIeiYwo9kxrwOjAMysAvgSsCqLry3AzJmw7bZhBpiIZOy5\nJwwdCg8+GDtJSWqz3N19EzAFmA08D9zj7svNbLKZTc4c9kPgSDNbCswDLnP3dwsVOm1mzoSRI8PU\nXhFp4vjj4c9/hnXrYicpOVnNc3f3We6+j7vv7e4/zjx3o7vfmPnzm+4+xt2HuPsB7n5XIUOnyUsv\nwYoVYVqviDRz/PFhrvvs2bGTlBzdoRpZ4784Ve4iW3HYYdCnj4ZmcqByj2zWrLB2e//+sZOIJFDn\nzjB2LDz8sO5WbSeVe0QNDfDoo+F+DRFpwfjx8N57sGBB7CQlReUe0fz5YVXTceNiJxFJsDFjwmbC\nDz0UO0lJUblHNGsWdO8OI0bETiKSYDvvDMOGqdzbSeUeiXuYADBqVNifQERaMW4cLFwI72qGdbZU\n7pG89BK8+mp4rUhE2jB2bLgi0t6qWVO5R9I4bffYY+PmECkJlZVheObhh2MnKRkq90hmz4Z99oEB\nA2InESkBnTuH9TnmzAlX8NImlXsEn3wCdXVhEoCIZGnMGHjrLVi6NHaSkqByj+CJJ+DDD7VQmEi7\nNF4NzZkTN0eJULlHMGcOdOkCNTWxk4iUkD32CLdzq9yzonKPYO5cGD4cevSInUSkxIweDY89Bh99\nFDtJ4qnci+y99+DppzUkI5KT0aPh44/DMsDSKpV7kc2fH17sV7mL5KCqKoxpPvJI7CSJp3Ivsrlz\nw96/h7a9A6KINNe9exjT1M1MbVK5F9kjj4RN3bt0iZ1EpEQdcwwsXhzGOKVFKvciWrUKXnklfG+K\nSI6OOSaMbdbWxk6SaCr3Ipo3L7wfNSpuDpGSdthhYapZ4w+UbJXKvYjmzYPdd4d9942dRKSEbbNN\neGFV5d4qlXuRbNkSZsqMHAlmsdOIlLhRo8LO8m+8ETtJYqnci2T5cnjnnVDuItJBjWObunpvkcq9\nSObPD+9V7iJ5sP/+sMsuelG1FSr3Ipk/H/beG/baK3YSkRTo1CksztR4V6B8jsq9CDZtgvp6XbWL\n5NXIkbB6dRh7l89RuRfB4sWwfr3KXSSvGpdV1dDMVqnci6BxvF1L/Irk0aBBYRngxh8w+QyVexHU\n1oZlqCsqYicRSRGzcMVUW6tx961QuRfYxo3w+OO6ahcpiJqaMMf4uediJ0kclXuBLVgAH3ygchcp\nCI27t0jlXmCN33NHHx03h0gqDRgQ5her3D8nq3I3s7Fm9qKZrTSzy1s4ptrMnjGz5WZWn9+Ypauu\nDoYMCfdbiEgB1NSEucZbtsROkihtlruZdQauB8YBg4GzzGxws2N6Ab8FTnD3/YHTCpC15GzYEHYD\nq66OnUQkxaqrw9ruy5fHTpIo2Vy5DwNWuvsqd98ATAUmNDtmIjDN3V8HcPe1+Y1ZmhYsCPv4arxd\npICqqsJ7Dc18hnkbU4jM7FRgrLtfmHl8DnC4u09pcsw1wDbA/kAP4NfufsdWvtYkYBJARUVF5dSp\nU3MK3dDQQPfu3XP63EJqnuuuu/bkttu+yPTpj9Oz56aIyUrnnCWFcrVP7FyHn3UWDYMGsfwHP/jM\n87FztaQjuWpqaha5e9sbdbp7q2/AqcCtTR6fA1zX7JjrgCeBbkAfYAWwT2tft7Ky0nNVW1ub8+cW\nUvNcxxzjfuCBcbI0VyrnLCmUq32i5zr/fPfevd03b/7M09FztaAjuYCF3kZvu3tWwzJrgH5NHvfN\nPNfUamC2u3/g7u8CjwJDs/jaqaXxdpEiqq6Gdetg2bLYSRIjm3JfAAwyswFm1hU4E5jR7JgHgBFm\n1sXMdgAOB57Pb9TSsnBhGG9vHA4UkQJqvIqq10S9Rm2Wu7tvAqYAswmFfY+7LzezyWY2OXPM88DD\nwBLgr4RhnLL+FVpXF95rfrtIEey1V3hr/METumRzkLvPAmY1e+7GZo+vBq7OX7TS1ji/vU+f2ElE\nykR1NTz4YJjv3kn3Z+oMFMDGjfDEExqSESmqqip49114vqxHhD+lci+ARYvCejIqd5EiavyB07g7\noHIvCI23i0QwYAD07atx9wyVewHU18N++8Guu8ZOIlJGzMLVe3291ndH5Z53mzaF9ds1v10kgupq\nWLsWXnwxdpLoVO55tngxNDRovF0kCo27f0rlnmeN31MabxeJYOBA2G03lTsq97yrr4d99gnfXyJS\nZBp3/5TKPY82b4bHHtOQjEhUVVXw5pvw8suxk0Slcs+jVau6s369hmREotK4O6Byz6tnn+0J6Mpd\nJKp99w37WqrcJV+WLOnFgAHQr1/bx4pIgZiFfz6r3CUf3GHJkp4akhFJgqoqeP11tnvrrdhJolG5\n58lzz8H69V01JCOSBJkfxJ7PPhs5SDwq9zxp/Begyl0kAQ44AHbaiV5LlsROEo3KPU/q66FPn08Y\nMCB2EhGhUyf48pd15S4d4x7KfejQ9zGLnUZEADj6aHZYsybMeS9DKvc8WLEC3n47lLuIJETjGOmj\nj8bNEYnKPQ8ax9tV7iIJctBBbNphh7KdEqlyz4P6eqiogH79PoodRUQadenC+iFDVO6Sm8bx9qoq\nNN4ukjDvDx0a9lRduzZ2lKJTuXfQK6/A6tWaAimSROuHDg1/KMNxd5V7B2l+u0hy/XOffaBbt7Ic\nmlG5d1CY3w6DB8dOIiLNeZcucOSRKndpv/r6sEaRxttFEqqqCpYuhXffjZ2kqFTuHfD66/DqqxqS\nEUm0xh/Qxx6Lm6PIVO4doPF2kRJw2GGw/fZlNzSjcu+A+nrYaScYMiR2EhFp0bbbwvDhKnfJXl0d\nfPnLYY0iEUmwqip49ln4+99jJyka1VKOVq8O++9WV8dOIiJtqqoKdxw+/njsJEWTVbmb2Vgze9HM\nVprZ5a0cd5iZbTKzU/MXMZk03i5SQg4/PAzP1NXFTlI0bZa7mXUGrgfGAYOBs8zsc7O6M8f9FJiT\n75BJVFcHPXtC4w1wIpJg220HRxyhcm9mGLDS3Ve5+wZgKjBhK8f9K3AfUBaLODTOb+/cOXYSEclK\ndTUsXgzvl8fqrdmU+x7AG00er8489ykz2wM4Cbghf9GSa82asIa7xttFSkh1dRh3L5P57l3y9HWu\nAS5z9y3Wyq2aZjYJmARQUVFBXY7/RGpoaMj5c/Nh7txdgcF067aQurqGxORqTVKzKVf7KFf7NM3V\nacMGRmyzDWvuvJOXe/RITK6CcfdW34DhwOwmj68Armh2zCvAq5m3BsLQzImtfd3KykrPVW1tbc6f\nmw8XXujeq5f7pk2ffT52rtYkNZtytY9ytc/nclVVuR9ySIwon9GR8wUs9DZ6292zGpZZAAwyswFm\n1hU4E5jR7BfEAHfv7+79gXuBi9x9eod/8yRUXZ3G20VKUk1N2Yy7t1nu7r4JmALMBp4H7nH35WY2\n2cwmFzpg0qxZAytXarxdpCQ1jruXwfruWY25u/ssYFaz525s4djzOx4ruWprw/uamrg5RCQHRxwR\npkXW1cEJJ8ROU1C6Q7WdamvDejIHHhg7iYi027bbhvXdG6/SUkzl3k61teGuVK0nI1KiamrCOjPv\nvRc7SUGpotrh1VfDnqkjR8ZOIiI5q6n5/zvbp5jKvR003i6SAocdBjvskPqhGZV7O8yfD7vsAvvv\nHzuJiOSsa9ewVve8ebGTFJTKPUvuodxHjtR+qSIlb+RIeP55eOut2EkKRuWepRUr4M03NSQjkgqN\nL5wlcMmEfFG5Z6nxX3B6MVUkBQ4+GHr1SvXQjMo9S/PmQb9+MHBg7CQi0mGdO4e7VVXu5W3LlvDC\nusbbRVJk5Mgwt/mVV2InKQiVexaeeQbWrYNRo2InEZG8afyBTunVu8o9C43/71XuIimy336w224q\n93I2bx4MHgy77x47iYjkjVm4Yps3L4y9pozKvQ2ffBJWB9VVu0gKjRoF77wDS5fGTpJ3Kvc2PPEE\nfPQRjB4dO4mI5N0xx4T3jzwSN0cBqNzbMHdumDVVVRU7iYjkXd++sO++4Qc9ZVTubZg7N6zvv+OO\nsZOISEGMHh3GXj/+OHaSvFK5t+K992DRIg3JiKTa6NFh7PWJJ2InySuVeyvmzQsLho0ZEzuJiBRM\ndTV06QJz5sROklcq91bMmQM9e4bln0UkpXr0CFvvqdzLg3v4fz1qVPilLiIpNmYMLF4Ma9fGTpI3\nKvcWvPACvPEGHHts7CQiUnCNP+gpunpXubfg4YfDe5W7SBk45BDo0wdmz46dJG9U7i2YPTtMf91r\nr9hJRKTgOnUKQzNz5qRmKQKV+1Z8+GHYoEVX7SJlZOzYMOa+eHHsJHmhct+Kurqwpsy4cbGTiEjR\nNM55fuihuDnyROW+FQ89BNtvryUHRMpKRQVUVqrc08odZs0KG2Fvt13sNCJSVOPHw5NPhtvTS5zK\nvZmXXoJVq+C442InEZGiGz8+vKCagimRKvdmHnwwvFe5i5Shww4LUyJnzoydpMNU7s3MnAn7768p\nkCJlqXPncPX+8MOwaVPsNB2SVbmb2Vgze9HMVprZ5Vv5+P8xsyVmttTMnjCzofmPWnjvvw+PPQb/\n8i+xk4hINMcfD+vWhbH3EtZmuZtZZ+B6YBwwGDjLzAY3O+wVoMrdhwA/BG7Od9BimDMn/LLWkIxI\nGRszJiwo9ac/xU7SIdlcuQ8DVrr7KnffAEwFJjQ9wN2fcPe/Zx4+CfTNb8zimDEjDLcNHx47iYhE\n07NnWAZ4xozYSTrE3L31A8xOBca6+4WZx+cAh7v7lBaO/w6wb+PxzT42CZgEUFFRUTl16tScQjc0\nNNC9e/ecPrclmzcbJ510JEce+R6XX/5CYnLlS1KzKVf7KFf75Jprj2nTGPSb3/DUHXfwUb9+ickF\nUFNTs8jdD23zQHdv9Q04Fbi1yeNzgOtaOLYGeB7Yua2vW1lZ6bmqra3N+XNbMn++O7jfd1/uX6MQ\nufIlqdmUq32Uq31yzvXqq6EQrr46r3kadeR8AQu9jX5196yGZdYATX919c089xlmdiBwKzDB3Uvu\nDoDp08NNS1pPRkTYay846KBQDCUqm3JfAAwyswFm1hU4E/jMYJSZ7QlMA85x95fyH7Ow3OH++8Pr\nKN26xU4jIolw4olhX9W3346dJCdtlru7bwKmALMJQy73uPtyM5tsZpMzh10J7Az81syeMbOFBUtc\nAIsWhY05Tj45dhIRSYyTTw5Xfg88EDtJTrLaQM7dZwGzmj13Y5M/Xwh87gXUUjFtWrh3QfPbReRT\nBxwAAwfCfffBpEmx07Rb2d+h6g733hsWCuvdO3YaEUkMMzjlFJg/P9zUVGLKvtyXLoUVK+C002In\nEZHEOfXUcGdjCQ7NlH25/+//hh22TjwxdhIRSZzKSujfPxRFiSnrcneH//mfMCSz666x04hI4pjB\n6afD3Lklt8Z7WZf74sVhSOb002MnEZHEOuOMMDQzbVrsJO1S1uV+992wzTbhNRMRka06+GAYNCgU\nRgkp23LfvDn8vzr2WNh559hpRCSxzGDiRKirgzWfuzk/scq23Ovrw/+ns8+OnUREEm/ixPAiXQld\nvZdtud95J/TooRuXRCQL++wDw4aF4igRZVnuDQ1hZtNpp8EOO8ROIyIl4bzzYMkSeOaZ2EmyUpbl\nPm0afPABfOUrsZOISMk480zo2hV+//vYSbJSluV+223hxe+jjoqdRERKRu/e4W7Hu+6CTz6JnaZN\nZVfuL70Ejz4KF1wQXgQXEcnahReGdWbuvz92kjaVXbnffHPY+/a882InEZGSM2oUDBgAN90UO0mb\nyqrcP/ooDJedeCJ84Qux04hIyenUKSz/W1cHL+S213KxlFW5T50a/kV10UWxk4hIyfrqV8MLq9df\nHztJq8qm3N3h2mth8GCoro6dRkRK1q67hvVm/vAHWL8+dpoWlU2519eH6anf/KZeSBWRDrrkknDD\nzG23xU7SorIp95//HPr00XIDIpIHlZVQVQW//jVs3Bg7zVaVRbkvXQoPPgjf+AZsv33sNCKSCpde\nCq+/ntj1Zsqi3H/8Y+jeHS6+OHYSEUmNceNg6FC46qqwzGzCpL7cly2De+6BKVO0AbaI5JEZfPe7\n8OKLYSpewqS+3L/73bD643e+EzuJiKTOSSeFq/f//E/YsCF2ms9Idbn/+c8wfXoodm3IISJ516lT\nGPd9+eXE3bWa2nLfsiVMe9x9d/jWt2KnEZHUGj8eRo6E730vUZtop7bcb7kFFi6En/0MunWLnUZE\nUssMrrkm3NB0xRWx03wqleW+ejVcdln4ZTpxYuw0IpJ6Q4aEoYJbbgnrziRA6sp9y5aw9MPGjWEI\nTHejikhRfP/7MHAgnH9+IpYlSF25X3UVzJ0Lv/pVOM8iIkXRrVvYY3XNmnCFuWVL1DipKvf77gtT\nHydOhK99LXYaESk7RxwRXuibNi1Mj4woq3I3s7Fm9qKZrTSzy7fycTOzazMfX2Jmh+Q/ausefDCU\n+vDhcOutGo4RkUi++c2w1duPfhSKPpI2y93MOgPXA+OAwcBZZja42WHjgEGZt0nADXnO2aLGpXwn\nTAivacycqfVjRCQiM7jxxrCh9mWXhXVPPv646DGyuXIfBqx091XuvgGYCkxodswE4A4PngR6mdlu\nec76GZs3h7H16uqw+ua4cTB/vpYYEJEE6NIlbKT97W/Db38b7mK9+2748MOiRcim3PcA3mjyeHXm\nufYekxcPPQTnnjuMnj1hzJiw09VNN8GMGbDjjoX4G0VEctC5c1hrfPbs8OeJE6FXLxg4kL733FPw\nv75Lwf+GJsxsEmHYhoqKCupymA/6yis7sueeX2DYsHUMGbKe4cPfo2vXLdTX5zlsDhoaGnL6byqG\npGZTrvZRrvZJRK6uXeG66+i1eDG9Fy1i27ff5h877FD4XO7e6hswHJjd5PEVwBXNjrkJOKvJ4xeB\n3Vr7upWVlZ6r2tranD+3kJKayz252ZSrfZSrfdKYC1jobfS2u2c1LLMAGGRmA8ysK3AmMKPZMTOA\nczOzZo4A1rv73zr+q0dERHLR5rCMu28ysynAbKAz8Dt3X25mkzMfvxGYBYwHVgIfAl8pXGQREWlL\nVmPu7j6LUOBNn7uxyZ8d0D5HIiIJkao7VEVEJFC5i4ikkMpdRCSFVO4iIimkchcRSSELE10i/MVm\n7wCv5fjpfYB38xgnX5KaC5KbTbnaR7naJ4259nL3Xdo6KFq5d4SZLXT3Q2PnaC6puSC52ZSrfZSr\nfco5l4ZlRERSSOUuIpJCpVruN8cO0IKk5oLkZlOu9lGu9inbXCU55i4iIq0r1St3ERFpRcmVe1ub\ndcdiZq+a2VIze8bMFkbM8TszW2tmy5o819vM5prZisz7nRKS63tmtiZzzp4xs/ERcvUzs1oze87M\nlpvZJZnno56zVnJFPWdmtp2Z/dXMns3k+n7m+djnq6Vc0b/HMjk6m9liM5uZeVzw81VSwzKZzbpf\nAkYTtvJbQNgk5LmowQjlDhzq7lHn1JrZ0UADYU/bAzLP/QxY5+4/yfxC3MndL0tAru8BDe7+82Jm\naZZrN8LGMk+bWQ9gEXAicD4Rz1kruU4n4jkzMwO6uXuDmW0DPA5cApxM3PPVUq6xRP4ey+T7FnAo\nsKO7H1+Mn8lSu3LPZrPusubujwLrmj09Abg98+fbCSVRVC3kis7d/+buT2f+/E/gecL+v1HPWSu5\nospsBtSQebhN5s2Jf75ayhWdmfUFjgNubfJ0wc9XqZV70TbizoEDj5jZosxesUlS0WRnrLeAiphh\nmvlXM1uSGbYp+nBRU2bWHzgYeIoEnbNmuSDyOcsMMTwDrAXmunsizlcLuSD+99g1wKXAlibPFfx8\nlVq5J9kIdz8IGAdcnBmGSJzMxiqJuKIBbgC+CBwE/A34RawgZtYduA/4prv/o+nHYp6zreSKfs7c\nfXPme70vMMzMDmj28Sjnq4VcUc+XmR0PrHX3RS0dU6jzVWrlvgbo1+Rx38xz0bn7msz7tcD9hCGk\npHg7M4bbOJa7NnIeANz97cwP5BbgFiKds8wY7X3AH919Wubp6Odsa7mScs4yWd4Hagnj2tHP19Zy\nJeB8HQVDQgb1AAABGklEQVSckHlNbiow0szuogjnq9TKPZvNuovOzLplXvTCzLoBY4BlrX9WUc0A\nzsv8+TzggYhZPtX4zZ1xEhHOWeaFuNuA5939l00+FPWctZQr9jkzs13MrFfmz9sTJje8QPzztdVc\nsc+Xu1/h7n3dvT+hr+a7+9kU43y5e0m9ETbifgl4Gfi/sfNkMn0ReDbztjxmLuBuwj8/NxJek7gA\n2BmYB6wAHgF6JyTXncBSYEnmm323CLlGEP5JvAR4JvM2PvY5ayVX1HMGHAgszvz9y4ArM8/HPl8t\n5Yr+PdYkYzUws1jnq6SmQoqISHZKbVhGRESyoHIXEUkhlbuISAqp3EVEUkjlLiKSQip3EZEUUrmL\niKSQyl1EJIX+H+DqpabomPJ8AAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "tau_00 = 10\n", "start_of_decrease = 30\n", "x = np.linspace(0,tau_00,100)\n", "y = 0.5*(1-np.cos(x*np.pi/tau_00))\n", "plt.plot(x,y, 'b-')\n", "x2 = np.linspace(start_of_decrease,start_of_decrease+tau_00,100)\n", "y2 = -0.5*(1-np.cos((x2-start_of_decrease)*np.pi/ tau_00)) + 1\n", "plt.plot(x2,y2,'r-')\n", "plt.grid('on')" ] }, { "cell_type": "raw", "metadata": {}, "source": [] }, { "cell_type": "raw", "metadata": { "collapsed": true }, "source": [ "IF( kt <= nn_tau000 ) THEN ! Increase the stress to its nominal value \n", " ! ! during the first nn_tau000 time-steps\n", " zfact = 0.5 * ( 1. - COS( rpi * REAL( kt, wp ) / REAL( nn_tau000, wp ) ) )\n", " zcoef = 1. / ( zrhoa * zcdrag )\n", " ztx = zfact * rn_utau0\n", " zty = zfact * rn_vtau0\n", " DO jj = 1, jpj\n", " DO ji = 1, jpi\n", " vtau (ji,jj) = yslopev * (REAL( jj+njmpp-1 - y0v, wp )) / y0v\n", " utau (ji,jj) = yslopeu * (REAL( jj+njmpp-1 - y0u, wp )) / y0u\n", " vtau(ji,jj) = vtau(ji,jj) * zty\n", " utau(ji,jj) = utau(ji,jj) * ztx\n", " taum(ji,jj) =SQRT( utau(ji,jj) * utau(ji,jj) + vtau(ji,jj) * vtau(ji,jj) )\n", " wndm(ji,jj) =SQRT( taum(ji,jj)*zcoef)\n", " END DO\n", " END DO\n", " ENDIF\n", "\n", "IF((( nn_tau002 + nn_tau001) >= kt) .AND.(kt >= nn_tau002)) THEN ! Decrease the stress to its nominal value \n", " ! ! during the first nn_tau000 time-steps\n", " zfact = -0.5 * ( 1. - COS( rpi * (REAL( (kt-nn_tau002), wp )) / REAL( nn_tau001, wp ) ) ) + 1\n", " zcoef = 1. / ( zrhoa * zcdrag )\n", " ztx = zfact * rn_utau0\n", " zty = zfact * rn_vtau0\n", " DO jj = 1, jpj\n", " DO ji = 1, jpi\n", " vtau (ji,jj) = yslopev * (REAL( jj+njmpp-1 - y0v, wp )) / y0v\n", " utau (ji,jj) = yslopeu * (REAL( jj+njmpp-1 - y0u, wp )) / y0u\n", " vtau(ji,jj) = vtau(ji,jj)\n", " utau(ji,jj) = utau(ji,jj)\n", " taum(ji,jj) =SQRT( utau(ji,jj) * utau(ji,jj) + vtau(ji,jj) * vtau(ji,jj) )\n", " wndm(ji,jj) =SQRT( taum(ji,jj)*zcoef)\n", " END DO\n", " END DO\n", " ENDIF\n", "\n", " IF (kt > (nn_tau002+nn_tau001)) THEN\n", " vtau(:,:) = 0.\n", " utau(:,:) = 0.\n", " ENDIF" ] }, { "cell_type": "code", "execution_count": 65, "metadata": { "collapsed": false }, "outputs": [], "source": [ "jpj = 564\n", "jpi = 284\n", "zcdrag = 0.001\n", "zrhoa = 1\n", "njmpp = 1\n", "nn_tau000=100\n", "nn_tau001=100\n", "nn_tau002=200\n", "rn_utau0 = 0\n", "rn_vtau0=0.169\n", "y0 = -50\n", "y_vtau = 400" ] }, { "cell_type": "code", "execution_count": 66, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "0.169\n", "0.0563333333333\n" ] } ], "source": [ "time = 400\n", "vtau = np.zeros((time,jpi, jpj))\n", "utau = np.zeros((time,jpi, jpj))\n", "taum = np.zeros((time,jpi, jpj))\n", "wndm = np.zeros((time,jpi, jpj))\n", "m = rn_vtau0/(y_vtau - y0)\n", "b = -m*y0\n", "for kt in range(time):\n", " if kt <= nn_tau000:\n", " zfact = 0.5 * ( 1. - np.cos( np.pi * kt / nn_tau000 ) )\n", " zcoef = 1. / ( zrhoa * zcdrag )\n", " #ztx = zfact * rn_utau0\n", " zty = zfact * rn_vtau0\n", " for jj in range(jpj):\n", " for ji in range(jpi):\n", " vtau[kt,ji,jj]=m*(jj+njmpp-1)+b\n", " #utau[kt,ji,jj] = yslopeu * ( jj+njmpp-1 - y0u) / y0u\n", " vtau[kt,...] = vtau[kt,...] * zty\n", " #utau[kt,ji,jj] = utau[kt,ji,jj] * ztx\n", " #taum[kt,...] = np.sqrt( utau[kt,ji,jj] * utau[kt,ji,jj] +\n", " #vtau[kt,...] * vtau[kt,ji,jj] )\n", " wndm[kt,...] = np.sqrt( np.abs(vtau[kt,...])*zcoef)\n", " elif kt>nn_tau000 and kt < nn_tau002:\n", " for jj in range(jpj):\n", " for ji in range(jpi):\n", " vtau[kt,ji,jj]=m*(jj+njmpp-1)+b\n", " #utau[kt,ji,jj] = yslopeu * ( jj+njmpp-1 - y0u) / y0u\n", " #utau[kt,ji,jj] = utau[kt,ji,jj] * ztx\n", " #taum[kt,...] = np.sqrt( utau[kt,ji,jj] * utau[kt,ji,jj] +\n", " #vtau[kt,...] * vtau[kt,ji,jj] )\n", " wndm[kt,...] = np.sqrt( np.abs(vtau[kt,...])*zcoef)\n", " elif((( nn_tau002 + nn_tau001) >= kt) and(kt >= nn_tau002)): \n", " zfact = -0.5 * ( 1. - np.cos( np.pi * (kt-nn_tau002) / nn_tau001 ) ) + 1\n", " zcoef = 1. / ( zrhoa * zcdrag )\n", " ztx = zfact * rn_utau0\n", " zty = zfact * rn_vtau0\n", " for jj in range( jpj):\n", " for ji in range ( jpi):\n", " vtau[kt,ji,jj]=m*(jj+njmpp-1)+b\n", " #utau[kt,ji,jj] = yslopeu * ( jj+njmpp-1 - y0u ) / y0u\n", " vtau[kt,...] = vtau[kt,...] * zty\n", " #utau[kt,ji,jj] = utau[kt,ji,jj] * ztx\n", " #taum[kt,...] = np.sqrt( utau[kt,ji,jj] * utau[kt,ji,jj] +\n", " #vtau[kt,...] * vtau[kt,ji,jj] )\n", " wndm[kt,...] = np.sqrt( np.abs(vtau[kt,...])*zcoef)\n", " elif (kt > (nn_tau002+nn_tau001)) :\n", " vtau[kt,:,:] = 0.\n", " #utau[kt,:,:] = 0. \n", "print(vtau[150,100,400])\n", "print(vtau[150,100,100])" ] }, { "cell_type": "code", "execution_count": 114, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 114, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAW0AAAD8CAYAAAC8TPVwAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAGA9JREFUeJzt3X+QXlWd5/H3x4aAIgoYjdkEl1jbNU5wBNlUZFfG0WHU\nhLEm7v5BhR01ZeGmqAoqW2NNhbF2nNkpqpgt1xqtArO9mJm4/shkFYouzRAxo2tNKZCoMRAC0sMP\nSUyICIgja0LCZ/+4J3ptuvu53Z3m6ZP+vKq6nvucc+5zvnlS9cnN6ftDtomIiDq8qN8FREREdwnt\niIiKJLQjIiqS0I6IqEhCOyKiIgntiIiKJLQjIqZB0gpJ90sakbR+jP4/lrRb0t2Svi3pglH9A5K+\nL+krXeZLaEdETJGkAeAGYCWwFLhC0tJRwx4Cfs/27wB/BQyN6v8wsLfrnAntiIipWw6M2H7Q9hFg\nM7CqPcD2t20/Wd7eASw+3idpMfCHwE1dJzxl2iWfAGef8yIvWjwrSomIWW7P3c8+bvuV0/mMS956\nup964rmu8+0BftlqGrJ9/Gh5EfBoq28f8KYJPu5K4B9a7/8G+FPgzE7FMEtCe9HiU9jy1Wn9HUTE\nHHH+a378yHQ/46knnuucOee/5se/tL1sunNKehtNaF9S3r8LOGT7u5Le2vVzZkVoR0RUaj9wbuv9\n4tL2GyS9gWYJZKXtn5bmNwN/JOky4HTgZZI+Z/s9E02YNe2IiKnbAQxKWiJpHrAaGG4PkPQa4Gbg\nvbZ/eLzd9rW2F9s+r+z3j70CG3KkHRExZbaPSroa2AYMABtt75F0VenfAPw58ArgRkkAR6ez3JLQ\njoiYBttbga2j2ja0tj8AfKDHZ3wT+GaX+bI8EhFRkYR2RERFEtoRERVJaEdEVCShHRFRkYR2RERF\nEtoRERVJaEdEVCShHRFRkYR2RERFEtoRERVJaEdEVCShHRFRkYR2RERFEtoRERWZFffTfvLYS7jl\n6Tf2u4yIqMKP+11AX3U60pb0sKS7Je2StLO0nSPpdkkPlNezW+OvlTQi6X5J75yp4iMi5prJLI+8\nzfaFrcfkrAe22x4Etpf3SFpK87yz84EVNI/YGTiBNUdEzFnTWdNeBWwq25uAd7faN9s+bPshYARY\nPo15IiJmLUkryqrCiKT1Y/T/saTdZbXi25IuaPVtlHRI0j1d5+sa2ga+Lum7ktaWtgW2D5Ttg8CC\nsr0IeLS1777SNvoPslbSTkk7f/Hkka71RkTMGmUV4QZgJbAUuKKsNrQ9BPye7d8B/goYavX9Hc2K\nRGddfxF5ie39kl4F3C7pvnanbUvyZCa2PUQpftH5Z01q34iIWWI5MGL7QQBJm2lWG+49PsD2t1vj\n7wAWt/q+Jem8yUzY6Ujb9v7yegi4pRT6mKSFpdCFwKEyfD9wbmv3xaUtIuJk02lloeVK4B+mM2HP\nI21JZwAvsv3zsv0O4L8Bw8Aa4PryemvZZRj4gqRPAP8KGATumk6REREn0uROM/7x/ONnzRVDZaVg\nUiS9jSa0L5nsvm1dlkcWALdIOj7+C7Zvk7QD2CLpSuAR4HIA23skbaH578FRYJ3tY9MpMiKijx5v\nnTU3WqeVBUlvAG4CVtr+6XSK6RnaZa3mgjHafwpcOs4+1wHXTaewiIgK7AAGJS2hCevVwH9qD5D0\nGuBm4L22fzjdCXMZe0TEFNk+ClwNbAP2AlvKasNVkq4qw/4ceAXNNSu72kstkr4IfAf4LUn7ysrF\nhGbFZewREbWyvRXYOqptQ2v7A8AHxtn3isnOlyPtiIiKJLQjIiqS0I6IqEhCOyKiIgntiIiKJLQj\nIiqS0I6IqEhCOyKiIgntiIiKJLQjIiqS0I6IqEhCOyKiIgntiIiKJLQjIiqS0I6IqEhCOyKiIgnt\niIiKJLQjIiqS0I6IqMiseEbk00dPZ9vB3+53GRFRha/2u4DfIGkF8ElgALjJ9vWj+l8H/C1wEfBR\n2x9v9f0XmudHGrgbeL/tX040X460IyKmSNIAcAOwElgKXCFp6ahhTwAfAj4+at9FpX2Z7dfThP7q\nXnMmtCMipm45MGL7QdtHgM3AqvYA24ds7wCeHWP/U4AXSzoFeAnw414TzorlkYiIF9LklmS/Ol/S\nzlbDkO2hsr0IeLTVtw94U5dPtb1f0seBHwH/D/ia7a/12i+hHRExscdtLzvRHyrpbJqj8iXAU8D/\nkfQe25+baL8sj0RETN1+4NzW+8WlrYs/AB6y/RPbzwI3A/++104J7YiIqdsBDEpaImkezS8Shzvu\n+yPgYkkvkSTgUmBvr52yPBIRMUW2j0q6GthGc/bHRtt7JF1V+jdIejWwE3gZ8Jyka4Cltu+U9CXg\ne8BR4PvA0JgTtXQO7XJqy05gv+13SToH+HvgPOBh4HLbT5ax1wJXAseAD9ne1nWeiIia2N4KbB3V\ntqG1fZBm2WSsfT8GfGwy801meeTD/Oah+3pgu+1BYHt5TzlHcTVwPrACuLEEfkRETFOn0Ja0GPhD\n4KZW8ypgU9neBLy71b7Z9mHbDwEjNOcyRkTENHU90v4b4E+B51ptC2wfKNsHgQVle6zzFheN/kBJ\nayXtlLTz2aeemVzVERFzVM/QlvQu4JDt7443xrZprp3vzPaQ7WW2l5161ksms2tExJzV5ReRbwb+\nSNJlwOnAyyR9DnhM0kLbByQtBA6V8dM5bzEiIibQ80jb9rW2F9s+j+YXjP9o+z005yKuKcPWALeW\n7WFgtaTTJC0BBoG7TnjlERFz0HTO074e2CLpSuAR4HKAco7iFuBemnMP19k+Nu1KIyJicqFt+5vA\nN8v2T2mu4Blr3HXAddOsLSIiRsll7BERFUloR0RUJKEdEVGRhHZEREUS2hERFUloR0RUJKEdEVGR\nhHZEREUS2hERFUloR0RUJKEdETENklZIul/SiKT1Y/S/TtJ3JB2W9JFRfQ9LulvSLkk7u8yXB/tG\nRExReZTiDcDbaR74skPSsO17W8OeAD7Er5/uNdrbbD/edc4caUdETN1yYMT2g7aPAJtpHrn4K7YP\n2d4BPHsiJsyRdkTMOUeOnMKP9s3vOnz+qKWLIdtDZXusxyu+aRKlGPi6pGPA/2x97rhmRWhP8guM\niHghPW572Qx99iW290t6FXC7pPtsf2uiHbI8EhExddN6vKLt/eX1EHALzXLLhBLaERFTtwMYlLRE\n0jyaRzIOd9lR0hmSzjy+DbwDuKfXfrNieSQioka2j0q6GtgGDAAbyyMXryr9GyS9GtgJvAx4TtI1\nwFJgPnCLJGiy+Au2b+s1Z0I7ImIabG8Fto5q29DaPkizbDLa08AFk50vyyMRERVJaEdEVCShHRFR\nkYR2RERFEtoRERVJaEdEVCShHRFRkYR2RERFEtoRERVJaEdEVCShHRFRkZ6hLel0SXdJ+oGkPZL+\nsrSfI+l2SQ+U17Nb+1xbnpd2v6R3zuQfICJiLulypH0Y+H3bFwAXAiskXQysB7bbHgS2l/dIWkpz\ne8LzgRXAjeU5ahERMU09Q9uNfylvTy0/pnkO2qbSvolfP7RyFbDZ9mHbDwEjdLixd0RE9NZpTVvS\ngKRdwCHgdtt3AgtsHyhDDgILyvZYz0xbNMZnrpW0U9LOYz//xZT/ABERc0mn0LZ9zPaFNPeEXS7p\n9aP6TXP03ZntIdvLbC8bOPOMyewaETFnTersEdtPAd+gWat+TNJCgPJ6qAyb1jPTIiJifF3OHnml\npLPK9ouBtwP30TwHbU0Ztga4tWwPA6slnSZpCTAI3HWiC4+ImIu6HGkvBL4haTfNQyxvt/0V4Hrg\n7ZIeAP6gvMf2HmALcC9wG7DO9rGZKD4iot8krSinN49IWj9G/+skfUfSYUkfabWfK+kbku4tp1N/\nuMt8PZ8RaXs38MYx2n8KXDrOPtcB13UpICKiVuV05htoViD2ATskDdu+tzXsCeBD/PoMu+OOAn9i\n+3vlqezflXT7qH2fJ1dERkRM3XJgxPaDto8Am2lOe/4V24ds7wCeHdV+wPb3yvbPgb2McabdaAnt\niIiJzT9+enL5Wdvq63SKcy+SzqNZ0biz19ieyyMREScbHRGn/Whe1+GP2142Y7VILwW+DFxj++le\n43OkHRExddM6xVnSqTSB/XnbN3fZJ6EdETF1O4BBSUskzaO579Jwlx0lCfgMsNf2J7pOmOWRiIgp\nsn1U0tXANmAA2Gh7j6SrSv8GSa8GdgIvA56TdA2wFHgD8F7g7nKbEIA/s711ojlnRWhPcn0pImLW\nKCG7dVTbhtb2QZplk9H+CdBk58vySERERRLaEREVSWhHRFQkoR0RUZGEdkRERRLaEREVSWhHRFQk\noR0RUZGEdkRERRLaEREVSWhHRFQkoR0RUZGEdkRERRLaEREVSWhHRFQkoR0RUZGEdkRERRLaEREV\nSWhHREyDpBWS7pc0Imn9GP2S9KnSv1vSRa2+D0u6R9Ke8uzInhLaERFTJGkAuAFYSfOw3iskLR01\nbCUwWH7WAp8u+74e+M/AcuAC4F2S/k2vORPaERFTtxwYsf2g7SPAZmDVqDGrgM+6cQdwlqSFwG8D\nd9p+xvZR4P8C/7HXhAntiIiJzZe0s/WzttW3CHi09X5faaPDmHuA35X0CkkvAS4Dzu1VzCm9Bkg6\nF/gssAAwMGT7k5LOAf4eOA94GLjc9pNln2uBK4FjwIdsb+s1T0TEC2XgCJz5iLsOf9z2shNdg+29\nkv4a+BrwC2AXTWZOqMuR9lHgT2wvBS4G1pU1m/XAdtuDwPbyntK3GjgfWAHcWNZ9IiJONvv5zaPj\nxaWt0xjbn7H9b22/BXgS+GGvCXuGtu0Dtr9Xtn8O7KU5tF8FbCrDNgHvLturgM22D9t+CBihWfeJ\niDjZ7AAGJS2RNI/mgHV41Jhh4H3lLJKLgZ/ZPgAg6VXl9TU069lf6DVhz+WRNknnAW8E7gQWHJ8Y\nOEizfAJNoN/R2m2sNZ6IiOrZPirpamAbMABstL1H0lWlfwOwlWa9egR4Bnh/6yO+LOkVwLPAOttP\n9Zqzc2hLeinwZeAa209LahduSZ0XiMrnraU5/YVTXn72ZHaNiJg1bG+lCeZ224bWtoF14+z7u5Od\nr9PZI5JOpQnsz9u+uTQ/Vk5bobweKu1d1niwPWR7me1lA2ecMdm6IyLmpJ6hreaQ+jPAXtufaHUN\nA2vK9hrg1lb7akmnSVpCc0L5XSeu5IiIuavL8sibgfcCd0vaVdr+DLge2CLpSuAR4HKAsp6zBbiX\n5syTdbZ7nsYSERG99Qxt2/8EaJzuS8fZ5zrguq5FTPKcyYiIOStXREZEVCShHRFRkYR2RERFEtoR\nERVJaEdEVCShHRFRkYR2RERFEtoRERVJaEdEVCShHRFRkYR2RERFEtoRERVJaEdEVCShHRFRkYR2\nRMQ0SFoh6X5JI5LWj9EvSZ8q/bslXdTqO0vSlyTdJ2mvpH/Xa76EdkTEFEkaAG4AVgJLgSskLR01\nbCXNE7wGaZ6L++lW3yeB22y/DrgA2NtrzoR2RMTULQdGbD9o+wiwGVg1aswq4LNu3AGcJWmhpJcD\nb6F5nCO2j3R5GntCOyJiYvMl7Wz9rG31LQIebb3fV9roMGYJ8BPgbyV9X9JNkno+5bzLMyIjIk4q\nA780L//nw12HP2572QyUcQpwEfBB23dK+iSwHvivE+2UI+2IiKnbD5zber+4tHUZsw/YZ/vO0v4l\nmhCfUEI7ImLqdgCDkpZImgesBoZHjRkG3lfOIrkY+JntA7YPAo9K+q0y7lLg3l4TZnkkImKKbB+V\ndDWwDRgANtreI+mq0r8B2ApcBowAzwDvb33EB4HPl8B/cFTfmBLaERHTYHsrTTC32za0tg2sG2ff\nXcCk1suzPBIRUZGEdkRERRLaEREVSWhHRFQkoR0RUZGEdkRERWbFKX+TvKQ0ImLO6nmkLWmjpEOS\n7mm1nSPpdkkPlNezW33XlvvG3i/pnTNVeETEXNRleeTvgBWj2tYD220PAtvLe8p9ZFcD55d9biz3\nm42IiBOgZ2jb/hbwxKjmVcCmsr0JeHerfbPtw7Yforlsc/kJqjUiYs6b6i8iF9g+ULYPAgvKdpd7\nywIgae3x+9MeefYXUywjImJumfbZI+W6ek9hvyHby2wvm3dqz/t+R0QEUw/txyQtBCivh0p7l3vL\nRkTEFE01tIeBNWV7DXBrq321pNMkLaF5kOVd0ysxIiKO63metqQvAm+leU7aPuBjwPXAFklXAo8A\nlwOU+8huobmR91Fgne1jM1R7RMSc0zO0bV8xTtel44y/DrhuOkVFRMTYchl7RERFEtoRERVJaEdE\nTIOkFeW2HSOS1o/RL0mfKv27JV1U2k+XdJekH0jaI+kvu8yX0I6ImKJym44bgJXAUuCKcjuPtpU0\nZ9INAmuBT5f2w8Dv274AuBBYUZ7WPqGEdkTE1C0HRmw/aPsIsJnmdh5tq4DPunEHcJakheX9v5Qx\np5afnhcqJrQjIqauy607xh0jaUDSLpoLFG+3fWevCWfF/bQjIl5I+uUR5t23r+vw+ZJ2tt4P2R46\nEXWU61gulHQWcIuk19u+Z6J9EtoRERN73Paycfq63Lqj5xjbT0n6Bs0trScM7SyPRERM3Q5gUNIS\nSfNonicwPGrMMPC+chbJxcDPbB+Q9MpyhI2kFwNvB+7rNWGOtCMipsj2UUlXA9uAAWBjuZ3HVaV/\nA7AVuIzm+QLPAO8vuy8ENpUzUF4EbLH9lV5zJrQjIqbB9laaYG63bWhtG1g3xn67gTdOdr4sj0RE\nVCShHRFRkYR2RERFEtoRERVJaEdEVCShHRFRkYR2RERFEtoRERVJaEdEVCShHRFRkVlxGfskb5MY\nETFn5Ug7IqIiCe2IiIoktCMiKpLQjoioSEI7IqIiCe2IiIoktCMiKpLQjoioyIyFtqQVku6XNCJp\n/UzNExHRT72yrjyF/VOlf7eki7ruO5YZCe3ydOEbgJXAUuAKSUtnYq6IiH7pmHUrgcHysxb49CT2\nfZ6ZOtJeDozYftD2EWAzsGqG5oqI6JcuWbcK+KwbdwBnSVrYcd/nmal7jywCHm293we8qT1A0lqa\nf3UADt928MZ7ZqiWF8J84PF+FzFFqb1/aq6/n7X/6+l+wNNHf7LttoM3zu84/HRJO1vvh2wPle2e\nWTfOmEUd932evt0wqvyhhwAk7bS9rF+1TFfN9af2/qm5/pprB7C9ot81TNVMhfZ+4NzW+8WlLSLi\nZNIl68Ybc2qHfZ9npta0dwCDkpZImgesBoZnaK6IiH7pknXDwPvKWSQXAz+zfaDjvs8zI0fato9K\nuhrYBgwAG23vmWCXoQn6alBz/am9f2quv+baT5jxsk7SVaV/A7AVuAwYAZ4B3j/Rvr3mlO0Z+cNE\nRMSJlysiIyIqktCOiKhI30O7tsvdJT0s6W5Ju46fuynpHEm3S3qgvJ7d7zqPk7RR0iFJ97Taxq1X\n0rXl7+J+Se/sT9W/qmWs2v9C0v7y/e+SdFmrbzbVfq6kb0i6V9IeSR8u7bP+u5+g9iq++5Oe7b79\n0Cy+/zPwWmAe8ANgaT9r6lDzw8D8UW3/HVhfttcDf93vOlu1vQW4CLinV700l9L+ADgNWFL+bgZm\nWe1/AXxkjLGzrfaFwEVl+0zgh6XGWf/dT1B7Fd/9yf7T7yPtk+Vy91XAprK9CXh3H2v5Dba/BTwx\nqnm8elcBm20ftv0QzW+7l78ghY5hnNrHM9tqP2D7e2X758BemivgZv13P0Ht45k1tc8F/Q7t8S7v\nnM0MfF3Sd8ul+AAL3Jx3CXAQWNCf0jobr95a/j4+WO6WtrG1vDBra5d0HvBG4E4q++5H1Q6Vffcn\no36Hdo0usX0hzZ251kl6S7vTzf8XqzmPsrZ6ae6Q9lrgQuAA8D/6W87EJL0U+DJwje2n232z/bsf\no/aqvvuTVb9Du7rL3W3vL6+HgFto/hv4WLlrF+X1UP8q7GS8emf934ftx2wfs/0c8L/49X/DZ13t\nkk6lCb3P2765NFfx3Y9Ve03f/cms36Fd1eXuks6QdObxbeAdwD00Na8pw9YAt/anws7Gq3cYWC3p\nNElLaO7/e1cf6hvX8cAr/gPN9w+zrHZJAj4D7LX9iVbXrP/ux6u9lu/+pNfv34TSXN75Q5rfOH+0\n3/X0qPW1NL8l/wGw53i9wCuA7cADwNeBc/pda6vmL9L8V/ZZmrXGKyeqF/ho+bu4H1g5C2v/38Dd\nwG6asFg4S2u/hGbpYzewq/xcVsN3P0HtVXz3J/tPLmOPiKhIv5dHIiJiEhLaEREVSWhHRFQkoR0R\nUZGEdkRERRLaEREVSWhHRFTk/wOR+Gx/W+FPXQAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.contourf(vtau[150,:,:].T\n", " )\n", "plt.colorbar()" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "0.4" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "Q0 = 10**3\n", "C = 10**(-3)\n", "Qa = 1\n", "\n", "windspeed = 20\n", "windspeed**2 * Qa * C / Q0 * Q0 " ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.6.0" } }, "nbformat": 4, "nbformat_minor": 2 }