Changeset 1541 for trunk/LMDZ.MARS/libf/aeronomars/moldiff_red.F90

Timestamp:

Apr 17, 2016, 10:23:21 PM (9 years ago)

Author:

emillour

Message:

Mars GCM:

• fix for 1D in writediagfi to enable writing at "ecritphy" rate.
• removed iniprint.h from phymars/dyn1d since it is in "misc"
• Some code cleanup in anticipation of future updates:
  • changed variable names in comgeomphy.F90: give them more explicit names: rlond => longitude , rlatd => latitude, airephy => cell_area, cuphy => dx , cvphy => dy
  • removed long(), lati() and area() from comgeomfi_h.F90, use longitude(), latitude() and cell_are() from comgeomphy.F90 instead

EM

File:

• trunk/LMDZ.MARS/libf/aeronomars/moldiff_red.F90 (modified) (2 diffs)

trunk/LMDZ.MARS/libf/aeronomars/moldiff_red.F90

@@ -3 +3 @@
 
 use tracer_mod, only: noms, mmol
-use comgeomfi_h, only: area
+use comgeomphy, only: cell_area
 
 implicit none
@@ -591 +591 @@
 ! the trend only at the end
 
-        PhiEscH=PhiEscH+wi(i_h)*Nrafk(nlraf,i_h)*area(ig) ! in s-1
-        PhiEscH2=PhiEscH2+wi(i_h2)*Nrafk(nlraf,i_h2)*area(ig) ! in s-1 (U in m/s, aire in m2, Nrafk in m-3)
-!       print*,'test',ig,wi(i_h),Nrafk(nlraf,i_h),wi(i_h2),Nrafk(nlraf,i_h2),area(ig),PhiEscH,PhiEscH2,i_h,i_h2
+        PhiEscH=PhiEscH+wi(i_h)*Nrafk(nlraf,i_h)*cell_area(ig) ! in s-1
+        PhiEscH2=PhiEscH2+wi(i_h2)*Nrafk(nlraf,i_h2)*cell_area(ig) ! in s-1 (U in m/s, aire in m2, Nrafk in m-3)
+!       print*,'test',ig,wi(i_h),Nrafk(nlraf,i_h),wi(i_h2),Nrafk(nlraf,i_h2),cell_area(ig),PhiEscH,PhiEscH2,i_h,i_h2
 !       stop