|
Last change
on this file since 1157 was
875,
checked in by jleconte, 12 years ago
|
11/02/2013 == JL
- Updated moist convection scheme to handle situations with a large water vapor content
- Added a keyword to enable ocean runoff in callphys.def (activerunoff)
|
|
File size:
769 bytes
|
| Line | |
|---|
| 1 | subroutine su_watercycle |
|---|
| 2 | |
|---|
| 3 | use watercommon_h |
|---|
| 4 | implicit none |
|---|
| 5 | #include "comcstfi.h" |
|---|
| 6 | #include "callkeys.h" |
|---|
| 7 | |
|---|
| 8 | |
|---|
| 9 | !================================================================== |
|---|
| 10 | ! |
|---|
| 11 | ! Purpose |
|---|
| 12 | ! ------- |
|---|
| 13 | ! Set up relevant constants and parameters for the water cycle, and water cloud properties |
|---|
| 14 | ! |
|---|
| 15 | ! Authors |
|---|
| 16 | ! ------- |
|---|
| 17 | ! Robin Wordsworth (2010) |
|---|
| 18 | ! Jeremy Leconte (2012) |
|---|
| 19 | ! |
|---|
| 20 | !================================================================== |
|---|
| 21 | |
|---|
| 22 | epsi = mH2O / mugaz |
|---|
| 23 | RCPD = cpp |
|---|
| 24 | |
|---|
| 25 | !RV = 1000.*R/mH2O |
|---|
| 26 | RV = 1000.*8.314/mH2O ! caution! R is R/mugaz already! |
|---|
| 27 | |
|---|
| 28 | RCPV = 1.88e3 ! specific heat capacity of water vapor at 350K |
|---|
| 29 | |
|---|
| 30 | RVTMP2 = RCPV/RCPD-1. ! not currently used... |
|---|
| 31 | |
|---|
| 32 | end subroutine su_watercycle |
|---|
Note: See
TracBrowser
for help on using the repository browser.
