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

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

@@ -39 +39 @@
 #include "dimensions.h"
       integer, parameter :: ngridmx = (2+(jjm-1)*iim - 1/jjm) 
-#include "dimphys.h"
-!#include "surfdat.h"
-!#include "comsoil.h"
-!#include "dimradmars.h"
 #include "paramet.h"
 #include "comconst.h"
 #include "comvert.h"
 #include "comgeom2.h"
-!#include "control.h"
 #include "logic.h"
 #include "description.h"
 #include "ener.h"
 #include "temps.h"
-!#include "lmdstd.h"
 #include "comdissnew.h"
 #include "clesph0.h"
 #include "serre.h"
 #include "netcdf.inc"
-!#include"advtrac.h"
-!#include"tracer.h"
 #include "datafile.h"
 c=======================================================================
@@ -105 +97 @@
 c variable physique
 c------------------
-!      REAL tsurf(ngridmx)      ! surface temperature 
-!      REAL tsoil(ngridmx,nsoilmx) ! soil temperature
-!      REAL co2ice(ngridmx)     ! CO2 ice layer
-!      REAL emis(ngridmx)       ! surface emissivity
-!      REAL tauscaling(ngridmx) ! dust conversion factor
       REAL tauscadyn(iip1,jjp1) ! dust conversion factor on the dynamics grid
-!      REAL,ALLOCATABLE :: qsurf(:,:)
-!      REAL q2(ngridmx,nlayermx+1)
-!      REAL rnaturfi(ngridmx)
       real alb(iip1,jjp1),albfi(ngridmx) ! albedos
       real ith(iip1,jjp1,nsoilmx),ithfi(ngridmx,nsoilmx) ! thermal inertia (3D)
@@ -203 +187 @@
 ! allocate arrays
       allocate(q(iip1,jjp1,llm,nqtot))
-!      allocate(qsurf(ngridmx,nqtot)) ! done in ini_surfdat_h
       allocate(coefvmr(nqtot))