1 | SUBROUTINE SUGFL |
---|
2 | |
---|
3 | USE PARKIND1 ,ONLY : JPIM ,JPRB |
---|
4 | USE YOMHOOK ,ONLY : LHOOK, DR_HOOK |
---|
5 | |
---|
6 | USE YOMLUN , ONLY : NULOUT |
---|
7 | USE YOM_YGFL , ONLY : YGFL,YGFLC,YL, YI, YA, YS, YR, YG, YTKE,& |
---|
8 | & YO3, YSRC, YEXT, YQ, YCVGQ, JPGFL,YQVA, & |
---|
9 | & YL_NL, YI_NL, YA_NL, YO3_NL, YSRC_NL, YEXT_NL, YQ_NL,YQVA_NL, & |
---|
10 | & YS_NL, YR_NL, YG_NL, YTKE_NL,YCVGQ_NL,NGFL_EXT,& |
---|
11 | & YCPF , YSPF ,YCPF_NL, YSPF_NL, YAERO, YAERO_NL, NAERO, YGRG, YGRG_NL, NGRG,& |
---|
12 | & YSDSAT, YCVV, YSDSAT_NL, YCVV_NL, YFORC, YFORC_NL, NGFL_FORC,& |
---|
13 | & YEZDIAG, YEZDIAG_NL, NGFL_EZDIAG, & |
---|
14 | & YGHG, YGHG_NL, NGHG, YTRAC, YTRAC_NL, NTRAC, & |
---|
15 | & YUOM ,YUAL ,YDOM ,YDAL ,YUEN ,YUNEBH, & |
---|
16 | & YUOM_NL ,YUAL_NL ,YDOM_NL ,YDAL_NL ,YUEN_NL ,YUNEBH_NL |
---|
17 | USE GFL_SUBS ,ONLY : DEFINE_GFL_COMP |
---|
18 | USE YOMDIM , ONLY : LADER |
---|
19 | USE YOMCST , ONLY : RV ,RCPV ,RCW ,RCS |
---|
20 | USE YOMGRB , ONLY : NGRBQ, NGRBO3, NGRBCLWC, NGRBCIWC,& |
---|
21 | & NGRBCC, NGRB149 |
---|
22 | USE YOMCT0 , ONLY : LFPART2,NCONF |
---|
23 | USE YOMFA , ONLY : YFAQ ,YFAL ,YFAI ,YFAO3 ,& |
---|
24 | & YFAS ,YFAR ,YFAG ,YFATKE ,& |
---|
25 | & YFACPF, YFASPF, YFACLF,YFASRC, YFACVGQ, YFASDSAT, YFACVV ,& |
---|
26 | & YFAQVA, & |
---|
27 | & YFAUOM, YFAUAL, YFADOM, YFADAL, YFAUEN, YFAUNEBH |
---|
28 | |
---|
29 | !**** *SUGFL* - Initialize definition of unified_treatment grid_point fields |
---|
30 | |
---|
31 | ! Purpose. |
---|
32 | ! -------- |
---|
33 | ! Initialize definition of unified_treatment fields (GFL) |
---|
34 | ! The GFL ordering in this routine is the ordering for "t" and |
---|
35 | ! "t+dt" values of the GFL in grid-point arrays (GFL and GFLT1), |
---|
36 | ! and also the ordering in the spectral array SPGFL. |
---|
37 | |
---|
38 | !** Interface. |
---|
39 | ! ---------- |
---|
40 | ! *CALL* *SUGFL |
---|
41 | |
---|
42 | ! Explicit arguments : |
---|
43 | ! -------------------- |
---|
44 | |
---|
45 | ! Implicit arguments : |
---|
46 | ! -------------------- |
---|
47 | ! MODULE YOM_YGFL |
---|
48 | |
---|
49 | ! Method. |
---|
50 | ! ------- |
---|
51 | ! See documentation |
---|
52 | |
---|
53 | ! Externals. |
---|
54 | ! ---------- |
---|
55 | ! None |
---|
56 | |
---|
57 | ! Reference. |
---|
58 | ! ---------- |
---|
59 | ! ECMWF Research Department documentation of the IFS |
---|
60 | |
---|
61 | ! Author. |
---|
62 | ! ------- |
---|
63 | ! Mats Hamrud *ECMWF* |
---|
64 | |
---|
65 | ! Modifications. |
---|
66 | ! -------------- |
---|
67 | ! Original : 2002-03-11 |
---|
68 | ! M.Hamrud : 2002-08-01 Extensive mods |
---|
69 | ! R. El Khatib : 2003-08-19 ARPEGE field names |
---|
70 | ! M.Hamrud 01-Oct-2003 CY28 Cleaning |
---|
71 | ! M. Tudor : 2003-10-17 introduce YCPF and YSPF in GFL |
---|
72 | ! A. Untch : 2005-03-11 introduce aerosols |
---|
73 | ! J. Flemming : 2005-04-11 replace aerosols with GRG |
---|
74 | ! Y. Seity : 2004-11-16 AROME GFL fields |
---|
75 | ! C. Fischer: 23-May-2005 No Ozone in the Meteo-France control variable |
---|
76 | ! Y. Bouteloup: 13-Oct-2005 introduce YCVGQ (Moisture convergence for French physics) |
---|
77 | ! J. Haseler : 2005-10-11 introduce LTRAJIO |
---|
78 | ! JJMorcrette 20060512 GEMS variables in cy31 |
---|
79 | ! S. Serrar 20060907 tracers used for diagnostics only (in GFL between GHG and GRG) |
---|
80 | ! A. Alias : 2006-10-13 introduction of YSDSAT and YCVV |
---|
81 | ! B. Sass : HIRLAM pseudo-prognostic field YQVA |
---|
82 | ! M. Bellus : 27-Sep-2006 introduce YUOM, YUAL, YDOM, YDAL, YUEN and |
---|
83 | ! YUNEBH (ALARO-0 prognostic convection) + spotted/corrected bug |
---|
84 | ! in index incrementation after YO3 |
---|
85 | ! ------------------------------------------------------------------ |
---|
86 | |
---|
87 | IMPLICIT NONE |
---|
88 | |
---|
89 | INTEGER(KIND=JPIM) :: JGFL,IGFLPTR,INCR |
---|
90 | LOGICAL :: LLDERS,LLT1,LL5 |
---|
91 | REAL(KIND=JPRB) :: ZHOOK_HANDLE |
---|
92 | |
---|
93 | !------------------------------------------------------------------------- |
---|
94 | |
---|
95 | #include "abor1.intfb.h" |
---|
96 | |
---|
97 | !------------------------------------------------------------------------- |
---|
98 | |
---|
99 | IF (LHOOK) CALL DR_HOOK('SUGFL',0,ZHOOK_HANDLE) |
---|
100 | |
---|
101 | !------------------------------------------------------------------------- |
---|
102 | |
---|
103 | IF( .NOT.LFPART2) THEN |
---|
104 | |
---|
105 | ! 1. CASE LFPART2=F |
---|
106 | |
---|
107 | ! 1.1 Initial settings. |
---|
108 | |
---|
109 | YGFL%NUMFLDS=0 |
---|
110 | |
---|
111 | LL5=NCONF/100 == 1.OR.NCONF == 801.OR.NCONF == 601.OR. & |
---|
112 | & NCONF == 401.OR.NCONF == 501 |
---|
113 | LLT1=.TRUE. |
---|
114 | ! All gridpoint fields have to be set up before the spectral ones |
---|
115 | ! (i.e. part 1.2 SHOULD be done before part 1.3) |
---|
116 | LLDERS=.FALSE. |
---|
117 | ! * Simple GFL variables: |
---|
118 | ! (order should be the same in parts 1.1, 1.2 and 1.3). |
---|
119 | INCR=0 |
---|
120 | YQ => YGFLC(JPGFL-INCR) |
---|
121 | INCR=INCR+1 |
---|
122 | YL => YGFLC(JPGFL-INCR) |
---|
123 | INCR=INCR+1 |
---|
124 | YI => YGFLC(JPGFL-INCR) |
---|
125 | INCR=INCR+1 |
---|
126 | YS => YGFLC(JPGFL-INCR) |
---|
127 | INCR=INCR+1 |
---|
128 | YR => YGFLC(JPGFL-INCR) |
---|
129 | INCR=INCR+1 |
---|
130 | YG => YGFLC(JPGFL-INCR) |
---|
131 | INCR=INCR+1 |
---|
132 | YTKE => YGFLC(JPGFL-INCR) |
---|
133 | INCR=INCR+1 |
---|
134 | YA => YGFLC(JPGFL-INCR) |
---|
135 | INCR=INCR+1 |
---|
136 | YO3 => YGFLC(JPGFL-INCR) |
---|
137 | INCR=INCR+1 |
---|
138 | YCPF => YGFLC(JPGFL-INCR) |
---|
139 | INCR=INCR+1 |
---|
140 | YSPF => YGFLC(JPGFL-INCR) |
---|
141 | INCR=INCR+1 |
---|
142 | YSRC => YGFLC(JPGFL-INCR) |
---|
143 | INCR=INCR+1 |
---|
144 | YCVGQ => YGFLC(JPGFL-INCR) |
---|
145 | INCR=INCR+1 |
---|
146 | YSDSAT => YGFLC(JPGFL-INCR) |
---|
147 | INCR=INCR+1 |
---|
148 | YCVV => YGFLC(JPGFL-INCR) |
---|
149 | INCR=INCR+1 |
---|
150 | YQVA => YGFLC(JPGFL-INCR) |
---|
151 | INCR=INCR+1 |
---|
152 | YUOM => YGFLC(JPGFL-INCR) |
---|
153 | INCR=INCR+1 |
---|
154 | YUAL => YGFLC(JPGFL-INCR) |
---|
155 | INCR=INCR+1 |
---|
156 | YDOM => YGFLC(JPGFL-INCR) |
---|
157 | INCR=INCR+1 |
---|
158 | YDAL => YGFLC(JPGFL-INCR) |
---|
159 | INCR=INCR+1 |
---|
160 | YUEN => YGFLC(JPGFL-INCR) |
---|
161 | INCR=INCR+1 |
---|
162 | YUNEBH => YGFLC(JPGFL-INCR) |
---|
163 | INCR=INCR+1 |
---|
164 | ! * Groups of multiple GFL variables: |
---|
165 | ! (order should be the same in parts 1.1, 1.2 and 1.3). |
---|
166 | IF (NGFL_EXT == 0) THEN |
---|
167 | YEXT => YGFLC(JPGFL-INCR:JPGFL-INCR) |
---|
168 | INCR=INCR+1 |
---|
169 | ENDIF |
---|
170 | IF (NGFL_FORC == 0) THEN |
---|
171 | YFORC => YGFLC(JPGFL-INCR:JPGFL-INCR) |
---|
172 | INCR=INCR+1 |
---|
173 | ENDIF |
---|
174 | IF (NGFL_EZDIAG == 0) THEN |
---|
175 | YEZDIAG => YGFLC(JPGFL-INCR:JPGFL-INCR) |
---|
176 | INCR=INCR+1 |
---|
177 | ENDIF |
---|
178 | IF (NGHG == 0) THEN |
---|
179 | YGHG => YGFLC(JPGFL-INCR:JPGFL-INCR) |
---|
180 | INCR=INCR+1 |
---|
181 | ENDIF |
---|
182 | IF (NTRAC == 0) THEN |
---|
183 | YTRAC => YGFLC(JPGFL-INCR:JPGFL-INCR) |
---|
184 | INCR=INCR+1 |
---|
185 | ENDIF |
---|
186 | IF (NGRG == 0) THEN |
---|
187 | YGRG => YGFLC(JPGFL-INCR:JPGFL-INCR) |
---|
188 | INCR=INCR+1 |
---|
189 | ENDIF |
---|
190 | IF (NAERO == 0) THEN |
---|
191 | YAERO => YGFLC(JPGFL-INCR:JPGFL-INCR) |
---|
192 | INCR=INCR+1 |
---|
193 | ENDIF |
---|
194 | |
---|
195 | ! 1.2 Grid-point GFL. |
---|
196 | |
---|
197 | IF(YQ_NL%LGP) THEN |
---|
198 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
199 | YQ=>YGFLC(IGFLPTR) |
---|
200 | CALL DEFINE_GFL_COMP(YDGFLC=YQ ,CDNAME=YFAQ%CLNAME,KGRIB=NGRBQ,LDGP=.TRUE.,& |
---|
201 | & KREQIN=YQ_NL%NREQIN, & |
---|
202 | & PREFVALI=YQ_NL%REFVALI,LDREQOUT=YQ_NL%LREQOUT,LDERS=LLDERS,LD5=LL5,& |
---|
203 | & LDT1=LLT1,LDGPINGP=YQ_NL%LGPINGP, & |
---|
204 | & LDTRAJIO=YQ_NL%LTRAJIO,LDTHERMACT=.TRUE.,PR=RV,PRCP=RCPV) |
---|
205 | ENDIF |
---|
206 | IF(YL_NL%LGP) THEN |
---|
207 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
208 | YL=>YGFLC(IGFLPTR) |
---|
209 | CALL DEFINE_GFL_COMP(YDGFLC=YL ,CDNAME=YFAL%CLNAME,KGRIB=NGRBCLWC, & |
---|
210 | & LDGP=.TRUE.,KREQIN=YL_NL%NREQIN, & |
---|
211 | & PREFVALI=YL_NL%REFVALI,LDREQOUT=YL_NL%LREQOUT,LDERS=LLDERS,LD5=LL5, & |
---|
212 | & LDT1=LLT1,LDTHERMACT=.TRUE.,PR=0.0_JPRB,PRCP=RCW) |
---|
213 | ENDIF |
---|
214 | IF(YI_NL%LGP) THEN |
---|
215 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
216 | YI=>YGFLC(IGFLPTR) |
---|
217 | CALL DEFINE_GFL_COMP(YDGFLC=YI ,CDNAME=YFAI%CLNAME,KGRIB=NGRBCIWC,& |
---|
218 | & LDGP=.TRUE.,KREQIN=YI_NL%NREQIN, & |
---|
219 | & PREFVALI=YI_NL%REFVALI,LDREQOUT=YI_NL%LREQOUT,LDERS=LLDERS,LD5=LL5, & |
---|
220 | & LDT1=LLT1,LDTHERMACT=.TRUE.,PR=0.0_JPRB,PRCP=RCS) |
---|
221 | ENDIF |
---|
222 | IF(YS_NL%LGP) THEN |
---|
223 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
224 | YS=>YGFLC(IGFLPTR) |
---|
225 | CALL DEFINE_GFL_COMP(YDGFLC=YS ,CDNAME=YFAS%CLNAME,KGRIB=NGRB149, & |
---|
226 | & LDGP=.TRUE.,KREQIN=YS_NL%NREQIN, & |
---|
227 | & PREFVALI=YS_NL%REFVALI,LDREQOUT=YS_NL%LREQOUT,LDERS=LLDERS,LD5=LL5, & |
---|
228 | & LDT1=LLT1,LDTHERMACT=.TRUE.,PR=0.0_JPRB,PRCP=RCS) |
---|
229 | ENDIF |
---|
230 | IF(YR_NL%LGP) THEN |
---|
231 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
232 | YR=>YGFLC(IGFLPTR) |
---|
233 | CALL DEFINE_GFL_COMP(YDGFLC=YR ,CDNAME=YFAR%CLNAME,KGRIB=NGRB149,& |
---|
234 | & LDGP=.TRUE.,KREQIN=YR_NL%NREQIN, & |
---|
235 | & PREFVALI=YR_NL%REFVALI,LDREQOUT=YR_NL%LREQOUT,LDERS=LLDERS,LD5=LL5, & |
---|
236 | & LDT1=LLT1,LDTHERMACT=.TRUE.,PR=0.0_JPRB,PRCP=RCW) |
---|
237 | ENDIF |
---|
238 | IF(YG_NL%LGP) THEN |
---|
239 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
240 | YG=>YGFLC(IGFLPTR) |
---|
241 | CALL DEFINE_GFL_COMP(YDGFLC=YG ,CDNAME=YFAG%CLNAME,KGRIB=NGRB149,& |
---|
242 | & LDGP=.TRUE.,KREQIN=YG_NL%NREQIN, & |
---|
243 | & PREFVALI=YG_NL%REFVALI,LDREQOUT=YG_NL%LREQOUT,LDERS=LLDERS,LD5=LL5, & |
---|
244 | & LDT1=LLT1,LDTHERMACT=.TRUE.,PR=0.0_JPRB,PRCP=RCS) |
---|
245 | ENDIF |
---|
246 | IF(YTKE_NL%LGP) THEN |
---|
247 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
248 | YTKE=>YGFLC(IGFLPTR) |
---|
249 | CALL DEFINE_GFL_COMP(YDGFLC=YTKE ,CDNAME=YFATKE%CLNAME,KGRIB=NGRB149,& |
---|
250 | & LDGP=.TRUE.,KREQIN=YTKE_NL%NREQIN, & |
---|
251 | & PREFVALI=YTKE_NL%REFVALI,LDREQOUT=YTKE_NL%LREQOUT,LDERS=LLDERS,LD5=LL5,& |
---|
252 | & LDT1=LLT1) |
---|
253 | ENDIF |
---|
254 | IF(YA_NL%LGP)THEN |
---|
255 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
256 | YA=>YGFLC(IGFLPTR) |
---|
257 | CALL DEFINE_GFL_COMP(YDGFLC=YA ,CDNAME=YFACLF%CLNAME,KGRIB=NGRBCC,& |
---|
258 | & LDGP=.TRUE.,KREQIN=YA_NL%NREQIN, & |
---|
259 | & PREFVALI=YA_NL%REFVALI,LDREQOUT=YA_NL%LREQOUT,LDERS=LLDERS,LD5=LL5,& |
---|
260 | & LDT1=LLT1) |
---|
261 | ENDIF |
---|
262 | IF(YO3_NL%LGP)THEN |
---|
263 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
264 | YO3=>YGFLC(IGFLPTR) |
---|
265 | CALL DEFINE_GFL_COMP(YDGFLC=YO3,CDNAME=YFAO3%CLNAME,KGRIB=NGRBO3,& |
---|
266 | & LDGP=.TRUE.,KREQIN=YO3_NL%NREQIN, & |
---|
267 | & PREFVALI=YO3_NL%REFVALI,LDREQOUT=YO3_NL%LREQOUT,LDERS=LLDERS,LD5=LL5,& |
---|
268 | & LDT1=LLT1,LDGPINGP=YO3_NL%LGPINGP, & |
---|
269 | & LDTRAJIO=YO3_NL%LTRAJIO) |
---|
270 | ENDIF |
---|
271 | IF(YCPF_NL%LGP)THEN |
---|
272 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
273 | YCPF=>YGFLC(IGFLPTR) |
---|
274 | CALL DEFINE_GFL_COMP(YDGFLC=YCPF,CDNAME=YFACPF%CLNAME,& |
---|
275 | & KGRIB=YCPF_NL%IGRBCODE,LDGP=.TRUE., & |
---|
276 | & KREQIN=YCPF_NL%NREQIN,PREFVALI= YCPF_NL%REFVALI,& |
---|
277 | & LDREQOUT=YCPF_NL%LREQOUT,LDERS=.FALSE.,LD5=LL5, & |
---|
278 | & LDT1=LLT1,LDGPINGP=.TRUE.) |
---|
279 | ENDIF |
---|
280 | IF(YSPF_NL%LGP)THEN |
---|
281 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
282 | YSPF=>YGFLC(IGFLPTR) |
---|
283 | CALL DEFINE_GFL_COMP(YDGFLC=YSPF,CDNAME=YFASPF%CLNAME,& |
---|
284 | & KGRIB=YSPF_NL%IGRBCODE,LDGP=.TRUE., & |
---|
285 | & KREQIN=YSPF_NL%NREQIN, PREFVALI=YSPF_NL%REFVALI,& |
---|
286 | & LDREQOUT=YSPF_NL%LREQOUT,LDERS=.FALSE.,LD5=LL5, & |
---|
287 | & LDT1=LLT1,LDGPINGP=.TRUE.) |
---|
288 | ENDIF |
---|
289 | IF(YSRC_NL%LGP) THEN |
---|
290 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
291 | YSRC=>YGFLC(IGFLPTR) |
---|
292 | CALL DEFINE_GFL_COMP(YDGFLC=YSRC ,CDNAME=YFASRC%CLNAME,KGRIB=NGRB149,& |
---|
293 | & LDGP=.TRUE.,KREQIN=YSRC_NL%NREQIN, & |
---|
294 | & PREFVALI=YSRC_NL%REFVALI,LDREQOUT=YSRC_NL%LREQOUT,LDERS=LLDERS,& |
---|
295 | & LD5=LL5,LDT1=LLT1) |
---|
296 | ENDIF |
---|
297 | IF(YCVGQ_NL%LGP) THEN |
---|
298 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
299 | YCVGQ=>YGFLC(IGFLPTR) |
---|
300 | CALL DEFINE_GFL_COMP(YDGFLC=YCVGQ ,CDNAME=YFACVGQ%CLNAME,KGRIB=NGRB149,& |
---|
301 | & LDGP=.TRUE.,KREQIN=YCVGQ_NL%NREQIN, & |
---|
302 | & PREFVALI=YCVGQ_NL%REFVALI, LDREQOUT=YCVGQ_NL%LREQOUT,LDERS=LLDERS,& |
---|
303 | & LD5=LL5,LDT1=LLT1) |
---|
304 | ENDIF |
---|
305 | IF(YSDSAT_NL%LGP) THEN |
---|
306 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
307 | YSDSAT=>YGFLC(IGFLPTR) |
---|
308 | CALL DEFINE_GFL_COMP(YDGFLC=YSDSAT,CDNAME=YFASDSAT%CLNAME,KGRIB=NGRB149,& |
---|
309 | & LDGP=.TRUE.,KREQIN=YSDSAT_NL%NREQIN, & |
---|
310 | & PREFVALI=YSDSAT_NL%REFVALI, LDREQOUT=YSDSAT_NL%LREQOUT,LDERS=LLDERS,& |
---|
311 | & LD5=LL5,LDT1=LLT1) |
---|
312 | ENDIF |
---|
313 | IF(YCVV_NL%LGP) THEN |
---|
314 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
315 | YCVV=>YGFLC(IGFLPTR) |
---|
316 | CALL DEFINE_GFL_COMP(YDGFLC=YCVV,CDNAME=YFACVV%CLNAME,KGRIB=NGRB149,& |
---|
317 | & LDGP=.TRUE.,KREQIN=YCVV_NL%NREQIN, & |
---|
318 | & PREFVALI=YCVV_NL%REFVALI, LDREQOUT=YCVV_NL%LREQOUT,LDERS=LLDERS,& |
---|
319 | & LD5=LL5,LDT1=LLT1) |
---|
320 | ENDIF |
---|
321 | IF(YQVA_NL%LGP) THEN |
---|
322 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
323 | YQVA=>YGFLC(IGFLPTR) |
---|
324 | CALL DEFINE_GFL_COMP(YDGFLC=YQVA,CDNAME=YFAQVA%CLNAME,& |
---|
325 | & KGRIB=YQVA_NL%IGRBCODE,LDGP=.TRUE., & |
---|
326 | & KREQIN=YQVA_NL%NREQIN, PREFVALI=YQVA_NL%REFVALI,& |
---|
327 | & LDREQOUT=YQVA_NL%LREQOUT,LDERS=.FALSE.,LD5=LL5, & |
---|
328 | & LDT1=LLT1,LDGPINGP=.TRUE.) |
---|
329 | ENDIF |
---|
330 | IF(YUOM_NL%LGP)THEN |
---|
331 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
332 | YUOM=>YGFLC(IGFLPTR) |
---|
333 | CALL DEFINE_GFL_COMP(YDGFLC=YUOM,CDNAME=YFAUOM%CLNAME,& |
---|
334 | & KGRIB=YUOM_NL%IGRBCODE,LDGP=.TRUE., & |
---|
335 | & KREQIN=YUOM_NL%NREQIN,& |
---|
336 | & PREFVALI=YUOM_NL%REFVALI,LDREQOUT=YUOM_NL%LREQOUT,& |
---|
337 | & LDERS=.FALSE.,LD5=LL5,LDT1=LLT1,LDGPINGP=.TRUE.) |
---|
338 | ENDIF |
---|
339 | IF(YUAL_NL%LGP)THEN |
---|
340 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
341 | YUAL=>YGFLC(IGFLPTR) |
---|
342 | CALL DEFINE_GFL_COMP(YDGFLC=YUAL,CDNAME=YFAUAL%CLNAME,& |
---|
343 | & KGRIB=YUAL_NL%IGRBCODE,LDGP=.TRUE., & |
---|
344 | & KREQIN=YUAL_NL%NREQIN,& |
---|
345 | & PREFVALI=YUAL_NL%REFVALI,LDREQOUT=YUAL_NL%LREQOUT,& |
---|
346 | & LDERS=.FALSE.,LD5=LL5,LDT1=LLT1,LDGPINGP=.TRUE.) |
---|
347 | ENDIF |
---|
348 | IF(YDOM_NL%LGP)THEN |
---|
349 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
350 | YDOM=>YGFLC(IGFLPTR) |
---|
351 | CALL DEFINE_GFL_COMP(YDGFLC=YDOM,CDNAME=YFADOM%CLNAME,& |
---|
352 | & KGRIB=YDOM_NL%IGRBCODE,LDGP=.TRUE., & |
---|
353 | & KREQIN=YDOM_NL%NREQIN,& |
---|
354 | & PREFVALI=YDOM_NL%REFVALI,LDREQOUT=YDOM_NL%LREQOUT,& |
---|
355 | & LDERS=.FALSE.,LD5=LL5,LDT1=LLT1,LDGPINGP=.TRUE.) |
---|
356 | ENDIF |
---|
357 | IF(YDAL_NL%LGP)THEN |
---|
358 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
359 | YDAL=>YGFLC(IGFLPTR) |
---|
360 | CALL DEFINE_GFL_COMP(YDGFLC=YDAL,CDNAME=YFADAL%CLNAME,& |
---|
361 | & KGRIB=YDAL_NL%IGRBCODE,LDGP=.TRUE., & |
---|
362 | & KREQIN=YDAL_NL%NREQIN,& |
---|
363 | & PREFVALI=YDAL_NL%REFVALI,LDREQOUT=YDAL_NL%LREQOUT,& |
---|
364 | & LDERS=.FALSE.,LD5=LL5,LDT1=LLT1,LDGPINGP=.TRUE.) |
---|
365 | ENDIF |
---|
366 | IF(YUEN_NL%LGP)THEN |
---|
367 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
368 | YUEN=>YGFLC(IGFLPTR) |
---|
369 | CALL DEFINE_GFL_COMP(YDGFLC=YUEN,CDNAME=YFAUEN%CLNAME,& |
---|
370 | & KGRIB=YUEN_NL%IGRBCODE,LDGP=.TRUE., & |
---|
371 | & KREQIN=YUEN_NL%NREQIN,& |
---|
372 | & PREFVALI=YUEN_NL%REFVALI,LDREQOUT=YUEN_NL%LREQOUT,& |
---|
373 | & LDERS=.FALSE.,LD5=LL5,LDT1=LLT1,LDGPINGP=.TRUE.) |
---|
374 | ENDIF |
---|
375 | IF(YUNEBH_NL%LGP)THEN |
---|
376 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
377 | YUNEBH=>YGFLC(IGFLPTR) |
---|
378 | CALL DEFINE_GFL_COMP(YDGFLC=YUNEBH,CDNAME=YFAUNEBH%CLNAME,& |
---|
379 | & KGRIB=YUNEBH_NL%IGRBCODE,LDGP=.TRUE., & |
---|
380 | & KREQIN=YUNEBH_NL%NREQIN,& |
---|
381 | & PREFVALI=YUNEBH_NL%REFVALI,LDREQOUT=YUNEBH_NL%LREQOUT,& |
---|
382 | & LDERS=.FALSE.,LD5=LL5,LDT1=LLT1,LDGPINGP=.TRUE.) |
---|
383 | ENDIF |
---|
384 | |
---|
385 | ! Extra GFL variables |
---|
386 | ! - do a preliminary test to see if, when there are extra-GFL and when at |
---|
387 | ! least one of them is a grid-point one, there is no extra-GFL which is |
---|
388 | ! a spectral one (not allowed): |
---|
389 | ! in this case all the extra-GFL must be grid-point ones. |
---|
390 | IF (NGFL_EXT > 1) THEN |
---|
391 | IF(ANY(YEXT_NL(1:NGFL_EXT)%LGP)) THEN |
---|
392 | DO JGFL=1,NGFL_EXT |
---|
393 | IF (YEXT_NL(JGFL)%LSP) CALL ABOR1& |
---|
394 | & ('SUGFL:EXTRA GFL FIELDS HAVE TO BE ALL EITHER SPECTRAL OR GRIDPOINT') |
---|
395 | ENDDO |
---|
396 | ENDIF |
---|
397 | ENDIF |
---|
398 | ! - define the GFL components (all the extra-GFL are grid-point ones). |
---|
399 | DO JGFL=1,NGFL_EXT |
---|
400 | IF(YEXT_NL(JGFL)%LGP) THEN |
---|
401 | IF(JGFL == 1) THEN |
---|
402 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
403 | YEXT=>YGFLC(IGFLPTR:IGFLPTR+NGFL_EXT-1) |
---|
404 | ENDIF |
---|
405 | CALL DEFINE_GFL_COMP(YDGFLC=YEXT(JGFL),CDNAME=YEXT_NL(JGFL)%CNAME ,& |
---|
406 | & KGRIB=YEXT_NL(JGFL)%IGRBCODE,LDGP=.TRUE.,KREQIN=YEXT_NL(JGFL)%NREQIN,& |
---|
407 | & PREFVALI=YEXT_NL(JGFL)%REFVALI, LDREQOUT=YEXT_NL(JGFL)%LREQOUT ,& |
---|
408 | & LDERS=LLDERS, LD5=LL5, LDT1=LLT1,LDGPINGP=YEXT_NL(JGFL)%LGPINGP, & |
---|
409 | & LDTRAJIO=YEXT_NL(JGFL)%LTRAJIO) |
---|
410 | ENDIF |
---|
411 | ENDDO |
---|
412 | |
---|
413 | ! Forcings |
---|
414 | DO JGFL=1,NGFL_FORC |
---|
415 | IF(YFORC_NL(JGFL)%LGP) THEN |
---|
416 | IF(JGFL == 1) THEN |
---|
417 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
418 | YFORC=>YGFLC(IGFLPTR:IGFLPTR+NGFL_FORC-1) |
---|
419 | ENDIF |
---|
420 | CALL DEFINE_GFL_COMP(YDGFLC=YFORC(JGFL),CDNAME=YFORC_NL(JGFL)%CNAME ,& |
---|
421 | & KGRIB=YFORC_NL(JGFL)%IGRBCODE,LDGP=.TRUE.,KREQIN=YFORC_NL(JGFL)%NREQIN,& |
---|
422 | & PREFVALI=YFORC_NL(JGFL)%REFVALI, LDREQOUT=YFORC_NL(JGFL)%LREQOUT ,& |
---|
423 | & LDERS=LLDERS, LD5=LL5, LDT1=LLT1,LDGPINGP=YFORC_NL(JGFL)%LGPINGP, & |
---|
424 | & LDTRAJIO=YFORC_NL(JGFL)%LTRAJIO) |
---|
425 | ENDIF |
---|
426 | ENDDO |
---|
427 | |
---|
428 | ! Easy diagnostics |
---|
429 | DO JGFL=1,NGFL_EZDIAG |
---|
430 | IF(YEZDIAG_NL(JGFL)%LGP) THEN |
---|
431 | IF(JGFL == 1) THEN |
---|
432 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
433 | YEZDIAG=>YGFLC(IGFLPTR:IGFLPTR+NGFL_EZDIAG-1) |
---|
434 | ENDIF |
---|
435 | CALL DEFINE_GFL_COMP(YDGFLC=YEZDIAG(JGFL),CDNAME=YEZDIAG_NL(JGFL)%CNAME ,& |
---|
436 | & KGRIB=YEZDIAG_NL(JGFL)%IGRBCODE,LDGP=.TRUE.,KREQIN=YEZDIAG_NL(JGFL)%NREQIN,& |
---|
437 | & PREFVALI=YEZDIAG_NL(JGFL)%REFVALI, LDREQOUT=YEZDIAG_NL(JGFL)%LREQOUT ,& |
---|
438 | & LDERS=LLDERS, LD5=LL5, LDT1=LLT1,LDGPINGP=YEZDIAG_NL(JGFL)%LGPINGP, & |
---|
439 | & LDTRAJIO=YEZDIAG_NL(JGFL)%LTRAJIO) |
---|
440 | ENDIF |
---|
441 | ENDDO |
---|
442 | |
---|
443 | !-- Greenhouse Gases |
---|
444 | DO JGFL=1,NGHG |
---|
445 | IF(YGHG_NL(JGFL)%LGP) THEN |
---|
446 | IF(JGFL == 1) THEN |
---|
447 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
448 | YGHG=>YGFLC(IGFLPTR:IGFLPTR+NGHG-1) |
---|
449 | ENDIF |
---|
450 | CALL DEFINE_GFL_COMP(YDGFLC=YGHG(JGFL),CDNAME=YGHG_NL(JGFL)%CNAME ,& |
---|
451 | & KGRIB=YGHG_NL(JGFL)%IGRBCODE, LDGP=.TRUE., KREQIN=YGHG_NL(JGFL)%NREQIN ,& |
---|
452 | & PREFVALI=YGHG_NL(JGFL)%REFVALI, LDREQOUT=YGHG_NL(JGFL)%LREQOUT ,& |
---|
453 | & LDERS=LLDERS, LD5=LL5, LDT1=LLT1,LDGPINGP=YGHG_NL(JGFL)%LGPINGP, & |
---|
454 | & LDTRAJIO=YGHG_NL(JGFL)%LTRAJIO) |
---|
455 | ENDIF |
---|
456 | ENDDO |
---|
457 | |
---|
458 | !-- Diagnostic tracers |
---|
459 | DO JGFL=1,NTRAC |
---|
460 | IF(YTRAC_NL(JGFL)%LGP) THEN |
---|
461 | IF(JGFL == 1) THEN |
---|
462 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
463 | YTRAC=>YGFLC(IGFLPTR:IGFLPTR+NTRAC-1) |
---|
464 | ENDIF |
---|
465 | CALL DEFINE_GFL_COMP(YDGFLC=YTRAC(JGFL),CDNAME=YTRAC_NL(JGFL)%CNAME ,& |
---|
466 | & KGRIB=YTRAC_NL(JGFL)%IGRBCODE, LDGP=.TRUE., KREQIN=YTRAC_NL(JGFL)%NREQIN ,& |
---|
467 | & PREFVALI=YTRAC_NL(JGFL)%REFVALI, LDREQOUT=YTRAC_NL(JGFL)%LREQOUT ,& |
---|
468 | & LDERS=LLDERS, LD5=LL5, LDT1=LLT1,LDGPINGP=YTRAC_NL(JGFL)%LGPINGP, & |
---|
469 | & LDTRAJIO=YTRAC_NL(JGFL)%LTRAJIO) |
---|
470 | ENDIF |
---|
471 | ENDDO |
---|
472 | |
---|
473 | !-- Reactive Gases |
---|
474 | DO JGFL=1,NGRG |
---|
475 | IF(YGRG_NL(JGFL)%LGP) THEN |
---|
476 | IF(JGFL == 1) THEN |
---|
477 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
478 | YGRG=>YGFLC(IGFLPTR:IGFLPTR+NGRG-1) |
---|
479 | ENDIF |
---|
480 | CALL DEFINE_GFL_COMP(YDGFLC=YGRG(JGFL),CDNAME=YGRG_NL(JGFL)%CNAME ,& |
---|
481 | & KGRIB=YGRG_NL(JGFL)%IGRBCODE, LDGP=.TRUE., KREQIN=YGRG_NL(JGFL)%NREQIN ,& |
---|
482 | & PREFVALI=YGRG_NL(JGFL)%REFVALI, LDREQOUT=YGRG_NL(JGFL)%LREQOUT ,& |
---|
483 | & LDERS=LLDERS, LD5=LL5, LDT1=LLT1,LDGPINGP=YGRG_NL(JGFL)%LGPINGP, & |
---|
484 | & LDTRAJIO=YGRG_NL(JGFL)%LTRAJIO) |
---|
485 | ENDIF |
---|
486 | ENDDO |
---|
487 | |
---|
488 | !-- Aerosols |
---|
489 | DO JGFL=1,NAERO |
---|
490 | IF(YAERO_NL(JGFL)%LGP) THEN |
---|
491 | IF(JGFL == 1) THEN |
---|
492 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
493 | YAERO=>YGFLC(IGFLPTR:IGFLPTR+NAERO-1) |
---|
494 | ENDIF |
---|
495 | CALL DEFINE_GFL_COMP(YDGFLC=YAERO(JGFL),CDNAME=YAERO_NL(JGFL)%CNAME ,& |
---|
496 | & KGRIB=YAERO_NL(JGFL)%IGRBCODE, LDGP=.TRUE., KREQIN=YAERO_NL(JGFL)%NREQIN ,& |
---|
497 | & PREFVALI=YAERO_NL(JGFL)%REFVALI, LDREQOUT=YAERO_NL(JGFL)%LREQOUT ,& |
---|
498 | & LDERS=LLDERS, LD5=LL5, LDT1=LLT1,LDGPINGP=YAERO_NL(JGFL)%LGPINGP, & |
---|
499 | & LDTRAJIO=YAERO_NL(JGFL)%LTRAJIO) |
---|
500 | ENDIF |
---|
501 | ENDDO |
---|
502 | |
---|
503 | ! 1.3 Spectral GFL. |
---|
504 | |
---|
505 | IF(YQ_NL%LSP) THEN |
---|
506 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
507 | YQ=>YGFLC(IGFLPTR) |
---|
508 | CALL DEFINE_GFL_COMP(YDGFLC=YQ ,CDNAME=YFAQ%CLNAME,KGRIB=NGRBQ,& |
---|
509 | & LDGP=.FALSE.,KREQIN=YQ_NL%NREQIN, & |
---|
510 | & PREFVALI=YQ_NL%REFVALI, LDREQOUT=YQ_NL%LREQOUT, LDERS=LADER,LD5=LL5, & |
---|
511 | & LDT1=LLT1,LDTHERMACT=.TRUE.,PR=RV,PRCP=RCPV) |
---|
512 | ENDIF |
---|
513 | IF(YL_NL%LSP) THEN |
---|
514 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
515 | YL=>YGFLC(IGFLPTR) |
---|
516 | CALL DEFINE_GFL_COMP(YDGFLC=YL ,CDNAME=YFAL%CLNAME,KGRIB=NGRBCLWC,& |
---|
517 | & LDGP=.FALSE.,KREQIN=YL_NL%NREQIN, & |
---|
518 | & PREFVALI=YL_NL%REFVALI, LDREQOUT=YL_NL%LREQOUT, LDERS=LADER,LD5=LL5, & |
---|
519 | & LDT1=LLT1,LDTHERMACT=.TRUE.,PR=0.0_JPRB,PRCP=RCW) |
---|
520 | ENDIF |
---|
521 | IF(YI_NL%LSP)THEN |
---|
522 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
523 | YI=>YGFLC(IGFLPTR) |
---|
524 | CALL DEFINE_GFL_COMP(YDGFLC=YI ,CDNAME=YFAI%CLNAME,KGRIB=NGRBCIWC,& |
---|
525 | & LDGP=.FALSE.,KREQIN=YI_NL%NREQIN, & |
---|
526 | & PREFVALI=YI_NL%REFVALI, LDREQOUT=YI_NL%LREQOUT, LDERS=LADER,LD5=LL5, & |
---|
527 | & LDT1=LLT1,LDTHERMACT=.TRUE.,PR=0.0_JPRB,PRCP=RCS) |
---|
528 | ENDIF |
---|
529 | IF(YS_NL%LSP)THEN |
---|
530 | CALL ABOR1('SUGFL: spectral snow not coded ') |
---|
531 | ENDIF |
---|
532 | IF(YR_NL%LSP)THEN |
---|
533 | CALL ABOR1('SUGFL: spectral rain not coded ') |
---|
534 | ENDIF |
---|
535 | IF(YG_NL%LSP)THEN |
---|
536 | CALL ABOR1('SUGFL: spectral graupels not coded ') |
---|
537 | ENDIF |
---|
538 | IF(YTKE_NL%LSP)THEN |
---|
539 | CALL ABOR1('SUGFL: spectral TKE not coded ') |
---|
540 | ENDIF |
---|
541 | IF(YA_NL%LSP)THEN |
---|
542 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
543 | YA=>YGFLC(IGFLPTR) |
---|
544 | CALL DEFINE_GFL_COMP(YDGFLC=YA ,CDNAME=YFACLF%CLNAME,KGRIB=NGRBCC,& |
---|
545 | & LDGP=.FALSE.,KREQIN=YA_NL%NREQIN, & |
---|
546 | & PREFVALI=YA_NL%REFVALI, LDREQOUT=YA_NL%LREQOUT, LDERS=LADER,LD5=LL5,& |
---|
547 | & LDT1=LLT1) |
---|
548 | ENDIF |
---|
549 | IF(YO3_NL%LSP)THEN |
---|
550 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
551 | YO3=>YGFLC(IGFLPTR) |
---|
552 | CALL DEFINE_GFL_COMP(YDGFLC=YO3,CDNAME=YFAO3%CLNAME,KGRIB=NGRBO3,& |
---|
553 | & LDGP=.FALSE.,KREQIN=YO3_NL%NREQIN, & |
---|
554 | & PREFVALI=YO3_NL%REFVALI, LDREQOUT=YO3_NL%LREQOUT, LDERS=LADER,LD5=LL5,& |
---|
555 | & LDT1=LLT1) |
---|
556 | ENDIF |
---|
557 | IF(YCPF_NL%LSP)THEN |
---|
558 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
559 | YCPF=>YGFLC(IGFLPTR) |
---|
560 | CALL DEFINE_GFL_COMP(YDGFLC=YCPF,CDNAME=YFACPF%CLNAME,& |
---|
561 | & KGRIB=YCPF_NL%IGRBCODE, LDGP=.FALSE., & |
---|
562 | & KREQIN=YCPF_NL%NREQIN,PREFVALI=YCPF_NL%REFVALI,LDREQOUT=YCPF_NL%LREQOUT,& |
---|
563 | & LDERS=.FALSE.,LD5=.FALSE., & |
---|
564 | & LDT1=LLT1) |
---|
565 | ENDIF |
---|
566 | IF(YSPF_NL%LSP)THEN |
---|
567 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
568 | YSPF=>YGFLC(IGFLPTR) |
---|
569 | CALL DEFINE_GFL_COMP(YDGFLC=YSPF,CDNAME=YFASPF%CLNAME,& |
---|
570 | & KGRIB=YSPF_NL%IGRBCODE,LDGP= .FALSE., & |
---|
571 | & KREQIN=YSPF_NL%NREQIN,PREFVALI=YSPF_NL%REFVALI,LDREQOUT=YSPF_NL%LREQOUT,& |
---|
572 | & LDERS=.FALSE.,LD5=.FALSE., & |
---|
573 | & LDT1=LLT1) |
---|
574 | ENDIF |
---|
575 | IF(YSRC_NL%LSP)THEN |
---|
576 | CALL ABOR1('SUGFL: spectral SRC not coded ') |
---|
577 | ENDIF |
---|
578 | IF(YCVGQ_NL%LSP)THEN |
---|
579 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
580 | YCVGQ=>YGFLC(IGFLPTR) |
---|
581 | CALL DEFINE_GFL_COMP(YDGFLC=YCVGQ,CDNAME=YFACVGQ%CLNAME,& |
---|
582 | & KGRIB=YCVGQ_NL%IGRBCODE,LDGP= .FALSE., & |
---|
583 | & KREQIN=YCVGQ_NL%NREQIN,PREFVALI=YCVGQ_NL%REFVALI,& |
---|
584 | & LDREQOUT=YCVGQ_NL%LREQOUT,LDERS= .TRUE.,LD5=.FALSE., & |
---|
585 | & LDT1=LLT1) |
---|
586 | ENDIF |
---|
587 | IF(YSDSAT_NL%LSP)THEN |
---|
588 | CALL ABOR1('SUGFL: spectral SDSAT not coded ') |
---|
589 | ENDIF |
---|
590 | IF(YCVV_NL%LSP)THEN |
---|
591 | CALL ABOR1('SUGFL: spectral CVV not coded ') |
---|
592 | ENDIF |
---|
593 | IF(YQVA_NL%LSP) THEN |
---|
594 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
595 | YQVA=>YGFLC(IGFLPTR) |
---|
596 | CALL DEFINE_GFL_COMP(YDGFLC=YQVA,CDNAME=YFAQVA%CLNAME,& |
---|
597 | & KGRIB=YQVA_NL%IGRBCODE, LDGP=.FALSE., & |
---|
598 | & KREQIN=YQVA_NL%NREQIN,PREFVALI=YQVA_NL%REFVALI,LDREQOUT=YQVA_NL%LREQOUT,& |
---|
599 | & LDERS=.FALSE.,LD5=.FALSE., & |
---|
600 | & LDT1=LLT1) |
---|
601 | ENDIF |
---|
602 | IF(YUOM_NL%LSP)THEN |
---|
603 | CALL ABOR1('SUGFL: spectral UOM not coded ') |
---|
604 | ! ---------------------------------------------------------------------- |
---|
605 | ! what is below must be checked when UOM will be coded! |
---|
606 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
607 | YUOM=>YGFLC(IGFLPTR) |
---|
608 | CALL DEFINE_GFL_COMP(YDGFLC=YUOM,CDNAME=YFAUOM%CLNAME,& |
---|
609 | & KGRIB=YUOM_NL%IGRBCODE,LDGP=.FALSE., & |
---|
610 | & KREQIN=YUOM_NL%NREQIN,& |
---|
611 | & PREFVALI=YUOM_NL%REFVALI,LDREQOUT=YUOM_NL%LREQOUT,& |
---|
612 | & LDERS=.FALSE.,LD5=.FALSE.,LDT1=LLT1) |
---|
613 | ! ---------------------------------------------------------------------- |
---|
614 | ENDIF |
---|
615 | IF(YUAL_NL%LSP)THEN |
---|
616 | CALL ABOR1('SUGFL: spectral UAL not coded ') |
---|
617 | ! ---------------------------------------------------------------------- |
---|
618 | ! what is below must be checked when UAL will be coded! |
---|
619 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
620 | YUAL=>YGFLC(IGFLPTR) |
---|
621 | CALL DEFINE_GFL_COMP(YDGFLC=YUAL,CDNAME=YFAUAL%CLNAME,& |
---|
622 | & KGRIB=YUAL_NL%IGRBCODE,LDGP=.FALSE., & |
---|
623 | & KREQIN=YUAL_NL%NREQIN,& |
---|
624 | & PREFVALI=YUAL_NL%REFVALI,LDREQOUT=YUAL_NL%LREQOUT,& |
---|
625 | & LDERS=.FALSE.,LD5=.FALSE.,LDT1=LLT1) |
---|
626 | ! ---------------------------------------------------------------------- |
---|
627 | ENDIF |
---|
628 | IF(YDOM_NL%LSP)THEN |
---|
629 | CALL ABOR1('SUGFL: spectral DOM not coded ') |
---|
630 | ! ---------------------------------------------------------------------- |
---|
631 | ! what is below must be checked when DOM will be coded! |
---|
632 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
633 | YDOM=>YGFLC(IGFLPTR) |
---|
634 | CALL DEFINE_GFL_COMP(YDGFLC=YDOM,CDNAME=YFADOM%CLNAME,& |
---|
635 | & KGRIB=YDOM_NL%IGRBCODE,LDGP=.FALSE., & |
---|
636 | & KREQIN=YDOM_NL%NREQIN,& |
---|
637 | & PREFVALI=YDOM_NL%REFVALI,LDREQOUT=YDOM_NL%LREQOUT,& |
---|
638 | & LDERS=.FALSE.,LD5=.FALSE.,LDT1=LLT1) |
---|
639 | ! ---------------------------------------------------------------------- |
---|
640 | ENDIF |
---|
641 | IF(YDAL_NL%LSP)THEN |
---|
642 | CALL ABOR1('SUGFL: spectral DAL not coded ') |
---|
643 | ! ---------------------------------------------------------------------- |
---|
644 | ! what is below must be checked when DAL will be coded! |
---|
645 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
646 | YDAL=>YGFLC(IGFLPTR) |
---|
647 | CALL DEFINE_GFL_COMP(YDGFLC=YDAL,CDNAME=YFADAL%CLNAME,& |
---|
648 | & KGRIB=YDAL_NL%IGRBCODE,LDGP=.FALSE., & |
---|
649 | & KREQIN=YDAL_NL%NREQIN,& |
---|
650 | & PREFVALI=YDAL_NL%REFVALI,LDREQOUT=YDOM_NL%LREQOUT,& |
---|
651 | & LDERS=.FALSE.,LD5=.FALSE.,LDT1=LLT1) |
---|
652 | ! ---------------------------------------------------------------------- |
---|
653 | ENDIF |
---|
654 | IF(YUEN_NL%LSP)THEN |
---|
655 | CALL ABOR1('SUGFL: spectral UEN not coded ') |
---|
656 | ! ---------------------------------------------------------------------- |
---|
657 | ! what is below must be checked when UEN will be coded! |
---|
658 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
659 | YUEN=>YGFLC(IGFLPTR) |
---|
660 | CALL DEFINE_GFL_COMP(YDGFLC=YUEN,CDNAME=YFAUEN%CLNAME,& |
---|
661 | & KGRIB=YUEN_NL%IGRBCODE,LDGP=.FALSE., & |
---|
662 | & KREQIN=YUEN_NL%NREQIN,& |
---|
663 | & PREFVALI=YUEN_NL%REFVALI,LDREQOUT=YUEN_NL%LREQOUT,& |
---|
664 | & LDERS=.FALSE.,LD5=.FALSE.,LDT1=LLT1) |
---|
665 | ! ---------------------------------------------------------------------- |
---|
666 | ENDIF |
---|
667 | IF(YUNEBH_NL%LSP)THEN |
---|
668 | CALL ABOR1('SUGFL: spectral UNEBH not coded ') |
---|
669 | ! ---------------------------------------------------------------------- |
---|
670 | ! what is below must be checked when UNEBH will be coded! |
---|
671 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
672 | YUNEBH=>YGFLC(IGFLPTR) |
---|
673 | CALL DEFINE_GFL_COMP(YDGFLC=YUNEBH,CDNAME=YFAUNEBH%CLNAME,& |
---|
674 | & KGRIB=YUNEBH_NL%IGRBCODE,LDGP=.FALSE., & |
---|
675 | & KREQIN=YUNEBH_NL%NREQIN,& |
---|
676 | & PREFVALI=YUNEBH_NL%REFVALI,LDREQOUT=YUNEBH_NL%LREQOUT,& |
---|
677 | & LDERS=.FALSE.,LD5=.FALSE.,LDT1=LLT1) |
---|
678 | ! ---------------------------------------------------------------------- |
---|
679 | ENDIF |
---|
680 | |
---|
681 | ! Extra GFL variables (assumed to be all spectral in this piece of code) |
---|
682 | DO JGFL=1,NGFL_EXT |
---|
683 | IF(YEXT_NL(JGFL)%LSP) THEN |
---|
684 | IF(JGFL == 1) THEN |
---|
685 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
686 | YEXT=>YGFLC(IGFLPTR:IGFLPTR+NGFL_EXT-1) |
---|
687 | ENDIF |
---|
688 | CALL DEFINE_GFL_COMP(YDGFLC=YEXT(JGFL),CDNAME=YEXT_NL(JGFL)%CNAME ,& |
---|
689 | & KGRIB=YEXT_NL(JGFL)%IGRBCODE,LDGP=.FALSE.,KREQIN=YEXT_NL(JGFL)%NREQIN,& |
---|
690 | & PREFVALI=YEXT_NL(JGFL)%REFVALI, LDREQOUT=YEXT_NL(JGFL)%LREQOUT ,& |
---|
691 | & LDERS=LADER, LD5=LL5, LDT1=LLT1) |
---|
692 | ENDIF |
---|
693 | ENDDO |
---|
694 | |
---|
695 | ! Forcings |
---|
696 | DO JGFL=1,NGFL_FORC |
---|
697 | IF(YFORC_NL(JGFL)%LSP) THEN |
---|
698 | CALL ABOR1('SUGFL: FORCINGS HAVE TO BE GRIDPOINT FIELDS') |
---|
699 | ! ---------------------------------------------------------------------- |
---|
700 | ! what is below must be checked when spectral forcings will be coded! |
---|
701 | IF(JGFL == 1) THEN |
---|
702 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
703 | YFORC=>YGFLC(IGFLPTR:IGFLPTR+NGFL_FORC-1) |
---|
704 | ENDIF |
---|
705 | CALL DEFINE_GFL_COMP(YDGFLC=YFORC(JGFL),CDNAME=YFORC_NL(JGFL)%CNAME,& |
---|
706 | & KGRIB=YFORC_NL(JGFL)%IGRBCODE,LDGP=.FALSE.,KREQIN=YFORC_NL(JGFL)%NREQIN,& |
---|
707 | & PREFVALI=YFORC_NL(JGFL)%REFVALI, LDREQOUT=YFORC_NL(JGFL)%LREQOUT, & |
---|
708 | & LDERS=.FALSE.,LD5=.FALSE.,LDT1=.FALSE.) |
---|
709 | ! ---------------------------------------------------------------------- |
---|
710 | ENDIF |
---|
711 | ENDDO |
---|
712 | |
---|
713 | ! Easy diagnostics |
---|
714 | DO JGFL=1,NGFL_EZDIAG |
---|
715 | IF(YEZDIAG_NL(JGFL)%LSP) THEN |
---|
716 | CALL ABOR1('SUGFL: EZDIAG HAVE TO BE GRIDPOINT FIELDS') |
---|
717 | ! ---------------------------------------------------------------------- |
---|
718 | ! what is below must be checked when spectral easy diag. will be coded! |
---|
719 | IF(JGFL == 1) THEN |
---|
720 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
721 | YEZDIAG=>YGFLC(IGFLPTR:IGFLPTR+NGFL_EZDIAG-1) |
---|
722 | ENDIF |
---|
723 | CALL DEFINE_GFL_COMP(YDGFLC=YEZDIAG(JGFL),CDNAME=YEZDIAG_NL(JGFL)%CNAME,& |
---|
724 | & KGRIB=YEZDIAG_NL(JGFL)%IGRBCODE,LDGP=.FALSE.,& |
---|
725 | & KREQIN=YEZDIAG_NL(JGFL)%NREQIN,& |
---|
726 | & PREFVALI=YEZDIAG_NL(JGFL)%REFVALI, LDREQOUT=YEZDIAG_NL(JGFL)%LREQOUT, & |
---|
727 | & LDERS=.FALSE.,LD5=.FALSE.,LDT1=.FALSE.) |
---|
728 | ! ---------------------------------------------------------------------- |
---|
729 | ENDIF |
---|
730 | ENDDO |
---|
731 | |
---|
732 | !Greenhouse Gases |
---|
733 | DO JGFL=1,NGHG |
---|
734 | IF(YGHG_NL(JGFL)%LSP) THEN |
---|
735 | IF(JGFL == 1) THEN |
---|
736 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
737 | YGHG=>YGFLC(IGFLPTR:IGFLPTR+NGHG-1) |
---|
738 | ENDIF |
---|
739 | CALL DEFINE_GFL_COMP(YDGFLC=YGHG(JGFL),CDNAME=YGHG_NL(JGFL)%CNAME ,& |
---|
740 | & KGRIB=YGHG_NL(JGFL)%IGRBCODE, LDGP=.FALSE., KREQIN=YGHG_NL(JGFL)%NREQIN ,& |
---|
741 | & PREFVALI=YGHG_NL(JGFL)%REFVALI, LDREQOUT=YGHG_NL(JGFL)%LREQOUT ,& |
---|
742 | & LDERS=LADER, LD5=LL5, LDT1=LLT1) |
---|
743 | ENDIF |
---|
744 | ENDDO |
---|
745 | |
---|
746 | !Diagnostic tracers |
---|
747 | DO JGFL=1,NTRAC |
---|
748 | IF(YTRAC_NL(JGFL)%LSP) THEN |
---|
749 | IF(JGFL == 1) THEN |
---|
750 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
751 | YTRAC=>YGFLC(IGFLPTR:IGFLPTR+NTRAC-1) |
---|
752 | ENDIF |
---|
753 | CALL DEFINE_GFL_COMP(YDGFLC=YTRAC(JGFL),CDNAME=YTRAC_NL(JGFL)%CNAME ,& |
---|
754 | & KGRIB=YTRAC_NL(JGFL)%IGRBCODE, LDGP=.FALSE., KREQIN=YTRAC_NL(JGFL)%NREQIN ,& |
---|
755 | & PREFVALI=YTRAC_NL(JGFL)%REFVALI, LDREQOUT=YTRAC_NL(JGFL)%LREQOUT ,& |
---|
756 | & LDERS=LADER, LD5=LL5, LDT1=LLT1) |
---|
757 | ENDIF |
---|
758 | ENDDO |
---|
759 | |
---|
760 | !Reactive Gases |
---|
761 | DO JGFL=1,NGRG |
---|
762 | IF(YGRG_NL(JGFL)%LSP) THEN |
---|
763 | IF(JGFL == 1) THEN |
---|
764 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
765 | YGRG=>YGFLC(IGFLPTR:IGFLPTR+NGRG-1) |
---|
766 | ENDIF |
---|
767 | CALL DEFINE_GFL_COMP(YDGFLC=YGRG(JGFL),CDNAME=YGRG_NL(JGFL)%CNAME ,& |
---|
768 | & KGRIB=YGRG_NL(JGFL)%IGRBCODE, LDGP=.FALSE., KREQIN=YGRG_NL(JGFL)%NREQIN ,& |
---|
769 | & PREFVALI=YGRG_NL(JGFL)%REFVALI, LDREQOUT=YGRG_NL(JGFL)%LREQOUT ,& |
---|
770 | & LDERS=LADER, LD5=LL5, LDT1=LLT1) |
---|
771 | ENDIF |
---|
772 | ENDDO |
---|
773 | |
---|
774 | !Aerosols |
---|
775 | DO JGFL=1,NAERO |
---|
776 | IF(YAERO_NL(JGFL)%LSP) THEN |
---|
777 | IF(JGFL == 1) THEN |
---|
778 | IGFLPTR=YGFL%NUMFLDS+1 |
---|
779 | YAERO=>YGFLC(IGFLPTR:IGFLPTR+NAERO-1) |
---|
780 | ENDIF |
---|
781 | CALL DEFINE_GFL_COMP(YDGFLC=YAERO(JGFL),CDNAME=YAERO_NL(JGFL)%CNAME ,& |
---|
782 | & KGRIB=YAERO_NL(JGFL)%IGRBCODE, LDGP=.FALSE., KREQIN=YAERO_NL(JGFL)%NREQIN ,& |
---|
783 | & PREFVALI=YAERO_NL(JGFL)%REFVALI, LDREQOUT=YAERO_NL(JGFL)%LREQOUT ,& |
---|
784 | & LDERS=LADER, LD5=LL5, LDT1=LLT1) |
---|
785 | ENDIF |
---|
786 | ENDDO |
---|
787 | |
---|
788 | ELSE |
---|
789 | |
---|
790 | ! 2. CASE LFPART2=T |
---|
791 | |
---|
792 | WRITE(NULOUT,*) ' SUGFL, LFPART2, YGFL%NUMFLDS=',YGFL%NUMFLDS |
---|
793 | ENDIF |
---|
794 | |
---|
795 | !------------------------------------------------------------------------- |
---|
796 | |
---|
797 | IF (LHOOK) CALL DR_HOOK('SUGFL',1,ZHOOK_HANDLE) |
---|
798 | END SUBROUTINE SUGFL |
---|