Attempt to write out coarsened convec fields

Gradient test with coarsened mass fluxes:

     alpha               J_bg              J_obs              J_tot                DJ1                DJ2         DJ2/DJ1       1-DJ2/DJ1
         0  4.09800000000e+03  1.62243225768e+03  5.72043225768e+03  3.34908136381e+04  0.00000000000e+00               0  1.00000000e+00
         1  1.25921665110e+04  2.11797385352e+06  2.13056602003e+06  3.34908136381e+04 -2.12484558777e+06      -63.445625  6.44456245e+01
       0.1  3.44033003739e+03  2.05143848727e+04  2.39547149101e+04  3.34908136381e+04 -1.82342826524e+05      -5.4445625  6.44456246e+00
      0.01  4.01716213760e+03  1.58419562407e+03  5.60135776167e+03  3.34908136381e+04  1.19074496006e+04      0.35554375  6.44456246e-01
     0.001  4.08976550510e+03  1.59933427536e+03  5.68909978046e+03  3.34908136381e+04  3.13324772131e+04      0.93555437  6.44456252e-02
    0.0001  4.09717504342e+03  1.61992971625e+03  5.71710475968e+03  3.34908136381e+04  3.32749799744e+04      0.99355544  6.44456316e-03
     1e-05  4.09791748927e+03  1.62218007610e+03  5.72009756537e+03  3.34908136381e+04  3.34692302501e+04      0.99935554  6.44456961e-04
     1e-06  4.09799174878e+03  1.62240702024e+03  5.72039876902e+03  3.34908136381e+04  3.34886552801e+04      0.99993555  6.44462705e-05
     1e-07  4.09799917488e+03  1.62242973374e+03  5.72042890862e+03  3.34908136381e+04  3.34905977434e+04      0.99999355  6.44638630e-06
     1e-08  4.09799991749e+03  1.62243200528e+03  5.72043192277e+03  3.34908136381e+04  3.34907919751e+04      0.99999935  6.46835094e-07
     1e-09  4.09799999175e+03  1.62243223244e+03  5.72043222418e+03  3.34908136381e+04  3.34908127115e+04      0.99999997  2.76659264e-08
     1e-10  4.09799999917e+03  1.62243225515e+03  5.72043225433e+03  3.34908136381e+04  3.34907781507e+04      0.99999894  1.05961454e-06

Gradient test with original fluxes:

     alpha               J_bg              J_obs              J_tot                DJ1                DJ2         DJ2/DJ1       1-DJ2/DJ1
         0  4.09800000000e+03  1.62243226210e+03  5.72043226210e+03  3.34908154610e+04  0.00000000000e+00               0  1.00000000e+00
         1  1.25921674144e+04  2.11797416251e+06  2.13056632992e+06  3.34908154610e+04 -2.12484589766e+06       -63.44563  6.44456303e+01
       0.1  3.44033004570e+03  2.05143878036e+04  2.39547178493e+04  3.34908154610e+04 -1.82342855872e+05       -5.444563  6.44456304e+00
      0.01  4.01716213761e+03  1.58419564143e+03  5.60135777904e+03  3.34908154610e+04  1.19074483064e+04       0.3555437  6.44456304e-01
     0.001  4.08976550509e+03  1.59933427829e+03  5.68909978338e+03  3.34908154610e+04  3.13324787243e+04      0.93555437  6.44456310e-02
    0.0001  4.09717504342e+03  1.61992972050e+03  5.71710476393e+03  3.34908154610e+04  3.32749817661e+04      0.99355544  6.44456374e-03
     1e-05  4.09791748927e+03  1.62218008051e+03  5.72009756978e+03  3.34908154610e+04  3.34692320697e+04      0.99935554  6.44457023e-04
     1e-06  4.09799174878e+03  1.62240702467e+03  5.72039877345e+03  3.34908154610e+04  3.34886571027e+04      0.99993555  6.44462743e-05
     1e-07  4.09799917488e+03  1.62242973817e+03  5.72042891304e+03  3.34908154610e+04  3.34905995987e+04      0.99999355  6.44541558e-06
     1e-08  4.09799991749e+03  1.62243200971e+03  5.72043192720e+03  3.34908154610e+04  3.34907937940e+04      0.99999935  6.46950952e-07
     1e-09  4.09799999175e+03  1.62243223686e+03  5.72043222861e+03  3.34908154610e+04  3.34908099830e+04      0.99999984  1.63564523e-07
     1e-10  4.09799999917e+03  1.62243225958e+03  5.72043225875e+03  3.34908154610e+04  3.34907872457e+04      0.99999916  8.42478047e-07

They are close, but different. I can only imagine that this is due to the mass fluxes being float32.

Now I will eliminate the block

if ( t2 > t_max ) then
   ! default t1,tmid,t2 would be : 21,00,03
   ! reset to : 21,00,00
   t2 = t_max
end if

from read_wind and see if the gradient test stll works.

     alpha               J_bg              J_obs              J_tot                DJ1                DJ2         DJ2/DJ1       1-DJ2/DJ1
         0  4.09800000000e+03  1.62243250586e+03  5.72043250586e+03  3.34908136342e+04  0.00000000000e+00               0  1.00000000e+00
         1  1.25921664977e+04  2.11797385498e+06  2.13056602148e+06  3.34908136342e+04 -2.12484558898e+06      -63.445625  6.44456246e+01
       0.1  3.44033003623e+03  2.05143851344e+04  2.39547151707e+04  3.34908136342e+04 -1.82342826648e+05      -5.4445625  6.44456246e+00
      0.01  4.01716213749e+03  1.58419587253e+03  5.60135801002e+03  3.34908136342e+04  1.19074495847e+04      0.35554375  6.44456247e-01
     0.001  4.08976550509e+03  1.59933452357e+03  5.68910002866e+03  3.34908136342e+04  3.13324772080e+04      0.93555437  6.44456253e-02
    0.0001  4.09717504342e+03  1.61992996444e+03  5.71710500787e+03  3.34908136342e+04  3.32749799703e+04      0.99355544  6.44456316e-03
     1e-05  4.09791748927e+03  1.62218032429e+03  5.72009781356e+03  3.34908136342e+04  3.34692302461e+04      0.99935554  6.44456964e-04
     1e-06  4.09799174878e+03  1.62240726843e+03  5.72039901721e+03  3.34908136342e+04  3.34886552764e+04      0.99993555  6.44462624e-05
     1e-07  4.09799917488e+03  1.62242998193e+03  5.72042915680e+03  3.34908136342e+04  3.34905977525e+04      0.99999355  6.44599798e-06
     1e-08  4.09799991749e+03  1.62243225347e+03  5.72043217096e+03  3.34908136342e+04  3.34907920660e+04      0.99999936  6.44002683e-07
     1e-09  4.09799999175e+03  1.62243248062e+03  5.72043247237e+03  3.34908136342e+04  3.34908118020e+04      0.99999995  5.47057122e-08
     1e-10  4.09799999917e+03  1.62243250334e+03  5.72043250251e+03  3.34908136342e+04  3.34907781507e+04      0.99999894  1.05949778e-06

It does! But... is the mass now continuous?