{"id":402,"date":"2022-06-26T08:59:43","date_gmt":"2022-06-25T23:59:43","guid":{"rendered":"https:\/\/www2.jpgu.org\/a\/?p=402"},"modified":"2022-07-05T12:52:25","modified_gmt":"2022-07-05T03:52:25","slug":"matlab-octave_flux","status":"publish","type":"post","link":"https:\/\/www2.jpgu.org\/a\/2022\/06\/26\/matlab-octave_flux\/","title":{"rendered":"MATLAB\/Octave\u3067\u6d77\u9762\u30d5\u30e9\u30c3\u30af\u30b9\u8a08\u7b97"},"content":{"rendered":"<p><\/p>\n<p style=\"text-align: right\">2022\u5e746\u6708\u3000\u5ddd\u5408\u7fa9\u7f8e<\/p>\n<p><span style=\"font-size: 14pt\">MATLAB\u307e\u305f\u306fGNU Octave\u3067\u6d77\u9762\u30d5\u30e9\u30c3\u30af\u30b9\u3092\u8a08\u7b97\u3057\u307e\u3059\u3002COARE3\u306e\u516c\u958b\u30d7\u30ed\u30b0\u30e9\u30e0\u3092\u4f7f\u7528\u3057\u307e\u3059\u3002GNU Octave\u3067\u3082\u52d5\u304d\u307e\u3059\u3002<\/span><\/p>\n<p>\u8ad6\u6587\u7b49\u3092\u57f7\u7b46\u3059\u308b\u306e\u306b\u4f7f\u7528\u3057\u305f\u5834\u5408\u306f\u53c2\u8003\u6587\u732e\u3084\u30d7\u30ed\u30b0\u30e9\u30e0\u306e\u30c0\u30a6\u30f3\u30ed\u30fc\u30c9\u5143\u3092\u660e\u8a18\u3057\u8b1d\u610f\u3092\u8868\u3057\u307e\u3057\u3087\u3046\u3002\u53c2\u8003\u6587\u732e\u3084\u30d0\u30fc\u30b8\u30e7\u30f3\u9593\u306e\u9055\u3044\u306f\u3054\u81ea\u5206\u3067\u8abf\u3079\u3066\u304f\u3060\u3055\u3044\u3002<\/p>\n<p>&nbsp;<\/p>\n<p><span style=\"font-size: 14pt\">\u4ee5\u4e0b\u306eFTP\u30b5\u30a4\u30c8\u304b\u3089\u3001cor3_0 \u307e\u305f\u306fcor3_5 \u307e\u305f\u306fcor3_6\u3092\u30c7\u30a3\u30ec\u30af\u30c8\u30ea\u3054\u3068\u30c0\u30a6\u30f3\u30ed\u30fc\u30c9\u3057\u307e\u3059\u3002<\/span><\/p>\n<p><span style=\"font-size: 14pt\">ftp:\/\/ftp1.esrl.noaa.gov\/BLO\/Air-Sea\/bulkalg\/<\/span><\/p>\n<p>&nbsp;<\/p>\n<p>\uff11\uff09coare3.0\u306e\u5834\u5408<\/p>\n<p>\u5165\u529b\u30c7\u30fc\u30bf\u306e\u884c\u5217\uff08\u4f8b\uff1adat\uff09\u3092\u4f5c\u6210\u3057\u3001\u95a2\u6570cor3_0ah.m\u306b\u6e21\u3057\u3001\u51fa\u529b\u30c7\u30fc\u30bf\u306e\u884c\u5217\uff08\u4f8b\uff1aflx\uff09\u3092\u5f97\u307e\u3059\u3002<\/p>\n<p>&gt;&gt; flx = cor3_0ah(dat);<\/p>\n<p>\u30d5\u30a1\u30a4\u30ebcor3_0ah.m\u306e\u4e2d\u3067\u3001\u98a8\u901f\u8a08\u3001\u6c17\u6e29\u8a08\u3001\u6e7f\u5ea6\u8a08\u306e\u6d77\u9762\u4e0b\u304b\u3089\u306e\u9ad8\u5ea6\uff08zu, zt, zq\uff09\u3001\u53ca\u3073\u30d0\u30eb\u30af\u6c34\u6e29\u306e\u6df1\u5ea6\uff08ts_depth\uff09\u3092\u6307\u5b9a\u3059\u308b\u5fc5\u8981\u304c\u3042\u308a\u307e\u3059\uff0884\uff0d87\u884c\u76ee\uff09\u3002<\/p>\n<p>Skin layer\u306e\u5f71\u97ff\u3092\u8003\u616e\u3057\u305f\u3044\u5834\u5408\u306f jcool = 1\uff08\u30c7\u30d5\u30a9\u30eb\u30c8\uff09, \u8003\u616e\u3057\u305f\u304f\u306a\u3044\u5834\u5408\u306f0\u306b\u3057\u3066\u304f\u3060\u3055\u3044\u3002<\/p>\n<p>Warm layer\u306e\u5f71\u97ff\u3092\u8003\u616e\u3057\u305f\u3044\u5834\u5408\u306f jwarm = 1\uff08\u30c7\u30d5\u30a9\u30eb\u30c8\uff09, \u8003\u616e\u3057\u305f\u304f\u306a\u3044\u5834\u5408\u306f0\u306b\u3057\u3066\u304f\u3060\u3055\u3044\u3002<\/p>\n<p>\u5165\u529b\u30c7\u30fc\u30bf\u306e\u884c\u5217\u306f\u4ee5\u4e0b\u306e\u901a\u308a\u3067\u3059\u3002\u5168\u3066\u306e\u5217\u306f\u540c\u3058\u9577\u3055\u306e\u30d9\u30af\u30c8\u30eb\u3067\u3059\u3002<\/p>\n<p>1\u5217\u76ee\/ \u65e5\u6642\uff1aYYYYMMDDhhmmss \u306e\u5f62\u3067\u4e0e\u3048\u308b\u3002\u4f8b\uff092022\u5e746\u670819\u65e512\u664230\u520600\u79d2\uff08\u56fd\u969b\u6a19\u6e96\u6642\uff09\u306e\u5834\u5408\u306f\u300120220619123000 \u3068\u3059\u308b\uff08\u6574\u6570\uff09\u3002<\/p>\n<p>2\u5217\u76ee\/ \u98a8\u901f (m\/s)<\/p>\n<p>3\u5217\u76ee\/ \u30b5\u30d6\u30b9\u30ad\u30f3\u6c34\u6e29\u3000\u203b\u4f7f\u7528\u3057\u306a\u3044\u306e\u3067\u30010\u306a\u3069\u9069\u5f53\u306a\u30c0\u30df\u30fc\u30c7\u30fc\u30bf\u3067\u53ef<\/p>\n<p>4\u5217\u76ee\/ \u6c17\u6e29 (\u00b0C)<\/p>\n<p>5\u5217\u76ee\/ \u6bd4\u6e7f (g\/kg)\u3000\u203b\u5358\u4f4d\u306b\u6ce8\u610f<\/p>\n<p>6\u5217\u76ee\/ \u4e0b\u5411\u304d\u77ed\u6ce2\u653e\u5c04 (W\/m2)<\/p>\n<p>7\u5217\u76ee\/ \u4e0b\u5411\u304d\u9577\u6ce2\u653e\u5c04 (W\/m2)<\/p>\n<p>8\u5217\u76ee\/ \u964d\u6c34\u91cf (mm\/hr)<\/p>\n<p>9\u5217\u76ee\/ \u7def\u5ea6 (\u00b0N)<\/p>\n<p>10\u5217\u76ee\/ \u7d4c\u5ea6 (\u00b0E)<\/p>\n<p>11\u5217\u76ee\/ \u30d0\u30eb\u30af\u6c34\u6e29 (\u00b0C)<\/p>\n<p>\u51fa\u529b\u30c7\u30fc\u30bf\u306e3, 4, 6\u5217\u76ee\u304c\u305d\u308c\u305e\u308chsb (\u9855\u71b1\u30d5\u30e9\u30c3\u30af\u30b9), hlb (\u6f5c\u71b1\u30d5\u30e9\u30c3\u30af\u30b9), taub (\u98a8\u5fdc\u529b)\u3067\u3059\u3002\u305d\u308c\u4ee5\u5916\u306e\u8aac\u660e\u306f\u5272\u611b\u3057\u307e\u3059\uff08\u30d7\u30ed\u30b0\u30e9\u30e0\u30d5\u30a1\u30a4\u30eb\u3092\u8aad\u3093\u3067\u304f\u3060\u3055\u3044\uff09\u3002<\/p>\n<p>&nbsp;<\/p>\n<p>2) coare3.5\u306e\u5834\u5408<\/p>\n<p>Warm layer\u3092\u8003\u616e\u3057\u306a\u3044\u5834\u5408\u306b\u306f\u3001<\/p>\n<p>&gt;&gt; flx = coare35vn(u, zu, t, zt, rh, zq, P, ts, Rs, Rl, lat, zi, rain, cp, sigH);<\/p>\n<p>\u3053\u3053\u3067\u3001u\uff08\u98a8\u901f\uff1am\/s\uff09, t\uff08\u6c17\u6e29\uff1a\u00b0C\uff09, rh\uff08\u76f8\u5bfe\u6e7f\u5ea6\uff1a%\uff09, ts\uff08\u6c34\u6e29\uff1a\u00b0C\uff09, P\uff08\u6c17\u5727\uff1ahPa\uff09, ts\uff08\u6c34\u6e29\uff1a\u00b0C\uff09, Rs\uff08\u4e0b\u5411\u304d\u77ed\u6ce2\u653e\u5c04\uff1aW\/m2\uff09, Rl\uff08\u4e0b\u5411\u304d\u9577\u6ce2\u653e\u5c04\uff1aW\/m2\uff09, rain\uff08\u964d\u6c34\u91cf\uff1amm\/hr\uff09, cp\uff08\u5353\u8d8a\u6ce2\u9ad8\u306e\u4f4d\u76f8\u901f\u5ea6\uff1am\/s\uff09, sigH\uff08\u6709\u7fa9\u6ce2\u9ad8\uff1am\uff09\u306f\u30d9\u30af\u30c8\u30eb\u3067\u4e0e\u3048\u307e\u3059\u3002Warm layer\u3092\u8003\u616e\u3057\u306a\u3044\u306e\u3067\u3001ts\u306b\u30d0\u30eb\u30af\u6c34\u6e29\u3092\u4e0e\u3048\u305f\u5834\u5408\u3067\u3082\u30b5\u30d6\u30b9\u30ad\u30f3\u6c34\u6e29\u3068\u898b\u306a\u3055\u308c\u307e\u3059\u3002<\/p>\n<p>zu\uff08\u98a8\u901f\u8a08\u306e\u9ad8\u5ea6\uff09, zt\uff08\u6c17\u6e29\u8a08\u306e\u9ad8\u5ea6\uff09, zq\uff08\u6e7f\u5ea6\u8a08\u306e\u9ad8\u5ea6\uff09, lat\uff08\u7def\u5ea6\uff09, zi\uff08\u5883\u754c\u5c64\u9ad8\u5ea6\uff09\u306f\u30b9\u30ab\u30e9\u30fc\u3067\u4e0e\u3048\u307e\u3059\u3002<\/p>\n<p>P, Rs, Rl, lat, zi\u306e\u30c7\u30fc\u30bf\u304c\u306a\u3044\u5834\u5408\u306b\u306f\u3001NaN\u3092\u5165\u529b\u3057\u307e\u3059\u3002\u305d\u306e\u5834\u5408\u3001\u30c7\u30d5\u30a9\u30eb\u30c8\u5024\u304c\u4f7f\u308f\u308c\u307e\u3059\uff08P=1015 hPa, Rs=150 W\/m2, Rl=370 W\/m2, lat=45\u00b0N, zi=600 m\uff09\u3002<\/p>\n<p>cp, sigH\u306e\u30c7\u30fc\u30bf\u304c\u306a\u3044\u5834\u5408\u306b\u306f\u3001NaN\u3092\u5165\u529b\u3057\u307e\u3059\u3002\u3053\u306e\u5834\u5408\u3001\u30d9\u30af\u30c8\u30eb\u3067\u3082\u30b9\u30ab\u30e9\u30fc\u3067\u3082\u304b\u307e\u3044\u307e\u305b\u3093\u3002<\/p>\n<p>Skin layer\u306e\u5f71\u97ff\u3092\u8003\u616e\u3057\u305f\u3044\u5834\u5408\u306f jcool = 1\uff08\u30c7\u30d5\u30a9\u30eb\u30c8\uff09, \u8003\u616e\u3057\u306a\u3044\u5834\u5408\u306f0\u306b\u3057\u3066\u304f\u3060\u3055\u3044\u3002<\/p>\n<p>\u51fa\u529b\u5024\u306e\u884c\u5217flx\u306e\u4e2d\u8eab\u306f\u3001\u524d\u304b\u3089usr (\u6469\u64e6\u901f\u5ea6), tau (\u98a8\u5fdc\u529b), hsb (\u9855\u71b1\u30d5\u30e9\u30c3\u30af\u30b9), hlb (\u6f5c\u71b1\u30d5\u30e9\u30c3\u30af\u30b9), hbb (\u6d6e\u529b\u30d5\u30e9\u30c3\u30af\u30b9), &#8230;\u3000\u3067\u3059\u3002coare3.0\u3068\u5c11\u3057\u7570\u306a\u308b\u306e\u3067\u6ce8\u610f\u3057\u307e\u3057\u3087\u3046\u3002<\/p>\n<p>&nbsp;<\/p>\n<p>Warm layer\u3092\u8003\u616e\u3057\u305f\u3044\u5834\u5408\u306b\u306f\u3001<\/p>\n<p>&gt;&gt; flx = coare35vnWarm(yday, ur, zu, t, zt, rh, zq, P, ts, Rs, Rl, lat, lon, zi, rain , ts_depth);<\/p>\n<p>yday: 1\u67081\u65e5\u304b\u3089\u306e\u65e5\u6570\uff08\u56fd\u969b\u6a19\u6e96\u6642\u30671\u67081\u65e50\u6642\u3092\u57fa\u6e96\u3068\u3059\u308b\uff09<\/p>\n<p>ur: \u76f8\u5bfe\u98a8\u901f\uff08m\/s\uff09 \u203b\u98a8\u901f\uff0d\u6d77\u9762\u6d41\u901f\u3001\u6d41\u901f\u304c\u308f\u304b\u3089\u306a\u3044\u5834\u5408\u306f\u98a8\u901f\u3067\u4ee3\u7528<\/p>\n<p>ts: \u6df1\u5ea6ts_depth\u306b\u304a\u3051\u308b\u6c34\u6e29\uff08\u00b0C\uff09<\/p>\n<p>lon: \u7d4c\u5ea6\uff08\u00b0E\uff09<\/p>\n<p>ts_depth: \u6c34\u6e29\u3092\u6e2c\u5b9a\u3057\u305f\u6df1\u5ea6\uff08m\uff09<\/p>\n<p>\u3053\u3053\u3067\u306flat\u3068lon\u306f\u30d9\u30af\u30c8\u30eb\u3067\u4e0e\u3048\u307e\u3059\u3002ts_depth\u306f\u30b9\u30ab\u30e9\u30fc\u3067\u3059\u3002<\/p>\n<p>!!! \u6ce8\u610f !!!\u3000\u305d\u306e\u307e\u307e\u3060\u3068\u3001grv\u3068\u3044\u3046\u5909\u6570\u304c\u306a\u3044\u3068\u3044\u3046\u30a8\u30e9\u30fc\u304c\u51fa\u307e\u3059\u3002\u30c7\u30a3\u30ec\u30af\u30c8\u30eacor3_6\u306b\u3042\u308bgrv.m\u3092cor3_5\u306b\u30b3\u30d4\u30fc\u3057\u3066\u7f6e\u304f\u3068\u89e3\u6c7a\u3057\u307e\u3059\u3002<\/p>\n<p>&nbsp;<\/p>\n<p>3) coare3.6\u306e\u5834\u5408<\/p>\n<p>Warm layer\u3092\u8003\u616e\u3057\u306a\u3044\u5834\u5408\u306b\u306f\u3001<\/p>\n<p>&gt;&gt; flx = coare36vn_zrf(ur, zu, t, zt, rh, zq, P, ts, Rs, Rl, lat, zi, rain, Ss, cp, sigH, zrf_u, zrf_t, zrf_q);<\/p>\n<p>\u3053\u3053\u3067\u3001ur\u306f\u76f8\u5bfe\u98a8\u901f\uff08\u98a8\u901f\uff0d\u6d77\u9762\u6d41\u901f\u3001m\/s\uff09\u3001Ss\u306f\u6d77\u9762\u306e\u5b9f\u7528\u5869\u5206\uff08psu\uff09\u3067\u3059\u3002<\/p>\n<p>zrf_u, zrf_t, zrf_q\u306f\u3001\u305d\u308c\u305e\u308c\u98a8\u901f\u3001\u6c17\u6e29\u3001\u6e7f\u5ea6\u306e\u57fa\u6e96\u9ad8\u5ea6\uff08m\uff09\u3067\u3059\u3002\u3053\u306e\u9ad8\u5ea6\u306b\u304a\u3051\u308b\u98a8\u901f\uff08Urf\uff09\u3001\u6c17\u6e29\uff08Trf\uff09\u3001\u6e7f\u5ea6\uff08Qrf\uff09\u306e\u63a8\u5b9a\u5024\u304c\u51fa\u529b\u3055\u308c\u307e\u3059\u3002\u5fc5\u8981\u306a\u3044\u5834\u5408\u306f\u9069\u5f53\u306a\u5024\uff0810\u3068\u304b\uff09\u3092\u5165\u308c\u3066\u304a\u3051\u3070\u3044\u3044\u3067\u3057\u3087\u3046\u3002<\/p>\n<p>zu, zt, zq, lat, zi, zrf_u, zrf_t, zrf_q\u306f\u30b9\u30ab\u30e9\u30fc\u3001\u305d\u308c\u4ee5\u5916\u306f\u30d9\u30af\u30c8\u30eb\u3067\u4e0e\u3048\u307e\u3059\u3002cp, sigH\u306fNaN\u3067\u3082\u30d9\u30af\u30c8\u30eb\u3067\u4e0e\u3048\u308b\u5fc5\u8981\u304c\u3042\u308a\u307e\u3059\u3002Warm layer\u3092\u8003\u616e\u3057\u306a\u3044\u306e\u3067\u3001ts\u306b\u30d0\u30eb\u30af\u6c34\u6e29\u3092\u4e0e\u3048\u305f\u5834\u5408\u3067\u3082\u30b5\u30d6\u30b9\u30ad\u30f3\u6c34\u6e29\u3068\u898b\u306a\u3055\u308c\u307e\u3059\u3002<\/p>\n<p>\u51fa\u529b\u5024flx\u306e\u524d\u304b\u30895\u5217\u306f\u3001coare3.5\u3068\u540c\u3058\u3067\u3059\u3002<\/p>\n<p>&nbsp;<\/p>\n<p>Warm layer\u3092\u8003\u616e\u3057\u305f\u3044\u5834\u5408\u306b\u306f\u3001<\/p>\n<p>&gt;&gt; flx = coare36vnWarm(yday, ur, zu, t, zt, rh, zq, P, ts, Rs, Rl, lat, lon, zi, rain , ts_depth, Ss, cp, sigH, zrf_u, zrf_t, zrf_q);<\/p>\n<p>\u3053\u3053\u3067\u306f\u3001zu, zt, zq, lat, lon, cp, sigH, zrf_u, zrf_t, zrf_q\u3082\u30d9\u30af\u30c8\u30eb\u3067\u4e0e\u3048\u308b\u5fc5\u8981\u304c\u3042\u308b\u3088\u3046\u306a\u306e\u3067\u3001\u6ce8\u610f\u3057\u3066\u304f\u3060\u3055\u3044\u3002zi\u3068ts_depth\u306f\u30b9\u30ab\u30e9\u30fc\u3067\u3088\u3044\u3088\u3046\u3067\u3059\u3002<\/p>\n<p>&nbsp;<\/p>","protected":false},"excerpt":{"rendered":"<p>2022\u5e746\u6708\u3000\u5ddd\u5408\u7fa9\u7f8e MATLAB\u307e\u305f\u306fGNU Octave\u3067\u6d77\u9762\u30d5\u30e9\u30c3\u30af\u30b9\u3092\u8a08\u7b97\u3057\u307e\u3059\u3002COARE3\u306e\u516c\u958b\u30d7\u30ed\u30b0\u30e9\u30e0\u3092\u4f7f\u7528\u3057\u307e\u3059\u3002GNU Octave\u3067\u3082\u52d5\u304d\u307e\u3059\u3002 \u8ad6\u6587\u7b49\u3092\u57f7\u7b46\u3059\u308b\u306e\u306b\u4f7f\u7528\u3057\u305f\u5834\u5408\u306f\u53c2\u8003\u6587\u732e\u3084\u30d7\u30ed [&hellip;]<\/p>\n","protected":false},"author":19,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":[],"categories":[5],"tags":[],"_links":{"self":[{"href":"https:\/\/www2.jpgu.org\/a\/wp-json\/wp\/v2\/posts\/402"}],"collection":[{"href":"https:\/\/www2.jpgu.org\/a\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www2.jpgu.org\/a\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www2.jpgu.org\/a\/wp-json\/wp\/v2\/users\/19"}],"replies":[{"embeddable":true,"href":"https:\/\/www2.jpgu.org\/a\/wp-json\/wp\/v2\/comments?post=402"}],"version-history":[{"count":13,"href":"https:\/\/www2.jpgu.org\/a\/wp-json\/wp\/v2\/posts\/402\/revisions"}],"predecessor-version":[{"id":455,"href":"https:\/\/www2.jpgu.org\/a\/wp-json\/wp\/v2\/posts\/402\/revisions\/455"}],"wp:attachment":[{"href":"https:\/\/www2.jpgu.org\/a\/wp-json\/wp\/v2\/media?parent=402"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www2.jpgu.org\/a\/wp-json\/wp\/v2\/categories?post=402"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www2.jpgu.org\/a\/wp-json\/wp\/v2\/tags?post=402"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}