Changeset 2429

Timestamp:

Nov 5, 2020, 4:28:35 PM (4 years ago)

Author:

emillour

Message:

Mars GCM:
Bug fix in jthermcalc_e107.F: correctly account for the actual Sun-Mars
distance.
FGG

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/aeronomars/jthermcalc_e107.F (modified) (3 diffs)

trunk/LMDZ.MARS/README

@@ -3193 +3193 @@
 to avoid unwanted temporal interpolation in read_dust_scenario, the time
 axis there should simply contain integers.
+
+== 05/11/2020 == FGG
+Bug fix in jthermcalc_e107.F: correctly account for the actual Sun-Mars
+distance.
+

trunk/LMDZ.MARS/libf/aeronomars/jthermcalc_e107.F

@@ -18 +18 @@
      .    nabs,e107,date_e107,e107_tab,
      .    coefit0,coefit1,coefit2,coefit3,coefit4
+      use comsaison_h, only: dist_sol
 
       implicit none
@@ -25 +26 @@
 c     input and output variables
 
-      integer    ig,nlayer
-      integer    chemthermod
-      integer    nesptherm                      !Number of species considered
-      real       rm(nlayer,nesptherm)         !Densities (cm-3)
-      real       tx(nlayer)                   !temperature
-      real       zenit                          !SZA
-      real       iz(nlayer)                   !Local altitude
-      real       zday                           !Martian day after Ls=0
-
+      integer,intent(in) :: ig,nlayer
+      integer,intent(in) :: chemthermod
+      integer,intent(in) :: nesptherm          !Number of species considered
+      real,intent(in) :: rm(nlayer,nesptherm)  !Densities (cm-3)
+      real,intent(in) :: tx(nlayer)            !temperature
+      real,intent(in) :: zenit         !SZA
+      real,intent(in) :: iz(nlayer)    !Local altitude
+      real,intent(in) :: zday          !Martian day after Ls=0
+
+! NB: no output variable! (computed jfotsout() is in module param_v4_h)
 
 c    local parameters and variables
@@ -1115 +1117 @@
 c     End of interpolation to obtain photoabsorption rates
 
-
-      return
+c     Correction to the actual Sun-Mars distance
+
+      jfotsout(:,:,:)=jfotsout(:,:,:)*(1.52/dist_sol)**2
 
       end