Changeset 1442 for trunk/LMDZ.VENUS/libf/phyvenus/nlte_tcool.F

Timestamp:

Jun 4, 2015, 4:23:32 PM (10 years ago)

Author:

slebonnois

Message:

SL: update of the Venus GCM, + corrections on routines used for newstart/start2archive for Titan and Venus, + some modifications on tools

File:

• trunk/LMDZ.VENUS/libf/phyvenus/nlte_tcool.F (modified) (6 diffs)

trunk/LMDZ.VENUS/libf/phyvenus/nlte_tcool.F

@@ -33 +33 @@
       implicit none
 
-#include "dimensions.h"
-#include "nlte_paramdef.h"
-#include "nlte_commons.h"
-#include "YOMCST.h"
+      include "dimensions.h"
+      include "nlte_paramdef.h"
+      include "nlte_commons.h"
+      include "YOMCST.h"
 
 c     Arguments
@@ -48 +48 @@
       real q15umco2_gcm(nlon,nlev) ! is in K/RDAY (see hrkday_convert) 
                                    ! but converted to K/s (see CONVERSION_KDAY_Ksec )
-!     real auxgcm(nlev)
       real*8 auxgcmd(nlev), aux2gcmd(nlev)
       real zmin_gcm
       integer ierr
       real*8 varerr
-
-
 
 c     local variables and constants
@@ -66 +63 @@
       real co2_ig(nlev),n2_ig(nlev),co_ig(nlev),o3p_ig(nlev)
       real mmean_ig(nlev),cpnew_ig(nlev)
-!!!!!
-c      real cpnew(nlon,nlev)
-c      real rnew(nlon,nlev)
-c      real mmean(nlon,nlev)
-
 
 
@@ -90 +82 @@
             mmean_ig(l)=mmean(ig,l)
             cpnew_ig(l)=cpnew(ig,l)
+
          enddo 
 
@@ -318 +311 @@
 !      print*, zl
 
-                    ! Creamos el perfil del NLTE modelo completo interpolando
-
-      call interhunt (    pl,zl,nl,      p_gcm,z_gcm,nlev, 2) ! [atm]
+! Creamos el perfil del NLTE modelo completo interpolando
+
+      call interhunt (    pl,zl,nl,  p_gcm,z_gcm,nlev, 2) ! [atm]
       call interhunt5veces
      $     ( t, co2vmr, n2vmr, covmr, o3pvmr,
@@ -354 +347 @@
             write (*,*) ' i, t(i), pl(i) =', i, t(i), pl(i)
          endif
-         if (t(i) .lt. 20.0) then
+         if (t(i) .lt. 50.0) then
             write (*,*) '!!!! WARNING    Temp lower than Histogram.'
             write (*,*) ' Histogram will be extrapolated. '