159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
183 use, INTRINSIC :: iso_c_binding, only: c_bool
184
185
186
187
188#include "implicit_f.inc"
189
190
191
192
193
194
195
196 TYPE (UNIT_TYPE_),INTENT(IN) ::UNITAB
197 CHARACTER*(*),INTENT(IN)::NAME
198 TYPE(SUBMODEL_DATA),INTENT(IN)::LSUBMODEL(*)
199
201 . INTENT(OUT)::dim_fac
202 LOGICAL,INTENT(OUT)::IS_AVAILABLE
203
204
205
206 INTEGER I,J,VALUE_TYPE,SUB_ID,IFLAGUNIT,UID
208 . fac_l,fac_m,fac_t,fac
209 real*8 dval,length_dim,mass_dim,time_dim
210 LOGICAL(KIND=C_BOOL) :: C_IS_AVAILABLE
211
212
213 c_is_available = .false.
214 length_dim = zero
215 mass_dim = zero
216 time_dim = zero
217 fac = one
218
219 CALL cpp_get_floatv_floatd(name(1:len_trim(name)),len_trim(name),dval,c_is_available,
220 . length_dim,mass_dim,time_dim,uid,sub_id)
221 is_available = c_is_available
222
223
224
225 IF(sub_id /= 0 .AND. uid == 0)THEN
226 IF(lsubmodel(sub_id)%UID /= 0)THEN
227 uid = lsubmodel(sub_id)%UID
228 ENDIF
229 ENDIF
230
231
232
233 iflagunit = 0
234 fac_m = zero
235 fac_l = zero
236 fac_t = zero
237 DO j=1,unitab%NUNITS
238 IF (unitab%UNIT_ID(j) == uid) THEN
239 fac_m = unitab%FAC_M(j)
240 fac_l = unitab%FAC_L(j)
241 fac_t = unitab%FAC_T(j)
242 iflagunit = 1
243 EXIT
244 ENDIF
245 ENDDO
246 IF (fac_m /= zero) fac = fac * (fac_m ** mass_dim )
247 IF (fac_l /= zero) fac = fac * (fac_l ** length_dim)
248 IF (fac_t /= zero) fac = fac * (fac_t ** time_dim )
249
250 dim_fac = fac
251
252
253 RETURN
254