Changeset 1379 for trunk

Timestamp:

Feb 20, 2015, 11:59:50 AM (10 years ago)

Author:

slebonnois

Message:

SL: correction of a bug in Titan's chemistry

File:

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

trunk/LMDZ.TITAN/libf/phytitan/calchim.F

@@ -132 +132 @@
          print*,l,rmil(l)-RA/1000.,press_c(l),nb(l),temp_c(l)
         ENDDO
+        rinter(klev+1)=(zlevmoy(klev+1)+RA)/1000.
 
 c au-dessus du GCM, dr regulier et rinter(NLEV)=1290+2575 km.
-       do l=klev+1,NLEV
-         rinter(l) = rinter(klev)
-     &          + (l-klev)*(3865.-rinter(klev))/(NLEV-klev)
-         rmil(l)   = rmil(klev)
-     &          + (l-klev)*(3865.-rinter(klev))/(NLEV-klev)
-       enddo        
+       do l=klev+2,NLEV
+         rinter(l) = rinter(klev+1)
+     &          + (l-klev-1)*(3865.-rinter(klev+1))/(NLEV-klev-1)
+         rmil(l-1) = (rinter(l-1)+rinter(l))/2.
+       enddo
+       rmil(NLEV) = rinter(NLEV)+(rinter(NLEV)-rinter(NLEV-1))/2.
 
 c lecture de tcp.ver, une seule fois
@@ -358 +359 @@
 c        print*,rinter(l)
        enddo
+       rinter(klev+1)=(RA+czlev(j,klev+1))/1000.
 
 c au-dessus du GCM, dr regulier et rinter(NLEV)=1290+2575 km.
-       do l=klev+1,NLEV
-         rinter(l) = rinter(klev)
-     &          + (l-klev)*(3865.-rinter(klev))/(NLEV-klev)
-         rmil(l)   = rmil(klev)
-     &          + (l-klev)*(3865.-rinter(klev))/(NLEV-klev)
-       enddo
-        
+       do l=klev+2,NLEV
+         rinter(l) = rinter(klev+1)
+     &          + (l-klev-1)*(3865.-rinter(klev+1))/(NLEV-klev-1)
+         rmil(l-1) = (rinter(l-1)+rinter(l))/2.
+       enddo
+       rmil(NLEV) = rinter(NLEV)+(rinter(NLEV)-rinter(NLEV-1))/2.
+
 c-----------------------------------------------------------------------
 c