Some tests on various versions of the WTG CRM. Test001: RCE calculation using dyn3d010, wind = 0, sst = 300 K. Test002: WTG calculation using reference profile from test001 with wind 5 m/s sst 300 K and no modifications to profiles. Set trelax = drelax = mrelax = 0.01 s^-1 and madvect = 0. Note that this is not like the Sobel model, as vertical advection is ignored. Test003: WTG calculation using mrelax = 0 and madvect = 1. Otherwise like test002. Test004: WTG calculation like test003 except trelax = 0.001 s^-1. Limit to 1000000 s. Go to dyn3d011 with smoothing in sthwtg. Test005: WTG calculation like test003 except dyn3d011. Go to dyn4d011x with smoothing removed from vzwtg. Test006: WTG calculation like test005 except dyn3d011x. Go to dyn4d012, which is like 11x except code simplification implemented; should be no changes in physics. (Section on forced heating removed.) Test007: WTG calculation like test006 except dyn3d012. ------------------------------------------------------------------ Make some tests on spec* and entropy* Test008: RCE calculation like test001 except spec008. Test009: RCE calculation like test008 except entropy004. ***Test 001 and test 9 show very similar RCE profiles of th, the, thes. Test010: Like test003 except entropy004, spectral = 0, wtg = 1, sobel = 0, hscale = 4775 m. This comes from h = ztrop*brunt/(PI*trelax) with ztrop = 15000 m, brunt = 0.01 s^-1, and trelax = 0.01 s^-1. Test011: Like test010 except spectral = 15. Test012: Like test010 except spec009. (Madvect = 1, Mrelax = 0) Test013: Like test012 with spec009 except spectral = 15. Compare with test011. The statistics on WTG vertical velocity for such weak winds is poor, so we cannot conclude much on the comparison of the 5 m/s WTG results between dyn3d010, spec009, and entropy004. Let's compare spec009 and entropy005 for higher wind speeds. Test014: WTG calc like test012 except spec009 and windspeed = 20 m/s. Reference profile from test008, averaged over 3000000 5000000 s. Test015: WTG calc like test010 (spectral = 0) except entropy005 and windspeed = 20 m/s. Reference profile from test009, averaged over 3000000-5000000 s. The results of these two tests are very similar. Test016: WTG calc like test015 except spectral = 15. (Entropy005, windspeed = 20 m/s. Same reference profile as test015.) Test017: RCE test of entropy006x to sort out checkpoint problem. Like test009. Test018: RCE test with extropy008x to test redo of radiation. Like test017 otherwise. Fixed radiation. Test019: RCE like test018 except real radiation. Test020: WTG calc like test016 except profile file is from test019. ------------------------------------------------------- Tests with entropy017 and entropy017x -- the latter mirrors vertical fluxes so that the flux is zero at the actual upper and lower boundaries. Entropy017y eliminates the term in the continuity equation that keeps the mean surface pressure at a specified value. Test021: RCE with entropy017. Note that drelax = 0. Test022: RCE with entropy017x. Note that drelax = 0. Test023: RCE with entropy017y. Note that drelax = 0. NOTE: EVERYTHING UP TO THIS POINT DELEATED. ------------------------------------------------------- Entropy017y has become the committed version entropy018. (Change to entropy019, as entropy018 is flawed.) Let's do some RCE tests. Test024: RCE with entropy017, drelax = 0.01 s^-1, vx = 0. Run for 10000000 s. Ran OK. Mean calculated for 8000000-10000000 s. The stratosphere develops very bad drift with this version, producing an adiabatic layer at the very top and increased theta values there. Test025: RCE with entropy019, drelax = 0.01 s^-1, vx = 0. Run for 10000000 s. Ran OK. Mean calculated for 8000000-10000000 s. The stratosphere problems of test024 are gone. Test026: RCE with entropy019, drelax = 0.01 s^-1, vx = 5. Run for 10000000 s. Ran OK. Mean calculated for 8000000-10000000 s. Test027: RCE with entropy019, drelax = 0.01 s^-1, vx = 10. Run for 10000000 s. Ran OK. Mean calculated for 8000000-10000000 s. Test028: RCE with entropy019, like test026 except mu = 0.01 rather than 0.1. Ran OK. However, there is still a lot of drift in the lower stratosphere for this case. Test029: RCE with entropy019, like test026 except mu = 0.001 rather than 0.1. Died after about 6000000 s, probably due to excessive stratospheric drift. OK, let's do some WTG comparisons between entropy017 and entropy019 starting with reference profiles from test024 and test025. Test030: WTG test with spec = 60 and hscale = 150 km. Entropy017, starting with test024. Wind = 10 m/s, entrain on, sobel off. Test031: WTG test with spec = 60 and hscale = 150 km. Entropy019, starting with test025. Wind = 10 m/s, entrain on, sobel off. Test032: RCE test like test026 except vy = 5 m/s. Initial file from test026. Entropy019. Finished. Let's do a standard array of WTG tests with RCE test026 as a reference profile and using entropy019. Note that an error was made in setting up test026 -- the 5 m/s wind was in the x rather than the z direction. Will compare later with a proper reference profile (test032). Test033: WTG with vy = 3 m/s. Finished. Test034: WTG with vy = 5 m/s. Finished. Test035: WTG with vy = 7 m/s. Finished. Test036: WTG with vy = 10 m/s. Finished. Similar WTG tests with RCE test032, where the error noted above has been corrected. Test037: WTG with vy = 0 m/s. Redone with entropy022. Test038: WTG with vy = 3 m/s. Redone with entropy022. Test039: WTG with vy = 5 m/s. Redone with entropy022. Test040: WTG with vy = 7 m/s. Redone with entropy022. Test041: WTG with vy = 10 m/s. Redone with entropy022. Test042: WTG with vy = 15 m/s. Redone with entropy022. Test043: WTG with vy = 20 m/s. Redone with entropy022. Entropy019 and earlier has a bug in the interactive radiation code -- the radiation is half a grid cell off in its proper point of application. Entropy020 tentatively fixes this. Let's do some comparison tests. Test044: RCE with interactive radiation and entropy019 using test032 output as a starting profile. Interrupted. Test045: RCE like test044 except entropy020. Interrupted. The differences between test044 and test045 are small. Meanwhile, entropy021 introduces a new diagostic variable, srtrain, which is the net source of rt due to rain processes. In principle the vertically integrated srtrain*rho should equal the rainflux. Test046: RCE test like test032, except entropy021. Test047: RCE test like test046, with entropy021, but with lambda = 0.01. rtcum and raincum balance quite well (though not perfectly) in test044 and test045. Rainfall rate balances well with the new variable srtrain in both test046 and test047. This justifies shutting down these four simulations prematurely. The RCE run test032 and the WTG runs test037-test043 are valid since they didn't use interactive radiation -- so, let's proceed with this sequence, with stablized reference profiles. Test048: WTG with wind = 0, thpert = 1, thmode = 2, entropy021. Redone with entropy022. Test049: WTG with wind = 3, thpert = 1, thmode = 2, entropy021. Redone with entropy022. Test050: WTG with wind = 5, thpert = 1, thmode = 2, entropy021. Redone with entropy022. Test051: WTG with wind = 7, thpert = 1, thmode = 2, entropy021. Redone with entropy022. Test052: WTG with wind = 10, thpert = 1, thmode = 2, entropy021. Redone with entropy022. Test053: WTG with wind = 15, thpert = 1, thmode = 2, entropy021. Redone with entropy022. Test054: WTG with wind = 20, thpert = 1, thmode = 2, entropy021. Redone with entropy022. WTG runs with destabilized reference profiles. Test055: WTG with wind = 0, thpert = -1, thmode = 2, entropy021. Not needed. Test056: WTG with wind = 3, thpert = -1, thmode = 2, entropy021. Not needed. Test057: WTG with wind = 5, thpert = -1, thmode = 2, entropy021. Not needed. Test058: WTG with wind = 7, thpert = -1, thmode = 2, entropy021. Redone with entropy022. Test059: WTG with wind = 10, thpert = -1, thmode = 2, entropy021. Redone with entropy022. Test060: WTG with wind = 15, thpert = -1, thmode = 2, entropy021. Redone with entropy022. Test061: WTG with wind = 20, thpert = -1, thmode = 2, entropy021. Redone with entropy022. WTG runs with no moisture or temperature perturbations, but SST perturbations. Test062: WTG with wind = 7 m/s, SST = 301 K. Finished. Test063: WTG with wind = 7 m/s, SST = 299 K. Finished. We still have some issues with non-conservation. Do some tests. Test064: RCE like test032 except entropy022x, which uses surface fluxes as lower BCs to vertical fluxes. CRASHED. Test065: RCE like test032 except entropy022y, which in addition imposes zero vertical flux at the model top. CRASHED. Both of these versions have severe problems with instability in the vertical thermodynamic structure. Let's take a different approach -- zero vertical fluxes at the lower boundary, approximate zero vertical flux at the upper gridpoint, put the lower boundary at z = 0, so that the grid points are at z = (iz + 0.5)*dz. Fix also the extrapolation of pi and th to the surface for the surface flux calculations. The fluxes are dumped at the lowest grid point. This makes the interactive radiation model's grid consistent with the convective model grid. Test066: RCE like test032 except starting from scratch, run for 15000000 s. STOPPED. Didn't work. Go back to entropy021 and add more cumulative versions of source terms to make diagnostics more accurate. Redo test066 with this new entropy022. Test067: WTG like test059 with entropy022. STOPPED. Both test066 and test067 were stopped, as they demonstrated success. RCE test066 demonstrated perfect budgets for entropy and water vapor and WTG test067 produced much more consistent measures of NGMS. They were not perfectly consistent, as a steady state was not reached, but for such a touchy case, they were much closer than previously. Redo a bunch of tests with entropy022: 037-043, 048-054, 058-061. Run some tests with the reference humidity profile modified. Use entropy023, as this has the same physics as entropy022 (aside from a trivial error in radiation calculation at initial time in 022). Test068: WTG run test032 as initial reference profile, sst = 300 K, wind = 0, no theta perturb, rhpert = +0.1, rhmode = 0. FINISHED. Test069: WTG run test032 as initial reference profile, sst = 300 K, wind = 3, no theta perturb, rhpert = +0.1, rhmode = 0. FINISHED. Test070: WTG run test032 as initial reference profile, sst = 300 K, wind = 5, no theta perturb, rhpert = +0.1, rhmode = 0. FINISHED. Test071: WTG run test032 as initial reference profile, sst = 300 K, wind = 7, no theta perturb, rhpert = +0.1, rhmode = 0. FINISHED. Test072: WTG run test032 as initial reference profile, sst = 300 K, wind = 10, no theta perturb, rhpert = +0.1, rhmode = 0. FINISHED. Test073: WTG run test032 as initial reference profile, sst = 300 K, wind = 15, no theta perturb, rhpert = +0.1, rhmode = 0. FINISHED. Test074: WTG run test032 as initial reference profile, sst = 300 K, wind = 20, no theta perturb, rhpert = +0.1, rhmode = 0. FINISHED. Test075: WTG run test032 as initial reference profile, sst = 300 K, wind = 7, no theta perturb, rhpert = -0.1, rhmode = 0. FINISHED. Test076: WTG run test032 as initial reference profile, sst = 300 K, wind = 10, no theta perturb, rhpert = -0.1, rhmode = 0. FINISHED. Test077: WTG run test032 as initial reference profile, sst = 300 K, wind = 15, no theta perturb, rhpert = -0.1, rhmode = 0. FINISHED. Test078: WTG run test032 as initial reference profile, sst = 300 K, wind = 20, no theta perturb, rhpert = -0.1, rhmode = 0. FINISHED. Test079: WTG run test032 as initial reference profile, sst = 299 K, wind = 7, no theta or rh perturbs. FINISHED. Test080: WTG run test032 as initial reference profile, sst = 301 K, wind = 7, no theta or rh perturbs. FINISHED. Test081: WTG run test032 as initial reference profile, sst = 302 K, wind = 7, no theta or rh perturbs. FINISHED. Test082: WTG run test032 as initial reference profile, sst = 303 K, wind = 7, no theta or rh perturbs. FINISHED. Let's do some tests with hscale = 3e4 m rather than 1.5e5 m. Reproduce tests 037-043 with the new value. Test083: Like test037 except hscale = 3e4 km, vy = 0 m/s, and entropy023. FINISHED. Test084: Like test037 except hscale = 3e4 km, vy = 3 m/s, and entropy023. FINISHED. Test085: Like test037 except hscale = 3e4 km, vy = 5 m/s, and entropy023. FINISHED. Test086: Like test037 except hscale = 3e4 km, vy = 7 m/s, and entropy023. FINISHED. Test087: Like test037 except hscale = 3e4 km, vy = 10 m/s, and entropy023. FINISHED. Test088: Like test037 except hscale = 3e4 km, vy = 15 m/s, and entropy023. FINISHED. Test089: Like test037 except hscale = 3e4 km, vy = 20 m/s, and entropy023. FINISHED. Let's run some RCE tests with strong radiative cooling a la Cohen and Craig. Test090: RCE like test032 except entropy023 and radiative cooling of 4 K/day, drag coeff = 0.002. Test091: RCE like test032 except entropy023 and radiative cooling of 8 K/day, drag coeff = 0.004. Test092: RCE like test032 except entropy023 and radiative cooling of 12 K/day, drag coeff = 0.006. Test093: RCE like test032 except entropy023 and radiative cooling of 16 K/day, drag coeff = 0.008. Test094: RCE repeat of test032 except entropy023. Test095: RCE like test032 except entropy023 and radiative cooling of 2 K/day, drag coeff = 0.001. +++++EVERYTHING UP TO THIS POINT DELETED. +++++Switch to entropy024, which has some significant changes. Test096: RCE like test094 except entropy024, nx = 192, nt = 1500. No profile file. +++++Test entropy023 vs. entropy024 interactive radiation. Test097: RCE like test094 except interactive radiation with entropy023. Test098: RCE like test094 except interactive radiation with entropy024. +++++Let's do some WTG runs with test096 as a reference profile. Use spectral WTG with hscale = 1.5e5 m, detain = 1, sobel = 0. First, omit th and rt perturbations and explore wind dependence. Test099: WTG using test096 as a reference, wind = 0. Entropy024. Test100: WTG using test096 as a reference, wind = 3. Entropy024. Test101: WTG using test096 as a reference, wind = 5. Entropy024. Test102: WTG using test096 as a reference, wind = 7. Entropy024. Test103: WTG using test096 as a reference, wind = 10. Entropy024. Test104: WTG using test096 as a reference, wind = 15. Entropy024. Test105: WTG using test096 as a reference, wind = 20. Entropy024. +++++Next, explore dependence on th perturbations at particular wind = 10 m/s. Test106: WTG using test096 as a reference, wind = 10, thm2 = -1. Entropy024. Test107: WTG using test096 as a reference, wind = 10, thm2 = -0.5. Entropy024. Test108: WTG using test096 as a reference, wind = 10, thm2 = 0.5. Entropy024. Test109: WTG using test096 as a reference, wind = 10, thm2 = 1. Entropy024. +++++Next, explore the effect of upper tropospheric perturbations in potential temperature. Test110: WTG using test096 as a reference, wind = 5, stemp = 300, like test101 except dthtrop = -1 K, pendepth = 5000 m. The rain just went to zero in this case! Test111: WTG using test096 as a reference, wind = 7, stemp = 300, like test102 except dthtrop = -3 K, pendepth = 10000 m. Test112: WTG using test096 as a reference, wind = 7, stemp = 300, like test102 except dthtrop = -1 K, pendepth = 10000 m. +++++The 7 m/s wind isn't enough to avoid a crash to zero rain, so let's increase to 10 m/s. Compare to test103. Test113: WTG using test096 as a reference, wind = 10, stemp = 300, like test103 except dthtrop = -1 K, pendepth = 10000 m. Entropy024. Test113x: Like test113 except using entropy025x with lt = 10 (effectively). Test113y: Like test113 except using entropy025 with lt = 100. Test113z: Like test113y except using entropy026 as a test. Test113a: Like test113y except using entropy027x as a test. (Trapezoid replaces simpson integration. Reduces nasty oscillations in thetaj.) Test113b: Like test113y except using entropy027y as a test. This has a major rewrite of the Fourier transform code in spectral WTG. Test114: WTG using test096 as a refernce, wind = 10, stemp = 300, like test103 except dthtrop = +1 K, pendepth = 10000 m. Entropy025 with lt = 10. Test115: WTG using test096 as a refernce, wind = 10, stemp = 300, like test103 except dthtrop = -2 K, pendepth = 5000 m. Entropy025 with lt = 10. Test116: WTG using test096 as a refernce, wind = 10, stemp = 300, like test103 except dthtrop = +2 K, pendepth = 5000 m. Entropy024. +++++Results: Test dthtrop pendep rain sfrac0 sfrac ii0 ii 103 0 0 11.0 0.66 0.73 20.5 20.4 113 -1 10000 13.6 0.67 0.73 19.9 20.2 113x -1 10000 15.2 0.67 0.73 19.9 20.1 113y -1 10000 13.1 0.67 0.73 19.9 20.1 114 +1 10000 7.4 0.65 0.70 21.0 20.9 115 -2 5000 11.9 0.67 0.72 20.7 20.6 116 +2 5000 9.4 0.65 0.72 20.2 20.2 +++++The potential temperature perturbation used above is non-zero at the surface. Use instead crmlinref.sh, which produces a linearly decreasing perturbation from the tropopause downward rather than an exponentially decreasing profile. Also, use entropy027. Finally, try a wind speed of 5 m/s. Test117: WTG using test096 as a reference, wind = 5, stemp = 300, dthtrop = 0 K, pendepth 10000 km. This repeats test101, but with entropy027. Kill. Test118: WTG like test117 except dthtrop = -2 K. Kill. +++++Something is odd with the test099 - test116 results: perhaps using an entropy024 RCE profile with later versions of entropy has problems. Let's generate a fresh RCE profile using 096 as a starting point. Test119: RCE starting with test096 but using entropy027, running for 10000000 seconds. DONE. +++++Examination suggests a systematic lowering of theta-e and theta-es of about 1 K at all altitudes between test096 and test119. This could explain the odd WTG results of tests 101 and 102, in which the precipitation rate crashed to zero. Let's now do some WTG tests like tests 099-105, but using a mean of 2000000-10000000 s of test119 as a reference profile. Use SST = 300 K, hscale = 150 km, spectral WTG. Use entropy027 except where noted. Note that a bug exists in the entropy027 output that makes the cumulative *wtg variables too small by a factor of 100. +++++Rerunning 120-126 using entropy029. DONE. Test120: WTG run using test119 as a reference, wind = 0 m/s. DONE. Test121: WTG run using test119 as a reference, wind = 3 m/s. DONE. Test122: WTG run using test119 as a reference, wind = 5 m/s. DONE. Test123: WTG run using test119 as a reference, wind = 7 m/s. DONE. Test124: WTG run using test119 as a reference, wind = 10 m/s. DONE. Test124x: Like test124 except using entropy029, which outputs vertical convective fluxes. TERMINATED AT HALFWAY POINT. Test125: WTG run using test119 as a reference, wind = 15 m/s. DONE. Test126: WTG run using test119 as a reference, wind = 20 m/s. DONE. +++++Let's explore the effect of thm2 = 1 on same range of winds and SST. Use entropy029 and test119 as a reference. Test127: WTG thm2 = 1, wind = 0 m/s. DONE. Test128: WTG thm2 = 1, wind = 3 m/s. DONE. Test129: WTG thm2 = 1, wind = 5 m/s. DONE. Test130: WTG thm2 = 1, wind = 7 m/s. DONE. Test131: WTG thm2 = 1, wind = 10 m/s. DONE. Test132: WTG thm2 = 1, wind = 15 m/s. DONE. Test133: WTG thm2 = 1, wind = 20 m/s. DONE. +++++Now for thm2 = -1 -- likely to have convection only for strong winds. Test134: WTG thm2 = -1, wind 7 m/s. DONE. Test135: WTG thm2 = -1, wind 10 m/s. DONE. Test136: WTG thm2 = -1, wind 15 m/s. DONE. Test137: WTG thm2 = -1, wind 20 m/s. DONE. +++++Let's make some relative humidity perturbations. Test138: WTG rhm0 = 0.05, wind = 0 m/s. DONE. Test139: WTG rhm0 = 0.05, wind = 3 m/s. DONE. Test140: WTG rhm0 = 0.05, wind = 5 m/s. DONE. Test141: WTG rhm0 = 0.05, wind = 7 m/s. DONE. Test142: WTG rhm0 = 0.05, wind = 10 m/s. DONE. Test143: WTG rhm0 = 0.05, wind = 15 m/s. DONE. Test144: WTG rhm0 = 0.05, wind = 20 m/s. DONE. Test145: WTG rhm0 = -0.05, wind = 7 m/s. DONE. Test146: WTG rhm0 = -0.05, wind = 10 m/s. DONE. Test147: WTG rhm0 = -0.05, wind = 15 m/s. DONE. Test148: WTG rhm0 = -0.05, wind = 20 m/s. DONE. +++++Let's now try sea surface temperature perturbations. First, +1 K. Test149: WTG SST = 301 K, wind = 0 m/s. DONE. Test150: WTG SST = 301 K, wind = 3 m/s. DONE. Test151: WTG SST = 301 K, wind = 5 m/s. DONE. Test152: WTG SST = 301 K, wind = 7 m/s. DONE. Test153: WTG SST = 301 K, wind = 10 m/s. DONE. Test154: WTG SST = 301 K, wind = 15 m/s. DONE. Test155: WTG SST = 301 K, wind = 20 m/s. DONE. +++++Next, -1 K. Test156: WTG SST = 299 K, wind = 7 m/s. DONE. Test157: WTG SST = 299 K, wind = 10 m/s. DONE. Test158: WTG SST = 299 K, wind = 15 m/s. DONE. Test159: WTG SST = 299 K, wind = 20 m/s. DONE. +++++Let's fill in our tests with thm0 runs. Test160: WTG thm0 = 1, wind = 0 m/s. DONE. Test161: WTG thm0 = 1, wind = 3 m/s. DONE. Test162: WTG thm0 = 1, wind = 5 m/s. DONE. Test163: WTG thm0 = 1, wind = 7 m/s. DONE. Test164: WTG thm0 = 1, wind = 10 m/s. DONE. Test165: WTG thm0 = 1, wind = 15 m/s. DONE. Test166: WTG thm0 = 1, wind = 20 m/s. DONE. Test167: WTG thm0 = -1, wind = 0 m/s. DONE. Test168: WTG thm0 = -1, wind = 3 m/s. DONE. Test169: WTG thm0 = -1, wind = 5 m/s. DONE. Test170: WTG thm0 = -1, wind = 7 m/s. DONE. Test171: WTG thm0 = -1, wind = 10 m/s. DONE. Test172: WTG thm0 = -1, wind = 15 m/s. DONE. Test173: WTG thm0 = -1, wind = 20 m/s. DONE. +++++Now do rhm2 cases. (Had to redo due to error.) Test174: WTG rhm2 = 0.05, wind = 0 m/s. DONE. Test175: WTG rhm2 = 0.05, wind = 3 m/s. DONE. Test176: WTG rhm2 = 0.05, wind = 5 m/s. DONE. Test177: WTG rhm2 = 0.05, wind = 7 m/s. DONE. Test178: WTG rhm2 = 0.05, wind = 10 m/s. DONE. Test179: WTG rhm2 = 0.05, wind = 15 m/s. DONE. Test180: WTG rhm2 = 0.05, wind = 20 m/s. DONE. Test181: WTG rhm2 = -0.05, wind = 0 m/s. DONE. Test182: WTG rhm2 = -0.05, wind = 3 m/s. DONE. Test183: WTG rhm2 = -0.05, wind = 5 m/s. DONE. Test184: WTG rhm2 = -0.05, wind = 7 m/s. DONE. Test185: WTG rhm2 = -0.05, wind = 10 m/s. DONE. Test186: WTG rhm2 = -0.05, wind = 15 m/s. DONE. Test187: WTG rhm2 = -0.05, wind = 20 m/s. DONE. +++++Let's try radiative cooling with same integrated cooling, different +++++shape. Moving the break point up or down leads to the following +++++equation: b = *(d_0 + h)/(d + h) where h is the tropopause height, +++++d_0 is the initial break height, d is the changed break height, +++++b_0 is the initial cooling rate, b is the final cooling rate +++++below the break point. So far, h = 15 km, d_0 = 12 km, b_0 = 2 K/day. +++++For d = 0 km, b = 2.8 K/day. Test188: A 1 time step RCE run with no profile file -- used to generate sentrad profile with d = 0 km, b = 2.8 K/day. Test189: Like test120 except newrad; vwind = 0 m/s. DELETE. Test190: Like test120 except newrad; vwind = 3 m/s. DELETE. Test191: Like test120 except newrad; vwind = 5 m/s. DELETE. Test192: Like test120 except newrad; vwind = 7 m/s. DELETE. Test193: Like test120 except newrad; vwind = 10 m/s. DELETE. Test194: Like test120 except newrad; vwind = 15 m/s. DELETE. Test195: Like test120 except newrad; vwind = 20 m/s. DELETE. +++++These results aren't too useful, since the vertically integrated +++++entropy sources don't match with the earlier calculations. Also, +++++the change in radiative cooling is very extreme and probably not +++++relevant. Delete. +++++Let's try some tests with high vertical wavenumber sine perturbations +++++in theta and rt. Test196: Starting with mean.119, using makeinput_sine.sh, with wavenumber = 10, thmval = 1 K, rhmval = 0, vwind = 10 m/s, sst = 300 K. Compare with test124. The potential temperature perturbation actually kills the convection in this case! DELETE. Test197: Starting with mean.119, using makeinput_sine.sh, with wavenumber = 10, thmval = 0.5 K, rhmval = 0, vwind = 10 m/s, sst = 300 K. Compare with test124. Test198: Starting with mean.119, using makeinput_sine.sh, with wavenumber = 10, thmval = -0.5 K, rhmval = 0, vwind = 10 m/s, sst = 300 K. Compare with test124. +++++Now do some tests with max theta perturbations at the tropopause, +++++decaying exponentially downward. Test199: Starting with mean.119, using makeinput_trop.sh, with dthtrop = 1 K, pendepth = 4000 m, vwind = 10 m/s, sst = 300 K. Compare with test 124. Test200: Starting with mean.119, using makeinput_trop.sh, with dthtrop = -1 K, pendepth = 4000 m, vwind = 10 m/s, sst = 300 K. Compare with test 124. +++++Let's finish off the initial sequence with runs with thm4 and rhm4. Test201: WTG with thm4 = 0.5, wind = 0 m/s, sst = 300 K. DONE. Test202: WTG with thm4 = 0.5, wind = 3 m/s, sst = 300 K. DONE. Test203: WTG with thm4 = 0.5, wind = 5 m/s, sst = 300 K. DONE. Test204: WTG with thm4 = 0.5, wind = 7 m/s, sst = 300 K. DONE. Test205: WTG with thm4 = 0.5, wind = 10 m/s, sst = 300 K. DONE. Test206: WTG with thm4 = 0.5, wind = 15 m/s, sst = 300 K. DONE. Test207: WTG with thm4 = 0.5, wind = 20 m/s, sst = 300 K. DONE. Test208: WTG with thm4 = -0.5, wind = 0 m/s, sst = 300 K. DONE. Test209: WTG with thm4 = -0.5, wind = 3 m/s, sst = 300 K. DONE. Test210: WTG with thm4 = -0.5, wind = 5 m/s, sst = 300 K. DONE. Test211: WTG with thm4 = -0.5, wind = 7 m/s, sst = 300 K. DONE. Test212: WTG with thm4 = -0.5, wind = 10 m/s, sst = 300 K. DONE. Test213: WTG with thm4 = -0.5, wind = 15 m/s, sst = 300 K. DONE. Test214: WTG with thm4 = -0.5, wind = 20 m/s, sst = 300 K. DONE. Test215: WTG with rhm4 = 0.03, wind = 0 m/s, sst = 300 K. DONE. Test216: WTG with rhm4 = 0.03, wind = 3 m/s, sst = 300 K. DONE. Test217: WTG with rhm4 = 0.03, wind = 5 m/s, sst = 300 K. DONE. Test218: WTG with rhm4 = 0.03, wind = 7 m/s, sst = 300 K. DONE. Test219: WTG with rhm4 = 0.03, wind = 10 m/s, sst = 300 K. DONE. Test220: WTG with rhm4 = 0.03, wind = 15 m/s, sst = 300 K. DONE. Test221: WTG with rhm4 = 0.03, wind = 20 m/s, sst = 300 K. DONE. Test222: WTG with rhm4 = -0.03, wind = 0 m/s, sst = 300 K. DONE. Test223: WTG with rhm4 = -0.03, wind = 3 m/s, sst = 300 K. DONE. Test224: WTG with rhm4 = -0.03, wind = 5 m/s, sst = 300 K. DONE. Test225: WTG with rhm4 = -0.03, wind = 7 m/s, sst = 300 K. DONE. Test226: WTG with rhm4 = -0.03, wind = 10 m/s, sst = 300 K. DONE. Test227: WTG with rhm4 = -0.03, wind = 15 m/s, sst = 300 K. DONE. Test228: WTG with rhm4 = -0.03, wind = 20 m/s, sst = 300 K. DONE. +++++Let's do a 3-D test for comparison with WRF. Test229: RCE like test119 except 3-D with 64 by 64 horizontal grid points. (Use 096 as a reference profile). +++++We need to fill in some non-precipitating cases so that the +++++fitting routines will have something to work on. Test230: WTG thm0 = -1, wind 0 m/s. DONE. Test231: WTG thm0 = -1, wind 3 m/s. DONE. Test232: WTG thm0 = -1, wind 5 m/s. DONE. Test233: WTG rhm0 = -0.05, wind = 0 m/s. DONE. Test234: WTG rhm0 = -0.05, wind = 3 m/s. DONE. Test235: WTG rhm0 = -0.05, wind = 5 m/s. DONE. Test236: WTG SST = 299 K, wind = 0 m/s. DONE. Test237: WTG SST = 299 K, wind = 3 m/s. DONE. Test238: WTG SST = 299 K, wind = 5 m/s. DONE. +++++We need to look at some cases with stronger positive moisture +++++perturbations. Start with mean119. Test239: WTG rhm0 = +0.1 wind 0 m/s. DONE. Test240: WTG rhm0 = +0.1 wind 3 m/s. DONE. Test241: WTG rhm0 = +0.1 wind 5 m/s. DONE. Test242: WTG rhm0 = +0.1 wind 7 m/s. DONE. Test243: WTG rhm0 = +0.1 wind 10 m/s. DONE. Test244: WTG rhm0 = +0.1 wind 15 m/s. DONE. Test245: WTG rhm0 = +0.1 wind 20 m/s. DONE. +++++We also need stronger potential temperature perturbations. Test246: WTG thm2 = +2 wind 0 m/s. DONE. Test247: WTG thm2 = +2 wind 3 m/s. DONE. Test248: WTG thm2 = +2 wind 5 m/s. DONE. Test249: WTG thm2 = +2 wind 7 m/s. DONE. Test250: WTG thm2 = +2 wind 10 m/s. DONE. Test251: WTG thm2 = +2 wind 15 m/s. DONE. Test252: WTG thm2 = +2 wind 20 m/s. DONE. +++++Let's try some strong humidity perturbations that keep the surface relative humidity unchanged but dry mid-levels. Test253: WTG rhm0 = -0.15, rhm2 = -0.2, wind = 0 m/s. DONE. Test254: WTG rhm0 = -0.15, rhm2 = -0.2, wind = 3 m/s. DONE. Test255: WTG rhm0 = -0.15, rhm2 = -0.2, wind = 5 m/s. DONE. Test256: WTG rhm0 = -0.15, rhm2 = -0.2, wind = 7 m/s. DONE. Test257: WTG rhm0 = -0.15, rhm2 = -0.2, wind = 10 m/s. DONE. Test258: WTG rhm0 = -0.15, rhm2 = -0.2, wind = 15 m/s. DONE. Test259: WTG rhm0 = -0.15, rhm2 = -0.2, wind = 20 m/s. DONE. +++++Try to resolve differences with Saska's runs. Areas of +++++difference are: +++++ 1. Use of interactive rather than fixed radiation. (Saska finds +++++ little difference. +++++ 2. Different initial reference profile. Shouldn't be important +++++ if satfrac0 and instabindex0 are the same. +++++ 3. Saska has wtermi = 1 m/s -- I have 5 m/s. +++++ 4. Saska has C_d = 0.0012, I have C_d = 0.001. +++++ 5. Saska has hscale = 4.77e4, I have hscale = 1.5e5. +++++ 6. Saska has drelax = 0.001, I have 0.01. Shouldn't be an issue. +++++ 7. Saska has 256 x grid pts, I have 192. Probably not important. Test260: WTG rhm* = 0, thm* = 0, wind = 10 m/s, fixed radiation, SST 300 K. Like test 124 except wtermi = 1. Test261: Like test 260 except C_d = 0.0012. Test262: Like test 260 except hscale = 4.77e4. Test263: Like test 260 except interactive radiation. Table of results: run params entflux iindex0 sfrac0 rain rain-fit 124 control 1.040 18.11 0.689 21.0 21.9 260 wtermi = 1 0.926 18.11 0.688 7.89 17.5 261 C_d = 0.0012 1.301 18.11 0.689 31.4 32.0 262 hscale = 4.77e4 1.003 18.11 0.689 17.1 20.5 263 interactive rad 0.980 18.11 0.689 22.0 19.6 +++++We need to fill in some more non-precipitating cases so that the +++++fitting routines will have something to work on. Test264: WTG thm2 = -1, wind 0 m/s. Test265: WTG thm2 = -1, wind 3 m/s. Test266: WTG thm2 = -1, wind 5 m/s.