418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439 use, INTRINSIC :: iso_c_binding, only: c_bool
443
444
445
446#include "implicit_f.inc"
447
448
449
450
451
452
453
454 TYPE (UNIT_TYPE_),INTENT(IN) :: UNITAB
455 CHARACTER*(*),INTENT(IN) :: NAME
456 TYPE(SUBMODEL_DATA),INTENT(IN) :: LSUBMODEL(NSUBMOD)
457
459 LOGICAL,INTENT(OUT) :: IS_AVAILABLE
460
461
462
463 INTEGER :: J,SUB_ID,IFLAGUNIT,UID
464 my_real :: fac_l,fac_m,fac_t,fac
465 real*8 :: dval,length_dim,mass_dim,time_dim
466 LOGICAL(KIND=C_BOOL) :: C_IS_AVAILABLE
467
468 c_is_available = .false.
469 length_dim = zero
470 mass_dim = zero
471 time_dim = zero
472 fac = one
473
474 CALL cpp_get_floatv_floatd(name(1:len_trim(name)),len_trim(name),dval,c_is_available,
475 . length_dim,mass_dim,time_dim,uid,sub_id)
476 is_available = c_is_available
477
478
479
480 IF(sub_id /= 0 .AND. uid == 0)THEN
481 IF(lsubmodel(sub_id)%UID /= 0)THEN
482 uid = lsubmodel(sub_id)%UID
483 ENDIF
484 ENDIF
485
486
487
488 iflagunit = 0
489 fac_m = zero
490 fac_l = zero
491 fac_t = zero
492 DO j=1,unitab%NUNITS
493 IF (unitab%UNIT_ID(j) == uid) THEN
494 fac_m = unitab%FAC_M(j)
495 fac_l = unitab%FAC_L(j)
496 fac_t = unitab%FAC_T(j)
497 iflagunit = 1
498 EXIT
499 ENDIF
500 ENDDO
501 IF (fac_m /= zero) fac = fac * (fac_m ** mass_dim )
502 IF (fac_l /= zero) fac = fac * (fac_l ** length_dim)
503 IF (fac_t /= zero) fac = fac * (fac_t ** time_dim )
504
505 rval = dval * fac
506
507
508 RETURN
509