Changeset 1684

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

Files:

• README (modified) (1 diff)
• libf/aeronomars/chemthermos.F90 (modified) (1 diff)
• libf/aeronomars/euvheat.F90 (modified) (1 diff)
• libf/aeronomars/hrtherm.F (modified) (2 diffs)
• libf/aeronomars/jthermcalc_e107.F (modified) (3 diffs)
• libf/phymars/callkeys.h (modified) (2 diffs)
• libf/phymars/conf_phys.F (modified) (1 diff)
• libf/phymars/physiq_mod.F (modified) (1 diff)

trunk/LMDZ.MARS/README

@@ -2420 +2420 @@
 - print_control_mod.F90 : make initialization occur via init_print_control_mod
                           to avoid circular module dependencies
+
+== 31/03/2017 == EM + FGG
+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.

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

trunk/LMDZ.MARS/libf/phymars/callkeys.h

@@ -19 +19 @@
      &   ,dustbin,nltemodel,nircorr,solvarmod,solvaryear
      
-      COMMON/callkeys_r/topdustref,solarcondate,semi,alphan,euveff,     &
-     &   tke_heat_flux,dustrefir
+      COMMON/callkeys_r/topdustref,semi,alphan,euveff,                  &
+     &   tke_heat_flux,dustrefir,fixed_euv_value
      
       LOGICAL callrad,calldifv,calladj,callcond,callsoil,               &
@@ -38 +38 @@
       real semi
       real alphan
-      real solarcondate
+      real fixed_euv_value
       real euveff
       real tke_heat_flux

trunk/LMDZ.MARS/libf/phymars/conf_phys.F

@@ -704 +704 @@
 
          write(*,*) "Method to include solar variability"
-         write(*,*) "0-> old method (using solarcondate); ",
-     &                  "1-> variability wit E10.7"
+         write(*,*) "0-> fixed value of E10.7 (fixed_euv_value); ",
+     &          "1-> daily evolution of E10.7 (for given solvaryear)"
          solvarmod=1
          call getin("solvarmod",solvarmod)
          write(*,*) " solvarmod = ",solvarmod
 
-         write(*,*) "date for solar flux calculation:",
-     &   " (1985 < date < 2002)",
-     $   " (Only used if solvarmod=0)"
-         write(*,*) "(Solar min=1996.4 ave=1993.4 max=1990.6)"
-         solarcondate=1993.4 ! default value
-         call getin("solarcondate",solarcondate)
-         write(*,*) " solarcondate = ",solarcondate
+         write(*,*) "Fixed euv (for solvarmod==0) 10.7 value?"
+         write(*,*) " (min=80 , ave=140, max=320)"
+         fixed_euv_value=140 ! default value
+         call getin("fixed_euv_value",fixed_euv_value)
+         write(*,*) " fixed_euv_value = ",fixed_euv_value
          
          write(*,*) "Solar variability as observed for MY: "

trunk/LMDZ.MARS/libf/phymars/physiq_mod.F

@@ -484 +484 @@
             call allocate_param_thermos(nlayer)
             call allocate_param_iono(nlayer,nreact)
-            if(solvarmod.eq.0) call param_read
-            if(solvarmod.eq.1) call param_read_e107
+            call param_read_e107
          endif
 #endif