Changeset 1420 for trunk/LMDZ.GENERIC

Timestamp:

Apr 17, 2015, 5:40:46 PM (10 years ago)

Author:

sglmd

Message:

some clean-up: hard coded parameters removed

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• dsolver.F (modified) (1 diff)
• dtridgl.F (modified) (1 diff)
• gfluxi.F (modified) (3 diffs)
• gfluxv.F (modified) (2 diffs)

trunk/LMDZ.GENERIC/libf/phystd/dsolver.F

@@ -6 +6 @@
 C DOUBLE PRECISION VERSION OF SOLVER
 
-      PARAMETER (NMAX=201)
+!!      PARAMETER (NMAX=201)
       IMPLICIT REAL*8  (A-H,O-Z)
       DIMENSION GAMA(NL),CP(NL),CM(NL),CPM1(NL),CMM1(NL),XK1(NL),
      *          XK2(NL),E1(NL),E2(NL),E3(NL),E4(NL)
-      DIMENSION AF(NMAX),BF(NMAX),CF(NMAX),DF(NMAX),XK(NMAX)
+      DIMENSION AF(2*NL),BF(2*NL),CF(2*NL),DF(2*NL),XK(2*NL)
 C*********************************************************
 C* THIS SUBROUTINE SOLVES FOR THE COEFFICIENTS OF THE    *

trunk/LMDZ.GENERIC/libf/phystd/dtridgl.F

@@ -5 +5 @@
 C     DOUBLE PRECISION VERSION OF TRIDGL
 
-      PARAMETER (NMAX=201)
+!!      PARAMETER (NMAX=201)
       IMPLICIT REAL*8  (A-H,O-Z)
       DIMENSION AF(L),BF(L),CF(L),DF(L),XK(L)
-      DIMENSION AS(NMAX),DS(NMAX)
+      DIMENSION AS(2*L),DS(2*L)
 
 C*    THIS SUBROUTINE SOLVES A SYSTEM OF TRIDIAGIONAL MATRIX

trunk/LMDZ.GENERIC/libf/phystd/gfluxi.F

@@ -36 +36 @@
       implicit none
 
-      INTEGER NLP
-      PARAMETER (NLP=101) ! MUST BE LARGER THAN NLEVEL 
 
       INTEGER NLL, NLAYER, L, NW, NT, NT2
@@ -48 +46 @@
       REAL*8  WAVEN, DW, UBARI, RSF
       REAL*8  BTOP, BSURF, FMIDP(L_NLAYRAD), FMIDM(L_NLAYRAD)
-      REAL*8  B0(NLP),B1(NLP),ALPHA(NLP),LAMDA(NLP),XK1(NLP),XK2(NLP)
-      REAL*8  GAMA(NLP),CP(NLP),CM(NLP),CPM1(NLP),CMM1(NLP),E1(NLP)
-      REAL*8  E2(NLP),E3(NLP),E4(NLP)
+      REAL*8  B0(L_NLAYRAD),B1(L_NLAYRAD),ALPHA(L_NLAYRAD)
+      REAL*8  LAMDA(L_NLAYRAD),XK1(L_NLAYRAD),XK2(L_NLAYRAD)
+      REAL*8  GAMA(L_NLAYRAD),CP(L_NLAYRAD),CM(L_NLAYRAD)
+      REAL*8  CPM1(L_NLAYRAD),CMM1(L_NLAYRAD),E1(L_NLAYRAD)
+      REAL*8  E2(L_NLAYRAD),E3(L_NLAYRAD),E4(L_NLAYRAD)
 
       REAL*8  FTOPUP, FLUXUP, FLUXDN
@@ -60 +60 @@
 C     WE GO WITH THE HEMISPHERIC CONSTANT APPROACH IN THE INFRARED
  
-
-      IF (NLL .GT. NLP) STOP 'PARAMETER NL TOO SMALL IN GFLUXI'
 
       NLAYER = L_NLAYRAD

trunk/LMDZ.GENERIC/libf/phystd/gfluxv.F

@@ -42 +42 @@
       implicit none
 
-      INTEGER NLP
-      PARAMETER (NLP=101) ! MUST BE LARGER THAN NLEVEL
+!!      INTEGER NLP
+!!      PARAMETER (NLP=101) ! MUST BE LARGER THAN NLEVEL
 
       REAL*8 EM, EP
@@ -50 +50 @@
       REAL*8 DTDEL(L_NLAYRAD), TDEL(L_NLEVRAD)
       REAL*8 FMIDP(L_NLAYRAD), FMIDM(L_NLAYRAD)
-      REAL*8 LAMDA(NLP), ALPHA(NLP), XK1(NLP), XK2(NLP)
-      REAL*8 G1(NLP), G2(NLP), G3(NLP), GAMA(NLP), CP(NLP), CM(NLP)
-      REAL*8 CPM1(NLP)
-      REAL*8 CMM1(NLP), E1(NLP), E2(NLP), E3(NLP), E4(NLP), EXPTRM(NLP)
+      REAL*8 LAMDA(L_NLAYRAD), ALPHA(L_NLAYRAD), XK1(L_NLAYRAD)
+      REAL*8 XK2(L_NLAYRAD),G1(L_NLAYRAD), G2(L_NLAYRAD)
+      REAL*8 G3(L_NLAYRAD), GAMA(L_NLAYRAD),CP(L_NLAYRAD),CM(L_NLAYRAD)
+      REAL*8 CPM1(L_NLAYRAD),CMM1(L_NLAYRAD), E1(L_NLAYRAD)
+      REAL*8 E2(L_NLAYRAD),E3(L_NLAYRAD),E4(L_NLAYRAD),EXPTRM(L_NLAYRAD)
       REAL*8 FLUXUP, FLUXDN
       REAL*8 FACTOR, TAUCUMIN(L_LEVELS), TAUCUM(L_LEVELS)