Ticket #619: depresys_control.mf77

File depresys_control.mf77, 2.6 KB (added by mx020105, 10 years ago)
Line 
1*/ Amanda mod: remember this is ppmm not ppmv
2*/
3*/ this is an ID for a mod
4*ID DEPRESYS_CONT
5*/
6*/
7*DECLARE S_PHYSCS
8*/
9*/ Pass latitude info through to R2_SET_GAS_MIX_RATIO routine
10*/
11*I S_PHYSCS.1732
12     &, sinlat
13*/
14*I S_PHYSCS.1948
15     &, sinlat
16*/
17*DECLARE RAD_CTL1
18*/
19*I ADB2F404.995
20     &, SIN_TRUE_LATITUDE(FIRST_POINT)
21*/
22*I ADB2F404.1025
23     &, SIN_TRUE_LATITUDE(FP_LOCAL(I))
24*/
25*DECLARE LWRAD3A
26*/
27*I ADB1F402.491
28     &, SIN_LATITUDE
29
30*I ADB1F402.500
31     &  , SIN_LATITUDE(NPD_FIELD)
32*/
33*/ Switch order of routines to provide pressure date input
34*/ to R2_SET_GAS_MIX_RATIO routine called from LWRAD3A
35*/
36*I LWRAD3A.434
37!     CALCULATE PRESSURES AND TEMPERATURES.
38      CALL R2_SET_THERMODYNAMIC(N_PROFILE, NLEVS, I_GATHER, .TRUE.
39     &   , PSTAR, TSTAR, AB, BB, AC, BC, PEXNER, TAC
40     &   , P, T, T_BDY, T_SURFACE, D_MASS
41     &   , NPD_FIELD, NPD_PROFILE, NPD_LAYER
42     &   )
43*/
44*/ Add extra input data latitude, p and tropopause level
45*/ to R2_SET_GAS_MIX_RATIO
46*/
47*I ADB1F405.353
48     &   , SIN_LATITUDE, P, TRINDX
49*/
50*D ADB1F405.354
51*D LWRAD3A.451
52*D ADB1F401.556
53*D ADB1F401.557
54*D LWRAD3A.454
55*D LWRAD3A.455
56*/
57*DECLARE SWRAD3A
58*/
59*/ Repeat same calls as for LWRAD3A
60*/
61*I ADB1F402.702
62     &   , SIN_LATITUDE
63*/
64*I ADB2F404.1529
65     &   , SIN_LATITUDE(NPD_FIELD)
66*/
67*I SWRAD3A.548
68!     SET THE THERMODYNAMIC PROPERTIES OF THE ATMOSPHERE.
69      CALL R2_SET_THERMODYNAMIC(NLIT, NLEVS, LIST, .FALSE.
70     &   , PSTAR, DUMMY, AB, BB, AC, BC
71     &   , DUMMY, TAC
72     &   , P, T, DUMMY, DUMMY, D_MASS
73     &   , NPD_FIELD, NPD_PROFILE, NPD_LAYER
74     &   )
75*/                                                             
76*I ADB1F405.953
77     &   , SIN_LATITUDE, P, TRINDX
78*/
79*D SWRAD3A.562
80*D SWRAD3A.563
81*D ADB2F404.1571
82*D ADB2F404.1572
83*D ADB2F404.1573
84*D SWRAD3A.567
85*D SWRAD3A.568
86*/
87*/
88*/ R2_SET_GAS_MIX_RATIO routine
89*/
90*DECLARE FILL3A
91*/
92*I ADB1F405.18
93     &   , SIN_LATITUDE, P, TRINDX
94*/
95*I FILL3A.79
96     &   , TRINDX(NPD_FIELD)
97*/
98*I  ADB1F405.34
99     &   , SIN_LATITUDE(NPD_FIELD)
100     &   , P(NPD_PROFILE, 0: NPD_LAYER)
101*/
102*I FILL3A.145
103     &   , GG !Tropopause level
104*/
105*D H3DBQLIM.1
106*/
107*I FILL3A.205
108        GG=TRINDX(LG)
109*/ If not in stratosphere leave H2O alone
110*/
111        IF(I .GT. NLEVS+1-GG) THEN
112             GAS_MIX_RATIO(L, I, IUMP_H2O)=H2O(LG, NLEVS-I+1)
113        ELSE
114*/ Set fixed SWV concentration of ~4ppmv.
115             H2OLMN=2.5E-6
116             GAS_MIX_RATIO(L, I, IUMP_H2O)=
117     &                        H2OLMN
118        ENDIF
119*/
120*D  H3DBQLIM.2
121*D  H3DBQLIM.3
122*/   
123*/  ************** end of modset ******************************