Climate under Cover
Digital Dynamic Simulation in Plant Bio-Engineering

Preface for 1st edition

Tokyo, October 1992 | Tadashi Takakura

The first edition framed mulching, row covers, floating mulch, and greenhouses as one family of covering techniques, emphasizing soil and boundary-air processes, nonlinear dynamic behavior, and computer simulation as a research and teaching tool.

Preface for 2nd edition

Nagasaki and Taipei, January 2002 | Tadashi Takakura and Wei Fang

The second edition explained the shift from CSMP to MATLAB in response to modern Windows environments, while keeping the goal of making all models accessible for readers to run, study, and extend.

Web edition

Taipei, April 2026 | Wei Fang

The web edition translates the MATLAB models into browser-based tools so the simulations can be opened, studied, and shared more easily in a modern online environment.

New Independent Project

CUC Soil Temperature Simulator

This browser tool is adapted from Cuc01.m and soil01.m to simulate how temperatures in five soil layers evolve over time under a surface-temperature boundary condition.

The simulator preserves the four original Matlab trial settings and also lets you enter custom parameters for comparison.

Source logic: Cuc01.m + soil01.m | Web version implemented with offline JavaScript for static hosting

Inputs

Case Selection and Parameters

Smaller dt gives a finer time grid for the plotted curves.

Results

Temperature Curves

Summary Table

Final and Peak Temperatures

Layer Final Temp (簞C) Minimum (簞C) Maximum (簞C) Peak Time (hr)

Model Notes

How the Matlab logic was mapped

  • The four original trial settings are preserved exactly; case 1 to case 4 only change ks and cs.
  • The governing factor matches soil01.m: val = ks / z^2 / cs.
  • The bottom boundary remains fixed at TBL, while the surface boundary changes sinusoidally over time.
  • The browser version preserves the same five-layer structure and boundary conditions as the Matlab model.
  • TF = T0 + TU * sin(ω(t - 8)), where ω = 2π / 24, t is time in hr, and T0 and TU are temperatures in °C.
  • val has units of 1/hr; ks is soil thermal conductivity, cs is volumetric heat capacity, and z is layer thickness in m.

Equation Set

Five-layer ODE structure

  • dT1/dt = [2(TF - T1) + (T2 - T1)] * val
  • dT2/dt = [(T1 - T2) + (T3 - T2)] * val
  • dT3/dt = [(T2 - T3) + (T4 - T3)] * val
  • dT4/dt = [(T3 - T4) + (T5 - T4)] * val
  • dT5/dt = [(T4 - T5) + 2(TBL - T5)] * val
  • dTi/dt is the temperature rate of change in °C/hr, while T1 to T5, TF, and TBL are temperatures in °C.
  • The five states represent soil nodes from surface to deep layer; the doubled top and bottom terms enforce the surface and bottom boundary conditions.

Soil Layer Diagram

Soil Layer Schematic

This schematic summarizes the five-layer soil structure used in Cuc01.m and soil01.m.

The top boundary is the surface temperature TF, the internal nodes are T1 to T5, and the bottom boundary is fixed at TBL.

Each layer thickness is represented by z, matching the spacing used in the original Matlab formulation.

TF T1 T2 T3 T4 T5 TBL * * * * * * * z z z z z z z

Recommended start file: cucindex.html