Changeset 459 for trunk/LMDZ.MARS/libf/aeronomars

Timestamp:

Dec 8, 2011, 10:04:00 AM (13 years ago)

Author:

emillour

Message:

Mars GCM: Include changes and updates for photochemistry by FL:

• aeronomars/calchim.F : change in units of surface density.
• aeronomars/surfacearea.F : new routine to compute ice and dust surface area

(m2/m3) available for heterogeneous reactions.

• phymars/initracer.F : bug correction: initialize igcm_ch4 and change loop

bounds (when guessing tracer names/properties with
old input files).

• phymars/watercloud.F : cleanup.
• phymars/physiq.F : add call to surfacearea; photochemistry is now called

after sedimentation to take into acount updated rdust
and rice.

EM

Location:

trunk/LMDZ.MARS/libf/aeronomars

Files:

• calchim.F (modified) (4 diffs)
• surfacearea.F (added)

trunk/LMDZ.MARS/libf/aeronomars/calchim.F

@@ -21 +21 @@
 c    update sept. 2008 identify tracers by their names (Ehouarn Millour)
 c    update 17/03/2011 synchronize with latest version of chemistry (Franck Lefevre)
+c    update 05/12/2011 synchronize with latest version of chemistry (Franck Lefevre)
 c
 c   Arguments:
@@ -42 +43 @@
 c    tauref(ngridmx)            Optical depth at 7 hPa
 c    co2ice(ngridmx)            co2 ice surface layer (kg.m-2)
-c    surfdust(ngridmx,nlayermx) dust surface area (micron2/cm3)
-c    surfice(ngridmx,nlayermx)  ice surface area (micron2/cm3)
+c    surfdust(ngridmx,nlayermx) dust surface area (m2/m3)
+c    surfice(ngridmx,nlayermx)  ice surface area (m2/m3)
 c
 c  Output:
@@ -87 +88 @@
       real    tauref(ngridmx)            ! optical depth at 7 hPa
       real    co2ice(ngridmx)            ! co2 ice surface layer (kg.m-2)
-      real    surfdust(ngridmx,nlayermx) ! dust surface area (micron2/cm3)
-      real    surfice(ngridmx,nlayermx)  !  ice surface area (micron2/cm3)
+      real    surfdust(ngridmx,nlayermx) ! dust surface area (m2/m3)
+      real    surfice(ngridmx,nlayermx)  !  ice surface area (m2/m3)
 
 c   outputs:
@@ -356 +357 @@
             zlocal(l) = zzlay(ig,l)/1000.
 c
-c           surfdust1d and surfice1d in cm2/cm3
-c
-            surfdust1d(l) = surfdust(ig,l)*1.e-8
-            surfice1d(l)  = surfice(ig,l)*1.e-8
+c           surfdust1d and surfice1d: conversion from m2/m3 to cm2/cm3
+c
+            surfdust1d(l) = surfdust(ig,l)*1.e-2
+            surfice1d(l)  = surfice(ig,l)*1.e-2
 c
 c           search for switch index between regions