| 1 | edhill | 1.1 | function experiment | 
| 2 |  |  |  | 
| 3 |  |  | global Nx Ny Nz | 
| 4 |  |  | global dx dy lat long | 
| 5 |  |  | global dz depth | 
| 6 |  |  | global delt iteration it | 
| 7 |  |  | global descriptor | 
| 8 |  |  | global this_path | 
| 9 |  |  | global g Cp rho_bar alpha | 
| 10 |  |  | global f0 B0 Q0 beta r_heat r_expt H | 
| 11 |  |  |  | 
| 12 |  |  | % prescribed parameters | 
| 13 |  |  | Nx = 124; Ny = 124; Nz = 30; | 
| 14 |  |  | dx = 0.01; dy = 0.01; dz = 0.005; | 
| 15 |  |  | delt = 0.1; | 
| 16 |  |  | iteration = ['0000250209'; ... | 
| 17 |  |  | '0000251210'; ... | 
| 18 |  |  | '0000253609'; ... | 
| 19 |  |  | '0000254610'; ... | 
| 20 |  |  | '0000255610'; ... | 
| 21 |  |  | '0000256610'; ... | 
| 22 |  |  | '0000257009'; ... | 
| 23 |  |  | '0000258010'; ... | 
| 24 |  |  | '0000259010'; ... | 
| 25 |  |  | '0000260010'; ... | 
| 26 |  |  | '0000260409'; ... | 
| 27 |  |  | '0000261410'; ... | 
| 28 |  |  | '0000262410'; ... | 
| 29 |  |  | '0000263410'; ... | 
| 30 |  |  | '0000264411'; ... | 
| 31 |  |  | '0000265411'; ... | 
| 32 |  |  | '0000266411'; ... | 
| 33 |  |  | '0000266810'; ... | 
| 34 |  |  | '0000267811'; ... | 
| 35 |  |  | '0000268811'; ... | 
| 36 |  |  | '0000269811'; ... | 
| 37 |  |  | '0000270210'; ... | 
| 38 |  |  | '0000274611'; ... | 
| 39 |  |  | '0000275611'; ... | 
| 40 |  |  | '0000276611'; ... | 
| 41 |  |  | '0000277010'; ... | 
| 42 |  |  | '0000278011'; ... | 
| 43 |  |  | '0000279011'; ... | 
| 44 |  |  | '0000280011'; ... | 
| 45 |  |  | '0000280410'; ... | 
| 46 |  |  | '0000281411'; ... | 
| 47 |  |  | '0000282411'; ... | 
| 48 |  |  | '0000283411'; ... | 
| 49 |  |  | '0000283910'; ... | 
| 50 |  |  | '0000284911'; ... | 
| 51 |  |  | '0000285911'; ... | 
| 52 |  |  | '0000286911'; ... | 
| 53 |  |  | '0000287410'; ... | 
| 54 |  |  | '0000288411'; ... | 
| 55 |  |  | '0000289411'; ... | 
| 56 |  |  | '0000290411'; ... | 
| 57 |  |  | '0000290910'; ... | 
| 58 |  |  | '0000291911'; ... | 
| 59 |  |  | '0000292911'; ... | 
| 60 |  |  | '0000293911'; ... | 
| 61 |  |  | '0000294410'; ... | 
| 62 |  |  | '0000295411'; ... | 
| 63 |  |  | '0000296411'; ... | 
| 64 |  |  | '0000297411'; ... | 
| 65 |  |  | '0000297910'; ... | 
| 66 |  |  | '0000298911'; ... | 
| 67 |  |  | '0000299911'; ... | 
| 68 |  |  | '0000300911'; ... | 
| 69 |  |  | '0000301410'; ... | 
| 70 |  |  | '0000302411'; ... | 
| 71 |  |  | '0000303411'; ... | 
| 72 |  |  | '0000304411'; ... | 
| 73 |  |  | '0000304810'; ... | 
| 74 |  |  | '0000305811'; ... | 
| 75 |  |  | '0000306811'; ... | 
| 76 |  |  | '0000307811'; ... | 
| 77 |  |  | '0000308310'; ... | 
| 78 |  |  | '0000310311'; ... | 
| 79 |  |  | '0000311810'; ... | 
| 80 |  |  | '0000312811'; ... | 
| 81 |  |  | '0000313811'; ... | 
| 82 |  |  | '0000314811'; ... | 
| 83 |  |  | '0000316211'; ... | 
| 84 |  |  | '0000317211'; ... | 
| 85 |  |  | '0000318211'; ... | 
| 86 |  |  | '0000318510'; ... | 
| 87 |  |  | '0000319511'; ... | 
| 88 |  |  | '0000320511'; ... | 
| 89 |  |  | '0000321511'; ... | 
| 90 |  |  | '0000321810'; ... | 
| 91 |  |  | '0000322811'; ... | 
| 92 |  |  | '0000323811'; ... | 
| 93 |  |  | '0000324811'; ... | 
| 94 |  |  | '0000325010'; ... | 
| 95 |  |  | '0000326011'; ... | 
| 96 |  |  | '0000327011'; ... | 
| 97 |  |  | '0000328011'; ... | 
| 98 |  |  | '0000329511'; ... | 
| 99 |  |  | '0000330511'; ... | 
| 100 |  |  | '0000331511'; ... | 
| 101 |  |  | '0000331910'; ... | 
| 102 |  |  | '0000332911'; ... | 
| 103 |  |  | '0000334911'; ... | 
| 104 |  |  | '0000335310'; ... | 
| 105 |  |  | '0000336311'; ... | 
| 106 |  |  | '0000337311'; ... | 
| 107 |  |  | '0000338610'; ... | 
| 108 |  |  | '0000339611']; | 
| 109 |  |  | descriptor = 'reference donut'; | 
| 110 |  |  | this_path  = ['/d04/helen/tank/R15donutprep03/dumps']; | 
| 111 |  |  | g = 9.81; Cp = 4000; rho_bar = 999.8; alpha = 2e-4; | 
| 112 |  |  | f0 = 0.75; Q0 = 500; beta = 0; r_heat = 0.2; r_expt = 0.48; | 
| 113 |  |  |  | 
| 114 |  |  | % derived  parameters | 
| 115 |  |  | lat = [1:1:Nx]; | 
| 116 |  |  | long = [1:1:Ny]; | 
| 117 |  |  | % depth = -cumsum(dz)+dz(1)/2; | 
| 118 |  |  | %b0 = g*alpha*Q0/rho_bar*Cp; | 
| 119 |  |  | % H = depth(Nz)-dz(Nz)/2; | 
| 120 |  |  |  | 
| 121 |  |  |  | 
| 122 |  |  |  | 
| 123 |  |  |  |