Changeset 1684 for trunk/LMDZ.MARS/libf/aeronomars

Timestamp:

Mar 31, 2017, 5:03:31 PM (8 years ago)

Author:

emillour

Message:

Mars GCM:
Add possibility to fix EUV input as E10.7 value and remove previous system
(which used parameter solarcondate). The E10.7 value is now set via
callphys.def by parameter "fixed_euv_value" which is only used if
solvarmod==0.
Guidelines for min/ave/max EUV input: fixed_euv_value=80/140/320.
EM + FGG

Location:

trunk/LMDZ.MARS/libf/aeronomars

Files:

• chemthermos.F90 (modified) (1 diff)
• euvheat.F90 (modified) (1 diff)
• hrtherm.F (modified) (2 diffs)
• jthermcalc_e107.F (modified) (3 diffs)

trunk/LMDZ.MARS/libf/aeronomars/chemthermos.F90

@@ -466 +466 @@
 
       !Photoabsorption coefficients
-      call flujo(solarcondate)
-      if(solvarmod.eq.0) then
-         call jthermcalc(ig,nlayer,chemthermod,rm,nesptherm,ztemp,zlocal,zenit)
-      else if(solvarmod.eq.1) then
-         call jthermcalc_e107(ig,nlayer,chemthermod,rm,nesptherm,ztemp,zlocal,zenit,zday)
-      endif
+      call jthermcalc_e107(ig,nlayer,chemthermod,rm,nesptherm,ztemp,zlocal,zenit,zday)
          
 

trunk/LMDZ.MARS/libf/aeronomars/euvheat.F90

@@ -359 +359 @@
       enddo
       
-      !Solar flux calculation
-      call flujo(solarcondate)
-
-                                         ! Not recommended for long runs 
-                                         ! (e.g. to build the MCD), as the
-                                         ! solar conditions at the end will 
-                                         ! be different to the ones initially
-                                         ! set
       
       do ig=1,ngrid

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

@@ -98 +98 @@
 
       !Calculation of photoabsortion coefficient
-      if(solvarmod.eq.0) then
-         call jthermcalc(ig,nlayer,euvmod,rm,nespeuv,tx,iz,zenit)
-      else if (solvarmod.eq.1) then
-         call jthermcalc_e107(ig,nlayer,euvmod,
+      call jthermcalc_e107(ig,nlayer,euvmod,
      .           rm,nespeuv,tx,iz,zenit,zday)
-      endif
 
       !Total photoabsorption coefficient
@@ -118 +114 @@
       end do
 
-
-      return
-
       end
 

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

@@ -20 +20 @@
 
       implicit none
+
+      include "callkeys.h"
 
 c     input and output variables
@@ -136 +138 @@
       end do
 
-      !E10.7 for the day: linear interpolation to tabulated values
-      realday=mod(zday,669.)
-      if(realday.lt.date_e107(1)) then
+      if (solvarmod==0) then
+        e107=fixed_euv_value
+      else
+        !E10.7 for the day: linear interpolation to tabulated values
+        realday=mod(zday,669.)
+        if(realday.lt.date_e107(1)) then
          e107=e107_tab(1)
-      else if(realday.ge.date_e107(669)) then
+        else if(realday.ge.date_e107(669)) then
          e107=e107_tab(669)   
-      else if(realday.ge.date_e107(1).and.
+        else if(realday.ge.date_e107(1).and.
      $        realday.lt.date_e107(669)) then
          do i=1,668
@@ -153 +158 @@
             endif
          enddo
-      endif
+        endif
+      endif ! of if (solvarmod==0)
 
       !Photoabsorption coefficients at TOA as a function of E10.7