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

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

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

Revision 1.3 - (hide 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 cnh 1.3 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 adcroft 1.2
4     #include "CPP_OPTIONS.h"
5    
6 cnh 1.3 SUBROUTINE INPHYS (FSG,HSG,RLAT,myThid)
7 adcroft 1.2
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 cnh 1.3 INTEGER NLON, NLAT, NLEV, NGP
27 adcroft 1.2 PARAMETER ( NLON=IX, NLAT=IL, NLEV=KX, NGP=NLON*NLAT )
28 cnh 1.3
29     INTEGER myThid
30 adcroft 1.2 C
31     C
32    
33     #include "GRID.h"
34     #include "EEPARAMS.h"
35     #include "PARAMS.h"
36 cnh 1.3
37     #include "Lev_def.h"
38 adcroft 1.2 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 cnh 1.3 INTEGER K
57 adcroft 1.2 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 cnh 1.3 IF ( myThid .EQ. 1 ) THEN
64 adcroft 1.2 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 cnh 1.3 WVI(K,2)=(LOG(SIGH(K))-SIGL(K))*WVI(K,1)
95 adcroft 1.2 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 cnh 1.3 ENDIF
211 adcroft 1.2
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