|
Last change
on this file since 867 was
716,
checked in by rwordsworth, 12 years ago
|
|
Mainly updates to radiative transfer and gas management scheme.
Most CIA data now read from standard HITRAN datafiles. For the H2O
continuum, two options have been added: the standard CKD continuum,
and the empirical formula in PPC (Pierrehumbert 2010). Use the toggle
'H2Ocont_simple' in callphys.def to choose.
Note to Martians: I've changed the default values of 'sedimentation' and
'co2cond' in inifis to false. Both these are defined in the standard deftank
callphys.def file, so there shouldn't be any compatibility problems.
|
|
File size:
460 bytes
|
| Rev | Line | |
|---|
| [305] | 1 | !------------------------------------------------------------------------- |
|---|
| [716] | 2 | subroutine bilinear(f,f11,f21,f12,f22,x,x1,x2,y,y1,y2) |
|---|
| 3 | ! Used for interpolation of continuum data |
|---|
| [305] | 4 | |
|---|
| [716] | 5 | implicit none |
|---|
| [305] | 6 | |
|---|
| [716] | 7 | real*8 x,y,x1,x2,y1,y2 |
|---|
| 8 | real*8 f,f11,f12,f21,f22,fA,fB |
|---|
| [305] | 9 | |
|---|
| [716] | 10 | ! 1st in x-direction |
|---|
| 11 | fA=f11*(x2-x)/(x2-x1)+f21*(x-x1)/(x2-x1) |
|---|
| 12 | fB=f12*(x2-x)/(x2-x1)+f22*(x-x1)/(x2-x1) |
|---|
| [305] | 13 | |
|---|
| [716] | 14 | ! then in y-direction |
|---|
| 15 | f=fA*(y2-y)/(y2-y1)+fB*(y-y1)/(y2-y1) |
|---|
| [305] | 16 | |
|---|
| [716] | 17 | return |
|---|
| 18 | end subroutine bilinear |
|---|
Note: See
TracBrowser
for help on using the repository browser.
