Changeset 1651

Timestamp:

Jan 30, 2017, 11:52:09 AM (8 years ago)

Author:

jaudouard

Message:

Modification of co2 clouds microtimestep and addition of a meteoritic flux of dust

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• deftank/callphys.def.co2clouds (modified) (2 diffs)
• libf/dyn3d/meteoriticflux_mod.F90 (added)
• libf/phymars/co2cloud.F (modified) (4 diffs)
• libf/phymars/conf_phys.F (modified) (2 diffs)
• libf/phymars/dyn1d/meteoriticflux_mod.F90 (added)

trunk/LMDZ.MARS/README

@@ -2374 +2374 @@
 -callphys.def.co2clouds was added in deftank. It contains the necessary flags and constant for including co2clouds
 -traceur.def.co2clouds was added in deftank. It supplies the necessary tracers for running the co2 clouds scheme
+
+== 30/01/2017 == JA
+-Modification of CO2 clouds microtimestep : now is set by imicro in callphys.def together with the water ice clouds microtimestep
+-Addition of a meteoritic flux of dust is now possible:
+          meteoriticflux_mod.F90 was added, it contains 3 variables that are provided in callphys.def (default values = no flux): meteo_flux_mass, meteo_flux_number and meteo_alt (in km)
+          meteo_flux_mass and meteo_flux_number are added to the dust tracers at the z-level corresponding to meteo_alt at every timestep inside co2cloud.F
+          conf_phys.F was modified to perform the initialization of the meteoritic flux with the callphys.def values (default values = no flux)
+          callphys.def.co2clouds in deftank has been accordingly modified
+          
+
+
+
+          

trunk/LMDZ.MARS/deftank/callphys.def.co2clouds

@@ -127 +127 @@
 naerkind = 1
 
+imicro=30
 #include co2 clouds scheme?
 co2clouds=.true.
@@ -164 +165 @@
 ## be used to compensate for low 15 um cooling)
 euveff = 0.21       
+#meteoritic flux of dust 
+meteo_flux_mass=2.21938e-14
+meteo_flux_number=  2.11935e+09
+meteo_alt= 90
 
 

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

@@ -12 +12 @@
       use tracer_mod, only: nqmx, igcm_co2, igcm_co2_ice,
      &     igcm_dust_mass, igcm_dust_number,
-     &     igcm_dust_mass, igcm_dust_number,
      &     igcm_ccnco2_mass, igcm_ccnco2_number,
      &     rho_dust, nuiceco2_sed, nuiceco2_ref,
      &     rho_ice_co2,r3n_q
+      USE meteoriticflux_mod, ONLY:meteo_flux_mass,meteo_flux_number
+     &     ,meteo_alt
+
       IMPLICIT NONE
 
@@ -126 +128 @@
       real masse (ngrid,nlay) ! Layer mass (kg.m-2)
     
-
+      double precision diff,diff0
+      integer meteo_lvl
       real tempo_traceur_t(ngrid,nlay)
       real tempo_traceurs(ngrid,nlay,nq)
@@ -152 +155 @@
 #endif
         call getin("imicro",imicro)
-             imicro=ptimestep/50.
-
         write(*,*)"imicro = ",imicro
         
@@ -159 +160 @@
         write(*,*)"Physical timestep is",ptimestep 
         write(*,*)"CO2 Microphysics timestep is",microtimestep 
-
+      
+              
+        write(*,*) "Warning ! Meteoritic flux of dust is turned on"
+        write(*,*) "Dust mass = ",meteo_flux_mass 
+        write(*,*) "Dust number = ",meteo_flux_number
+        write(*,*) "Are added at the z-level (km)",meteo_alt
+        write(*,*) "Every timestep in co2cloud.F"
+     
         firstcall=.false.
       ENDIF                     ! of IF (firstcall)
    
 c-----Initialization
+
+c add meteoritic flux of dust 
+    !convert meteo_alt (in km) to z-level 
+        !pzlay altitudes of the layers
+    
+      do ig=1, ngrid
+         diff0=1000
+         meteo_lvl=1
+         do  l=1,nlay
+            diff=pzlay(ig,l)/1000-meteo_alt
+            if (abs(diff) .lt. diff0) then
+               diff0=abs(diff)
+               meteo_lvl=l
+            endif
+         enddo
+c         write(*,*) "meteo_lvl=",meteo_lvl
+         pq(ig,meteo_lvl,igcm_dust_mass)=pq(ig,meteo_lvl,igcm_dust_mass)
+     &        +meteo_flux_mass
+         pq(ig,meteo_lvl,igcm_dust_number)=
+     &        pq(ig,meteo_lvl,igcm_dust_number)+meteo_flux_number
+      enddo
+
+       call WRITEDIAGFI(ngrid,"pzlay","altitude","km",3,
+     &        pzlay)
       beta=0.85
       subpdq(1:ngrid,1:nlay,1:nq) = 0

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

@@ -44 +44 @@
       use dimradmars_mod, only: naerkind, name_iaer,
      &                      ini_scatterers,tauvis
+      USE meteoriticflux_mod !meteo_flux_mass,meteo_flux_number,alt_meteo
 
       IMPLICIT NONE
@@ -711 +712 @@
          write(*,*) " euveff = ", euveff
 
+        !Meteoritic flux of dust : number, mass and altitude
+        !dust tracers are incremented by thoses values in co2cloud.F
+        ! Thus only apply if co2clouds=.true. in callphys.def
+        meteo_flux_mass=0.
+        call getin("meteo_flux_mass",meteo_flux_mass)
+        meteo_flux_number=0.
+        call getin("meteo_flux_number",meteo_flux_number)
+        meteo_alt=5.
+        call getin("meteo_alt",meteo_alt)
+        WRITE(*,*) ""
+        WRITE(*,*) "meteoritic flux of dust"
+        WRITE(*,*) "dust mass = ",meteo_flux_mass
+        WRITE(*,*) "dust number = ",meteo_flux_number
+        WRITE(*,*) "at altitude = ",meteo_alt
+         
+
          if (.not.callthermos) then
            if (thermoswater) then