Ignore:
Timestamp:
Jun 20, 2015, 9:22:53 AM (10 years ago)
Author:
emillour
Message:

Common dynamics: A couple of bug fixes

  • calfis[_p].F array boundaries must be explicitely specified when underlying arrays are of different sizes.
  • advect_new_p.F : missing initializations of intermediate variables at topmost layer.

EM

Location:
trunk/LMDZ.COMMON/libf/dynlonlat_phylonlat
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.COMMON/libf/dynlonlat_phylonlat/calfis.F

    r1422 r1459  
    9494c    ------------------
    9595
    96 #include "dimensions.h"
    97 #include "paramet.h"
     96      include "dimensions.h"
     97      include "paramet.h"
    9898
    9999      INTEGER ngridmx
    100100      PARAMETER( ngridmx = 2+(jjm-1)*iim - 1/jjm   )
    101101
    102 #include "comgeom2.h"
    103 #include "iniprint.h"
     102      include "comgeom2.h"
     103      include "iniprint.h"
    104104
    105105c    Arguments :
    … …  
    259259        call tpot2t(llm,tetamoy,tmoy,pkmoy)
    260260c SI ON TIENT COMPTE DE LA VARIATION DE G AVEC L'ALTITUDE:
    261         zlaymoy(:) = g*rad*rad/(g*rad-phimoy(:))-rad
     261        zlaymoy(1:llm) = g*rad*rad/(g*rad-phimoy(1:llm))-rad
    262262        zlevmoy(1) = phimoy(0)/g
    263263        DO l=2,llm

  • trunk/LMDZ.COMMON/libf/dynlonlat_phylonlat/calfis_p.F

    r1422 r1459  
    103103c    ------------------
    104104
    105 #include "dimensions.h"
    106 #include "paramet.h"
     105      include "dimensions.h"
     106      include "paramet.h"
    107107
    108108      INTEGER ngridmx
    109109      PARAMETER( ngridmx = 2+(jjm-1)*iim - 1/jjm   )
    110110
    111 #include "comgeom2.h"
    112 #include "iniprint.h"
     111      include "comgeom2.h"
     112      include "iniprint.h"
    113113#ifdef CPP_MPI
    114114      include 'mpif.h'
    … …  
    348348        call tpot2t_p(1,llm,tetamoy,tmoy,pkmoy)
    349349c SI ON TIENT COMPTE DE LA VARIATION DE G AVEC L'ALTITUDE:
    350         zlaymoy(:) = g*rad*rad/(g*rad-phimoy(:))-rad
     350        zlaymoy(1:llm) = g*rad*rad/(g*rad-phimoy(1:llm))-rad
    351351        zlevmoy(1) = phimoy(0)/g
    352352        DO l=2,llm
Note: See TracChangeset for help on using the changeset viewer.