105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127 use, INTRINSIC :: iso_c_binding, only: c_bool
131
132
133
134#include "implicit_f.inc"
135
136
137
138
139
140
141
142 TYPE (UNIT_TYPE_),INTENT(IN) :: UNITAB
143 INTEGER,INTENT(IN) :: INDEX
144 CHARACTER*(*),INTENT(IN) :: NAME
145 TYPE(SUBMODEL_DATA),INTENT(IN) :: LSUBMODEL(NSUBMOD)
146
148 LOGICAL,INTENT(OUT) :: IS_AVAILABLE
149
150
151
152 INTEGER J,SUB_ID,IFLAGUNIT,UID
153 my_real :: fac_l,fac_m,fac_t,fac
154 real*8 dval,length_dim,mass_dim,time_dim
155 LOGICAL(KIND=C_BOOL) C_IS_AVAILABLE
156
157 c_is_available
158 length_dim
159 mass_dim = zero
160 time_dim = zero
161 fac = one
162
163 CALL cpp_get_floatv_floatd_index(name(1:len_trim(name)),len_trim(name),dval,index,c_is_available,
164 . length_dim,mass_dim,time_dim,uid,sub_id)
165 is_available = c_is_available
166
167
168
169 IF(sub_id /= 0 .AND. uid == 0)THEN
170 IF(lsubmodel(sub_id)%UID /= 0)THEN
171 uid = lsubmodel(sub_id)%UID
172 ENDIF
173 ENDIF
174
175
176
177 iflagunit = 0
178 fac_m = zero
179 fac_l = zero
180 fac_t = zero
181 DO j=1,unitab%NUNITS
182 IF (unitab%UNIT_ID(j) == uid) THEN
183 fac_m = unitab%FAC_M(j)
184 fac_l = unitab%FAC_L(j)
185 fac_t = unitab%FAC_T(j)
186 iflagunit = 1
187 EXIT
188 ENDIF
189 ENDDO
190 IF (fac_m /= zero) fac = fac * (fac_m ** mass_dim )
191 IF (fac_l /= zero) fac = fac * (fac_l ** length_dim)
192 IF (fac_t /= zero) fac = fac * (fac_t ** time_dim )
193
194 rval = dval * fac
195
196 RETURN
197