Changeset 1266 for trunk/LMDZ.MARS/libf/phymars/gwprofil.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/gwprofil.F (modified) (3 diffs)

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

@@ -58 +58 @@
 C
 
-#include "dimensions.h"
-#include "dimphys.h"
-!#include "dimradmars.h"
       integer klon,klev,kidia,kfdia
 #include "yoegwd.h"
@@ -91 +88 @@
 
       integer ji,jk,jl,ilevh
-      REAL ZDZ2 (NDLO2,nlayermx) , ZNORM(NDLO2) , zoro(NDLO2)
-      REAL ZTAU (NDLO2,nlayermx+1)
+      REAL ZDZ2 (NDLO2,klev) , ZNORM(NDLO2) , zoro(NDLO2)
+      REAL ZTAU (NDLO2,klev+1)
 C
 C-----------------------------------------------------------------------
@@ -182 +179 @@
 c     write(*,*) 'ptau'
 c     write(*,99) ((ji,ilevh,ptau(ji,ilevh),ji=1,NDLO2),
-c    .                  ilevh=1,nlayermx+1)
+c    .                  ilevh=1,klev+1)
  99   FORMAT(i3,i3,f15.5)