-
Notifications
You must be signed in to change notification settings - Fork 0
/
type.f90
539 lines (459 loc) · 26.1 KB
/
type.f90
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
MODULE TYPES
! Definitions of various derived data types
USE PRECISION_PARAMETERS
USE GLOBAL_CONSTANTS, ONLY : NULL_BOUNDARY,NEUMANN
IMPLICIT NONE
CHARACTER(255), PARAMETER :: typeid='$Id: type.f90 9814 2012-01-11 17:49:04Z randy.mcdermott $'
CHARACTER(255), PARAMETER :: typerev='$Revision: 9814 $'
CHARACTER(255), PARAMETER :: typedate='$Date: 2012-01-11 09:49:04 -0800 (Wed, 11 Jan 2012) $'
TYPE PARTICLE_CLASS_TYPE
CHARACTER(30) :: ID,SPEC_ID,DEVC_ID='null',CTRL_ID='null',QUANTITIES(10),SMOKEVIEW_BAR_LABEL(10),SURF_ID='null',PROP_ID='null',&
RADIATIVE_PROPERTY_TABLE_ID='null',CNF_RAMP_ID='null',DISTRIBUTION='ROSIN-RAMMLER-LOGNORMAL',&
CHILD_DISTRIBUTION='ROSIN-RAMMLER-LOGNORMAL',CHILD_CNF_RAMP_ID='null'
CHARACTER(60) :: SMOKEVIEW_LABEL(10)
CHARACTER(25) :: VEG_DEGRADATION
REAL(EB) :: DENSITY,TMP_V=-1._EB,TMP_MELT=-1._EB,FTPR,HEAT_OF_COMBUSTION,ADJUST_EVAPORATION, &
LIFETIME,DIAMETER,MINIMUM_DIAMETER,MAXIMUM_DIAMETER,GAMMA,KILL_RADIUS, &
TMP_INITIAL,SIGMA,VERTICAL_VELOCITY,HORIZONTAL_VELOCITY, &
VEG_SV,VEG_MOISTURE,VEG_CHAR_FRACTION,VEG_ASH_FRACTION,VEG_DRAG_COEFFICIENT,VEG_BULK_DENSITY, &
VEG_DENSITY,VEG_BURNING_RATE_MAX,VEG_DEHYDRATION_RATE_MAX,VEG_INITIAL_TEMPERATURE, &
VEG_FUEL_MPV_MIN,VEG_MOIST_MPV_MIN,VEG_ASH_MPV_MAX,USER_DRAG_COEFFICIENT, &
SURFACE_TENSION,BREAKUP_CHILD_DIAMETER,BREAKUP_CHILD_GAMMA,BREAKUP_CHILD_SIGMA,DENSE_VOLUME_FRACTION, &
REAL_REFRACTIVE_INDEX,COMPLEX_REFRACTIVE_INDEX,TOL_INT,VEG_CHAR_MPV_MIN
REAL(EB) :: VEG_H_H2O,VEG_A_H2O,VEG_E_H2O,VEG_H_PYR,VEG_A_PYR,VEG_E_PYR,VEG_H_CHAR,VEG_A_CHAR,VEG_E_CHAR, &
VEG_BETA_CHAR,VEG_NU_O2_CHAR,VEG_CHAR_ENTHALPY_FRACTION
INTEGER :: VEG_NDT_SUBCYCLES
REAL(EB), POINTER, DIMENSION(:) :: R_CDF,CDF,CHILD_R_CDF,CHILD_CDF,W_CDF,R50,LAMBDA
REAL(EB), POINTER, DIMENSION(:,:) :: WQABS,WQSCA,TABLE_ROW
REAL(EB), POINTER, DIMENSION(:,:) :: ORIENTATION
INTEGER :: SAMPLING,N_QUANTITIES,QUANTITIES_INDEX(10),EVAP_INDEX=0,RGB(3),RADIATIVE_PROPERTY_INDEX=0, &
SURF_INDEX=-1,DRAG_LAW=1,DEVC_INDEX=0,CTRL_INDEX=0,PROP_INDEX=-1,N_SPLIT=1,SPLIT_PART_INDEX=0,Z_INDEX=-1,&
Y_INDEX=-1,CNF_RAMP_INDEX=-1,CHILD_CNF_RAMP_INDEX=-1
INTEGER, POINTER, DIMENSION(:) :: IL_CDF,IU_CDF
LOGICAL :: STATIC=.FALSE.,WATER=.FALSE.,FUEL=.FALSE.,MASSLESS,TREE=.FALSE.,MONODISPERSE=.FALSE.,EVAPORATE=.TRUE., &
BREAKUP=.FALSE.,CHECK_DISTRIBUTION=.FALSE.
LOGICAL VEG_REMOVE_CHARRED,VEG_STEM,VEG_CHAR_OXIDATION
END TYPE PARTICLE_CLASS_TYPE
TYPE (PARTICLE_CLASS_TYPE), DIMENSION(:), ALLOCATABLE, TARGET :: PARTICLE_CLASS
TYPE PARTICLE_TYPE
REAL(EB) :: X,Y,Z,TMP,U,V,W,R,PWT,ACCEL_X,ACCEL_Y,ACCEL_Z,T,RE=0._EB, &
VEG_KAPPA,VEG_EMISS,VEG_DIVQC,VEG_DIVQR,VEG_PACKING_RATIO,VEG_FUEL_MASS,VEG_MOIST_MASS,VEG_IGN_TON, &
VEG_IGN_TOFF,VEG_IGN_TRAMPON,VEG_IGN_TRAMPOFF,VEG_CHAR_MASS,VEG_ASH_MASS,VEG_SV,VEG_MLR,&
VEG_VOLFRACTION
LOGICAL :: IGNITOR,VEG_IGNITED
INTEGER :: IOR,CLASS,TAG,WALL_INDEX=0,VEG_N_TREE_OUTPUT,SPLIT_IOR,SURF_INDEX=0
LOGICAL :: SHOW,SPLAT=.FALSE.
!! Shift to use 1D structure
! REAL(EB), POINTER, DIMENSION(:) :: TMP_S, LAYER_THICKNESS,X_S
! REAL(EB), POINTER, DIMENSION(:,:) :: RHO_S
! INTEGER, POINTER, DIMENSION(:) :: N_LAYER_CELLS
! REAL(EB), POINTER, DIMENSION(:) :: AWM_AEROSOL,A_LP_MPUA,MASSFLUX_ACTUAL,MASSFLUX,SURFACE_DENSITY,LP_CPUA,LP_MPUA
! LOGICAL :: SHRINKING
END TYPE PARTICLE_TYPE
TYPE WALL_TYPE
REAL(EB), POINTER, DIMENSION(:) :: TMP_S, LAYER_THICKNESS,X_S
REAL(EB), POINTER, DIMENSION(:,:) :: ILW,RHO_S
INTEGER, POINTER, DIMENSION(:) :: N_LAYER_CELLS
LOGICAL :: SHRINKING, BURNAWAY,ALREADY_ALLOCATED=.FALSE.
REAL(EB), POINTER, DIMENSION(:) :: VEG_ASHMASS_L,VEG_CHARMASS_L,VEG_FUELMASS_L,VEG_MOISTMASS_L,VEG_TMP_L
REAL(EB) :: VEG_HEIGHT
REAL(EB) :: AREA_ADJUST=1._EB,AW=0._EB,DUWDT=0._EB,E_WALL,EW=0._EB,HEAT_TRANS_COEF=0._EB,KW=1._EB, &
QCONF=0._EB,QRADIN=0._EB,QRADOUT=0._EB,RAW=0._EB,RDN=1._EB,RHO_F=1._EB,TMP_B=0._EB,TMP_F=0._EB, &
TW=0._EB,U_TAU=0._EB,UVW_GHOST=0._EB,UW=0._EB,UWS=0._EB,UW0=0._EB,XW=0._EB,YW=0._EB,ZW=0._EB
REAL(EB), POINTER, DIMENSION(:) :: AWM_AEROSOL,A_LP_MPUA,MASSFLUX_ACTUAL,MASSFLUX,SURFACE_DENSITY,LP_CPUA,LP_MPUA,RHODW,ZZ_F
INTEGER :: BACK_INDEX=0,BOUNDARY_TYPE=NULL_BOUNDARY,SURF_INDEX_ORIG=0,OBST_INDEX=0,NPPCW=1,PRESSURE_BC_INDEX=NEUMANN,&
PRESSURE_ZONE_WALL=0,VENT_INDEX=0
INTEGER :: II=0,JJ=0,KK=0,IIG=0,JJG=0,KKG=0,SURF_INDEX=0,IOR=0,NOM=0,NOM_IB(6)=0
END TYPE WALL_TYPE
TYPE SPECIES_TYPE
REAL(EB) :: MW=0._EB,YY0=0._EB,RCON,MAXMASS,MASS_EXTINCTION_COEFFICIENT=0._EB,&
SPECIFIC_HEAT=-1._EB,REFERENCE_ENTHALPY=-1._EB,&
REFERENCE_TEMPERATURE,MU_USER=-1._EB,K_USER=-1._EB,D_USER=-1._EB,EPSK=-1._EB,SIG=-1._EB,&
FLD_LETHAL_DOSE=0._EB,FIC_CONCENTRATION=0._EB,&
SPECIFIC_HEAT_LIQUID=-1,DENSITY_LIQUID,VAPORIZATION_TEMPERATURE,HEAT_OF_VAPORIZATION=-1._EB,MELTING_TEMPERATURE,&
H_V_REFERENCE_TEMPERATURE=-1._EB,H_V_CORRECTOR=0._EB ,TMP_V=-1._EB,TMP_MELT=-1._EB,ATOMS(118)=0._EB,&
MEAN_DIAMETER=1.E-6_EB,CONDUCTIVITY_SOLID,DENSITY_SOLID
LOGICAL :: ABSORBING=.FALSE.,ISFUEL=.FALSE.,BACKGROUND=.FALSE.,SMIX_COMPONENT_ONLY=.FALSE.,LISTED=.FALSE.
CHARACTER(30) :: ID,RAMP_CP,RAMP_CP_L,RAMP_K,RAMP_MU,RAMP_D
CHARACTER(100) :: FORMULA
INTEGER :: MODE=2,RAMP_CP_INDEX=-1,RAMP_CP_L_INDEX=-1,RAMP_K_INDEX=-1,RAMP_MU_INDEX=-1,RAMP_D_INDEX=-1
REAL(EB),POINTER,DIMENSION(:) :: H_V,C_P_L,C_P_L_BAR,H_L
END TYPE SPECIES_TYPE
TYPE (SPECIES_TYPE), DIMENSION(:), ALLOCATABLE, TARGET :: SPECIES
TYPE SPECIES_MIXTURE_TYPE
REAL(EB), POINTER, DIMENSION(:) :: MASS_FRACTION,VOLUME_FRACTION
REAL(EB) :: MW , RCON, ZZ0=0._EB, MASS_EXTINCTION_COEFFICIENT=0._EB,ADJUST_NU=1._EB,ATOMS(118)=0._EB,MEAN_DIAMETER,&
DENSITY_SOLID,CONDUCTIVITY_SOLID
CHARACTER(30), POINTER, DIMENSION(:) :: SPEC_ID
CHARACTER(30) :: ID='null'
CHARACTER(100) :: FORMULA='null'
INTEGER :: AWM_INDEX = -1
LOGICAL :: DEPOSITING=.FALSE.,VALID_ATOMS=.TRUE.
END TYPE SPECIES_MIXTURE_TYPE
TYPE (SPECIES_MIXTURE_TYPE), DIMENSION(:), ALLOCATABLE, TARGET :: SPECIES_MIXTURE
TYPE REACTION_TYPE
CHARACTER(30) :: FUEL,OXIDIZER,PRODUCTS,ID
CHARACTER(30), DIMENSION(:), POINTER :: SMIX_ID,SMIX_ID_READ,SPEC_ID,SPEC_ID_READ
REAL(EB) :: C,H,N,O,EPUMO2,HEAT_OF_COMBUSTION,HOC_COMPLETE,A,A_IN,E,E_IN,TOL_INT, &
Y_O2_INFTY,Y_N2_INFTY=0._EB,MW_FUEL,MW_SOOT, &
CO_YIELD,SOOT_YIELD,H2_YIELD, SOOT_H_FRACTION,RHO_EXPONENT, CRIT_FLAME_TMP, &
NU_O2=0._EB,NU_N2=0._EB,NU_H2O=0._EB,NU_CO2=0._EB,NU_CO=0._EB,NU_SOOT=0._EB
REAL(EB) :: AUTO_IGNITION_TEMPERATURE=0._EB,THRESHOLD_TEMP=0._EB,N_T=0._EB
REAL(EB), DIMENSION(:), POINTER :: NU,NU_READ,NU_SPECIES,N_S,N_S_READ
INTEGER :: FUEL_SMIX_INDEX=-1,MODE,N_SMIX,N_SPEC
LOGICAL :: IDEAL,CHECK_ATOM_BALANCE,REVERSIBLE,NEW_MIX
CHARACTER(100) :: FYI='null'
CHARACTER(255) :: EQUATION
CHARACTER(100) :: FWD_ID
END TYPE REACTION_TYPE
TYPE (REACTION_TYPE), DIMENSION(:), ALLOCATABLE, TARGET :: REACTION
TYPE MATERIAL_TYPE
REAL(EB) :: K_S,C_S,RHO_S,EMISSIVITY,DIFFUSIVITY,KAPPA_S,MOISTURE_FRACTION,TMP_BOIL,INIT_VAPOR_FLUX
INTEGER :: PYROLYSIS_MODEL
CHARACTER(30) :: RAMP_K_S,RAMP_C_S
INTEGER :: N_REACTIONS,PROP_INDEX=-1,CABL_INDEX=-1
INTEGER, DIMENSION(MAX_REACTIONS) :: N_RESIDUE
INTEGER, DIMENSION(MAX_MATERIALS,MAX_REACTIONS) :: RESIDUE_MATL_INDEX
REAL(EB), DIMENSION(MAX_REACTIONS) :: TMP_REF,TMP_IGN,TMP_THR,RATE_REF,THR_SIGN
REAL(EB), DIMENSION(MAX_MATERIALS,MAX_REACTIONS) :: NU_RESIDUE
REAL(EB), DIMENSION(MAX_REACTIONS) :: A,E,H_R,N_S,N_T,HEATING_RATE,PYROLYSIS_RANGE,HEAT_OF_COMBUSTION,TOL_INT
REAL(EB), DIMENSION(MAX_SPECIES,MAX_REACTIONS) :: NU_SPEC,NU_GAS,ADJUST_BURN_RATE
LOGICAL, DIMENSION(MAX_REACTIONS) :: PCR = .FALSE.
CHARACTER(30), DIMENSION(MAX_MATERIALS,MAX_REACTIONS) :: RESIDUE_MATL_NAME
CHARACTER(30), DIMENSION(MAX_SPECIES,MAX_REACTIONS) :: SPEC_ID
CHARACTER(100) :: FYI='null'
LOGICAL :: USER_DEFINED=.TRUE.,CONDUIT=.FALSE.,AIR=.FALSE.
END TYPE MATERIAL_TYPE
TYPE (MATERIAL_TYPE), DIMENSION(:), ALLOCATABLE, TARGET :: MATERIAL
TYPE SURFACE_TYPE
REAL(EB) :: TMP_FRONT=-1._EB,VEL,PLE,Z0,CONVECTIVE_HEAT_FLUX,NET_HEAT_FLUX, &
VOLUME_FLUX,HRRPUA,MLRPUA,T_IGN,SURFACE_DENSITY,CELL_SIZE_FACTOR, &
E_COEFFICIENT,TEXTURE_WIDTH,TEXTURE_HEIGHT,THICKNESS,EXTERNAL_FLUX,TMP_BACK=-1._EB, &
DXF,DXB,MASS_FLUX_TOTAL,PARTICLE_MASS_FLUX,EMISSIVITY,MAX_PRESSURE,&
TMP_IGN,H_V,REGRID_FACTOR,LAYER_DIVIDE,ROUGHNESS,RADIUS=-1._EB,LENGTH=-1._EB,WIDTH=-1._EB, &
DT_INSERT,H_FIXED=-1._EB,EMISSIVITY_BACK,CONV_LENGTH,XYZ(3),FIRE_SPREAD_RATE, &
MINIMUM_LAYER_THICKNESS
REAL(EB), POINTER, DIMENSION(:) :: DX,RDX,RDXN,X_S,DX_WGT,MF_FRAC,PARTICLE_INSERT_CLOCK
REAL(EB), DIMENSION(0:MAX_SPECIES) :: MASS_FRACTION,MASS_FLUX
REAL(EB), DIMENSION(-5:MAX_SPECIES) :: TAU,ADJUST_BURN_RATE=1._EB
INTEGER, DIMENSION(-5:MAX_SPECIES) :: RAMP_INDEX=0
INTEGER, DIMENSION(3) :: RGB
REAL(EB) :: TRANSPARENCY
REAL(EB), DIMENSION(2) :: VEL_T
INTEGER, DIMENSION(2) :: LEAK_PATH,DUCT_PATH
INTEGER :: THERMAL_BC_INDEX,NPPC,SPECIES_BC_INDEX,VELOCITY_BC_INDEX,SURF_TYPE,N_CELLS,PART_INDEX,PROP_INDEX=-1,CABL_INDEX=-1
INTEGER :: PYROLYSIS_MODEL
INTEGER :: N_LAYERS,N_MATL
INTEGER, POINTER, DIMENSION(:) :: N_LAYER_CELLS,LAYER_INDEX,MATL_INDEX
INTEGER, DIMENSION(MAX_LAYERS,MAX_MATERIALS) :: LAYER_MATL_INDEX
INTEGER, DIMENSION(MAX_LAYERS) :: N_LAYER_MATL
INTEGER, POINTER, DIMENSION(:,:,:) :: RESIDUE_INDEX
REAL(EB), POINTER, DIMENSION(:) :: MIN_DIFFUSIVITY
REAL(EB), DIMENSION(MAX_LAYERS) :: LAYER_THICKNESS, LAYER_DENSITY, TMP_INNER, STRETCH_FACTOR
CHARACTER(30), POINTER, DIMENSION(:) :: MATL_NAME
CHARACTER(30), DIMENSION(MAX_LAYERS,MAX_MATERIALS) :: LAYER_MATL_NAME
REAL(EB), DIMENSION(MAX_LAYERS,MAX_MATERIALS) :: LAYER_MATL_FRAC
LOGICAL :: BURN_AWAY,ADIABATIC,THERMALLY_THICK,INTERNAL_RADIATION,SHRINK,USER_DEFINED=.TRUE., &
FREE_SLIP=.FALSE.,NO_SLIP=.FALSE.,SPECIFIED_NORMAL_VELOCITY=.FALSE.,SPECIFIED_TANGENTIAL_VELOCITY=.FALSE.
INTEGER :: GEOMETRY,BACKING,PROFILE
CHARACTER(30) :: PART_ID,RAMP_Q,RAMP_V,RAMP_T,RAMP_EF,RAMP_PART
CHARACTER(30), DIMENSION(0:MAX_SPECIES) :: RAMP_MF
CHARACTER(60) :: ID,TEXTURE_MAP
CHARACTER(100) :: FYI='null'
! Boundary vegetation
LOGICAL :: VEGETATION=.FALSE.,VEG_NO_BURN=.FALSE.,VEG_GROUND_ZERO_RAD=.TRUE., &
VEG_CHAR_OXIDATION,VEG_HCONV_CYLRE,VEG_HCONV_CYLLAM
INTEGER :: NVEG_L
REAL(EB) :: VEG_ASH_FRACTION,VEG_CHAR_FRACTION,VEG_DRAG_INI,VEG_ELEMENT_DENSITY,VEG_HEIGHT,VEG_INITIAL_TEMP,VEG_LOAD, &
VEG_MOISTURE,VEG_SV,VEG_PACKING,VEG_KAPPA,FIRELINE_MLR_MAX,VEG_GROUND_TEMP
REAL(EB) :: VEG_H_H2O,VEG_A_H2O,VEG_E_H2O,VEG_H_PYR,VEG_A_PYR,VEG_E_PYR,VEG_H_CHAR,VEG_A_CHAR,VEG_E_CHAR, &
VEG_BETA_CHAR,VEG_NU_O2_CHAR,VEG_CHAR_ENTHALPY_FRACTION
INTEGER, POINTER, DIMENSION(:) :: VEG_KGAS_L
REAL(EB), POINTER, DIMENSION(:) :: VEG_FUEL_FLUX_L,VEG_MOIST_FLUX_L,VEG_DIVQNET_L
REAL(EB), POINTER, DIMENSION(:) :: VEG_INCM_RADFCT_L,VEG_FINCM_RADFCT_L,VEG_FINCP_RADFCT_L !add index for mult veg
REAL(EB), POINTER, DIMENSION(:,:) :: VEG_SEMISSM_RADFCT_L,VEG_SEMISSP_RADFCT_L !add index for mult veg
! Level Set Firespread
CHARACTER(25) :: VEG_DEGRADATION,VEG_LSET_SURFACE_FIRE_HEAD_ROS_MODEL,VEG_LSET_CROWN_FIRE_HEAD_ROS_MODEL
LOGICAL :: VEG_LSET_CROWN_VEG,VEG_LSET_HEADWIDTH_DEPENDENCE,VEG_LSET_IGNITE,VEG_LSET_SPREAD,VEG_LSET_ELLIPSE, &
VEG_LSET_TAN2
REAL(EB) :: VEG_LSET_IGNITE_TIME,VEG_LSET_ROS_HEAD,VEG_LSET_ROTH_ZEROWINDSLOPE_ROS,VEG_LSET_QCON,VEG_LSET_ROS_FLANK, &
VEG_LSET_ROS_BACK, &
VEG_LSET_WIND_EXP,VEG_LSET_PHIDEPTH,VEG_LSET_HEAT_OF_COMBUSTION,VEG_LSET_SURF_LOAD,VEG_LSET_CHAR_FRACTION, &
VEG_LSET_SURF_EFFM,VEG_LSET_CANOPY_BULK_DENSITY,VEG_LSET_CANOPY_HEIGHT,VEG_LSET_FUEL_STRATA_GAP, &
VEG_LSET_CRUZ_PROB_PASSIVE,VEG_LSET_CRUZ_PROB_ACTIVE,VEG_LSET_CRUZ_PROB_CROWN,VEG_LSET_WAF_UNSHELTERED, &
VEG_LSET_WAF_SHELTERED
!Veg SAV and height (fuel depth) for Farsite emulation
REAL(EB) :: VEG_LSET_SIGMA,VEG_LSET_SURF_HEIGHT,VEG_LSET_BETA
! Unstructured geometry
CHARACTER(60) :: BC_FILENAME='null'
END TYPE SURFACE_TYPE
TYPE (SURFACE_TYPE), DIMENSION(:), ALLOCATABLE, TARGET :: SURFACE
TYPE WALL_ILW_TYPE
REAL(EB), POINTER, DIMENSION(:,:) :: ILW
END TYPE WALL_ILW_TYPE
TYPE OMESH_TYPE
REAL(EB), POINTER, DIMENSION(:,:,:) :: MU,RHO,RHOS,U,V,W,US,VS,WS,H,HS,FVX,FVY,FVZ,D,DS,KRES
REAL(EB), POINTER, DIMENSION(:,:,:,:) :: ZZ,ZZS
INTEGER, POINTER, DIMENSION(:,:) :: IJKW
INTEGER, POINTER, DIMENSION(:) :: BOUNDARY_TYPE
TYPE(WALL_ILW_TYPE), POINTER, DIMENSION(:) :: WALL_ILW
TYPE(PARTICLE_TYPE), POINTER, DIMENSION(:) :: PARTICLE
INTEGER :: N_DROP_ORPHANS,N_DROP_ORPHANS_DIM,N_DROP_ADOPT, &
I_MIN_R=-10,I_MAX_R=-10,J_MIN_R=-10,J_MAX_R=-10,K_MIN_R=-10,K_MAX_R=-10,NIC_R=0, &
I_MIN_S=-10,I_MAX_S=-10,J_MIN_S=-10,J_MAX_S=-10,K_MIN_S=-10,K_MAX_S=-10,NIC_S=0
REAL(EB), POINTER, DIMENSION(:) :: &
REAL_SEND_PKG1,REAL_SEND_PKG2,REAL_SEND_PKG3,REAL_SEND_PKG4,REAL_SEND_PKG5,REAL_SEND_PKG6,REAL_SEND_PKG7, &
REAL_RECV_PKG1,REAL_RECV_PKG2,REAL_RECV_PKG3,REAL_RECV_PKG4,REAL_RECV_PKG5,REAL_RECV_PKG6,REAL_RECV_PKG7
INTEGER , POINTER, DIMENSION(:) :: INTG_SEND_PKG1,INTG_RECV_PKG1
INTEGER , POINTER, DIMENSION(:,:) :: INTG_SEND_PKGW
LOGICAL , POINTER, DIMENSION(:) :: LOGI_SEND_PKG1,LOGI_RECV_PKG1
END TYPE OMESH_TYPE
TYPE OBSTRUCTION_TYPE
CHARACTER(30) :: DEVC_ID='null',CTRL_ID='null',PROP_ID='null'
INTEGER, DIMENSION(-3:3) :: SURF_INDEX=0
LOGICAL, DIMENSION(-3:3) :: SHOW_BNDF=.TRUE.
INTEGER, DIMENSION(3) :: RGB=(/0,0,0/)
INTEGER, DIMENSION(3) :: DIMENSIONS=0
REAL(EB) :: TRANSPARENCY=1._EB,BULK_DENSITY=-1._EB,VOLUME_ADJUST=1._EB
REAL(EB), DIMENSION(3) :: TEXTURE=0._EB
REAL(EB) :: X1=0._EB,X2=1._EB,Y1=0._EB,Y2=1._EB,Z1=0._EB,Z2=1._EB,MASS=1.E6_EB
REAL(EB), DIMENSION(3) :: FDS_AREA=-1._EB,INPUT_AREA=-1._EB
INTEGER :: I1=-1,I2=-1,J1=-1,J2=-1,K1=-1,K2=-1,COLOR_INDICATOR=-1,TYPE_INDICATOR=-1,ORDINAL=0
INTEGER :: DEVC_INDEX=-1,CTRL_INDEX=-1,PROP_INDEX=-1
LOGICAL :: SAWTOOTH=.TRUE.,HIDDEN=.FALSE.,PERMIT_HOLE=.TRUE.,ALLOW_VENT=.TRUE.,CONSUMABLE=.FALSE.,REMOVABLE=.FALSE., &
THIN=.FALSE.,HOLE_FILLER=.FALSE.,NOTERRAIN=.FALSE.
END TYPE OBSTRUCTION_TYPE
TYPE GEOMETRY_TYPE
CHARACTER(30) :: ID='geom',SMVOBJECT='sensor',TFILE='null'
CHARACTER(30) :: SHAPE='SPHERE',SURF_ID='null'
REAL(EB) :: X1=0._EB,X2=0._EB,Y1=0._EB,Y2=0._EB,Z1=0._EB,Z2=0._EB
REAL(EB) :: X0=0._EB,Y0=0._EB,Z0=0._EB,X=0._EB,Y=0._EB,Z=0._EB,XOR=0._EB,YOR=0._EB,ZOR=0._EB
REAL(EB) :: RADIUS=1.E6_EB
REAL(EB) :: U0=0._EB,V0=0._EB,W0=0._EB,U=0._EB,V=0._EB,W=0._EB
REAL(EB) :: OMEGA=0._EB ! rotation rate (rad/s) about orientation axis
REAL(EB) :: OMEGA_X=0._EB,OMEGA_Y=0._EB,OMEGA_Z=0._EB ! rotation rate about grid basis vectors
LOGICAL :: TRANSLATE=.FALSE.,ROTATE=.FALSE.,HOLE=.FALSE.,TWO_SIDED=.FALSE.
INTEGER :: ISHAPE
INTEGER, POINTER, DIMENSION(:) :: MIN_I,MAX_I,MIN_J,MAX_J,MIN_K,MAX_K
INTEGER, DIMENSION(3) :: RGB=(/192,192,192/)
REAL(EB), DIMENSION(3) :: NN=(/0._EB,0._EB,0._EB/),HL=(/0.5_EB,0.5_EB,0.5_EB/)
REAL(EB) :: PIXELS=1.0
REAL(EB) :: ROUGHNESS
END TYPE GEOMETRY_TYPE
TYPE(GEOMETRY_TYPE), POINTER, DIMENSION(:) :: GEOMETRY
TYPE VERTEX_TYPE
REAL(EB) :: X,Y,Z
END TYPE VERTEX_TYPE
TYPE(VERTEX_TYPE), POINTER, DIMENSION(:) :: VERTEX
! http://www.sdsc.edu/~tkaiser/f90.html#Linked lists
TYPE LINKED_LIST_TYPE
INTEGER :: INDEX ! data
TYPE(LINKED_LIST_TYPE), POINTER :: NEXT ! pointer to the next element
END TYPE LINKED_LIST_TYPE
TYPE CUTCELL_LINKED_LIST_TYPE
INTEGER :: INDEX ! data
TYPE(CUTCELL_LINKED_LIST_TYPE), POINTER :: NEXT ! pointer to the next element
REAL(EB) :: AREA ! cutcell area for index
END TYPE CUTCELL_LINKED_LIST_TYPE
TYPE FACET_TYPE
INTEGER :: VERTEX(3)=0,SURF_INDEX=0
CHARACTER(30) :: SURF_ID='null'
REAL(EB) :: NVEC(3)=0._EB,TMP_F=300._EB,TMP_G=300._EB,AREA=0._EB
TYPE(CUTCELL_LINKED_LIST_TYPE), POINTER :: CUTCELL_LIST
TYPE(SURFACE_TYPE), POINTER :: SURFACE_QUANTITY
END TYPE FACET_TYPE
TYPE(FACET_TYPE), POINTER, DIMENSION(:) :: FACET
TYPE VOLUME_TYPE
INTEGER :: VERTEX(4)=0
CHARACTER(30) :: MATL_ID='null'
END TYPE VOLUME_TYPE
TYPE(VOLUME_TYPE), POINTER, DIMENSION(:) :: VOLUME
TYPE CSVF_TYPE
CHARACTER(256) :: CSVFILE,UVWFILE
END TYPE CSVF_TYPE
TYPE(CSVF_TYPE), POINTER, DIMENSION(:) :: CSVFINFO
TYPE VENTS_TYPE
INTEGER :: I1=-1,I2=-1,J1=-1,J2=-1,K1=-1,K2=-1,BOUNDARY_TYPE=0,IOR=0,SURF_INDEX=0,DEVC_INDEX=-1,CTRL_INDEX=-1, &
COLOR_INDICATOR=99,TYPE_INDICATOR=0,ORDINAL=0,PRESSURE_RAMP_INDEX=0,NODE_INDEX=-1
INTEGER, DIMENSION(3) :: RGB=-1
REAL(EB) :: TRANSPARENCY = 1._EB
REAL(EB), DIMENSION(3) :: TEXTURE=0._EB
REAL(EB) :: X1=0._EB,X2=0._EB,Y1=0._EB,Y2=0._EB,Z1=0._EB,Z2=0._EB,FDS_AREA=0._EB,TOTAL_INPUT_AREA=0._EB, &
X0=-9.E6_EB,Y0=-9.E6_EB,Z0=-9.E6_EB,FIRE_SPREAD_RATE=0.05_EB,INPUT_AREA=0._EB,TMP_EXTERIOR=-1000._EB, &
DYNAMIC_PRESSURE=0._EB,UVW(3)=-1.E12_EB
LOGICAL :: ACTIVATED=.TRUE.,GHOST_CELLS_ONLY=.FALSE.
CHARACTER(30) :: DEVC_ID='null',CTRL_ID='null',ID='null'
! turbulent inflow (experimental)
INTEGER :: N_EDDY=0
REAL(EB) :: R_IJ(3,3)=0._EB,A_IJ(3,3),SIGMA_IJ(3,3),EDDY_BOX_VOLUME=0._EB, &
X_EDDY_MIN=0._EB,X_EDDY_MAX=0._EB, &
Y_EDDY_MIN=0._EB,Y_EDDY_MAX=0._EB, &
Z_EDDY_MIN=0._EB,Z_EDDY_MAX=0._EB
REAL(EB), POINTER, DIMENSION(:,:) :: U_EDDY,V_EDDY,W_EDDY
REAL(EB), POINTER, DIMENSION(:) :: X_EDDY,Y_EDDY,Z_EDDY,CU_EDDY,CV_EDDY,CW_EDDY
END TYPE VENTS_TYPE
TYPE TABLES_TYPE
INTEGER :: NUMBER_ROWS,NUMBER_COLUMNS
REAL(EB), POINTER, DIMENSION (:,:) :: TABLE_DATA
END TYPE TABLES_TYPE
TYPE (TABLES_TYPE), DIMENSION(:), ALLOCATABLE, TARGET :: TABLES
TYPE RAMPS_TYPE
REAL(EB) :: SPAN,DT,T_MIN,T_MAX
REAL(EB), POINTER, DIMENSION(:) :: INDEPENDENT_DATA,DEPENDENT_DATA,INTERPOLATED_DATA
INTEGER :: NUMBER_DATA_POINTS,NUMBER_INTERPOLATION_POINTS,DEVC_INDEX=-1,CTRL_INDEX=-1
CHARACTER(30) :: DEVC_ID='null',CTRL_ID='null'
END TYPE RAMPS_TYPE
TYPE (RAMPS_TYPE), DIMENSION(:), ALLOCATABLE, TARGET :: RAMPS
TYPE HUMAN_TYPE
CHARACTER(60) :: NODE_NAME='null'
CHARACTER(30) :: FFIELD_NAME='null'
REAL(EB) :: X=0._EB,Y=0._EB,Z=0._EB,U=0._EB,V=0._EB,W=0._EB,F_X=0._EB,F_Y=0._EB,&
X_old=0._EB,Y_old=0._EB,X_group=0._EB,Y_group=0._EB
REAL(EB) :: UBAR=0._EB, VBAR=0._EB, UBAR_Center=0._EB, VBAR_Center=0._EB
REAL(EB) :: Speed=1.25_EB, Radius=0.255_EB, Mass=80.0_EB, Tpre=1._EB, Tau=1._EB, &
Eta=0._EB, Ksi=0._EB, Tdet=0._EB
REAL(EB) :: r_torso=0.15_EB, r_shoulder=0.095_EB, d_shoulder=0.055_EB, angle=0._EB, &
torque=0._EB, m_iner=4._EB
REAL(EB) :: tau_iner=0.2_EB, angle_old=0._EB, omega=0._EB
REAL(EB) :: A=2000._EB, B=0.08_EB, C_Young=120000._EB, Gamma=16000._EB, Kappa=40000._EB, &
Lambda=0.5_EB, Commitment=0._EB
REAL(EB) :: SumForces=0._EB, IntDose=0._EB, DoseCrit1=0._EB, DoseCrit2=0._EB, SumForces2=0._EB
REAL(EB) :: TempMax1=0._EB, FluxMax1=0._EB, TempMax2=0._EB, FluxMax2=0._EB, Density=0._EB, DensityR=0._EB, DensityL=0._EB
REAL(EB) :: P_detect_tot=0._EB, v0_fac=1._EB, D_Walls=0._EB
REAL(EB) :: T_FallenDown=0._EB, F_FallDown=0._EB, Angle_FallenDown=0._EB, SizeFac_FallenDown=0._EB, T_CheckFallDown=0._EB
INTEGER :: IOR=-1, ILABEL=0, COLOR_INDEX=0, INODE=0, IMESH=-1, IPC=0, IEL=0, I_FFIELD=0, I_Target2=0
INTEGER :: GROUP_ID=0, DETECT1=0, GROUP_SIZE=0, I_Target=0, I_DoorAlgo=0, I_Door_Mode=0, STRS_Direction = 1
INTEGER :: STR_SUB_INDX, SKIP_WALL_FORCE_IOR
LOGICAL :: SHOW=.TRUE., NewRnd=.TRUE.
LOGICAL :: SeeDoorXB1=.FALSE., SeeDoorXB2=.FALSE., SeeDoorXYZ1=.FALSE., SeeDoorXYZ2=.FALSE.
END TYPE HUMAN_TYPE
TYPE HUMAN_GRID_TYPE
! (x,y,z) Centers of the grid cells in the main evacuation meshes
! SOOT_DENS: Smoke density at the center of the cell (mg/m3)
! FED_CO_CO2_O2: Purser's FED for co, co2, and o2
REAL(EB) :: X,Y,Z,SOOT_DENS,FED_CO_CO2_O2,TMP_G,RADFLUX
INTEGER :: N, N_old, IGRID, IHUMAN, ILABEL
! IMESH: (x,y,z) which fire mesh, if any
! II,JJ,KK: Fire mesh cell reference
INTEGER :: IMESH,II,JJ,KK
END TYPE HUMAN_GRID_TYPE
TYPE SLICE_TYPE
INTEGER :: I1,I2,J1,J2,K1,K2,INDEX,INDEX2=0,Z_INDEX=-999,Y_INDEX=-999,PART_INDEX=0,VELO_INDEX=0
REAL(FB), DIMENSION(2) :: MINMAX
REAL(EB):: SLICE_AGL
LOGICAL :: TERRAIN_SLICE=.FALSE.,CELL_CENTERED=.FALSE., FIRE_LINE=.FALSE.,LEVEL_SET_FIRE_LINE=.FALSE.
CHARACTER(60) :: SMOKEVIEW_LABEL
CHARACTER(30) :: SMOKEVIEW_BAR_LABEL,ID='null'
END TYPE SLICE_TYPE
TYPE BOUNDARY_FILE_TYPE
INTEGER :: INDEX,PROP_INDEX,Z_INDEX=-999,Y_INDEX=-999,PART_INDEX=0
CHARACTER(60) :: SMOKEVIEW_LABEL
CHARACTER(30) :: SMOKEVIEW_BAR_LABEL
LOGICAL :: CELL_CENTERED=.FALSE.
END TYPE BOUNDARY_FILE_TYPE
TYPE (BOUNDARY_FILE_TYPE), DIMENSION(:), ALLOCATABLE, TARGET :: BOUNDARY_FILE
TYPE BOUNDARY_ELEMENT_FILE_TYPE
INTEGER :: INDEX,PROP_INDEX,Z_INDEX=-999,Y_INDEX=-999,PART_INDEX=0
CHARACTER(60) :: SMOKEVIEW_LABEL
CHARACTER(30) :: SMOKEVIEW_BAR_LABEL
LOGICAL :: CELL_CENTERED=.TRUE.
END TYPE BOUNDARY_ELEMENT_FILE_TYPE
TYPE (BOUNDARY_ELEMENT_FILE_TYPE), DIMENSION(:), ALLOCATABLE, TARGET :: BOUNDARY_ELEMENT_FILE
TYPE ISOSURFACE_FILE_TYPE
INTEGER :: INDEX=1,INDEX2=0,REDUCE_TRIANGLES=1,N_VALUES=1,Y_INDEX=-999,Y_INDEX2=-999,Z_INDEX=-999,Z_INDEX2=-999,&
VELO_INDEX=0
REAL(FB) :: VALUE(10)
CHARACTER(60) :: SMOKEVIEW_LABEL,SMOKEVIEW_LABEL2
CHARACTER(30) :: SMOKEVIEW_BAR_LABEL,SMOKEVIEW_BAR_LABEL2
END TYPE ISOSURFACE_FILE_TYPE
TYPE (ISOSURFACE_FILE_TYPE), DIMENSION(:), ALLOCATABLE, TARGET :: ISOSURFACE_FILE
TYPE PROFILE_TYPE
REAL(EB) :: X,Y,Z
INTEGER :: IOR,IW,ORDINAL,MESH
CHARACTER(30) :: ID,QUANTITY
END TYPE PROFILE_TYPE
TYPE (PROFILE_TYPE), DIMENSION(:), ALLOCATABLE, TARGET :: PROFILE
TYPE INITIALIZATION_TYPE
REAL(EB) :: TEMPERATURE,DENSITY,MASS_FRACTION(MAX_SPECIES)=0._EB,X1,X2,Y1,Y2,Z1,Z2,MASS_PER_VOLUME,MASS_PER_TIME,DT_INSERT, &
X0,Y0,Z0,U0,V0,W0,RADIUS,VOLUME,HRRPUV=0._EB
REAL(EB), POINTER, DIMENSION(:) :: PARTICLE_INSERT_CLOCK
INTEGER, POINTER, DIMENSION(:) :: WALL_INDEX_START
INTEGER :: PART_INDEX=0,NUMBER_INITIAL_PARTICLES,LU_PARTICLE,PROF_INDEX=0,DEVC_INDEX=0,CTRL_INDEX=0
LOGICAL :: ADJUST_DENSITY=.FALSE.,ADJUST_TEMPERATURE=.FALSE.,POINTWISE_PARTICLE_INIT=.FALSE., &
SINGLE_INSERTION=.TRUE.,ALREADY_INSERTED=.FALSE.,UNIFORM=.FALSE.
CHARACTER(30) :: SHAPE,PROF_ID,DEVC_ID,CTRL_ID
END TYPE INITIALIZATION_TYPE
TYPE (INITIALIZATION_TYPE), DIMENSION(:), ALLOCATABLE, TARGET :: INITIALIZATION
TYPE P_ZONE_TYPE
REAL(EB) :: X1,X2,Y1,Y2,Z1,Z2,DPSTAR=0._EB
REAL(EB), POINTER, DIMENSION(:) :: LEAK_AREA=>NULL()
INTEGER :: N_DUCTNODES
INTEGER, POINTER, DIMENSION(:) :: NODE_INDEX=>NULL()
CHARACTER(30) :: ID
END TYPE P_ZONE_TYPE
TYPE (P_ZONE_TYPE), DIMENSION(:), ALLOCATABLE, TARGET :: P_ZONE
TYPE MULTIPLIER_TYPE
REAL(EB) :: DXB(6),DX0,DY0,DZ0
INTEGER :: I_LOWER,I_UPPER,J_LOWER,J_UPPER,K_LOWER,K_UPPER,N_LOWER,N_UPPER,N_COPIES
CHARACTER(30) :: ID
LOGICAL :: SEQUENTIAL=.FALSE.
END TYPE MULTIPLIER_TYPE
TYPE (MULTIPLIER_TYPE), DIMENSION(:), ALLOCATABLE, TARGET :: MULTIPLIER
TYPE FILTER_TYPE
INTEGER :: RAMP_INDEX = -1
CHARACTER(30) :: ID,TABLE_ID
REAL(EB), DIMENSION(:), POINTER :: EFFICIENCY,MULTIPLIER
REAL(EB) :: CLEAN_LOSS,LOADING_LOSS
END TYPE FILTER_TYPE
TYPE (FILTER_TYPE), DIMENSION(:), ALLOCATABLE, TARGET :: FILTER
TYPE DUCTNODE_TYPE
INTEGER :: FILTER_INDEX=-1, N_DUCTS,VENT_INDEX = -1, MESH_INDEX = -1,ZONE_INDEX=-1
INTEGER, DIMENSION(:), POINTER :: DUCT_INDEX
CHARACTER(30) :: ID,TABLE_ID,VENT_ID='null'
REAL(EB), DIMENSION(:,:), POINTER :: LOSS_ARRAY, FILTER_LOADING
REAL(EB), DIMENSION(:), POINTER :: ZZ, DIR, BRANCH_LOSS, ZZ_V
REAL(EB) :: LOSS, P,TMP,RHO,RSUM,CP,XYZ(3),FILTER_LOSS,TMP_V,RHO_V,RSUM_V,CP_V
LOGICAL :: UPDATED, FIXED, AMBIENT = .FALSE.,LEAKAGE=.FALSE.
END TYPE DUCTNODE_TYPE
TYPE (DUCTNODE_TYPE), DIMENSION(:), ALLOCATABLE, TARGET :: DUCTNODE
TYPE NODE_BC_TYPE
REAL(EB), DIMENSION(:), POINTER :: ZZ_V
REAL(EB) :: TMP_V,RHO_V,RSUM_V,CP_V, P
END TYPE NODE_BC_TYPE
TYPE DUCT_TYPE
INTEGER :: NODE_INDEX(2)=-1,DEVC_INDEX=-1,CTRL_INDEX=-1,FAN_INDEX=-1,AIRCOIL_INDEX=-1,RAMP_INDEX=0
REAL(EB) :: ROUGHNESS,LENGTH, DIAMETER, AREA_INITIAL, AREA, CP_D, RHO_D,TMP_D,LOSS(2)=0._EB,VEL(4)=0._EB,RSUM_D=0._EB,&
DP_FAN=0._EB,TOTAL_LOSS = 0._EB,COIL_Q=0._EB,VOLUME_FLOW=1.E6_EB,VOLUME_FLOW_INITIAL=1.E6_EB,&
X(2),Y(2),Z(2),TAU=-1._EB
REAL(EB), DIMENSION(:), POINTER :: ZZ
LOGICAL :: ROUND = .TRUE.,SQUARE = .FALSE.,DAMPER = .FALSE.,DAMPER_OPEN = .TRUE.,FAN_OPERATING=.TRUE.,COIL_OPERATING=.TRUE.,&
FIXED=.FALSE.,REVERSE=.FALSE.,UPDATED,LEAKAGE=.FALSE.
CHARACTER(30) :: ID,RAMP_ID
REAL(EB) :: FAN_ON_TIME = 1.E10_EB
END TYPE DUCT_TYPE
TYPE (DUCT_TYPE), DIMENSION(:), ALLOCATABLE, TARGET :: DUCT
TYPE FAN_TYPE
INTEGER :: FAN_TYPE,RAMP_INDEX,SPIN_INDEX=0
REAL(EB) :: VOL_FLOW,MAX_FLOW,MAX_PRES,OFF_LOSS=0._EB,TAU
CHARACTER(30) :: ID,FAN_RAMP
END TYPE FAN_TYPE
TYPE(FAN_TYPE), DIMENSION(:), ALLOCATABLE, TARGET :: FAN
TYPE AIRCOIL_TYPE
REAL(EB) :: COOLANT_TEMPERATURE=293.15_EB,COOLANT_CP=4186._EB, EFFICIENCY, COOLANT_MDOT=-9999._EB,FIXED_Q=-1.E10_EB
CHARACTER(30) :: ID
END TYPE AIRCOIL_TYPE
TYPE(AIRCOIL_TYPE), DIMENSION(:), ALLOCATABLE, TARGET :: AIRCOIL
TYPE NETWORK_TYPE
INTEGER :: N_DUCTS,N_DUCTNODES,N_MATRIX
INTEGER, DIMENSION(:), POINTER :: DUCT_INDEX,NODE_INDEX,MATRIX_INDEX
END TYPE NETWORK_TYPE
TYPE(NETWORK_TYPE), DIMENSION(:), ALLOCATABLE, TARGET :: NETWORK
CONTAINS
SUBROUTINE GET_REV_type(MODULE_REV,MODULE_DATE)
INTEGER,INTENT(INOUT) :: MODULE_REV
CHARACTER(255),INTENT(INOUT) :: MODULE_DATE
WRITE(MODULE_DATE,'(A)') typerev(INDEX(typerev,':')+1:LEN_TRIM(typerev)-2)
READ (MODULE_DATE,'(I5)') MODULE_REV
WRITE(MODULE_DATE,'(A)') typedate
END SUBROUTINE GET_REV_type
END MODULE TYPES