Update solid solver for flexible wings and SuperLU library

Makefile

@@ -121,7 +121,7 @@ ifdef NOHDF5
 PRE_FLAGS += -DNOHDF5
 endif
 FFLAGS += -module $(OBJDIR) # specify directory for modules.
-FFLAGS += -vec_report0
+#FFLAGS += -vec_report0
 # debug flags
 # FFLAGS+= -g -debug all -traceback -check all -CB
 endif

src/flexible_solver/save_wing.f90

@@ -22,8 +22,6 @@ subroutine SaveWingData( time, wings )
       write(ns1_string,'(i3)') np+1
       format_ns1 = '('//ns1_string//'(es15.8,1x))'
 
-      write(*,*) ns1_string, format_ns1
-
       !-- save control point data
       if (wings(i)%ControlPoint /= 0) then
       open  (14, file = 'wing_data'//wingstr//'_point'//ctrlpoint//'.t', status = 'unknown',position='append')

src/insects/wings_geometry.f90

@@ -360,7 +360,7 @@ subroutine draw_blade_fourier(xx0, ddx, mask, mask_color, us,Insect,color_wing,M
                           xte = Radius_Fourier(1.0d0-theta,Insect,wingID)
 
                           !-- amplitude
-                          R0 = xle-xte
+                          R0 = 0.5*(xle-xte)
 
                           !-- get smooth rectangular function
                           R = dabs ( x_wing(1) - 0.5*(xte+xle) )
@@ -1002,12 +1002,12 @@ subroutine draw_wing_bristled(xx0, ddx, mask, mask_color, us,Insect,color_wing,M
               -0.159718,-0.153299,-0.14512,-0.138701,-0.130518,-0.11999/)
   ! Bristle roots
   ybr_root = (/0.245961,0.275969,0.307151,0.553556,0.572956,0.59412,0.614697,0.634098,0.651735,0.66937,0.682306, &
-               0.694063,0.70347,0.711118,0.711121,0.70701,0.697018,0.687025,0.675857,0.664099,0.651754,0.639998, &
+               0.694063,0.70250,0.711118,0.711121,0.70701,0.697018,0.687025,0.675857,0.664099,0.651754,0.639998, &
                0.627068,0.617075,0.604732,0.592388,0.582396,0.568287,0.555942,0.541834,0.528309,0.514787,0.439526, &
                0.426,0.414239,0.397773,0.384247,0.369543,0.355426,0.338368,0.318955,0.296013,0.271301,0.247763, &
                0.225406,0.20599,0.190102,0.17539/)
   xbr_root = (/-0.067707,-0.07412,-0.078773,-0.02397,-0.014569,-0.003993,0.00541,0.014811,0.023624,0.03361,0.038314, &
-               0.04419,0.048304,0.047726,0.045382,0.039515,0.03247,0.026596,0.019549,0.013674,0.007212,0.000164, &
+               0.04419,0.046000,0.047726,0.045382,0.039515,0.03247,0.026596,0.019549,0.013674,0.007212,0.000164, &
                -0.00923,-0.01569,-0.02391,-0.031545,-0.038591,-0.045641,-0.052103,-0.060326,-0.063859,-0.070323, &
                -0.097961,-0.100907,-0.103852,-0.107975,-0.110921,-0.112111,-0.1133,-0.113906,-0.112756,-0.111611, &
                -0.106364,-0.099358,-0.095867,-0.091787,-0.087116,-0.081858/)

src/insects/wings_motion.f90

@@ -315,6 +315,26 @@ subroutine FlappingMotion(time, Insect, protocoll, phi, alpha, theta, phi_dt, &
     theta = 0.d0
     theta_dt = 0.d0
 
+  case ("revolving-set3")
+    ! revolving wing kinematics, pre-defined set. We fix alpha and increase
+    ! phi linearily with a quadratic startup transient until phi=pi/16,
+    ! which is reached at time equal to pi/8
+    ! we use PHI_DOT = 1 as normalization
+    ! Therefore, the acceleration is twice as fast as in "revolving-set2"
+    if (time < pi/8.0d0) then
+      phi = 4.0d0*time**2/pi
+      phi_dt = 8.0d0*time/pi
+    else
+      phi = time - pi/16.0d0
+      phi_dt = 1.0d0	
+    endif
+    ! feathering angle is constant
+    alpha = - Insect%init_alpha_phi_theta(1) ! Mind the "-" sign
+    alpha_dt = 0.d0
+    ! elevation angle is always zero
+    theta = 0.d0
+    theta_dt = 0.d0
+
   case ("Drosophila_hovering_fry")
     !---------------------------------------------------------------------------
     ! motion protocoll digitalized from Fry et al JEB 208, 2303-2318 (2005)