Changeset 2678

Timestamp:

May 17, 2022, 3:22:49 PM (3 years ago)

Author:

tpierron

Message:

Mars GCM :
Add incoming solar direct flux at surface in stats outputs for the MCD.
+ declare new variable "flux_1AU" for solar constant.
TP

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/callradite_mod.F (modified) (2 diffs)
• libf/phymars/dimradmars_mod.F90 (modified) (1 diff)
• libf/phymars/physiq_mod.F (modified) (5 diffs)

trunk/LMDZ.MARS/README

@@ -3706 +3706 @@
 arguments were wrong).
 Turned inichim_newstart.F90 into a module.
+
+== 17/05/2022 == TP
+Add incoming solar direct flux at surface in stats outputs for the MCD.
++ declare new variable "flux_1AU" for solar constant.
+

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

@@ -20 +20 @@
      &            iaer_dust_conrath,iaer_dust_doubleq,
      &            iaer_dust_submicron, iaer_h2o_ice, iaer_co2_ice,
-     &            iaer_stormdust_doubleq,iaer_topdust_doubleq
+     &            iaer_stormdust_doubleq,iaer_topdust_doubleq,flux_1AU
       use yomlw_h, only: gcp, nlaylte
       use comcstfi_h, only: g,cpp
@@ -559 +559 @@
 c       -------------------------------------------------------
 c          Mars solar constant (W m-2)
-c          1370 W.m-2 is the solar constant at 1 AU.
-           cste_mars=1370./(dist_sol*dist_sol)
+c          flux_1AU = 1370 W.m-2 is the solar constant at 1 AU.
+           cste_mars=flux_1AU/(dist_sol*dist_sol)
            zzdtsw(:,:) = 0.
            call swmain ( nd, nflev,

trunk/LMDZ.MARS/libf/phymars/dimradmars_mod.F90

@@ -107 +107 @@
   integer,parameter :: nabsmx=2
   integer,parameter :: nt_pademx=19
+  
+! Solar flux constant at 1AU for slope scheme
+!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+  real, parameter :: flux_1AU = 1370. ! solar constant (W.m-2)
 
 !!

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

@@ -61 +61 @@
      &                          dtrad, fluxrad_sky, fluxrad, albedo,
      &                          naerkind, iaer_dust_doubleq, 
-     &                          iaer_stormdust_doubleq
+     &                          iaer_stormdust_doubleq, iaer_h2o_ice,
+     &                          flux_1AU
       use dust_param_mod, only: doubleq, lifting, callddevil,
      &                          tauscaling, odpref, dustbin,
@@ -283 +284 @@
       REAL sl_ra, sl_di0
       REAL sky
+      REAL direct_flux(ngrid)    ! Incident direct solar flux on Mars at surface (W.m-2)
 
       REAL,PARAMETER :: stephan = 5.67e-08 ! Stephan Boltzman constant
@@ -1076 +1078 @@
                  sl_di0 = 0.
                  if (mu0(ig) .gt. 0.) then
-                  sl_di0 = mu0(ig)*(exp(-sl_tau/mu0(ig)))
-                  sl_di0 = sl_di0*1370./dist_sol/dist_sol
+                  sl_di0 = mu0(ig)*(exp(-sl_tau/mu0(ig)))
+                  sl_di0 = sl_di0*flux_1AU/dist_sol/dist_sol
                   sl_di0 = sl_di0/ztim1
                   sl_di0 = fluxsurf_sw(ig,l)*sl_di0
@@ -2446 +2448 @@
 c          diagfi files
 c        -------------------------------------------------------------------
+
+
+         do ig=1,ngrid
+           if(mu0(ig).le.0.01) then 
+            direct_flux(ig) = 0.
+           else 
+            direct_flux(ig) = flux_1AU/dist_sol/dist_sol*mu0(ig)*
+     &                    exp(-(tau(ig,iaer_dust_doubleq)+
+     &                          tau(ig,iaer_h2o_ice))/mu0(ig))
+           endif
+         enddo
 
          if (tracer) then
@@ -2672 +2685 @@
 c          call wstats(ngrid,"lw_htrt","lw heat.rate",
 c    &                 "W.m-2",3,zdtlw)
+          call wstats(ngrid,"direct_flux",
+     &                "Direct incoming flux at surface",
+     &                "W.m-2",2,direct_flux)     
 
           if (calltherm) then