0000000000000000000000000000000000000000 071d675d69f3e65d78451735d397a3ecddb24677 Emeline Bolmont <emeline@arguin> 1395231574 +0100	branch: Created from tides
071d675d69f3e65d78451735d397a3ecddb24677 bfbc0f78c1c423f5c3a9cd4e4fa1868a21a24ad7 Emeline Bolmont <emeline@arguin> 1395233772 +0100	commit: I added two functions: sigma (as a function of K2Deltat) and the dissipation timescale
bfbc0f78c1c423f5c3a9cd4e4fa1868a21a24ad7 89409bf172265a6d6cd151c12f667451558479ac Emeline Bolmont <emeline@arguin> 1395235184 +0100	commit: I created horb scal spin to have cos obliquity
89409bf172265a6d6cd151c12f667451558479ac 09d864d059e917982d086ef788976d414d87e111 Emeline Bolmont <emeline@arguin> 1395238059 +0100	commit: I added functions, including the one to calculate Edot_tides
09d864d059e917982d086ef788976d414d87e111 9e5cac1a14a49a6b50b33f6c026b2b281fc2e9b2 Emeline Bolmont <emeline@arguin> 1395331696 +0100	commit: I just added jupiter=4 for this new model for planets
9e5cac1a14a49a6b50b33f6c026b2b281fc2e9b2 dbcd56f8fb9dcefb27dbd015697a7d471906e462 Emeline Bolmont <emeline@arguin> 1395332296 +0100	commit: I added if jupiter.eq.4 in user_module, but loop empty ahora
dbcd56f8fb9dcefb27dbd015697a7d471906e462 ebb0edb20cac81242ee2da65067bcb5ba32b89df Emeline Bolmont <emeline@arguin> 1395396691 +0100	commit: defined several k2deltat, I do not know if it will be necessary
ebb0edb20cac81242ee2da65067bcb5ba32b89df 11d4e4374122d3740297acc87e46c2b8b5a6d27a Emeline Bolmont <emeline@arguin> 1395409300 +0100	commit: I added Gabriels files in directory models_planet/
11d4e4374122d3740297acc87e46c2b8b5a6d27a b5e0ba3ba998642d8ddf018e3abee5b0125fdafe Emeline Bolmont <emeline@arguin> 1395410537 +0100	commit: I forgot to add something: files of profile depth, radius, density, viscosity etc..
b5e0ba3ba998642d8ddf018e3abee5b0125fdafe 26426d7364f03e1876f51c637e91a6c2b8588e95 Emeline Bolmont <emeline@arguin> 1395410679 +0100	commit: I removed the comment on the mysterious INTERFACE thing
26426d7364f03e1876f51c637e91a6c2b8588e95 0a7341d6ed02caf45c2f9ea3df73b074fdd5603c Emeline Bolmont <emeline@arguin> 1395413033 +0100	commit: I added commented lines to create a output file guving k2, Dt and Ediss
0a7341d6ed02caf45c2f9ea3df73b074fdd5603c 1ba2100fd3459ffbba07254f721e311c1f210baa Emeline Bolmont <emeline@arguin> 1395845203 +0100	commit: modified just a tiny bit Forced_tide.f90
1ba2100fd3459ffbba07254f721e311c1f210baa 0b16e4e52ec1f089d04a280203af47da00840a54 Emeline Bolmont <emeline@arguin> 1396004401 +0100	commit: Somehow there were no close(1) each time I opened a file open(1,...). I added them.
0b16e4e52ec1f089d04a280203af47da00840a54 e94831122ff9929eebc1275d16de20af607b4878 Emeline Bolmont <emeline@arguin> 1396018675 +0100	commit: added stuff to take into account jupiter=4, load the corresponding file
e94831122ff9929eebc1275d16de20af607b4878 da05dc0721b09f6899a63fd252fc0d769c905460 Emeline Bolmont <emeline@arguin> 1398173832 +0200	commit: I modified the format numbers written in tidesGR.out
da05dc0721b09f6899a63fd252fc0d769c905460 c9c6a907bca278c08b463a5980b4c0aa79d93772 Emeline Bolmont <emeline@arguin> 1398183740 +0200	commit: I wrote stuff to read the files with frequency, love number and delta_t for the planets
c9c6a907bca278c08b463a5980b4c0aa79d93772 bf7ff9962bfdb63b89b71bbc53a71a2e4141fde2 Emeline Bolmont <emeline@arguin> 1398261340 +0200	commit: I saved the data from these: freq, k2, delta_t
bf7ff9962bfdb63b89b71bbc53a71a2e4141fde2 e130f221d549994ae778e858fcc722091826da09 Emeline Bolmont <emeline@arguin> 1401972165 +0200	commit: I added a few lines to compute the tidal forcing frequency
e130f221d549994ae778e858fcc722091826da09 6bcfc73a43770108d49fdbb7f92e5fb452d8be52 Emeline Bolmont <emeline@arguin> 1401972267 +0200	commit: I commented it because I did not define the new variables, and it would not compile
6bcfc73a43770108d49fdbb7f92e5fb452d8be52 bbe16de489ea3d03ef1bc692859c4bcbef84df93 Emeline Bolmont <emeline@arguin> 1403185189 +0200	commit: I put the definition of the forcing frequency, I interpolated k2 and Dt to have the real one. The thing compiles, but I guess there are mistakes. I need to take care of the difference between k2_tides and k2_fluid next.
bbe16de489ea3d03ef1bc692859c4bcbef84df93 0f5702fbbb2386ed3d85680f6eca798cf21cd50d Emeline Bolmont <emeline@arguin> 1403773525 +0200	commit: I added the files of Gabriel to compute k2_fluid
0f5702fbbb2386ed3d85680f6eca798cf21cd50d 15e9dc2e827e67068beb773627fea7335295652a Emeline Bolmont <emeline@arguin> 1403773836 +0200	commit: I modified forced_tides.f90.
15e9dc2e827e67068beb773627fea7335295652a 5d32ef9142337a2269906d8edb0454e41aec4105 Emeline Bolmont <emeline@arguin> 1403773883 +0200	commit (amend): I modified forced_tides.f90.
5d32ef9142337a2269906d8edb0454e41aec4105 94e0e722cd7d55af0ec06dee8509f57ad076cbbb Emeline Bolmont <emeline@arguin> 1403800651 +0200	commit: I finished the initiation for jupiter=4, I did not put the interpolation within the code yet. I put k2_fluid for the stuff concerning rotation flatening. I did not do yet the writing of the files with k2_tides, Delta_t and E_dissipation, but I will. Now I will test if what I did works
94e0e722cd7d55af0ec06dee8509f57ad076cbbb 38c02ca590e0d93bf512857228be785e25733a46 Emeline Bolmont <emeline@arguin> 1403802369 +0200	commit: I removed access=append from where I read PLANETi.in
38c02ca590e0d93bf512857228be785e25733a46 4ca848d922aaef080899f211889f394a127ec131 Emeline Bolmont <emeline@arguin> 1403813767 +0200	commit: I changed some if conditions that were screwing things up, I put for all calculations of sigmap the use of the subroutine and I fixed the fact that in PLANETi.in, the radius was in meters, not AU
4ca848d922aaef080899f211889f394a127ec131 5712acec8ca886cda9f9fe1645819bf4c31c61ff Emeline Bolmont <emeline@arguin> 1403814021 +0200	commit: there was a typo
5712acec8ca886cda9f9fe1645819bf4c31c61ff 5833c56e28b5f8224a1930164566920813ab32b9 Emeline Bolmont <emeline@arguin> 1403814522 +0200	commit: n_orbit was in day^-1, so I put it in seconds
5833c56e28b5f8224a1930164566920813ab32b9 892d2813b1c9b7b8a1acfc1290964227398e6e26 Emeline Bolmont <emeline@arguin> 1403815324 +0200	commit: actually no, n_orbit has to be in day-1
892d2813b1c9b7b8a1acfc1290964227398e6e26 07826ad20057057884783f04448b60566c641ab8 Emeline Bolmont <emeline@arguin> 1403850767 +0200	commit: I moved the thing to calculate sigmap, I was trying to do it before initializing the spin. I have to think about the frequency with which I will call that. But first I will do some tests
07826ad20057057884783f04448b60566c641ab8 a31d5e806ce9a2027d0a68b000e5b43053cecb81 Emeline Bolmont <emeline@arguin> 1403851541 +0200	commit: in file PLANETi.in, freq is in s-1, but forcing_frequency in user_module is in day-1. I fixed that in interpolation
a31d5e806ce9a2027d0a68b000e5b43053cecb81 01e615e78b8d7b86276832de073c3cfaadc6b3f9 Emeline Bolmont <emeline@arguin> 1403852028 +0200	commit: I forgot to save nptk2 which is the number of freq in fil PLANETi.in
01e615e78b8d7b86276832de073c3cfaadc6b3f9 fb3378b09d75656a8a80d2183b6567e880220a64 Emeline Bolmont <emeline@arguin> 1404141315 +0200	commit: I just put a better version of makespin.sh
fb3378b09d75656a8a80d2183b6567e880220a64 1199ef4e7883d9e8cc0bf20a15e215646ea5b861 Emeline Bolmont <emeline@arguin> 1404286748 +0200	commit: I added the thing that if forcing-frequency is outside the range, it assumes the min or max value of table of freq and I added forcing_freq in the output
1199ef4e7883d9e8cc0bf20a15e215646ea5b861 80e7d86fb81ba5859982b21fbb22ec85a94102e8 Emeline Bolmont <emeline@arguin> 1404288213 +0200	commit: I modified Forced_tide_new to have freq in day-1 and a range between freq=1e-10 and 1e-4 s-1, to avoid the resonances. I also modified user_module to take into account this
80e7d86fb81ba5859982b21fbb22ec85a94102e8 e002428dbd9a5a69b4ee37be653f53c5998de2b0 Emeline Bolmont <emeline@arguin> 1404288265 +0200	commit: I modified the example files PLANET1.in and PLANET2.in
e002428dbd9a5a69b4ee37be653f53c5998de2b0 752c05ee2bb0aa92a53f84a9247ce073adda40cf Emeline Bolmont <emeline@arguin> 1404288754 +0200	commit: I added something so that the dissipation stuff is calculated very ten timestep
752c05ee2bb0aa92a53f84a9247ce073adda40cf 36dd0c23990725f254e58097cff0ade729160403 Emeline Bolmont <emeline@arguin> 1404288941 +0200	commit: I put the real length of files PLANETi.in
36dd0c23990725f254e58097cff0ade729160403 7e275cc68cd3e69ec4dab55fe2aa7ecf34e7459d Emeline Bolmont <emeline@arguin> 1404290105 +0200	commit: I modified the ploting routines to incorporate Ediss and the files coef_tides
7e275cc68cd3e69ec4dab55fe2aa7ecf34e7459d b29c057a856bd648c71981467abc015ed22c2339 Emeline Bolmont <emeline@arguin> 1404295951 +0200	commit: I was using the same file for all planets... I fixed that
b29c057a856bd648c71981467abc015ed22c2339 66c9dd93c048f8ba87abe8112505246f7173d5af Emeline Bolmont <emeline@arguin> 1404309948 +0200	commit: it was not doing well the freq when it was outside the range
66c9dd93c048f8ba87abe8112505246f7173d5af b5c19707f0474588c2fe6ad0927e387db943b3bf Emeline Bolmont <emeline@arguin> 1404314746 +0200	commit: I modified slightly the script_plot thing
b5c19707f0474588c2fe6ad0927e387db943b3bf 92233d6aa389dd4273db05bfbd8a9d022ed4f323 Emeline Bolmont <emeline@arguin> 1404314797 +0200	commit: I added a file to plot forcing_frequency, k2
92233d6aa389dd4273db05bfbd8a9d022ed4f323 654114151eb0fa4cf365bd1578715eb4d37dc430 Emeline Bolmont <emeline@arguin> 1409736967 +0200	commit: I changed a line of place in tides_constant_GR
654114151eb0fa4cf365bd1578715eb4d37dc430 d2d7be94d248d05eee0a358ff75c18e4bd9f00c7 Emeline Bolmont <emeline@arguin> 1413463344 +0200	commit: I added a README file in the models_planet directory to be able to remember how to use the stuff there
d2d7be94d248d05eee0a358ff75c18e4bd9f00c7 536b2f001bdc35ad534d0fa1570a7e8ebcbb3a0d Emeline Bolmont <emeline@arguin> 1413468940 +0200	commit: There was a mistake in the expression of Ediss, I was writing cos(obliquity)=h.scal.spin, but I didn't normalize it, so it was wrong...
536b2f001bdc35ad534d0fa1570a7e8ebcbb3a0d 9fadab6cf8c0edbae6aeab0dfafb7b682e7b0e7c Emeline Bolmont <emeline@arguin> 1413469051 +0200	commit: And I put a 2 for the forcing frequency, coz apparently there is a 2. Francoise Remus told me that
9fadab6cf8c0edbae6aeab0dfafb7b682e7b0e7c ab54c23ac1a658e1786b82229a77b94c6676b899 Emeline Bolmont <emeline@arguin> 1413477023 +0200	commit: I was using an integer to read files bigger than 1177, and integer does not go that far, so I changed the numbers to read files, Jupiter, BD and all and dissipation factor
ab54c23ac1a658e1786b82229a77b94c6676b899 2ae44b461676323a9c09c5982411bdeb636021ce Emeline Bolmont <emeline@arguin> 1413477583 +0200	commit: nptmss should be integer(kind=8), long integer
