Changeset 1849

Timestamp:

Dec 12, 2017, 11:28:55 AM (7 years ago)

Author:

jvatant

Message:

Fixed (?) some weird dark magic incoherence
between dynetat0 and newstart ...

File:

• trunk/LMDZ.TITAN/libf/dynphy_lonlat/phytitan/newstart.F (modified) (2 diffs)

trunk/LMDZ.TITAN/libf/dynphy_lonlat/phytitan/newstart.F

@@ -368 +368 @@
 c-----------------------------------------------------------------------
 
-      kappa = tab_cntrl(9) 
-      etot0 = tab_cntrl(12)
-      ptot0 = tab_cntrl(13)
-      ztot0 = tab_cntrl(14)
-      stot0 = tab_cntrl(15)
-      ang0 = tab_cntrl(16)
+      ! JVO 2017 :  For this part it seems to need a shift of 1 in the 
+      ! tab_cntrl indexes compared to dynetat0 ... I introduced it to have
+      ! correct outputs of newstart but I don't really understand why ...
+
+      kappa = tab_cntrl(10) 
+      etot0 = tab_cntrl(13)
+      ptot0 = tab_cntrl(14)
+      ztot0 = tab_cntrl(15)
+      stot0 = tab_cntrl(16)
+      ang0 = tab_cntrl(17)
       write(*,*) "Newstart: kappa,etot0,ptot0,ztot0,stot0,ang0"
       write(*,*) kappa,etot0,ptot0,ztot0,stot0,ang0
 
       ! for vertical coordinate
-      preff=tab_cntrl(18) ! reference surface pressure
-      pa=tab_cntrl(17)  ! reference pressure at which coord is purely pressure
+      preff=tab_cntrl(19) ! reference surface pressure
+      pa=tab_cntrl(18)  ! reference pressure at which coord is purely pressure
       !NB: in start_archive files tab_cntrl(17)=tab_cntrl(18)=0
       write(*,*) "Newstart: preff=",preff," pa=",pa
@@ -389 +393 @@
       if (yes.eq.'y') then
         write(*,*)"New value of reference surface pressure preff?"
-        write(*,*)"   (for Mars, typically preff=610)"
+        write(*,*)"   (for Titan, typically preff=140000)"
         read(*,*) preff
         write(*,*)"New value of reference pressure pa for purely"
         write(*,*)"pressure levels (for hybrid coordinates)?"
-        write(*,*)"   (for Mars, typically pa=20)"
+        write(*,*)"   (for Titan, typically pa=10000)"
         read(*,*) pa
       endif