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

Timestamp:

Sep 11, 2019, 9:13:17 AM (5 years ago)

Author:

emillour

Message:

Mars GCM:

• Updated chemical core to include ionospheric chemistry

FGG

File:

• trunk/LMDZ.MARS/libf/phymars/physiq_mod.F (modified) (8 diffs)

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

@@ -31 +31 @@
      &                      nuice_ref, rho_ice, rho_dust, ref_r0,
      &                      igcm_he, igcm_stormdust_mass,
-     &                      igcm_stormdust_number 
+     &                      igcm_stormdust_number
       use comsoil_h, only: inertiedat, ! soil thermal inertia
      &                     tsoil, nsoilmx ! number of subsurface layers
@@ -576 +576 @@
          if(callnlte.and.nltemodel.eq.2) call nlte_setup
          if(callnirco2.and.nircorr.eq.1) call NIR_leedat
-         if(thermochem) call chemthermos_readini
+
 
         IF (tracer.AND.water.AND.(ngrid.NE.1)) THEN
@@ -1614 +1614 @@
 c        photochemistry :
 c        --------------
-         IF (photochem .or. thermochem) then
+         IF (photochem) then
 
 !           dust and ice surface area
@@ -1665 +1665 @@
            endif
 
-         END IF  ! of IF (photochem.or.thermochem)
+         END IF  ! of IF (photochem)
 #endif
 
@@ -2396 +2396 @@
            endif ! (dustbin.ne.0)
 
-           if (thermochem .or. photochem) then
+           if (photochem) then
               do iq=1,nq
                  if (noms(iq) .ne. "dust_mass" .and.
@@ -2428 +2428 @@
      &                           (1000*mmol(iq))
 
-!                   call wstats(ngrid,"rho_"//trim(noms(iq)),
-!    $                     "Number density","cm-3",3,rhopart)
-!                   call writediagfi(ngrid,"rho_"//trim(noms(iq)),
-!    $                     "Number density","cm-3",3,rhopart)
+                   call wstats(ngrid,"rho_"//trim(noms(iq)),
+     $                   "Number density","cm-3",3,rhopart)
+                   call writediagfi(ngrid,"rho_"//trim(noms(iq)),
+     $                  "Number density","cm-3",3,rhopart)
 
 !                   vertical column (molecule.cm-2)
@@ -2461 +2461 @@
                  end if ! of if (noms(iq) .ne. "dust_mass" ...)
               end do ! of do iq=1,nq
-           end if ! of if (thermochem .or. photochem)
+           end if ! of if (photochem)
 
            end if ! of if (tracer)
@@ -2939 +2939 @@
             call WRITEDIAGFI(ngrid,"qnir","NIR heating","K/s",
      $           3,zdtnirco2)
+
+!            call wstats(ngrid,"q15um","15 um cooling","K/s",
+!     $           3,zdtnlte)
+!            call wstats(ngrid,"quv","UV heating","K/s",
+!     $           3,zdteuv)
+!            call wstats(ngrid,"cond","Thermal conduction","K/s",
+!     $           3,zdtconduc)
+!            call wstats(ngrid,"qnir","NIR heating","K/s",
+!     $           3,zdtnirco2)
 
          endif  !(callthermos)