/[MITgcm]/MITgcm/pkg/aim/ini_inphys.F
ViewVC logotype

Contents of /MITgcm/pkg/aim/ini_inphys.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph

Revision 1.3 - (show annotations) (download)
Tue May 29 19:28:53 2001 UTC (23 years ago) by cnh
Branch: MAIN
CVS Tags: ecco_c44_e19, ecco_c44_e18, ecco_c44_e17, ecco_c44_e16, checkpoint40pre7, checkpoint43a-release1mods, checkpoint44e_post, release1_p12, release1_p13, release1_p10, release1_p11, release1_p16, release1_p17, release1_p14, release1_p15, checkpoint46g_pre, release1_p13_pre, chkpt44a_pre, ecco_c44_e23, ecco_c44_e20, ecco_c44_e21, ecco_c44_e26, ecco_c44_e27, ecco_c44_e24, ecco_c44_e25, icebear5, icebear4, icebear3, icebear2, checkpoint46f_post, checkpoint46d_pre, checkpoint46e_post, release1-branch_tutorials, checkpoint46c_post, checkpoint46b_post, ecco-branch-mod1, checkpoint46e_pre, ecco-branch-mod3, ecco-branch-mod4, ecco-branch-mod5, checkpoint45d_post, release1_beta1, checkpoint45b_post, checkpoint44g_post, release1-branch-end, release1_final_v1, checkpoint44b_pre, release1_p12_pre, checkpoint46c_pre, checkpoint43, checkpoint40, checkpoint41, checkpoint46, checkpoint44, checkpoint45, checkpoint44f_post, checkpoint40pre3, checkpoint40pre2, checkpoint40pre1, checkpoint44b_post, checkpoint40pre6, checkpoint40pre5, checkpoint46h_post, checkpoint40pre9, checkpoint40pre8, checkpoint40pre4, checkpoint44h_pre, release1_b1, chkpt44d_post, checkpoint46h_pre, release1_p8, release1_p9, release1_p1, release1_p2, release1_p3, release1_p4, release1_p5, release1_p6, release1_p7, checkpoint46a_post, chkpt44a_post, chkpt44c_post, checkpoint44f_pre, checkpoint46a_pre, checkpoint45c_post, checkpoint44e_pre, checkpoint46b_pre, ecco_ice2, ecco_ice1, checkpoint44h_post, checkpoint46d_post, checkpoint46g_post, checkpoint45a_post, ecco-branch-mod2, checkpoint42, release1-branch_branchpoint, ecco_c44_e22, release1_chkpt44d_post, chkpt44c_pre
Branch point for: release1_final, c24_e25_ice, ecco-branch, release1, release1_coupled, icebear, release1-branch, release1_50yr
Changes since 1.2: +12 -40 lines 
Updates for multi-threaded AIM with support for both latlon
and CS.
Needs compatible changes to verfication/
1 C $Header: /u/gcmpack/models/MITgcmUV/pkg/aim/ini_inphys.F,v 1.2 2001/02/02 21:36:29 adcroft Exp $
2 C $Name: $
3
4 #include "CPP_OPTIONS.h"
5
6 SUBROUTINE INPHYS (FSG,HSG,RLAT,myThid)
7
8 IMPLICIT rEAL*8 ( A-H,O-Z)
9
10 C--
11 C-- SUBROUTINE INPHYS (FSG,HSG,RLAT)
12 C--
13 C-- Purpose: Initialize common blocks for physical parametrization routines
14 C-- Input : FSG : sigma at full levels
15 C-- HSG : sigma at half levels
16 C-- RLAT : gaussian-grid latitudes
17 C-- Initialized common blocks: PHYCON, FSIGMU, VDICON,
18 C-- FORCON, SFLCON, CNVCON, LSCCON, RADCON
19 C--
20 C Resolution parameters
21 C
22 C
23 #include "atparam.h"
24 #include "atparam1.h"
25 C
26 INTEGER NLON, NLAT, NLEV, NGP
27 PARAMETER ( NLON=IX, NLAT=IL, NLEV=KX, NGP=NLON*NLAT )
28
29 INTEGER myThid
30 C
31 C
32
33 #include "GRID.h"
34 #include "EEPARAMS.h"
35 #include "PARAMS.h"
36
37 #include "Lev_def.h"
38 C
39 C Physical constants + functions of sigma and latitude
40 C
41 #include "com_physcon.h"
42 C
43 C Constants for sub-grid-scale physics
44 C
45 #include "com_forcon.h"
46 #include "com_sflcon.h"
47 #include "com_cnvcon.h"
48 #include "com_lsccon.h"
49 #include "com_radcon.h"
50 #include "com_vdicon.h"
51 C
52 C == Routine Arguments ==
53 REAL FSG(NLEV), HSG(0:NLEV), RLAT(NLAT)
54
55 C == Local Variables ==
56 INTEGER K
57 DATA Pground /1. _d +5/
58
59 C--- 1. Time independent parameters and arrays
60 C
61 C 1.1 Physical constants
62 C
63 IF ( myThid .EQ. 1 ) THEN
64 P0 = 1. _d +5
65 GG = 9.81 _d 0
66 RD = 287. _d 0
67 CP = 1004. _d 0
68 C Latent heat is in J/g for consistency with spec.hum. in g/Kg
69 ALHC = 2501. _d 0
70 SBC = 5.67 _d -8
71 C
72 C 1.2 Functions of sigma and latitude
73 C
74 SIGH(0)=HSG(0)
75 C
76 DO K=1,NLEV
77 SIG(K) = FSG(K)
78 Cg77 SIGL(K) = DLOG(FSG(K))
79 SIGL(K) = LOG(FSG(K))
80 SIGH(K) = HSG(K)
81 DSIG(K) = HSG(K)-HSG(K-1)
82 C POUT(K) = PRLEV(FSG(K))
83 GRDSIG(K) = GG/(DSIG(K)*P0)
84 GRDSCP(K) = GRDSIG(K)/CP
85 ENDDO
86 C
87 C Weights for vertical interpolation at half-levels(1,nlev) and surface
88 C Note that for phys.par. half-lev(k) is between full-lev k and k+1
89 C Fhalf(k) = Ffull(k)+WVI(K,2)*(Ffull(k+1)-Ffull(k))
90 C Fsurf = Ffull(nlev)+WVI(nlev,2)*(Ffull(nlev)-Ffull(nlev-1))
91 C
92 DO K=1,NLEV-1
93 WVI(K,1)=1./(SIGL(K+1)-SIGL(K))
94 WVI(K,2)=(LOG(SIGH(K))-SIGL(K))*WVI(K,1)
95 ENDDO
96 C
97 WVI(NLEV,1)=0.
98 WVI(NLEV,2)=-SIGL(NLEV)*WVI(NLEV-1,2)
99 cchdbg
100 WVI(NLEV,2)=-SIGL(NLEV)*WVI(NLEV-1,1)
101 cchdbg
102 C
103 C--- 2. Constants for boundary forcing (common FORCON):
104 C
105 SOLC = 342. _d 0
106
107 ALBSEA = 0.07 _d 0
108 ALBICE = 0.60 _d 0
109 ALBSN = 0.60 _d 0
110 SDALB = 60. _d 0
111 SWCAP = 75. _d 0
112 SWWIL = 0. _d 0
113 C
114 C--- 3. Constants for surface fluxes (common SFLCON):
115 C
116 FWIND0 = 0.6 _d 0
117 cchdbg *****************************************************
118 FWIND0 = 0.6 _d 0*0.7 _d 0
119 cchdbg ******************************************************
120 cchdbg FTEMP0 = 1.
121 cchdbg FHUM0 = 1.
122 FTEMP0 = 0. _d 0
123 FHUM0 = 0.5 _d 0
124 FHUM0 = 1. _d 0
125
126 c CDL = 3.0E-3
127 c CDS = 1.2E-3
128 c CHL = 1.2E-3
129 c CHS = 1.2E-3
130 CDL = 1.8 _d -3
131 CDS = 0.8 _d -3
132 CHL = 2.0 _d -3
133 CHS = 0.9 _d -3
134 C
135 VGUST = 4. _d 0
136 DTHETAF = 3. _d 0
137 FSTAB = 0.67 _d 0
138 C
139 C--- 4. Constants for convection (common CNVCON):
140 C
141 RHBL = 0.8 _d 0
142 TRCNV = 6. _d 0
143 ENTMAX = 0.5 _d 0
144 C
145 C--- 5. Constants for large-scale condensation (common LSCCON):
146 C
147 RHLSC = 0.9 _d 0
148 TRLSC = 4. _d 0
149 C
150 C--- 6. Constants for radiation (common RADCON):
151 C
152 ABSSW = 0.06 _d 0
153 cchdbg ********************
154 ABSSW = 0.04 _d 0
155 ABSSW = 0.06 _d 0
156 cchdbg *********************
157 ABSLW = 1.40 _d 0
158 cchdbg ****************************************************
159 cch ABSLW = 1.00
160 cchdbg ****************************************************
161
162 ABWSW = 0.06 _d 0
163 ABWLW = 0.55 _d 0
164 cchdbg ****************************************************
165 ABWLW = 0.60 _d 0
166 cchdbg ABWLW = 0.65
167 cchdbg ABWLW = 0.90
168 cchdbg ****************************************************
169
170 ABCSW = 0.0 _d 0
171 ABCLW = 0.15 _d 0
172 ABCLW = 0.15 _d 0
173 cchdbg *******************************************
174 c ABCLW = 3.0
175 c ABCLW = 2.0
176 cch ABCLW = 3.5
177 cch ABCLW = 2.5
178 cch ABCLW = 2.8
179 cchdbg ********************************************
180
181 EPSSW = 0.03 _d 0
182 EPSLW = 0.05 _d 0
183
184 ALBCL = 0.50 _d 0
185 cchdbg *******************************************
186 ALBCL = 0.49 _d 0
187 ALBCL = 0.47 _d 0
188 ALBCL = 0.40 _d 0
189 ALBCL = 0.35 _d 0
190 ALBCL = 0.33 _d 0
191 ALBCL = 0.36 _d 0
192 ALBCL = 0.34 _d 0
193 cchdbg ********************************************
194 cchdbg RHCL1 = 0.45
195 cchdbg RHCL2 = 0.85
196 cchdbg QACL = 1.00
197 cchdbg *******************************************
198 RHCL1 = 0.50 _d 0
199 RHCL2 = 0.90 _d 0
200 QACL = 0.50 _d 0
201 cchdbg *******************************************
202 C
203 C 2.6 Constants for vertical diffusion and sh. conv. (common VDICON):
204 C
205 TRVDI = 48. _d 0
206 TRSHC = 24. _d 0
207 cchdbg *******************************************
208 TRSHC = 12. _d 0
209 cchdbg *******************************************
210 ENDIF
211
212 CcnhDebugStarts
213 C write(0,*) 'MAXVAL(NLEVxy)=',MAXVAL(NLEVxy)
214 C write(0,*) 'MINVAL(NLEVxy)=',MINVAL(NLEVxy)
215 CcnhDebugEnds
216 C
217 RETURN
218 END
  ViewVC Help
Powered by ViewVC 1.1.22