Changeset 421 for trunk/LMDZ.MARS/libf

Timestamp:

Nov 24, 2011, 4:05:22 PM (14 years ago)

Author:

acolaitis

Message:

## deposition.F: corrected a bug with z0, replaced with z0(ngrid) with the correct call to surfdat.h ## makemeso: corrected syntax of call to compile to make sure that there are no conflicts with other instances of compile in the environment of the user.

File:

• trunk/LMDZ.MARS/libf/aeronomars/deposition.F (modified) (4 diffs)

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

@@ -14 +14 @@
 #include "chimiedata.h"
 #include "conc.h"
+#include "surfdat.h"
 c
 c     input
@@ -72 +73 @@
 c        friction velocity (m.s-1)
 c
-         ustar = ubar*karman/log(zzlay(ig,1)/z0)
+         ustar = ubar*karman/log(zzlay(ig,1)/z0(ig))
 c
 c        aerodynamic resistance (s.m-1)
 c
-         ra = 1./(karman*ustar)*log(zzlay(ig,1)/z0)
+         ra = 1./(karman*ustar)*log(zzlay(ig,1)/z0(ig))
 c
 c        molecular diffusivity of methane in *air*  (cm2.s-1)
@@ -127 +128 @@
          gravity = 3.7
          dp = pplev(ig,1) - pplev(ig,2)
-         cd = (karman/log(zzlay(ig,1)/z0))**2.
+         cd = (karman/log(zzlay(ig,1)/z0(ig)))**2.
 c
          gam = (4./nuch4)*dp
@@ -142 +143 @@
          if (ig .eq. ig_vl1) then
             print*,'**** deposition ****'
-            print*,'z0     = ', z0, 'm'
+            print*,'z0     = ', z0(ig), 'm'
             print*,'deltaz = ', deltaz, ' m'
             print*,'deltap = ', dp, 'pa'