Changeset 2464 for trunk/LMDZ.VENUS/libf/phyvenus/hrtherm.F

Timestamp:

Feb 22, 2021, 4:43:11 PM (4 years ago)

Author:

slebonnois

Message:

SL: major update related to extension to 250 km. New EUV heating as used in Martian GCM, debug of a few routines, adding He, clean up, updating mmean regularly, modification of the order of processes, add the possibility to use Hedin composition in rad transfer

File:

• trunk/LMDZ.VENUS/libf/phyvenus/hrtherm.F (modified) (5 diffs)

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

@@ -1 +1 @@
 c**********************************************************************
 
-      subroutine hrtherm(ig,euvmod,rm,nespeuv,tx,iz,zenit,zday,jtot)
+      subroutine hrtherm(ig,euvmod,rm,nespeuv,tx,iz,zenit,jtot)
 
 
@@ -9 +9 @@
 c**********************************************************************
       use dimphy
-      use conc
+      use param_v4_h, only: ninter,nabs,jfotsout,fluxtop,freccen
+
       implicit none
 
@@ -15 +16 @@
 
 
-#include "param.h"
-#include "param_v4.h"
 #include "clesphys.h"
 
@@ -38 +37 @@
       real       zenit
       real       iz(klev)
-      real       zday
 
       ! tracer indexes for the EUV heating:
@@ -95 +93 @@
 
       !Calculation of photoabsortion coefficient
-      if(solvarmod.eq.0) then
-         call jthermcalc(ig,euvmod,rm,nespeuv,tx,iz,zenit)
-      else if (solvarmod.eq.1) then
-         call jthermcalc_e107(ig,euvmod,rm,nespeuv,tx,iz,zenit,zday)
-         do indexint=1,ninter
-            fluxtop(indexint)=1.
-         enddo
-      endif
+      call jthermcalc_e107(ig,klev,euvmod,rm,nespeuv,tx,iz,zenit)
 
       !Total photoabsorption coefficient    !  erg/(s*cm3)