UPP (develop)
Loading...
Searching...
No Matches
cuparm.f
Go to the documentation of this file.
1
3 module cuparm_mod
4
5 implicit none
6
7 real, parameter :: h1=1.e0, &
8 h1d5=1.5e0, &
9 h2d5=2.5e0, &
10 h3000=3000.e0, &
11 h10e5=100000.e0, &
12 d00=0.e0, &
13 d125=.125e0, &
14 d50=.5e0, &
15 d608=.608e0, &
16 g=9.8e0, &
17 cp=1004.6e0, &
18 capa=0.28589641e0, &
19 rog=287.04/9.8, &
20 elwv=2.50e6, &
21 elivw=2.72e6, &
22 row=1.e3, &
23 epsq=2.e-12, &
24 a2=17.2693882e0, &
25 a3=273.16e0, &
26 a4=35.86e0, &
27 t0=273.16e0, &
28 t1=274.16e0, &
29 pq0=379.90516e0, &
30 stresh=1.10e0, &
31 stabs=1.0e0, &
32 stabd=.90e0, &
33 stabfc=1.00e0, &
34 dttop=0.0e0, &
35!---VVVVV
36 rhf=0.10, &
37 epsup=1.00, &
38 epsdn=1.05, &
39 epsth=0.0, &
40 pbm=13000., &
41 pqm=20000., &
42 pno=1000., &
43 pone=2500., &
44 zsh=2000., &
45 pfrz=15000., &
46 pshu=45000., &
47
48!AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
49!
50! AUGUST '91: SCHEME HAVING THE OPTION OF USING DIFFERENT FAST AND
51! SLOW PROFILES FOR SEA AND FOR LAND POINTS; AND ALSO THE "SEA" AND
52! THE "LAND" SCHEME EVERYWHERE. OVER LAND PROFILES DEPART FROM THE
53! FAST (DRY) PROFILES ONLY FOR PRECIPITATION/TIME STEP >
54! A PRESCRIBED VALUE (CURRENTLY, IN THE VERSION #3 DONE WEDNESDAY
55! 18 SEPTEMBER, 1/4 INCH/24 H). USE OF VARIOUS SWITCHES AS FOLLOWS.
56!
57! THE "OLD" ("HARD", =ZAVISA OCT.1990) LAND SCHEME WITH FIXED
58! LAND PROFILES IS RUN BY
59! * SETTING OCT90=.TRUE. IN THE FIRST EXECUTABLE LINE FOLLOWING
60! THESE COMMENTS (THIS REACTIVATES EFI=H1 OVER LAND IF
61! .NOT.UNIS, AND CLDEFI(K)=EFIMN AS THE LEFTOVER CLDEFI
62! VALUE AT SWAP POINTS);
63! * DEFINING FAST LAND PROFILES SAME AS FAST SEA PROFILES (OR BY
64! CHOOSING ANOTHER SET OF LAND PROFILES ZAVISA USED AT
65! EARLIER TIME);
66! * SETTING FSL=1.; AND
67! * DEFINING STEFI (STARTING EFI) EQUAL TO AVGEFI.
68! (THE LAST THREE POINTS ARE HANDLED BY SWITCHING AROUND THE
69! "CFM" COMMENTS AT TWO PLACES)
70!
71! THE "OLD,OLD" (APPR. ORIGINAL BETTS) SCHEME IS RUN BY
72! * SPECIFYING UNIL=.TRUE.;
73! * SETTING FSL=1.;
74! * SETTING OCT90=.TRUE.
75! (WITH THESE SETTINGS FAST LAND PROFILES ONLY ARE USED).
76! FM
77 fss=.85e0, &
78 efimn=.20e0, &
79 efmnt=.70e0, &
80 fcc=.50, &
81 fcp=h1-fcc, &
82!
83! IN THIS VERSION 3.5, OVER LAND AND FOR THE FAST PROFILES, DSPB
84! IS PRESCRIBED TO BE 25 PERCENT DRIER THAN THE FAST SEA VALUE
85! (IN ROUGH AGREEMENT WITH BINDER, QJ, IN PRESS) WHILE DSP0 AND
86! DSPT ARE EACH 20 PERCENT DRIER THAN THE CORRESPONDING FAST
87! SEA VALUES. WITH FSL=.875 THIS MAKES THE
88! AVERAGE OF THE FAST AND THE SLOW LAND PROFILES SOMEWHAT DRIER
89! THAN THE OCT90 FIXED LAND PROFILES. FM
90!
95!
96 dspbfl=-4843.75e0, &
97 dsp0fl=-7050.00e0, &
98 dsptfl=-2250.0e0, &
99 fsl=.850e0, &
100!*** ACTIVATE THE FOLLOWING LINE IF OCT90=.TRUE. (AND COMMENT OUT THE
101!*** PRECEDING LINE):
102! DSPBFL=-3875.E0, & !< Deficit Saturation Pressure (DSP) at cloud base - fast land ?
103! DSP0FL=-5875.E0, & !< Deficit Saturation Pressure (DSP) at freezing level - fast land ?
104! DSPTFL=-1875.E0, & !< Deficit Saturation Pressure (DSP) at cloud top - fast land ?
105! FSL=1.0E0, & !< Factor to obtain "slow" profiles over land ?
106 dspbfs=-3875.e0, &
107 dsp0fs=-5875.e0, &
108 dsptfs=-1875.e0, &
109 dspbsl=dspbfl*fsl, &
110 dsp0sl=dsp0fl*fsl, &
111 dsptsl=dsptfl*fsl, &
112 dspbss=dspbfs*fss, &
113 dsp0ss=dsp0fs*fss, &
114 dsptss=dsptfs*fss, &
115!*** NEW CONVECTION SCHEME WITH CROSSING DSP PROFILES ******************
116!+- &, UNIS=.FALSE.,EFIMN=.71E0,EFMNT=.71,FCC=0.5,FCP=H1-FCC
117!+- &, DSPBL=-3875.E0,DSP0L=-5875.E0,DSPTL=-1875.E0
118!+- &, DSPBS=-2875.E0,DSP0S=-5125.E0,DSPTS=-4875.E0
119!+- &, DSPBF=-4375.E0,DSP0F=-4375.E0,DSPTF=-1000.E0
120!*** BETTS CONVECTION SCHEME *******************************************
121! &, UNIS=.FALSE.,EFIMN=.9999E0,EFMNT=.9999E0,FCC=.50,FCP=H1-FCC
122! &, DSPBL=-3875.E0,DSP0L=-5875.E0,DSPTL=-1875.E0
123! &, DSPBF=-3875.E0,DSP0F=-5875.E0,DSPTF=-1875.E0
124! &, DSPBS=-3875.E0,DSP0S=-5875.E0,DSPTS=-1875.E0
125!***********************************************************************
126 trel=3000., &
127 epsntp=.0010e0, &
128 efifc=5.0e0, &
129 avgefi=(efimn+1.e0)*.5e0, &
130 dspc=-3000.e0, &
131 epsp=1.e-7, &
132 stefi=1.e0, &
133!*** ACTIVATE THE FOLLOWING LINE AND COMMENT OUT THE PRECEDING LINE IF
134!*** OCT90=.TRUE.
135! STEFI=AVGEFI, &
136 slopbl=(dspbfl-dspbsl)/(h1-efimn), &
137 slop0l=(dsp0fl-dsp0sl)/(h1-efimn), &
138 sloptl=(dsptfl-dsptsl)/(h1-efimn), &
139 slopbs=(dspbfs-dspbss)/(h1-efimn), &
140 slop0s=(dsp0fs-dsp0ss)/(h1-efimn), &
141 slopts=(dsptfs-dsptss)/(h1-efimn), &
142 slope=(h1 -efmnt)/(h1-efimn)
143 real, parameter :: &
144 a23m4l=a2*(a3-a4)*elwv, &
145 elocp=elivw/cp, &
146 cprlg=cp/(row*g*elwv), &
147 rcp=h1/cp
148 logical,parameter :: &
149 unis=.false., &
150 unil=.false., &
151 oct90=.false.
152 end module cuparm_mod
153