Changeset 2243

Timestamp:

Feb 24, 2020, 3:36:36 PM (5 years ago)

Author:

jvatant

Message:

Get rid of the old 'sparadrap' in iniorbit.F assuming Mkm instead of AU if periastre or apostre gt 25, (otherwise Neptune is a pb)
Everybody should be using AU, that's all.
--JVO

Location:

trunk/LMDZ.GENERIC

Files:

• README (modified) (1 diff)
• libf/phystd/iniorbit.F (modified) (1 diff)

trunk/LMDZ.GENERIC/README

@@ -1500 +1500 @@
 - In thermcell_plume, restore initial formula to compute the vertical speed and cleanup the file
 
-== 29/01/2019 == AB
+== 29/01/2020 == AB
 - clean up thermcell_dq
 
-== 30/01/2019 == AB
+== 30/01/2020 == AB
 - fix a bug in thermcell_env. Now zqt is correctly initialized when tracer h2o_vap is missing (consistency with flag water is assumed).
 - clean up thermcell_dv2 and use plume height (zmin-zmax) instead of maximal altitude (zmax) in computations
 - the thermal plume model is able to manage several plumes in the same column and work without the convective adjustment.
+
+== 24/02/2020 (r2243) == JVO 
+- Get rid of the old 'sparadrap' in iniorbit.F assuming Mkm instead of AU if periastre or apostre gt 25, (otherwise Neptune is a pb)
+  Everybody should be using AU, that's all.

trunk/LMDZ.GENERIC/libf/phystd/iniorbit.F

@@ -32 +32 @@
       peri_day=pperi_day
 
-
-      !!!! SPARADRAP TEMPORAIRE !!!!
-      !!!! SPARADRAP TEMPORAIRE !!!!
-      !!!! We hope that all cases are above 25 Mkm [OK with Gliese 581d]
-      IF ( apoastr .gt. 25.) THEN
-        PRINT*,'!!!!! WARNING !!!!!'
-        PRINT*,'!!!!! YOU ARE ABOUT TO WITNESS A DIRT HACK !!!!!'
-        PRINT*,'This must be an old start file.'
-        PRINT*,'The code changed 19/03/2012: we now use AU.'
-        PRINT*,'So I am assuming units are in Mkm here'
-        PRINT*,'and I am performing a conversion towards AU.'
-        periastr = periastr / 149.598 ! Mkm to AU
-        apoastr = apoastr / 149.598 ! Mkm to AU
-      ENDIF
-      !!!! SPARADRAP TEMPORAIRE !!!!
-      !!!! SPARADRAP TEMPORAIRE !!!!
-
- 
       PRINT*,'iniorbit: Periastron in AU  ',periastr
       PRINT*,'iniorbit: Apoastron in AU  ',apoastr