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 ****************************** |
---|