Changeset 1266 for trunk/LMDZ.MARS/libf/phymars/orodrag.F

Timestamp:

May 19, 2014, 10:05:08 AM (11 years ago)

Author:

aslmd

Message:

LMDZ.MARS
IMPORTANT CHANGE

• Remove all reference/use of nlayermx and dimphys.h
• Made use of automatic arrays whenever arrays are needed with dimension nlayer
• Remove lots of obsolete reference to dimensions.h
• Converted iono.h and param_v4.h into corresponding modules

(with embedded subroutine to allocate arrays)
(no arrays allocated if thermosphere not used)

• Deleted param.h and put contents into module param_v4_h
• Adapted testphys1d, newstart, etc...
• Made DATA arrays in param_read to be initialized by subroutine

fill_data_thermos in module param_v4_h

• Optimized computations in paramfoto_compact (twice less dlog10 calculations)
• Checked consistency before/after modification in debug mode
• Checked performance is not impacted (same as before)

File:

• trunk/LMDZ.MARS/libf/phymars/orodrag.F (modified) (3 diffs)

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

@@ -53 +53 @@
 C
 C     comcstfi.h
-C     dimphys.h
 C     yoegwd.h
 C
@@ -76 +75 @@
 C
 C
-#include "dimensions.h"
-#include "dimphys.h"
-!#include "dimradmars.h"
       integer klon,klev,kidia
       parameter(kidia=1)
@@ -123 +119 @@
       real zconb,zabsv,zzd1,ratio,zust,zvst,zdis,ztemp
 C
-      REAL   ZTAU(NDLO2,nlayermx+1),
-     *       ZSTAB(NDLO2,nlayermx+1),
-     *       ZVPH(NDLO2,nlayermx+1),
-     *       ZRHO(NDLO2,nlayermx+1),
-     *       ZRI(NDLO2,nlayermx+1),
-     *       ZpsI(NDLO2,nlayermx+1),
-     *       Zzdep(NDLO2,nlayermx)
+      REAL   ZTAU(NDLO2,klev+1),
+     *       ZSTAB(NDLO2,klev+1),
+     *       ZVPH(NDLO2,klev+1),
+     *       ZRHO(NDLO2,klev+1),
+     *       ZRI(NDLO2,klev+1),
+     *       ZpsI(NDLO2,klev+1),
+     *       Zzdep(NDLO2,klev)
       REAL   ZDUDT(NDLO2),
      *       ZDVDT(NDLO2),