MoDeNa  1.0
Software framework facilitating sequential multi-scale modelling
rheologyexactdummy.f90
1 @ingroup mod_src_rheology
2 
3 program rheologyexactdummy
4 ! use iso_c_binding
5 ! use fmodena
6 ! use tfem_m
7  implicit none
8  real(8) :: temp, shear, conv
9  real(8) :: mu, surfacetension, polymerviscosity
10  integer :: ret
11 
12 ! integer(c_size_t) :: surfaceTempPos
13 ! integer(c_size_t) :: viscosityTempPos, viscosityConvPos
14 ! type(c_ptr) :: surfaceModel, surfaceInputs, surfaceOutputs
15 ! type(c_ptr) :: viscosityModel, viscosityInputs, viscosityOutputs
16 
17 ! surfaceModel = modena_model_new (c_char_"SurfaceTension"//c_null_char);
18 ! surfaceInputs = modena_inputs_new (surfaceModel);
19 ! surfaceOutputs = modena_outputs_new (surfaceModel);
20 ! surfaceTempPos = modena_model_inputs_argPos(surfaceModel, &
21 ! c_char_"T"//c_null_char);
22 !
23 ! viscosityModel = modena_model_new (c_char_"polymerViscosity"//c_null_char);
24 ! viscosityInputs = modena_inputs_new (viscosityModel);
25 ! viscosityOutputs = modena_outputs_new (viscosityModel);
26 ! viscosityTempPos = modena_model_inputs_argPos(viscosityModel, &
27 ! c_char_"T"//c_null_char);
28 ! viscosityConvPos = modena_model_inputs_argPos(viscosityModel, &
29 ! c_char_"X"//c_null_char);
30 
31 ! call modena_model_argPos_check(surfaceModel)
32 
33 ! call modena_model_argPos_check(viscosityModel)
34 
35 
36  open(14, file='RheologyExact.in')
37  read(14, * ) temp, shear, conv
38  close(14)
39 
40 ! call modena_inputs_set( surfaceInputs, surfaceTempPos, temp )
41 
42 ! call modena_inputs_set( viscosityInputs, viscosityTempPos, temp )
43 ! call modena_inputs_set( viscosityInputs, viscosityConvPos, conv )
44 !
45 ! ret = modena_model_call(surfaceModel, surfaceInputs, surfaceOutputs)
46 ! if ( ret /= 0 ) then
47 ! call modena_inputs_destroy( surfaceInputs )
48 ! call modena_outputs_destroy( surfaceOutputs )
49 ! call modena_model_destroy( surfaceModel )
50 ! print*,ret
51 ! call exit(ret)
52 ! end if
53 
54 ! ret = modena_model_call(viscosityModel, viscosityInputs, viscosityOutputs)
55 ! if ( ret /= 0 ) then
56 ! call modena_inputs_destroy( viscosityInputs )
57 ! call modena_outputs_destroy( viscosityOutputs )
58 ! call modena_model_destroy( viscosityModel )
59 ! print*,ret
60 ! call exit(ret)
61 ! end if
62 
63 ! surfaceTension = modena_outputs_get(surfaceOutputs, 0_c_size_t);
64 ! polymerViscosity = modena_outputs_get(viscosityOutputs, 0_c_size_t);
65 
66 
67  mu = 1.0
68  open(15, file='RheologyExact.out')
69  write(15,*) mu
70  close(15)
71 
72 
73 end program rheologyexactdummy