NCEPLIBS-w3emc 2.12.0
Loading...
Searching...
No Matches
w3fi70.f
Go to the documentation of this file.
1C> @file
2C> @brief Computes scaling constants used by grdprt().
3C> @author Ralph Jones @date 1991-10-26
4C
5C> Computes the four scaling constants used by grdprt, w3fp03,
6C> or w3fp05 from the 28 byte (pds) product definition section of
7C> grib edition one.
8C>
9C> Program history log:
10C> - Ralph Jones 1991-10-26
11C> - Ralph Jones 1993-03-29 Add save statement
12C> - Ralph Jones 1993-08-08 Add 156 (cin), 158 (tke) to tables
13C> - Ralph Jones 1993-10-16 Changes for o.n. 388 ver. oct. 8,1993
14C>
15C> @param[in] PDS 28 byte (pds) grib product definition section.
16C> @param[out] CNST 4 constant's used by grdprt(), w3fp05(), or w3fp03().
17C> @param[out] IER 0 = normal return | 1 = .
18C>
19C> @author Ralph Jones @date 1991-10-26
20 SUBROUTINE w3fi70(PDS,CNST,IER)
21C
22C SET DEFAULT VALUES FOR NMC FIELDS GRID PRINTING
23C
24 REAL CNST(4)
25C
26 INTEGER ID(25)
27 INTEGER Q
28C
29 CHARACTER * 1 PDS(28)
30C
31 SAVE
32C
33C UNPACK 28 BYTE (PDS) INTO 25 INTEGER WORDS
34C
35 CALL w3fi69(pds,id)
36C
37 ier = 0
38C
39C INDICATOR OF PARAMETER AND UNITS
40C
41 q = id(8)
42C
43C INDICATOR OF LEVEL OR LAYERS
44C
45 itypes = id(9)
46 i9 = id(9)
47C
48C HEIGHTS, PRESSURE, ETC. OF THE LEVEL OR LAYER
49C
50 IF ((i9.GE.1.AND.i9.LE.100).OR.i9.EQ.102.OR.
51 & i9.EQ.103.OR.i9.EQ.105.OR.i9.EQ.107.OR.
52 & i9.EQ.109.OR.i9.EQ.111.OR.i9.EQ.113.OR.
53 & i9.EQ.125.OR.i9.EQ.160.OR.i9.EQ.200.OR.
54 & i9.EQ.201) THEN
55 ilvl = id(11)
56 ELSE
57 ilvl = id(10)
58 END IF
59
60 IF (q.EQ.1.OR.q.EQ.2.OR.q.EQ.26) THEN
61C
62C*** PRESSURE, PRESSURE REDUCED TO MSL, PRESSURE ANOMALY (Pa)
63C
64 cnst(1) = 0.0
65 cnst(2) = 0.01
66 cnst(3) = 4.0
67 cnst(4) = 0.0
68C*** IF SFC, TROPOPAUSE PRESSURE, SIGMA ..
69 IF (itypes.EQ.1.OR.itypes.EQ.6.OR.itypes.EQ.7)cnst(3)=25.0
70 IF (itypes.EQ.107) cnst(3) = 25.0
71C
72 ELSE IF (q.EQ.3) THEN
73C
74C*** PRESSURE TENDENCY (Pa/s)
75C
76 cnst(1) = 0.0
77 cnst(2) = 1.0
78 cnst(3) = 4.0
79 cnst(4) = 0.0
80C
81 ELSE IF (q.EQ.6) THEN
82C
83C*** GEOPOTENTIAL (m**2/s**2)
84C
85 cnst(1) = 0.0
86 cnst(2) = 1.0
87 cnst(3) = 4.0
88 cnst(4) = 0.0
89C
90 ELSE IF (q.EQ.7.OR.q.EQ.8.OR.q.EQ.27.OR.q.EQ.222) THEN
91C
92C*** GEOPOTENTIAL, GEOPOTENTIAL HEIGHT, ANOMALY
93C*** 5-WAVE GEOPOTENTIAL HEIGHT ............
94C
95 cnst(3) = 60.
96 IF (ilvl.LT.500) cnst(3) = 120.
97C*** IF SFC OR TROPOPAUSE PRESSURE ..
98 IF ((itypes.EQ.1) .OR. (itypes.EQ.7)) cnst(3) = 500.0
99 IF (itypes.EQ.107) cnst(3) = 500.0
100
101 cnst(1) = 0.0
102 cnst(2) = 1.0
103 cnst(4) = 0.0
104 IF (cnst(3) .EQ. 500.) cnst(4) = 2.0
105C
106 ELSE IF (q.EQ.11.OR.q.EQ.12.OR.q.EQ.13.OR.q.EQ.14.OR.
107 & q.EQ.15.OR.q.EQ.16.OR.q.EQ.17.OR.q.EQ.18.OR.
108 & q.EQ.25.OR.q.EQ.85) THEN
109C
110
111C*** TEMPERATURES (deg. K)
112C*** VIRTUAL TEMPERATURE (deg. K)
113C*** POTENTIAL TEMPERATURE (deg. K)
114C*** PSEUDO-ADIABATIC POTENTIAL TEMPERATURE (deg. K)
115C*** MAXIMUN TEMPERATURE (deg. K)
116C*** MINUMUN TEMPERATURE (deg. K)
117C*** DEW POINT TEMPERATURE (deg. K)
118C*** DEW POINT DEPRESSION (OR DEFICIT) (deg. K)
119C
120C*** TEMP (DEG K) CONVERT TO DEG C, EXCEPT POTENTIAL TEMPERATURE
121C
122C CNST(1) = -273.15
123 cnst(1) = 0.0
124 cnst(2) = 1.0
125 cnst(3) = 5.0
126 cnst(4) = 0.0
127 IF (q.EQ.13) cnst(1) = 0.0
128C
129 ELSE IF (q.EQ.19) THEN
130C
131C*** LAPSE RATE, deg. K/m ...............
132C
133 cnst(1) = 0.0
134 cnst(2) = 1.0
135 cnst(3) = 4.0
136 cnst(4) = 0.0
137C
138 ELSE IF (q.EQ.21.OR.q.EQ.22.OR.q.EQ.23) THEN
139C
140C*** RADAR SPECTRA (1), (2), (3) ...............
141C
142 cnst(1) = 0.0
143 cnst(2) = 1.0
144 cnst(3) = 10.0
145 cnst(4) = 0.0
146C
147 ELSE IF (q.EQ.28.OR.q.EQ.29.OR.q.EQ.30) THEN
148C
149C*** WAVE SPECTRA (1), (2), (3) ...............
150C
151 cnst(1) = 0.0
152 cnst(2) = 1.0
153 cnst(3) = 10.0
154 cnst(4) = 0.0
155C
156 ELSE IF (q.EQ.31) THEN
157C
158C*** WIND DIRECTION (deg. true)
159C
160 cnst(1) = 0.0
161 cnst(2) = 1.0
162 cnst(3) = 10.0
163 cnst(4) = 0.0
164C
165 ELSE IF (q.EQ.32.OR.q.EQ.33.OR.q.EQ.34) THEN
166C
167C*** WIND SPEED, U-COMPONENT OF WIND,
168C*** V-COMPONENT OF WIND m/s -------------------
169C
170 cnst(1) = 0.0
171 cnst(2) = 1.0
172 cnst(3) = 10.0
173 IF (itypes.EQ.1.AND.ilvl.EQ.0) cnst(3) = 3.0
174 IF (itypes.EQ.107) cnst(3) = 3.0
175 cnst(4) = 0.0
176C
177 ELSE IF (q.EQ.35.OR.q.EQ.36) THEN
178C
179C*** STREAM FUNCTION, VELOCITY POTENTIAL (m**2/s)
180C*** STREAM FUNCTION OR VELOCITY POTENTIAL (m**2/s) CONVERTED TO M.
181C*** CONVERT TO METERS. (M*M/SEC * FOG)
182C
183 cnst(1) = 0.
184 cnst(2) = 1.03125e-4 / 9.8
185 cnst(3) = 60.
186 cnst(4) = 0.
187C
188 ELSE IF (q.EQ.37) THEN
189C
190C*** MONTGOMERY STREAM FUNCTION (m**2/s**2)
191C
192 cnst(1) = 0.0
193 cnst(2) = 1.0
194 cnst(3) = 2.0
195 cnst(4) = 0.0
196C
197 ELSE IF (q.EQ.38) THEN
198C
199C*** SIGMA COORD. VERTICAL VELOCITY (/s) TO MICROBARS/SEC
200C
201 cnst(1) = 0.0
202 cnst(2) = 1.0
203 cnst(3) = 2.0
204 cnst(4) = 0.0
205C
206 ELSE IF (q.EQ.39) THEN
207C
208C*** VERTICAL VELOCITY (Pa/s) TO MICROBARS/SEC
209C*** SIGN CHANGED SUCH THAT POSITIVE VALUES INDICATE UPWARD MOTION.
210C
211 cnst(1) = 0.0
212 cnst(2) = -1.e1
213 cnst(3) = 2.0
214 cnst(4) = 0.0
215C
216 ELSE IF (q.EQ.40) THEN
217C
218C*** GEOMETRIC VERTICAL VELOCITY -DZDT- (m/s)
219C
220 cnst(1) = 0.0
221 cnst(2) = 1.0
222 cnst(3) = 10.0
223 cnst(4) = 0.0
224C
225 ELSE IF (q.EQ.41.OR.q.EQ.42.OR.q.EQ.43.OR.q.EQ.44.OR.
226 & q.EQ.45.OR.q.EQ.46) THEN
227C
228C*** ABSOLUTE VORTICITY -ABS-V (/s)
229C*** ABSOLUTE DIVERGENCE -ABS-V (/s)
230C*** RELATIVE VORTICITY -REL-V (/s)
231C*** RELATIVE DIVERGENCE -REL-D (/s)
232C*** VERTICAL U-COMPONENT SHEAR -VUCSH (/s)
233C*** VERTICAL V-COMPONENT SHEAR -VVCSH (/s)
234C
235 cnst(1) = 0.0
236 cnst(2) = 1.0e+6
237 cnst(3) = 40.0
238 cnst(4) = 0.0
239C
240 ELSE IF (q.EQ.47) THEN
241C
242C*** DIRECTION OF CURRENT -DIR-C (deg. true)
243C
244 cnst(1) = 0.0
245 cnst(2) = 1.0
246 cnst(3) = 10.0
247 cnst(4) = 0.0
248C
249 ELSE IF (q.EQ.48.OR.q.EQ.49.OR.q.EQ.50) THEN
250C
251C*** SPEED OF CURRENT (m/s)
252C*** U AND V COMPONENTS OF CURRENT (m/s)
253C
254 cnst(1) = 0.
255 cnst(2) = 1.
256 cnst(3) = 2.
257 cnst(4) = 0.
258C
259 ELSE IF (q.EQ.51.OR.q.EQ.53) THEN
260C
261C*** SPECIFIC HUMIDITY SPF H (kg/kg)
262C*** HUMIDITY MIXING RATIO MIXR (kg/kg)
263C
264 cnst(1) = 0.0
265 cnst(2) = 1.e+3
266 cnst(3) = 2.0
267 cnst(4) = 0.0
268C
269 ELSE IF (q.EQ.52) THEN
270C
271C*** RELATIVE HUMIDITY R H (%)
272C
273 cnst(1) = 0.0
274 cnst(2) = 1.0
275 cnst(3) = 20.0
276 cnst(4) = 0.0
277C
278 ELSE IF (q.EQ.54.OR.q.EQ.57.OR.q.EQ.58) THEN
279C
280C*** PRECIPITABLE WATER (kg/m**2) OR .1 GRAM/CM*CM OR MILLIMETERS/CM*CM
281C*** CHANGE TO CENTI-INCHES/CM*CM
282C*** EVAPERATION
283C*** CLOUD ICE (kg/m**2)
284C
285 cnst(1) = 0.0
286 cnst(2) = 3.937
287 cnst(3) = 10.0
288 cnst(4) = 0.0
289C
290 ELSE IF (q.EQ.55.OR.q.EQ.56) THEN
291C
292C*** VAPOR PRESSURE VAPP, SATURATION DEFICIT SAT D (Pa)
293C
294 cnst(1) = 0.0
295 cnst(2) = 1.0
296 cnst(3) = 10.0
297 cnst(4) = 0.0
298C
299 ELSE IF (q.EQ.59) THEN
300C
301C*** PRECIPITATION RATE (kg/m**2/s)
302C
303 cnst(1) = 0.0
304 cnst(2) = 1.0
305 cnst(3) = 20.0
306 cnst(4) = 0.0
307C
308 ELSE IF (q.EQ.60) THEN
309C
310C*** THUNDERSTORM PROBABILITY (%)
311C
312 cnst(1) = 0.0
313 cnst(2) = 1.0
314 cnst(3) = 20.0
315 cnst(4) = 0.0
316C
317 ELSE IF (q.EQ.61.OR.q.EQ.62.OR.q.EQ.63.OR.q.EQ.64.OR.
318 & q.EQ.65) THEN
319C
320C*** TOTAL PRECIPITATION A PCP (kg/m**2)
321C*** LARGE SCALE PRECIPITATION NCPCP (kg/m**2)
322C*** CONVECTIVE PRECIPITATION ACPCP (kg/m**2)
323C*** SNOWFALL RATE WATER EQUIVALENT SRWEQ (kg/m**2/s)
324C*** WATER EQUIV. OF ACCUM. SNOW DEPTH WEASD (kg/m**2)
325C
326 cnst(1) = 0.0
327 cnst(2) = 1.0
328 cnst(3) = 2.0
329 cnst(4) = 0.0
330
331 ELSE IF (q.EQ.66) THEN
332C
333C*** SNOW DEPTH (METERS) (1 or 0) for snow or no snow
334C
335 cnst(1) = 0.0
336 cnst(2) = 1.0
337 cnst(3) = 1.0
338 cnst(4) = 0.0
339C
340 ELSE IF (q.EQ.67.OR.q.EQ.68.OR.q.EQ.69.OR.q.EQ.70) THEN
341C
342C*** MIXING LAYER DEPTH MIXHT (m)
343C*** TRANSIENT THEMOCLINE DEPTH TTHDP (m)
344C*** MAIN THERMOCLINE DEPTH MTHCD (m)
345C*** MAIN THERMOCLINE ANOMALY MTHCA (m)
346C
347 cnst(1) = 0.0
348 cnst(2) = 39.37
349 cnst(3) = 06.0
350 cnst(4) = 0.0
351C
352 ELSE IF (q.EQ.120.OR.q.EQ.121) THEN
353C
354C*** WAVE COMPONENT OF GEOPOTENTIAL (GEOP M)
355C
356 cnst(1) = 0.0
357 cnst(2) = 1.0
358 cnst(3) = 10.0
359 cnst(4) = 0.0
360C
361 ELSE IF (q.EQ.71.OR.q.EQ.72.OR.q.EQ.73.OR.q.EQ.74.OR.
362 & q.EQ.75) THEN
363C
364C*** TOTAL CLOUD COVER T CDC (%)
365C*** CONVECTIVE CLOUD COVER CDCON (%)
366C*** LOW CLOUD COVER L CDC (%)
367C*** MEDIUM CLOUD COVER M CDC (%)
368C*** HIGH CLOUD COVER H CDC (%)
369C
370 cnst(1) = 0.0
371 cnst(2) = 1.0
372 cnst(3) = 10.0
373 cnst(4) = 0.0
374C
375 ELSE IF (q.EQ.76) THEN
376C
377C*** CLOUD WATER -C-WAT (kg/m**2)
378C
379 cnst(1) = 0.0
380 cnst(2) = 1.0
381 cnst(3) = 10.0
382 cnst(4) = 0.0
383C
384 ELSE IF (q.EQ.78) THEN
385C
386C*** CONVECTIVE SNOW -C-SNO (kg/m**2)
387C
388 cnst(1) = 0.0
389 cnst(2) = 1.0
390 cnst(3) = 10.0
391 cnst(4) = 0.0
392C
393 ELSE IF (q.EQ.79) THEN
394C
395C*** LARGE SCALE SNOW -LSSNO (kg/m**2)
396C
397 cnst(1) = 0.0
398 cnst(2) = 0.1
399 cnst(3) = 500.0
400 cnst(4) = 0.0
401C
402 ELSE IF (q.EQ.80) THEN
403C
404C*** WATER TEMPERAUTER -WTMP- (deg. K)
405C
406 cnst(1) = 0.0
407 cnst(2) = 1.0
408 cnst(3) = 2.0
409 cnst(4) = 0.0
410C
411 ELSE IF (q.EQ.81) THEN
412C
413C*** LAND/SEA (1=LAND; 0=SEA)
414C*** ICE CONCENTRATION (ICE=1; NO ICE=0)
415C
416 cnst(1) = 0.0
417 cnst(2) = 1.0
418 cnst(3) = 1.0
419 cnst(4) = 0.5
420C
421 ELSE IF (q.EQ.82.OR.q.EQ.83.OR.q.EQ.92.OR.q.EQ.97) THEN
422C
423C*** DEVIATION OF SEA LEVEL FROM MEAN (m)
424C*** SUFACE ROUGHNESS (m)
425C*** ICE THICKNESS (m)
426C*** ICE GROWTH (m)
427C
428 cnst(1) = 0.0
429 cnst(2) = 1.0
430 cnst(3) = 2.0
431 cnst(4) = 0.0
432C
433 ELSE IF (q.EQ.84) THEN
434C
435C*** ALBEDO (%)
436C
437 cnst(1) = 0.0
438 cnst(2) = 1.0
439 cnst(3) = 10.0
440 cnst(4) = 0.0
441C
442 ELSE IF (q.EQ.86) THEN
443C
444C*** SOIL MOISTURE CONTENT (kg/m**2) -SOILM
445C
446 cnst(1) = 0.0
447 cnst(2) = 1.0
448 cnst(3) = 10.0
449 cnst(4) = 0.0
450C
451 ELSE IF (q.EQ.87) THEN
452C
453C*** VEGETATION -VEG- (%)
454C
455 cnst(1) = 0.0
456 cnst(2) = 1.0
457 cnst(3) = 10.0
458 cnst(4) = 0.0
459C
460 ELSE IF (q.EQ.88) THEN
461C
462C*** SALINITY -SALTY- (kg/kg)
463C
464 cnst(1) = 0.0
465 cnst(2) = 1.0
466 cnst(3) = 10.0
467 cnst(4) = 0.0
468C
469 ELSE IF (q.EQ.89) THEN
470C
471C*** DENSITY -DEN-- (kg/m**3)
472C
473 cnst(1) = 0.0
474 cnst(2) = 1.0
475 cnst(3) = 10.0
476 cnst(4) = 0.0
477C
478 ELSE IF (q.EQ.90) THEN
479C
480C*** WATER RUNOFF -WAT-R (kg/m**2)
481C
482 cnst(1) = 0.0
483 cnst(2) = 1.0
484 cnst(3) = 10.0
485 cnst(4) = 0.0
486C
487 ELSE IF (q.EQ.93) THEN
488C
489C*** DIRECTION OF ICE DRIFT -DICED (deg. true)
490C
491 cnst(1) = 0.0
492 cnst(2) = 1.0
493 cnst(3) = 10.0
494 cnst(4) = 0.0
495C
496 ELSE IF (q.EQ.94.OR.q.EQ.95.OR.q.EQ.96) THEN
497C
498C*** SPEED OF ICE DRIFT -SICED (m/s)
499C*** U-COMPONENT OF ICE DRIFT -U-ICE (m/s)
500C*** V-COMPONENT OF ICE DRIFT -V-ICE (m/s)
501C
502 cnst(1) = 0.0
503 cnst(2) = 1.0
504 cnst(3) = 2.0
505 cnst(4) = 0.0
506C
507 ELSE IF (q.EQ.98) THEN
508C
509C*** ICE DIVERGENCE -ICE D (/s)
510C
511 cnst(1) = 0.0
512 cnst(2) = 1.0
513 cnst(3) = 10.0
514 cnst(4) = 0.0
515C
516 ELSE IF (q.EQ.99) THEN
517C
518C*** SNO MELT -SNO- M (kg/m**2)
519C
520 cnst(1) = 0.0
521 cnst(2) = 1.0
522 cnst(3) = 10.0
523 cnst(4) = 0.0
524C
525 ELSE IF (q.EQ.100.OR.q.EQ.102.OR.q.EQ.105) THEN
526C
527C*** HEIGHT OF WIND DRIVEN OCEAN WAVES, SEA SWELLS, OR COMBINATION
528C*** (m)
529C
530 cnst(1) = 0.0
531 cnst(2) = 1.0
532 cnst(3) = 1.0
533 cnst(4) = 0.0
534C
535 ELSE IF (q.EQ.101.OR.q.EQ.104.OR.q.EQ.107.OR.q.EQ.109) THEN
536C
537C*** DIRECTION OF WIND WAVES, SWELLS WAVES, PRIMARY WAVE, SECONDARY
538C*** WAVE (deg. true) --------------------
539C
540 cnst(1) = 0.0
541 cnst(2) = 1.0
542 cnst(3) = 20.0
543 cnst(4) = 0.0
544C
545 ELSE IF (q.EQ.103.OR.q.EQ.106.OR.q.EQ.108.OR.q.EQ.110) THEN
546C
547C*** MEAN PERIOD OF WIND WAVES, SWELLS WAVES, PRIMARY WAVE, SECONDARY
548C*** WAVE (s) --------------------
549C
550 cnst(1) = 0.0
551 cnst(2) = 1.0
552 cnst(3) = 2.0
553 cnst(4) = 0.0
554C
555 ELSE IF (q.EQ.111.OR.q.EQ.112.OR.q.EQ.113.OR.q.EQ.114.OR.
556 & q.EQ.115.OR.q.EQ.116.OR.q.EQ.117.OR.q.EQ.121.OR.
557 & q.EQ.122.OR.q.EQ.123) THEN
558C
559C*** NET SHORTWAVE RADITION (SURFACE) -NSWRS w/m **2
560C*** NET LONGWAVE RADITION (SURFACE) -SHTFL w/m**2
561C*** NET SHORTWAVE RADITION (TOP OF ATOMS.) -NSWRT w/m**2
562C*** NET LONGWAVE RADITION (TOP OF ATOMS.) -NLWRT w/m**2
563C*** LONG WAVE RADITION -LWAVR w/m**2
564C*** SHORT WAVE RADITION -SWAVE w/m**2
565C*** GLOBAL RADITION -G-RAD w/m**2
566C*** LATENT HEAT FLUX -LHTFL w/m**2
567C*** SENSIBLE HEAT FLUX -SHTFL w/m**2
568C*** BOUNDARY LAYER DISSIPATION -BLYDP w/m**2
569C
570 cnst(1) = 0.0
571 cnst(2) = 1.0
572 cnst(3) = 5.0
573 IF (q.EQ.114) cnst(3) = 20.0
574 cnst(4) = 0.0
575C
576 ELSE IF (q.EQ.127) THEN
577C
578C IMAGE DATA -IMG-D
579C
580 cnst(1) = 0.0
581 cnst(2) = 1.0
582 cnst(3) = 10.0
583 cnst(4) = 0.0
584C
585 ELSE IF (q.EQ.128) THEN
586C
587C Mean Sea Level Pressure -MSLSA (Pa)
588C (Standard Atmosphere Reduction)
589C
590 cnst(1) = 0.0
591 cnst(2) = 0.01
592 cnst(3) = 4.0
593 cnst(4) = 0.0
594C
595 ELSE IF (q.EQ.129) THEN
596C
597C Mean Sea Level Pressure -MSLMA (Pa)
598C (Maps System Reduction)
599C
600 cnst(1) = 0.0
601 cnst(2) = 0.01
602 cnst(3) = 4.0
603 cnst(4) = 0.0
604C
605 ELSE IF (q.EQ.130) THEN
606C
607C Mean Sea Level Pressure -MSLET (Pa)
608C (ETA Model Reduction)
609C
610 cnst(1) = 0.0
611 cnst(2) = 0.01
612 cnst(3) = 4.0
613 cnst(4) = 0.0
614C
615 ELSE IF (q.EQ.131.OR.q.EQ.132.OR.q.EQ.133.OR.q.EQ.134) THEN
616C
617C*** SURFACE LIFTED INDEX ..(DEG K)
618C*** BEST (4 LAYER) LIFTED INDEX ..(DEG K)
619C*** K INDEX ..(DEG K) TO DEG C.
620C*** SWEAT INDEX ..(DEG K) TO DEG C.
621C
622 IF (q.EQ.131.OR.q.EQ.132) THEN
623 cnst(1) = 0.0
624 ELSE
625 cnst(1) = -273.15
626 END IF
627 cnst(2) = 1.0
628 cnst(3) = 4.0
629 cnst(4) = 0.0
630C
631 ELSE IF (q.EQ.135) THEN
632C
633C*** HORIZONTIAL MOISTURE DIVERGENCE (KG/KG/S) -MCONV
634C
635 cnst(1) = 0.0
636 cnst(2) = 1.e+8
637 cnst(3) = 10.0
638 cnst(4) = 0.0
639C
640 ELSE IF (q.EQ.136) THEN
641C
642C*** VERTICAL SPEED SHEAR (1/SEC)... TO BE CONVERTED TO KNOTS/1000 FT
643C
644 cnst(1) = 0.0
645 cnst(2) = 592.086
646 cnst(3) = 2.0
647 cnst(4) = 0.0
648C
649 ELSE IF (q.EQ.137) THEN
650C
651C*** 3-hr pressure tendency (TSLSA) (Pa/s)
652C
653 cnst(1) = 0.0
654 cnst(2) = 1000.0
655 cnst(3) = 10.0
656 cnst(4) = 0.0
657C
658 ELSE IF (q.EQ.156) THEN
659C
660C*** CONVECTIVE INHIBITION -CIN-- (J/kg)
661C
662 cnst(1) = 0.0
663 cnst(2) = 1.0
664 cnst(3) = 10.0
665 cnst(4) = 0.0
666C
667 ELSE IF (q.EQ.157) THEN
668C
669C*** CONVECTIVE AVAILABLE POTENTIAL ENERGY -CAPE- (J/kg)
670C
671 cnst(1) = 0.0
672 cnst(2) = 1.0
673 cnst(3) = 500.0
674 cnst(4) = 0.0
675C
676 ELSE IF (q.EQ.158) THEN
677C
678C*** TURBULENT KINETIC ENERGY -TKE-- (J/kg)
679C
680 cnst(1) = 0.0
681 cnst(2) = 1.0
682 cnst(3) = 100.0
683 cnst(4) = 0.0
684C
685 ELSE IF (q.EQ.175) THEN
686C
687C*** MODEL LAYER NUMBER (FROM BOTTOM UP) -SGLYR (non-dim)
688C
689 cnst(1) = 0.0
690 cnst(2) = 1.0
691 cnst(3) = 1.0
692 cnst(4) = 0.0
693C
694 ELSE IF (q.EQ.176) THEN
695C
696C*** LATITUDE (-90 TO +90) -NLAT- (deg)
697C
698 cnst(1) = 0.0
699 cnst(2) = 1.0
700 cnst(3) = 10.0
701 cnst(4) = 0.0
702C
703 ELSE IF (q.EQ.177) THEN
704C
705C*** EAST LATITUDE (0-360) -ELON- (deg)
706C
707 cnst(1) = 0.0
708 cnst(2) = 1.0
709 cnst(3) = 10.0
710 cnst(4) = 0.0
711C
712 ELSE IF (q.EQ.201) THEN
713C
714C*** ICE-FREE WATER SURFACE -ICWAT (%)
715C
716 cnst(1) = 0.0
717 cnst(2) = 1.0
718 cnst(3) = 10.0
719 cnst(4) = 0.0
720C
721 ELSE IF (q.EQ.204) THEN
722C
723C*** DOWNWARD SHORT WAVE RAD. FLUX -DSWRF (W/m**2)
724C
725 cnst(1) = 0.0
726 cnst(2) = 1.0
727 cnst(3) = 10.0
728 cnst(4) = 0.0
729C
730 ELSE IF (q.EQ.205) THEN
731C
732C*** DOWNWARD LONG WAVE RAD. FLUX -DLWRF (W/m**2)
733C
734 cnst(1) = 0.0
735 cnst(2) = 1.0
736 cnst(3) = 10.0
737 cnst(4) = 0.0
738C
739 ELSE IF (q.EQ.207) THEN
740C
741C*** MOISTURE AVAILABILITY -MSTAV (%)
742C
743 cnst(1) = 0.0
744 cnst(2) = 1.0
745 cnst(3) = 10.0
746 cnst(4) = 0.0
747C
748 ELSE IF (q.EQ.208) THEN
749C
750C*** EXCHANGE COEFFICIENT -SFEXC (kg/m**3)(m/s)
751C
752 cnst(1) = 0.0
753 cnst(2) = 1.0
754 cnst(3) = 10.0
755 cnst(4) = 0.0
756CC
757 ELSE IF (q.EQ.209) THEN
758C
759C*** NO. OF MIXED LAYERS NEXT TO SURFACE -MIXLY (integer)
760C
761 cnst(1) = 0.0
762 cnst(2) = 1.0
763 cnst(3) = 10.0
764 cnst(4) = 0.0
765C
766 ELSE IF (q.EQ.211) THEN
767C
768C*** UPWARD SHORT WAVE RAD. FLUX -USWRF (W/m**2)
769C
770 cnst(1) = 0.0
771 cnst(2) = 1.0
772 cnst(3) = 10.0
773 cnst(4) = 0.0
774C
775 ELSE IF (q.EQ.212) THEN
776C
777C*** UPWARD LONG WAVE RAD. FLUX -ULWRF (W/m**2)
778C
779 cnst(1) = 0.0
780 cnst(2) = 1.0
781 cnst(3) = 10.0
782 cnst(4) = 0.0
783C
784 ELSE IF (q.EQ.213) THEN
785C
786C*** AMOUNT OF NON-CONVECTIVE CLOUD -CDLYR (%)
787C
788 cnst(1) = 0.0
789 cnst(2) = 1.0
790 cnst(3) = 10.0
791 cnst(4) = 0.0
792C
793 ELSE IF (q.EQ.216) THEN
794C
795C*** TEMPERATURE TENDENCY BY ALL RADIATION -TTRAD (Deg. K/s)
796C
797 cnst(1) = 0.0
798 cnst(2) = 1.0
799 cnst(3) = 10.0
800 cnst(4) = 0.0
801C
802 ELSE IF (q.EQ.218) THEN
803C
804C*** PRECIP. INDEX (0.0-1.00) -PREIX (note will look like %)
805C
806 cnst(1) = 0.0
807 cnst(2) = 100.0
808 cnst(3) = 10.0
809 cnst(4) = 0.0
810C
811 ELSE IF (q.EQ.220) THEN
812C
813C*** NATURAL LOG OF SURFACE PRESSURE -NLGSP ln(kPa)
814C
815 cnst(1) = 0.0
816 cnst(2) = 1.0
817 cnst(3) = 10.0
818 cnst(4) = 0.0
819C
820C*** NONE OF THE ABOVE ....
821C
822 ELSE
823C
824C SET DEFAULT VALUES
825C
826 cnst(1) = 0.0
827 cnst(2) = 1.0
828 cnst(3) = 5.0
829 cnst(4) = 0.0
830 ier = 1
831 END IF
832C
833 RETURN
834 END
w3fi69
subroutine w3fi69(pds, id)
Converts an edition 1 grib produce definition section (pds) to a 25, or 27 word integer array.
Definition w3fi69.f:29
w3fi70
subroutine w3fi70(pds, cnst, ier)
Computes the four scaling constants used by grdprt, w3fp03, or w3fp05 from the 28 byte (pds) product ...
Definition w3fi70.f:21