Changeset 1927

Timestamp:

May 29, 2018, 12:53:32 AM (6 years ago)

Author:

jvatant

Message:

Correct a mistake in eddy diffusion profiles in calchim
--JVO

Location:

trunk/LMDZ.TITAN

Files:

• README (modified) (1 diff)
• libf/phytitan/calchim.F90 (modified) (3 diffs)
• libf/phytitan/physiq_mod.F90 (modified) (1 diff)

trunk/LMDZ.TITAN/README

@@ -1444 +1444 @@
 + Add supplementary diagnostics outputs for microphysics ( precip, flux, rc ... ) ( new muphy_diag.F90 module )
 + Correct the outputs of microphys tracers to be in X/m-3 to be comparable to "standard values"
++ Correct outputs of chemical tracers to be in mol/mol
 
 + Also update the deftank callphys.def with past revs modifs of microphysics

trunk/LMDZ.TITAN/libf/phytitan/calchim.F90

@@ -315 +315 @@
      ! 1E6 at 300 km ~ 10-1 mbar
         IF     ( logp.ge.1.0 .and. logp.le.4.0 ) THEN 
-              kedd(l) = 10.**(6.0+1.1*(logp-1.0))
+              kedd(l) = 10.**(6.0+1.3*(logp-1.0)/3.0)
      ! 2E7 at 600 km ~ 10-4 mbar
         ELSEIF ( logp.gt.4.0 .and. logp.le.6.0 ) THEN
@@ -323 +323 @@
              kedd(l) = 1.e8
         ENDIF
+        print *, l, kedd(l)
      ENDDO
      
@@ -329 +330 @@
         temp1   = ( log10(press_c(l)/press_c(klev-15)) ) / ( log10(press_c(klev+1)/press_c(klev-15)) )
         kedd(l) = 10.**( 3.0 + log10(kedd(klev+1)/1.e3) * temp1 )
+        print *, l, kedd(l)
      ENDDO
      

trunk/LMDZ.TITAN/libf/phytitan/physiq_mod.F90

@@ -1594 +1594 @@
       if (callchim) then
         do iq=1,nkim
-          CALL send_xios_field(cnames(iq),zq(:,:,iq+nmicro)/rat_mmol(iq+nmicro))
+          CALL send_xios_field(cnames(iq),zq(:,:,iq+nmicro)/rat_mmol(iq+nmicro)) ! mol/mol
         enddo
       endif