| 1 | 101c101,103 |
|---|
| 2 | < |
|---|
| 3 | --- |
|---|
| 4 | > !****Mars |
|---|
| 5 | > REAL :: x_param,y_param,rho_param,dilat |
|---|
| 6 | > !****Mars |
|---|
| 7 | 171,172c173,180 |
|---|
| 8 | < hm = 100. |
|---|
| 9 | < xa = 5.0 |
|---|
| 10 | --- |
|---|
| 11 | > ! hm = 100. |
|---|
| 12 | > ! xa = 5.0 |
|---|
| 13 | > !****Mars |
|---|
| 14 | > !! hm = 4000. |
|---|
| 15 | > !! xa = 5.0 |
|---|
| 16 | > hm = 2000. |
|---|
| 17 | > xa = 6.0 |
|---|
| 18 | > !****Mars |
|---|
| 19 | 174a183,191 |
|---|
| 20 | > !****Mars |
|---|
| 21 | > jcm = jde/2 |
|---|
| 22 | > ! |
|---|
| 23 | > ! xa1 = 5000./500. |
|---|
| 24 | > ! xal1 = 4000./500. |
|---|
| 25 | > ! pii = 2.*asin(1.0) |
|---|
| 26 | > ! hm1 = 250. |
|---|
| 27 | > !! hm1 = 1000. |
|---|
| 28 | > !****Mars |
|---|
| 29 | 177,184c194,197 |
|---|
| 30 | < xa1 = 5000./500. |
|---|
| 31 | < xal1 = 4000./500. |
|---|
| 32 | < pii = 2.*asin(1.0) |
|---|
| 33 | < hm1 = 250. |
|---|
| 34 | < ! hm1 = 1000. |
|---|
| 35 | < |
|---|
| 36 | < |
|---|
| 37 | < stretch_grid = .true. |
|---|
| 38 | --- |
|---|
| 39 | > !****Mars |
|---|
| 40 | > ! stretch_grid = .true. |
|---|
| 41 | > stretch_grid = .false. |
|---|
| 42 | > !****Mars |
|---|
| 43 | 301a315,336 |
|---|
| 44 | > |
|---|
| 45 | > !****Mars |
|---|
| 46 | > !!3D hill |
|---|
| 47 | > ! grid%ht(i,j) = hm/(1.+(float(i-icm)/xa)**2+(float(j-jcm)/xa)**2) |
|---|
| 48 | > !!3D crater |
|---|
| 49 | > ! grid%ht(i,j) = hm - hm/(1.+(float(i-icm)/xa)**2+(float(j-jcm)/xa)**2) |
|---|
| 50 | > !!3D crater w/ rims |
|---|
| 51 | > x_param = float(i-icm) |
|---|
| 52 | > y_param = float(j-jcm) |
|---|
| 53 | > dilat = xa/2 |
|---|
| 54 | > rho_param = sqrt(x_param**2 + y_param**2) |
|---|
| 55 | > |
|---|
| 56 | > ! revolution surface ; seed is a fourth order polynom |
|---|
| 57 | > grid%ht(i,j) = (rho_param+6*dilat)*(rho_param+10*dilat) |
|---|
| 58 | > grid%ht(i,j) = (rho_param-6*dilat)*(rho_param-10*dilat)*grid%ht(i,j) |
|---|
| 59 | > ! flat terrain elsewhere - smooth gradient (no abrupt fall) |
|---|
| 60 | > grid%ht(i,j) = grid%ht(i,j)*(tanh(rho_param+7*dilat)/2 - tanh(rho_param-7*dilat)/2) |
|---|
| 61 | > grid%ht(i,j) = hm - (hm*.4/1500)*grid%ht(i,j)/(dilat**4) |
|---|
| 62 | > |
|---|
| 63 | > !if (rho_param .GE. dilat*10) ht(i,j) = hm |
|---|
| 64 | > |
|---|
| 65 | > |
|---|
| 66 | 580c615,618 |
|---|
| 67 | < parameter( debug = .false.) |
|---|
| 68 | --- |
|---|
| 69 | > |
|---|
| 70 | > ! parameter( debug = .false.) |
|---|
| 71 | > !****Mars |
|---|
| 72 | > parameter( debug = .true.) |
|---|
| 73 | 596c634,637 |
|---|
| 74 | < parameter (p1000mb = 1.e+05, r = 287, cp = 1003., cv = cp-r, cvpm = -cv/cp, g=9.81 ) |
|---|
| 75 | --- |
|---|
| 76 | > !****Mars |
|---|
| 77 | > ! parameter (p1000mb = 1.e+05, r = 287, cp = 1003., cv = cp-r, cvpm = -cv/cp, g=9.81 ) |
|---|
| 78 | > parameter (p1000mb = 610, r = 192, cp = 845., cv = cp-r, cvpm = -cv/cp,g=3.72) |
|---|
| 79 | > !****Mars |
|---|
