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

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/phymars

Files:

• callkeys.h (modified) (2 diffs)
• conf_phys.F (modified) (1 diff)
• physiq_mod.F (modified) (1 diff)

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