Changeset 1530


Ignore:
Timestamp:
Apr 6, 2016, 8:51:34 AM (9 years ago)
Author:
emillour
Message:

Venus and Titan GCMs:
Updates in the physics to keep up with updates in LMDZ5 (up to
LMDZ5 trunk, rev 2350) concerning dynamics/physics separation:

  • Adapted makelmdz and makelmdz_fcm script to stop if trying to compile 1d model or newstart or start2archive in parallel.
  • got rid of references to "dimensions.h" in physics. Within physics packages, use nbp_lon (=iim), nbp_lat (=jjmp1) and nbp_lev (=llm) from module mod_grid_phy_lmdz (in phy_common) instead. Only partially done for Titan, because of many hard-coded commons; a necessary first step will be to clean these up (using modules).

EM

Location:
trunk
Files:
6 deleted
98 edited
4 moved

Legend:

Unmodified
Added
Removed
  • trunk/DOC/chantiers/commit_importants.log

    r1525 r1530  
    16591659Ehouarn: Further adaptations to enforce dynamics/physics separation:
    16601660- get rid of references to "control_mod" from physics packages.
     1661
     1662**********************
     1663**** commit_v1530 ****
     1664**********************
     1665Ehouarn: Updates in the physics (Venus and Titan) to keep up with updates
     1666in LMDZ5 (up to LMDZ5 trunk, rev 2350) concerning dynamics/physics separation:
     1667- Adapted makelmdz and makelmdz_fcm script to stop if trying to compile
     1668  1d model or newstart or start2archive in parallel.
     1669- got rid of references to "dimensions.h" in physics.
     1670  Within physics packages, use nbp_lon (=iim), nbp_lat (=jjmp1) and
     1671  nbp_lev (=llm) from module mod_grid_phy_lmdz (in phy_common) instead.
     1672  Only partially done for Titan, because of many hard-coded commons; a
     1673  necessary first step will be to clean these up (using modules).
  • trunk/LMDZ.COMMON/libf/dynlonlat_phylonlat/phyvenus/readstartphy.F

    r1520 r1530  
    1515c======================================================================
    1616      use dimphy
     17      USE mod_grid_phy_lmdz, only: klon_glo
    1718      IMPLICIT none
    18 #include "dimensions.h"
    1919#include "netcdf.inc"
    2020#include "dimsoil.h"
    2121c======================================================================
    22       integer ngridmx
    23       parameter (ngridmx=(2+(jjm-1)*iim - 1/jjm))
    2422
    2523      CHARACTER*(*) fichnom
    26       REAL rlat(ngridmx), rlon(ngridmx)
    27       REAL tsol(ngridmx)
    28       REAL tsoil(ngridmx,nsoilmx)
    29       REAL albe(ngridmx)
    30       REAL radsol(ngridmx)
    31       REAL sollw(ngridmx)
    32       real solsw(ngridmx)
    33       real fder(ngridmx)
    34       real dlw(ngridmx)
    35       real sollwdown(ngridmx)
    36       REAL zmea(ngridmx), zstd(ngridmx)
    37       REAL zsig(ngridmx), zgam(ngridmx), zthe(ngridmx)
    38       REAL zpic(ngridmx), zval(ngridmx)
     24      REAL rlat(klon_glo), rlon(klon_glo)
     25      REAL tsol(klon_glo)
     26      REAL tsoil(klon_glo,nsoilmx)
     27      REAL albe(klon_glo)
     28      REAL radsol(klon_glo)
     29      REAL sollw(klon_glo)
     30      real solsw(klon_glo)
     31      real fder(klon_glo)
     32      real dlw(klon_glo)
     33      real sollwdown(klon_glo)
     34      REAL zmea(klon_glo), zstd(klon_glo)
     35      REAL zsig(klon_glo), zgam(klon_glo), zthe(klon_glo)
     36      REAL zpic(klon_glo), zval(klon_glo)
    3937      INTEGER length
    4038      PARAMETER (length=100)
     
    129127         xmin = 1.0E+20
    130128         xmax = -1.0E+20
    131          DO i = 1, ngridmx
     129         DO i = 1, klon_glo
    132130            xmin = MIN(tsol(i),xmin)
    133131            xmax = MAX(tsol(i),xmax)
     
    148146         PRINT*, "phyetat0: Le champ <Tsoil"//str2//"> est absent"
    149147         PRINT*, "          Il prend donc la valeur de surface"
    150          DO i=1, ngridmx
     148         DO i=1, klon_glo
    151149             tsoil(i,isoil)=tsol(i)
    152150         ENDDO
     
    185183         xmin = 1.0E+20
    186184         xmax = -1.0E+20
    187          DO i = 1, ngridmx
     185         DO i = 1, klon_glo
    188186            xmin = MIN(albe(i),xmin)
    189187            xmax = MAX(albe(i),xmax)
     
    213211      xmin = 1.0E+20
    214212      xmax = -1.0E+20
    215       DO i = 1, ngridmx
     213      DO i = 1, klon_glo
    216214         xmin = MIN(solsw(i),xmin)
    217215         xmax = MAX(solsw(i),xmax)
     
    239237      xmin = 1.0E+20
    240238      xmax = -1.0E+20
    241       DO i = 1, ngridmx
     239      DO i = 1, klon_glo
    242240         xmin = MIN(sollw(i),xmin)
    243241         xmax = MAX(sollw(i),xmax)
     
    266264      xmin = 1.0E+20
    267265      xmax = -1.0E+20
    268       DO i = 1, ngridmx
     266      DO i = 1, klon_glo
    269267         xmin = MIN(fder(i),xmin)
    270268         xmax = MAX(fder(i),xmax)
     
    293291      xmin = 1.0E+20
    294292      xmax = -1.0E+20
    295       DO i = 1, ngridmx
     293      DO i = 1, klon_glo
    296294         xmin = MIN(dlw(i),xmin)
    297295         xmax = MAX(dlw(i),xmax)
     
    320318      xmin = 1.0E+20
    321319      xmax = -1.0E+20
    322       DO i = 1, ngridmx
     320      DO i = 1, klon_glo
    323321         xmin = MIN(sollwdown(i),xmin)
    324322         xmax = MAX(sollwdown(i),xmax)
     
    345343      xmin = 1.0E+20
    346344      xmax = -1.0E+20
    347       DO i = 1, ngridmx
     345      DO i = 1, klon_glo
    348346         xmin = MIN(radsol(i),xmin)
    349347         xmax = MAX(radsol(i),xmax)
     
    372370      xmin = 1.0E+20
    373371      xmax = -1.0E+20
    374       DO i = 1, ngridmx
     372      DO i = 1, klon_glo
    375373         xmin = MIN(zmea(i),xmin)
    376374         xmax = MAX(zmea(i),xmax)
     
    396394      xmin = 1.0E+20
    397395      xmax = -1.0E+20
    398       DO i = 1, ngridmx
     396      DO i = 1, klon_glo
    399397         xmin = MIN(zstd(i),xmin)
    400398         xmax = MAX(zstd(i),xmax)
     
    420418      xmin = 1.0E+20
    421419      xmax = -1.0E+20
    422       DO i = 1, ngridmx
     420      DO i = 1, klon_glo
    423421         xmin = MIN(zsig(i),xmin)
    424422         xmax = MAX(zsig(i),xmax)
     
    444442      xmin = 1.0E+20
    445443      xmax = -1.0E+20
    446       DO i = 1, ngridmx
     444      DO i = 1, klon_glo
    447445         xmin = MIN(zgam(i),xmin)
    448446         xmax = MAX(zgam(i),xmax)
     
    468466      xmin = 1.0E+20
    469467      xmax = -1.0E+20
    470       DO i = 1, ngridmx
     468      DO i = 1, klon_glo
    471469         xmin = MIN(zthe(i),xmin)
    472470         xmax = MAX(zthe(i),xmax)
     
    492490      xmin = 1.0E+20
    493491      xmax = -1.0E+20
    494       DO i = 1, ngridmx
     492      DO i = 1, klon_glo
    495493         xmin = MIN(zpic(i),xmin)
    496494         xmax = MAX(zpic(i),xmax)
     
    516514      xmin = 1.0E+20
    517515      xmax = -1.0E+20
    518       DO i = 1, ngridmx
     516      DO i = 1, klon_glo
    519517         xmin = MIN(zval(i),xmin)
    520518         xmax = MAX(zval(i),xmax)
  • trunk/LMDZ.COMMON/makelmdz

    r1521 r1530  
    363363fi
    364364
     365if [[ "$parallel" != "none" && \
     366     ( "$code" == "newstart" || "$code" == "start2archive" ) ]]
     367then
     368    echo "newstart or start2archive utilities should be run in serial!"
     369    echo " Compile without -parallel [mpi|omp|mpi_omp] option!"
     370    exit
     371fi
     372
    365373if [[ "$parallel" == "none" ]]
    366374then
     
    580588   INCLUDE="$INCLUDE "'-I$(LIBF)/dyn2d'
    581589elif (( $dimc == 1 )) ; then
     590   ## Sanity check: 1D models should be used in serial
     591   if [[ $parallel != "none" ]]
     592   then
     593     echo "Error: a 1D model should not be compiled with "
     594     echo "       -parallel [mpi|omp|mpi_omp] option!"
     595     exit
     596   fi
    582597   echo pas de dynamique
    583598   src_dirs="$src_dirs phy${physique}/dyn1d"
  • trunk/LMDZ.COMMON/makelmdz_fcm

    r1521 r1530  
    296296fi
    297297
     298if [[ "$parallel" != "none" && \
     299     ( "$code" == "newstart" || "$code" == "start2archive" ) ]]
     300then
     301    echo "newstart or start2archive utilities should be run in serial!"
     302    echo " Compile without -parallel [mpi|omp|mpi_omp] option!"
     303    exit
     304fi
     305
    298306if [[ "$parallel" == "mpi" ]]
    299307then
     
    495503if (( $dimension == 1 ))
    496504then
     505## Sanity check: 1D models should be used in serial
     506  if [[ $parallel != "none" ]]
     507  then
     508    echo "Error: a 1D model should not be compiled with "
     509    echo "       -parallel [mpi|omp|mpi_omp] option!"
     510    exit
     511  fi
    497512  CPP_KEY="$CPP_KEY CPP_1D"
    498513##in 1D, add dyn3d to include path (because main prog is in physics)
  • trunk/LMDZ.TITAN/libf/phytitan/aaam_bud.F

    r102 r1530  
    77c
    88      use dimphy
     9      use mod_grid_phy_lmdz, only: nbp_lon, nbp_lat, klon_glo
    910      implicit none
    1011c======================================================================
     
    5253c ===================
    5354c
    54 c iim--common-I: Number of longitude intervals
    55 c jjm--common-I: Number of latitude intervals
     55c nbp_lon--common-I: Number of longitude intervals
     56c nbp_lat-1--common-I: Number of latitude intervals
    5657c klon-common-I: Number of points seen by the physics
    57 c                iim*(jjm-1)+2 for instance
     58c                nbp_lon*(nbp_lat-1-1)+2 for instance
    5859c klev-common-I: Number of vertical layers
    5960c======================================================================
     
    7879c======================================================================
    7980
    80 #include "dimensions.h"
    8181c
    8282c ARGUMENTS
     
    106106C  PUT AAM QUANTITIES AT ZERO:
    107107C
    108       if(iim+1.gt.801.or.jjm+1.gt.401)then
     108      if(nbp_lon+1.gt.801.or.nbp_lat.gt.401)then
    109109      print *,' Pb de dimension dans aaam_bud'
    110110      stop
     
    114114      hadley=1.e18
    115115      hadday=1.e18*1.e7
    116       dlat=xpi/float(jjm)
    117       dlon=2.*xpi/float(iim)
     116      IF (klon_glo.EQ.1) THEN
     117        dlat=xpi
     118      ELSE
     119        dlat=xpi/float(nbp_lat-1)
     120      ENDIF
     121      dlon=2.*xpi/float(nbp_lon)
    118122     
    119123      do iax=1,3
     
    140144          zlat(1)=plat(l)*xpi/180.
    141145
    142         do i=1,iim+1
     146        do i=1,nbp_lon+1
    143147
    144148          zs(i,1)=phis(l)/rg
     
    154158
    155159
    156       do j = 2,jjm
     160      do j = 2,nbp_lat-1
    157161
    158162C Values at Greenwich (Periodicity)
    159163
    160       zs(iim+1,j)=phis(l+1)/rg
    161       ps(iim+1,j)=p(l+1,1)
    162           ssou(iim+1,j)=dragu(l+1)+liftu(l+1)
    163           ssov(iim+1,j)=dragv(l+1)+liftv(l+1)
    164           blsu(iim+1,j)=clu(l+1)
    165           blsv(iim+1,j)=clv(l+1)
    166       zlon(iim+1)=-plon(l+1)*xpi/180.
     164      zs(nbp_lon+1,j)=phis(l+1)/rg
     165      ps(nbp_lon+1,j)=p(l+1,1)
     166          ssou(nbp_lon+1,j)=dragu(l+1)+liftu(l+1)
     167          ssov(nbp_lon+1,j)=dragv(l+1)+liftv(l+1)
     168          blsu(nbp_lon+1,j)=clu(l+1)
     169          blsv(nbp_lon+1,j)=clv(l+1)
     170      zlon(nbp_lon+1)=-plon(l+1)*xpi/180.
    167171      zlat(j)=plat(l+1)*xpi/180.
    168172
    169       ub(iim+1,j)=0.
    170       vb(iim+1,j)=0.
     173      ub(nbp_lon+1,j)=0.
     174      vb(nbp_lon+1,j)=0.
    171175         do k=1,nlev
    172          ub(iim+1,j)=ub(iim+1,j)+u(l+1,k)*(p(l+1,k)-p(l+1,k+1))/rg
    173          vb(iim+1,j)=vb(iim+1,j)+v(l+1,k)*(p(l+1,k)-p(l+1,k+1))/rg
     176         ub(nbp_lon+1,j)=ub(nbp_lon+1,j)+u(l+1,k)*
     177     &                                   (p(l+1,k)-p(l+1,k+1))/rg
     178         vb(nbp_lon+1,j)=vb(nbp_lon+1,j)+v(l+1,k)*
     179     &                                   (p(l+1,k)-p(l+1,k+1))/rg
    174180         enddo
    175181     
    176182
    177       do i=1,iim
     183      do i=1,nbp_lon
    178184
    179185      l=l+1
     
    201207
    202208      l=l+1
    203       ub(1,jjm+1)=0.
    204       vb(1,jjm+1)=0.
     209      ub(1,nbp_lat)=0.
     210      vb(1,nbp_lat)=0.
    205211      do k=1,nlev
    206          ub(1,jjm+1)=ub(1,jjm+1)+u(l,k)*(p(l,k)-p(l,k+1))/rg
    207          vb(1,jjm+1)=vb(1,jjm+1)+v(l,k)*(p(l,k)-p(l,k+1))/rg
    208       enddo
    209       zlat(jjm+1)=plat(l)*xpi/180.
    210 
    211       do i=1,iim+1
    212       zs(i,jjm+1)=phis(l)/rg
    213       ps(i,jjm+1)=p(l,1)
    214           ssou(i,jjm+1)=dragu(l)+liftu(l)
    215           ssov(i,jjm+1)=dragv(l)+liftv(l)
    216           blsu(i,jjm+1)=clu(l)
    217           blsv(i,jjm+1)=clv(l)
    218       ub(i,jjm+1)=ub(1,jjm+1)                               
    219       vb(i,jjm+1)=vb(1,jjm+1)                               
     212         ub(1,nbp_lat)=ub(1,nbp_lat)+u(l,k)*(p(l,k)-p(l,k+1))/rg
     213         vb(1,nbp_lat)=vb(1,nbp_lat)+v(l,k)*(p(l,k)-p(l,k+1))/rg
     214      enddo
     215      zlat(nbp_lat)=plat(l)*xpi/180.
     216
     217      do i=1,nbp_lon+1
     218      zs(i,nbp_lat)=phis(l)/rg
     219      ps(i,nbp_lat)=p(l,1)
     220          ssou(i,nbp_lat)=dragu(l)+liftu(l)
     221          ssov(i,nbp_lat)=dragv(l)+liftv(l)
     222          blsu(i,nbp_lat)=clu(l)
     223          blsv(i,nbp_lat)=clv(l)
     224      ub(i,nbp_lat)=ub(1,nbp_lat)                               
     225      vb(i,nbp_lat)=vb(1,nbp_lat)                               
    220226      enddo
    221227
     
    223229C  MOMENT ANGULAIRE
    224230C
    225         DO j=1,jjm   
    226         DO i=1,iim
     231        DO j=1,nbp_lat-1   
     232        DO i=1,nbp_lon
    227233
    228234           raam(1)=raam(1)-rea**3*dlon*dlat*0.5*
     
    261267C
    262268
    263         DO j=1,jjm
    264         DO i=1,iim
     269        DO j=1,nbp_lat-1
     270        DO i=1,nbp_lon
    265271           tmou(1)=tmou(1)-rea**2*dlon*0.5*sin(zlon(i))
    266272     c  *(zs(i,j)-zs(i,j+1))
     
    272278        ENDDO
    273279           
    274         DO j=2,jjm
    275         DO i=1,iim
     280        DO j=2,nbp_lat-1
     281        DO i=1,nbp_lon
    276282           tmou(1)=tmou(1)+rea**2*dlat*0.5*sin(zlat(j))
    277283     c  *(zs(i+1,j)-zs(i,j))
     
    288294C
    289295        l=1
    290         DO j=2,jjm
    291         DO i=1,iim
     296        DO j=2,nbp_lat-1
     297        DO i=1,nbp_lon
    292298        l=l+1
    293299           tsso(1)=tsso(1)-rea**3*cos(zlat(j))*dlon*dlat*
     
    335341c00   format(F12.5,5(1x,F12.5))
    336342
    337       write(iam+1,*)((zs(i,j),i=1,iim),j=1,jjm+1)
    338       write(iam+1,*)((ps(i,j),i=1,iim),j=1,jjm+1)
    339       write(iam+1,*)((ub(i,j),i=1,iim),j=1,jjm+1)
    340       write(iam+1,*)((vb(i,j),i=1,iim),j=1,jjm+1)
    341       write(iam+1,*)((ssou(i,j),i=1,iim),j=1,jjm+1)
    342       write(iam+1,*)((ssov(i,j),i=1,iim),j=1,jjm+1)
    343       write(iam+1,*)((blsu(i,j),i=1,iim),j=1,jjm+1)
    344       write(iam+1,*)((blsv(i,j),i=1,iim),j=1,jjm+1)
     343      write(iam+1,*)((zs(i,j),i=1,nbp_lon),j=1,nbp_lat)
     344      write(iam+1,*)((ps(i,j),i=1,nbp_lon),j=1,nbp_lat)
     345      write(iam+1,*)((ub(i,j),i=1,nbp_lon),j=1,nbp_lat)
     346      write(iam+1,*)((vb(i,j),i=1,nbp_lon),j=1,nbp_lat)
     347      write(iam+1,*)((ssou(i,j),i=1,nbp_lon),j=1,nbp_lat)
     348      write(iam+1,*)((ssov(i,j),i=1,nbp_lon),j=1,nbp_lat)
     349      write(iam+1,*)((blsu(i,j),i=1,nbp_lon),j=1,nbp_lat)
     350      write(iam+1,*)((blsv(i,j),i=1,nbp_lon),j=1,nbp_lat)
    345351
    346352      RETURN
  • trunk/LMDZ.TITAN/libf/phytitan/ajsec.F

    r1048 r1530  
    77
    88      use dimphy
     9      use mod_grid_phy_lmdz, only: nbp_lev
    910      use cpdet_mod, only: t2tpot, tpot2t
    1011      IMPLICIT none
     
    2728c d_qfi-----output-R-Incrementation des traceurs
    2829c======================================================================
    29 #include "dimensions.h"
    3030#include "YOMCST.h"
    3131      REAL paprs(klon,klev+1), pplay(klon,klev)
     
    6666c------------------------------------- passage en temperature potentielle
    6767! ADAPTATION GCM POUR CP(T)
    68       call t2tpot(klon*llm,tfi,zh,ppk)
     68      call t2tpot(klon*nbp_lev,tfi,zh,ppk)
    6969c
    7070      DO k = limbas, limhau
     
    168168c------------------------------------- et calcul du d_t
    169169! ADAPTATION GCM POUR CP(T)
    170       call tpot2t(klon*llm,zh,zt,ppk)
     170      call tpot2t(klon*nbp_lev,zh,zt,ppk)
    171171
    172172      DO k = limbas, limhau
  • trunk/LMDZ.TITAN/libf/phytitan/ballon.F

    r102 r1530  
    33
    44      use dimphy
     5      use mod_grid_phy_lmdz, only: nbp_lon, nbp_lat
    56      implicit none
    67
     
    9495c======================================================================
    9596
    96 #include "dimensions.h"
    9797#include "YOMCST.h"
    9898c
     
    110110      INTEGER jj,ii,ll
    111111
    112       REAL zlon(iim+1),zlat(jjm+1)
    113       save zlon,zlat
     112      REAL,SAVE,ALLOCATABLE :: zlon(:),zlat(:)
    114113
    115114      REAL time
     
    141140      print*,"BALLOONS ACTIVATED"
    142141
     142      allocate(zlon(nbp_lon+1))
     143      allocate(zlat(nbp_lat))
     144
    143145C Latitudes:
    144146      zlat(1)=plat(1)*RPI/180.
    145       do j = 2,jjm
    146          k=(j-2)*iim+2
     147      do j = 2,nbp_lat-1
     148         k=(j-2)*nbp_lon+2
    147149         zlat(j)=plat(k)*RPI/180.
    148150      enddo
    149       zlat(jjm+1)=plat(klon)*RPI/180.
     151      zlat(nbp_lat)=plat(klon)*RPI/180.
    150152
    151153C Longitudes:
    152       do i = 1,iim
     154      do i = 1,nbp_lon
    153155         k=i+1
    154156         zlon(i)=plon(k)*RPI/180.
    155157      enddo
    156       zlon(iim+1)=zlon(1)+2.*RPI
     158      zlon(nbp_lon+1)=zlon(1)+2.*RPI
    157159
    158160c verif init     lat de 90 à -90, lon de -180 à 180
     
    256258
    257259      use dimphy
     260      use mod_grid_phy_lmdz, only: nbp_lon, nbp_lat
    258261      implicit none
    259262
     
    283286c======================================================================
    284287
    285 #include "dimensions.h"
    286288#include "YOMCST.h"
    287289c
     
    289291c
    290292      real map_u(klon),map_v(klon),map_a(klon)
    291       real latit(jjm+1),longit(iim)
     293      real latit(nbp_lat),longit(nbp_lon)
    292294      real phi,lam,ubal,vbal,abal
    293295c
     
    296298      INTEGER i,j,k
    297299      INTEGER jj,ii
    298       REAL    ujj(iim+1),vjj(iim+1),ajj(iim+1)
     300      REAL    ujj(nbp_lon+1),vjj(nbp_lon+1),ajj(nbp_lon+1)
    299301      REAL    factlat,factlon
    300302
     
    302304c-------------------------------------------------
    303305        jj=1  ! POLE NORD
    304         do j=2,jjm
     306        do j=2,nbp_lat-1
    305307          if (phi.lt.latit(j)) jj=j
    306308        enddo
     
    309311c pole nord
    310312        if (jj.eq.1) then
    311          do i=1,iim
     313         do i=1,nbp_lon
    312314          ujj(i) = map_u(i+1)*factlat + map_u(1)*(1-factlat)
    313315          vjj(i) = map_v(i+1)*factlat + map_v(1)*(1-factlat)
     
    315317         enddo
    316318c pole sud
    317         elseif (jj.eq.jjm) then
    318          do i=1,iim
    319           k = (jj-2)*iim+1+i
     319        elseif (jj.eq.nbp_lat-1) then
     320         do i=1,nbp_lon
     321          k = (jj-2)*nbp_lon+1+i
    320322          ujj(i) = map_u(klon)*factlat + map_u(k)*(1-factlat)
    321323          vjj(i) = map_v(klon)*factlat + map_v(k)*(1-factlat)
     
    324326c autres latitudes
    325327        else
    326          do i=1,iim
    327           k = (jj-2)*iim+1+i
    328           ujj(i) = map_u(k+iim)*factlat + map_u(k)*(1-factlat)
    329           vjj(i) = map_v(k+iim)*factlat + map_v(k)*(1-factlat)
    330           ajj(i) = map_a(k+iim)*factlat + map_a(k)*(1-factlat)
     328         do i=1,nbp_lon
     329          k = (jj-2)*nbp_lon+1+i
     330          ujj(i) = map_u(k+nbp_lon)*factlat + map_u(k)*(1-factlat)
     331          vjj(i) = map_v(k+nbp_lon)*factlat + map_v(k)*(1-factlat)
     332          ajj(i) = map_a(k+nbp_lon)*factlat + map_a(k)*(1-factlat)
    331333         enddo
    332334        endif
    333         ujj(iim+1)=ujj(1)
    334         vjj(iim+1)=vjj(1)
    335         ajj(iim+1)=ajj(1)
     335        ujj(nbp_lon+1)=ujj(1)
     336        vjj(nbp_lon+1)=vjj(1)
     337        ajj(nbp_lon+1)=ajj(1)
    336338
    337339c Interpolation in longitudes
    338340c-------------------------------------------------
    339341        ii=1  ! lon=-180
    340         do i=2,iim
     342        do i=2,nbp_lon
    341343          if (lam.gt.longit(i)) ii=i
    342344        enddo
     
    370372c======================================================================
    371373
    372 #include "dimensions.h"
    373374#include "YOMCST.h"
    374375c
  • trunk/LMDZ.TITAN/libf/phytitan/calchim.F

    r1379 r1530  
    1919      USE comgeomphy,  only: rlatd
    2020      USE moyzon_mod, only: tmoy,playmoy,zlaymoy,zlevmoy,klat
     21      use mod_grid_phy_lmdz, only: nbp_lat
    2122      implicit none
    22 #include "dimensions.h"
    2323#include "clesphys.h"
    24 #include "paramet.h"
    2524#include "YOMCST.h"
    2625
     
    105104c ************************************
    106105
    107         allocate(krpd(15,ND+1,NLRT,jjp1),krate(NLEV,NR),md(NLEV,NC))
     106        allocate(krpd(15,ND+1,NLRT,nbp_lat),krate(NLEV,NR),md(NLEV,NC))
    108107
    109108c Verification du nombre de composes: coherence common_mod et nqmax-nmicro
     
    265264c taux de photodissociations:
    266265c --------------------------
    267         call disso(krpd,jjp1)
     266        call disso(krpd,nbp_lat)
    268267
    269268c reactions chimiques:
  • trunk/LMDZ.TITAN/libf/phytitan/clmain.F

    r1443 r1530  
    3535      USE interface_surf
    3636      use dimphy
     37      use mod_grid_phy_lmdz, only: nbp_lev
    3738      use cpdet_mod, only: t2tpot
    3839      IMPLICIT none
     
    6667cAA la premiere couche
    6768c======================================================================
    68 #include "dimensions.h"
    6969c$$$ PB ajout pour soil
    7070#include "dimsoil.h"
     
    324324
    325325! ADAPTATION GCM POUR CP(T)
    326          call t2tpot(knon*llm,yt,yteta,ppk)
     326         call t2tpot(knon*nbp_lev,yt,yteta,ppk)
    327327
    328328         yzlev(1:knon,1)=0.
     
    479479      USE interface_surf
    480480      use dimphy
     481      use mod_grid_phy_lmdz, only: nbp_lon, nbp_lat, nbp_lev
    481482      use cpdet_mod, only: t2tpot,tpot2t,cpdet
    482483
     
    486487c Objet: diffusion verticale de "h"
    487488c======================================================================
    488 #include "dimensions.h"
    489489#include "YOMCST.h"
    490490#include "dimsoil.h"
     
    594594
    595595c passage en enthalpie potentielle
    596       call t2tpot(knon*llm,t,local_h,ppk)
     596      call t2tpot(knon*nbp_lev,t,local_h,ppk)
    597597c     print*,"tpot en entree de clqh=",local_h(klon/2,:)
    598598
     
    685685! ADAPTATION GCM POUR CP(T)
    686686      CALL interfsurf_hq(itime, dtime, rmu0,
    687      e klon, iim, jjm, knon,
     687     e klon, nbp_lon, nbp_lat-1, knon,
    688688     e rlon, rlat, cufi, cvfi,
    689689     e debut, lafin, soil_model, nsoilmx,tsoil,
     
    733733      ENDDO
    734734      ENDDO
    735       call tpot2t(knon*llm,local_h,d_t,ppk)
     735      call tpot2t(knon*nbp_lev,local_h,d_t,ppk)
    736736
    737737c     print*,"tpot en sortie de clqh=",local_h(klon/2,:)
     
    781781c flux_v---output-R- (diagnostic) flux du vent: (kg m/s)/(m**2 s)
    782782c======================================================================
    783 #include "dimensions.h"
    784783#include "iniprint.h"
    785784      INTEGER knon
     
    908907c pcfh-----output-R- coefficients a calculer (chaleur et humidite)
    909908c======================================================================
    910 #include "dimensions.h"
    911909#include "YOMCST.h"
    912910#include "iniprint.h"
     
    11721170
    11731171      use dimphy
     1172      use mod_grid_phy_lmdz, only: nbp_lev
    11741173      use cpdet_mod, only: cpdet
    11751174      IMPLICIT none
     
    11881187c pcfh-----output-R- coefficients a calculer (chaleur et humidite)
    11891188c======================================================================
    1190 #include "dimensions.h"
    11911189#include "YOMCST.h"
    11921190#include "iniprint.h"
  • trunk/LMDZ.TITAN/libf/phytitan/cltrac.F

    r102 r1530  
    2727c flux_tr--output-R- flux de tr
    2828c======================================================================
    29 #include "dimensions.h"
    3029      REAL dtime
    3130      REAL coef(klon,klev)
  • trunk/LMDZ.TITAN/libf/phytitan/coefkzmin.F

    r1056 r1530  
    99      IMPLICIT NONE
    1010
    11 #include "dimensions.h"
    1211#include "YOMCST.h"
    1312
  • trunk/LMDZ.TITAN/libf/phytitan/common_mod.F90

    r1525 r1530  
    9090use dimphy
    9191IMPLICIT NONE
    92 #include "dimensions.h"
    9392#include "microtab.h"
    9493
  • trunk/LMDZ.TITAN/libf/phytitan/diagphy.F

    r1048 r1530  
    5151      implicit none
    5252
    53 #include "dimensions.h"
    5453#include "YOMCST.h"
    5554C
     
    212211      IMPLICIT NONE
    213212C
    214 #include "dimensions.h"
    215213#include "YOMCST.h"
    216214C
  • trunk/LMDZ.TITAN/libf/phytitan/drag_noro.F

    r1356 r1530  
    1111      use dimphy
    1212      IMPLICIT none
    13 
    14 #include "dimensions.h"
    15 #include "paramet.h"
    1613
    1714c======================================================================
     
    3027c paprs---input-R-Pressure in semi layers    (Pa)
    3128c pplay---input-R-Pressure model-layers      (Pa)
    32 c pgeop---input-R-Geopotential model layers (reference ground)
     29c pgeop---input-R-Geopotential model layers (reference to ground)
    3330c pn2-----input-R-Brunt-Vaisala freq.^2 at 1/2 layers
    3431c t-------input-R-temperature (K)
  • trunk/LMDZ.TITAN/libf/phytitan/grid_noro.F

    r1356 r1530  
    5252C=======================================================================
    5353
     54      use mod_grid_phy_lmdz, only: nbp_lon, nbp_lat
    5455      IMPLICIT none
    5556     
    56 #include "dimensions.h"
    5757#include "YOMCST.h"
    5858
     
    7979C INTERMEDIATE FIELDS  (CORRELATIONS OF OROGRAPHY GRADIENT)
    8080
    81       REAL ztz(iim+1,jjm+1),zxtzx(iim+1,jjm+1)
    82       REAL zytzy(iim+1,jjm+1),zxtzy(iim+1,jjm+1)
    83       REAL weight(iim+1,jjm+1)
     81      REAL ztz(nbp_lon+1,nbp_lat),zxtzx(nbp_lon+1,nbp_lat)
     82      REAL zytzy(nbp_lon+1,nbp_lat),zxtzy(nbp_lon+1,nbp_lat)
     83      REAL weight(nbp_lon+1,nbp_lat)
    8484
    8585C CORRELATIONS OF USN OROGRAPHY GRADIENTS
     
    107107c   
    108108c
    109       if(iim.ne.imar) STOP 'Problem dim. x'
    110       if(jjm.ne.jmar-1) STOP 'Problem dim. y'
     109      if(nbp_lon.ne.imar) STOP 'Problem dim. x'
     110      if(nbp_lat-1.ne.jmar-1) STOP 'Problem dim. y'
    111111      IF (imar.GT.2200 .OR. jmar.GT.1100) THEN
    112112         PRINT*, 'imar or jmar too big', imar, jmar
     
    114114      ENDIF
    115115
    116       IF(imar+1.ne.iim+1.or.jmar.ne.jjm+1)THEN
     116      IF(imar+1.ne.nbp_lon+1.or.jmar.ne.nbp_lat)THEN
    117117        print *,' imar or jmar bad dimensions:',imar,jmar
    118118        call abort
     
    297297C  FIRST FILTER, MOVING AVERAGE OVER 9 POINTS.
    298298
    299        CALL MVA9(zmea,iim+1,jjm+1)
    300        CALL MVA9(zstd,iim+1,jjm+1)
    301        CALL MVA9(zpic,iim+1,jjm+1)
    302        CALL MVA9(zval,iim+1,jjm+1)
    303        CALL MVA9(zxtzx,iim+1,jjm+1)
    304        CALL MVA9(zxtzy,iim+1,jjm+1)
    305        CALL MVA9(zytzy,iim+1,jjm+1)
     299       CALL MVA9(zmea,nbp_lon+1,nbp_lat)
     300       CALL MVA9(zstd,nbp_lon+1,nbp_lat)
     301       CALL MVA9(zpic,nbp_lon+1,nbp_lat)
     302       CALL MVA9(zval,nbp_lon+1,nbp_lat)
     303       CALL MVA9(zxtzx,nbp_lon+1,nbp_lat)
     304       CALL MVA9(zxtzy,nbp_lon+1,nbp_lat)
     305       CALL MVA9(zytzy,nbp_lon+1,nbp_lat)
    306306
    307307       DO ii = 1, imar
  • trunk/LMDZ.TITAN/libf/phytitan/gwprofil.F

    r1056 r1530  
    3838      use dimphy
    3939      IMPLICIT NONE
    40 
    41 #include "dimensions.h"
    42 #include "paramet.h"
    4340
    4441#include "YOMCST.h"
  • trunk/LMDZ.TITAN/libf/phytitan/gwstress.F

    r1056 r1530  
    5454      use dimphy
    5555      implicit none
    56 
    57 #include "dimensions.h"
    58 #include "paramet.h"
    5956
    6057#include "YOMCST.h"
  • trunk/LMDZ.TITAN/libf/phytitan/hgardfou.F

    r102 r1530  
    33!
    44      SUBROUTINE hgardfou (t,tsol,text)
     5
     6      use dimphy
     7      IMPLICIT none
    58c======================================================================
    69c Verifier la temperature
    710c======================================================================
    8       use dimphy
    9       IMPLICIT none
    10 #include "dimensions.h"
    1111#include "YOMCST.h"
    1212      REAL t(klon,klev), tsol(klon)
     
    8484           ok = .FALSE.
    8585           DO i = 1, jbad
    86              PRINT *,'i,temperature sol =',jadrs(i),zt(jadrs(i))
     86             PRINT *,'i,temperature =',jadrs(i),zt(jadrs(i))
    8787           ENDDO
    8888         ENDIF
     
    101101           ok = .FALSE.
    102102           DO i = 1, jbad
    103              PRINT *,'i,temperature sol =',jadrs(i),zt(jadrs(i))
     103             PRINT *,'i,temperature =',jadrs(i),zt(jadrs(i))
    104104           ENDDO
    105105         ENDIF
  • trunk/LMDZ.TITAN/libf/phytitan/ini_histday.h

    r1056 r1530  
    2222
    2323         CALL histdef(nid_day, "phis", "Surface geop. height", "-",
    24      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     24     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    2525     .                "once",  zsto,zout)
    2626
    2727         CALL histdef(nid_day, "aire", "Grid area", "-",
    28      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     28     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    2929     .                "once",  zsto,zout)
    3030
    3131ccccccc axe Ls
    3232         CALL histdef(nid_day, "ls", "Solar longitude", "degrees",
    33      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     33     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    3434     .                "ave(X)", zsto,zout)
    3535
     
    3737
    3838         CALL histdef(nid_day, "tsol", "Surface Temperature", "K",
    39      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     39     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    4040     .                "ave(X)", zsto,zout)
    4141
    4242         CALL histdef(nid_day, "psol", "Surface Pressure", "Pa",
    43      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     43     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    4444     .                "ave(X)", zsto,zout)
    4545
    4646c        CALL histdef(nid_day, "ue", "Zonal energy transport", "-",
    47 c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     47c    .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    4848c    .                "ave(X)", zsto,zout)
    4949
    5050c        CALL histdef(nid_day, "ve", "Merid energy transport", "-",
    51 c     .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     51c     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    5252c     .                "ave(X)", zsto,zout)
    5353
     
    6060
    6161         CALL histdef(nid_day, "temp", "Air temperature", "K",
    62      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     62     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    6363     .                "ave(X)", zsto,zout)
    6464
    6565         CALL histdef(nid_day, "pres", "Air pressure", "Pa",
    66      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     66     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    6767     .                "ave(X)", zsto,zout)
    6868
    6969         CALL histdef(nid_day, "geop", "Geopotential height", "m",
    70      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     70     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    7171     .                "ave(X)", zsto,zout)
    7272
    7373         CALL histdef(nid_day, "vitu", "Zonal wind", "m/s",
    74      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     74     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    7575     .                "ave(X)", zsto,zout)
    7676
    7777         CALL histdef(nid_day, "vitv", "Meridional wind", "m/s",
    78      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     78     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    7979     .                "ave(X)", zsto,zout)
    8080
    8181         CALL histdef(nid_day, "vitw", "Vertical wind", "Pa/s",
    82      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     82     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    8383     .                "ave(X)", zsto,zout)
    8484
    8585         CALL histdef(nid_day, "tops", "Solar rad. at TOA", "W/m2",
    86      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     86     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    8787     .                "ave(X)", zsto1,zout)
    8888
    8989         CALL histdef(nid_day, "duvdf", "Boundary-layer dU", "m/s2",
    90      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     90     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    9191     .                "ave(X)", zsto,zout)
    9292
    9393         CALL histdef(nid_day, "dudyn", "Dynamics dU", "m/s2",
    94      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     94     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    9595     .                "ave(X)", zsto,zout)
    9696
     
    101101c           DO iq=1,nmicro
    102102c             CALL histdef(nid_day, tname(iq), ttext(iq), "n/m2",
    103 c     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     103c     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    104104c     .                "ave(X)", zsto,zout)
    105105c           ENDDO
    106106             CALL histdef(nid_day, "qaer","nb tot aer" , "n/m2",
    107      .                    iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     107     .                    nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    108108     .                    "ave(X)", zsto,zout)
    109109
    110110            if (clouds.eq.1) then
    111111             CALL histdef(nid_day, "qnoy","nb tot noy" , "n/m2",
    112      .                    iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     112     .                    nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    113113     .                    "ave(X)", zsto,zout)
    114114             CALL histdef(nid_day, "qgl1","V tot gl1" , "m3/m2",
    115      .                    iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     115     .                    nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    116116     .                    "ave(X)", zsto,zout)
    117117             CALL histdef(nid_day, "qgl2","V tot gl2" , "m3/m2",
    118      .                    iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     118     .                    nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    119119     .                    "ave(X)", zsto,zout)
    120120             CALL histdef(nid_day, "qgl3","V tot gl3" , "m3/m2",
    121      .                    iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     121     .                    nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    122122     .                    "ave(X)", zsto,zout)
    123123c--------------
    124124c ----- SATURATION ESP NUAGES
    125125               CALL histdef(nid_day,"ch4sat", "saturation CH4", "--",
    126      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     126     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    127127     .                "ave(X)", zsto,zout)
    128128               CALL histdef(nid_day,"c2h6sat", "saturation C2H6", "--",
    129      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     129     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    130130     .                "ave(X)", zsto,zout)
    131131               CALL histdef(nid_day,"c2h2sat", "saturation C2H2", "--",
    132      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     132     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    133133     .                "ave(X)", zsto,zout)
    134134c --------------
    135135c ----- RESERVOIR DE SURFACE
    136136               CALL histdef(nid_day, "reserv", "Reservoir surface","m",
    137      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     137     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    138138     .                "ave(X)", zsto,zout)
    139139c --------------
    140140c ----- ECHANGE GAZ SURF/ATM (evaporation)
    141141               CALL histdef(nid_day, "evapch4", "Evaporation CH4","m",
    142      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     142     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    143143     .                "ave(X)", zsto,zout)
    144144c --------------
    145145c ----- PRECIPITATIONS (precipitations cumulatives)
    146146               CALL histdef(nid_day,"prech4","Precip CH4","m",
    147      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     147     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    148148     .                "ave(X)", zsto,zout)
    149149               CALL histdef(nid_day,"prec2h6","Precip C2H6",
    150      .                "m",iim,jj_nb,nhori, 1,1,1, nvert, 32,
     150     .                "m",nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    151151     .                "ave(X)", zsto,zout)
    152152               CALL histdef(nid_day,"prec2h2","Precip C2H2",
    153      .                "m",iim,jj_nb,nhori, 1,1,1, nvert, 32,
     153     .                "m",nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    154154     .                "ave(X)", zsto,zout)
    155155               CALL histdef(nid_day,"prenoy","Precip NOY",
    156      .                "um/s",iim,jj_nb,nhori, 1,1,1, nvert, 32,
     156     .                "um/s",nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    157157     .                "ave(X)", zsto,zout)
    158158               CALL histdef(nid_day,"preaer","Precip AER",
    159      .                "um/s",iim,jj_nb,nhori, 1,1,1, nvert, 32,
     159     .                "um/s",nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    160160     .                "ave(X)", zsto,zout)
    161161c --------------
    162162c ----- FLUX GLACE
    163163               CALL histdef(nid_day,"flxgl1", "flux gl CH4",
    164      .              "kg/m2/s",iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     164     .              "kg/m2/s",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    165165     .              "ave(X)", zsto,zout)
    166166               CALL histdef(nid_day,"flxgl2", "flux gl C2H6",
    167      .              "kg/m2/s",iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     167     .              "kg/m2/s",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    168168     .              "ave(X)", zsto,zout)
    169169               CALL histdef(nid_day,"flxgl3", "flux gl C2H2",
    170      .              "kg/m2/s",iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     170     .              "kg/m2/s",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    171171     .              "ave(X)", zsto,zout)
    172172c --------------
    173173c ----- RAYON DES GOUTTES
    174174               CALL histdef(nid_day,"rcldbar", "rayon moyen goutte",
    175      .                "m",iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     175     .                "m",nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    176176     .                "ave(X)", zsto,zout)
    177177            endif
     
    182182           DO iq=nmicro+1,nqmax
    183183         CALL histdef(nid_day, tname(iq), ttext(iq), "ppm",
    184      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     184     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    185185     .                "ave(X)", zsto,zout)
    186186           ENDDO
     
    198198
    199199         CALL histdef(nid_day, "topl", "IR rad. at TOA", "W/m2",
    200      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     200     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    201201     .                "ave(X)", zsto1,zout)
    202202
    203203         CALL histdef(nid_day, "sols", "Solar rad. at surf.", "W/m2",
    204      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     204     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    205205     .                "ave(X)", zsto1,zout)
    206206
    207207         CALL histdef(nid_day, "soll", "IR rad. at surface", "W/m2",
    208      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     208     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    209209     .                "ave(X)", zsto1,zout)
    210210
     
    212212
    213213         CALL histdef(nid_day, "SWnet", "Net SW flux","W/m2",
    214      .                iim,jj_nb,nhori, klev,1,klev,nvert,
     214     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
    215215     .                32, "ave(X)", zsto1,zout)
    216216
    217217         CALL histdef(nid_day, "LWnet", "Net LW flux","W/m2",
    218      .                iim,jj_nb,nhori, klev,1,klev,nvert,
     218     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
    219219     .                32, "ave(X)", zsto1,zout)
    220220
     
    224224           write(str2,'(i2.2)') k
    225225         CALL histdef(nid_day,"thv"//str2,"Haze Opa Vis",
    226      .                "--",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     226     .                "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    227227     .                "ave(X)",zsto1,zout)
    228228         ENDDO
     
    231231           write(str2,'(i2.2)') k
    232232         CALL histdef(nid_day,"thi"//str2,"Haze Opa IR",
    233      .                "--",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     233     .                "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    234234     .                "ave(X)",zsto1,zout)
    235235         ENDDO
     
    240240           write(str2,'(i2.2)') k
    241241         CALL histdef(nid_day,"khv"//str2,"Haze ext Vis ",
    242      .                "m-1",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     242     .                "m-1",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    243243     .                "ave(X)",zsto1,zout)
    244244         ENDDO
     
    247247           write(str2,'(i2.2)') k
    248248         CALL histdef(nid_day,"khi"//str2,"Haze ext IR ",
    249      .                "m-1",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     249     .                "m-1",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    250250     .                "ave(X)",zsto1,zout)
    251251         ENDDO
     
    256256           write(str2,'(i2.2)') k
    257257         CALL histdef(nid_day,"tgv"//str2,"Gas Opa Vis",
    258      .                "--",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     258     .                "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    259259     .                "ave(X)",zsto1,zout)
    260260         ENDDO
     
    263263           write(str2,'(i2.2)') k
    264264         CALL histdef(nid_day,"tgi"//str2,"Gas Opa IR",
    265      .                "--",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     265     .                "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    266266     .                "ave(X)",zsto1,zout)
    267267         ENDDO
     
    272272           write(str2,'(i2.2)') k
    273273         CALL histdef(nid_day,"kgv"//str2,"Gas ext Vis ",
    274      .                "m-1",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     274     .                "m-1",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    275275     .                "ave(X)",zsto1,zout)
    276276         ENDDO
     
    279279           write(str2,'(i2.2)') k
    280280         CALL histdef(nid_day,"kgi"//str2,"Gas ext IR ",
    281      .                "m-1",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     281     .                "m-1",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    282282     .                "ave(X)",zsto1,zout)
    283283         ENDDO
     
    287287         if (clouds.eq.1) then
    288288           CALL histdef(nid_day,"tcld","Cld Opa proxy",
    289      .                "--",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     289     .                "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    290290     .                "ave(X)",zsto,zout)
    291291
     
    293293c ----- EXTINCTION NUAGES
    294294           CALL histdef(nid_day,"kcld","Cld Ext proxy",
    295      .                "m-1",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     295     .                "m-1",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    296296     .                "ave(X)",zsto,zout)
    297297         endif
     
    303303
    304304         CALL histdef(nid_day, "dtdyn", "Dynamics dT", "K/s",
    305      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     305     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    306306     .                "ave(X)", zsto,zout)
    307307
    308308         CALL histdef(nid_day, "dtphy", "Physics dT", "K/s",
    309      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     309     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    310310     .                "ave(X)", zsto,zout)
    311311
    312312         CALL histdef(nid_day, "dtvdf", "Boundary-layer dT", "K/s",
    313      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     313     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    314314     .                "ave(X)", zsto,zout)
    315315
    316316         CALL histdef(nid_day, "dtajs", "Dry adjust. dT", "K/s",
    317      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     317     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    318318     .                "ave(X)", zsto,zout)
    319319
    320320         CALL histdef(nid_day, "dtswr", "SW radiation dT", "K/s",
    321      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     321     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    322322     .                "ave(X)", zsto,zout)
    323323
    324324         CALL histdef(nid_day, "dtlwr", "LW radiation dT", "K/s",
    325      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     325     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    326326     .                "ave(X)", zsto,zout)
    327327
    328328c        CALL histdef(nid_day, "dtec", "Cinetic dissip dT", "K/s",
    329 c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     329c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    330330c    .                "ave(X)", zsto,zout)
    331331
     
    338338c        call histdef(nid_day, "taux",
    339339c    $         "Zonal wind stress", "Pa", 
    340 c    $         iim,jj_nb,nhori, 1,1,1, nvert, 32,
     340c    $         nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    341341c    $         "ave(X)", zsto,zout)
    342342
    343343c        call histdef(nid_day, "tauy",
    344344c    $         "Meridional xind stress", "Pa", 
    345 c    $         iim,jj_nb,nhori, 1,1,1, nvert, 32,
     345c    $         nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    346346c    $         "ave(X)", zsto,zout)
    347347
    348348c        CALL histdef(nid_day, "cdrm", "Momentum drag coef.", "-",
    349 c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     349c    .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    350350c    .                "ave(X)", zsto,zout)
    351351
    352352c        CALL histdef(nid_day, "cdrh", "Heat drag coef.", "-",
    353 c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     353c    .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    354354c    .                "ave(X)", zsto,zout)
    355355
  • trunk/LMDZ.TITAN/libf/phytitan/ini_histins.h

    r1056 r1530  
    2121
    2222         CALL histdef(nid_ins, "phis", "Surface geop. height", "-",
    23      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     23     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    2424     .                "once",  zsto,zout)
    2525
    2626         CALL histdef(nid_ins, "aire", "Grid area", "-",
    27      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     27     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    2828     .                "once",  zsto,zout)
    2929
    3030ccccccc axe Ls
    3131         CALL histdef(nid_ins, "ls", "Solar longitude", "degrees",
    32      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     32     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    3333     .                "inst(X)", zsto,zout)
    3434
     
    3636
    3737         CALL histdef(nid_ins, "tsol", "Surface Temperature", "K",
    38      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     38     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    3939     .                "inst(X)", zsto,zout)
    4040
    4141         CALL histdef(nid_ins, "psol", "Surface Pressure", "Pa",
    42      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     42     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    4343     .                "inst(X)", zsto,zout)
    4444
    4545c        CALL histdef(nid_ins, "ue", "Zonal energy transport", "-",
    46 c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     46c    .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    4747c    .                "inst(X)", zsto,zout)
    4848
    4949c        CALL histdef(nid_ins, "ve", "Merid energy transport", "-",
    50 c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     50c    .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    5151c    .                "inst(X)", zsto,zout)
    5252
     
    5959
    6060         CALL histdef(nid_ins, "temp", "Air temperature", "K",
    61      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     61     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    6262     .                "inst(X)", zsto,zout)
    6363
    6464         CALL histdef(nid_ins, "pres", "Air pressure", "Pa",
    65      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     65     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    6666     .                "inst(X)", zsto,zout)
    6767
    6868         CALL histdef(nid_ins, "geop", "Geopotential height", "m",
    69      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     69     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    7070     .                "inst(X)", zsto,zout)
    7171
    7272         CALL histdef(nid_ins, "vitu", "Zonal wind", "m/s",
    73      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     73     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    7474     .                "inst(X)", zsto,zout)
    7575
    7676         CALL histdef(nid_ins, "vitv", "Meridional wind", "m/s",
    77      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     77     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    7878     .                "inst(X)", zsto,zout)
    7979
    8080         CALL histdef(nid_ins, "vitw", "Vertical wind", "Pa/s",
    81      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     81     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    8282     .                "inst(X)", zsto,zout)
    8383
    8484         CALL histdef(nid_ins, "tops", "Solar rad. at TOA", "W/m2",
    85      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     85     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    8686     .                "inst(X)", zsto,zout)
    8787
    8888c        CALL histdef(nid_ins, "duvdf", "Boundary-layer dU", "m/s2",
    89 c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     89c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    9090c    .                "inst(X)", zsto,zout)
    9191
    9292c        CALL histdef(nid_ins, "dudyn", "Dynamics dU", "m/s2",
    93 c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     93c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    9494c    .                "inst(X)", zsto,zout)
    9595
     
    105105           DO iq=1,nmicro
    106106         CALL histdef(nid_ins, tname(iq), ttext(iq), "n/m2",
    107      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     107     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    108108     .                "inst(X)", zsto,zout)
    109109           ENDDO
     
    112112           DO iq=nmicro+1,nqmax
    113113         CALL histdef(nid_ins, tname(iq), ttext(iq), "ppm",
    114      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     114     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    115115     .                "inst(X)", zsto,zout)
    116116           ENDDO
     
    123123
    124124         CALL histdef(nid_ins, "topl", "IR rad. at TOA", "W/m2",
    125      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     125     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    126126     .                "inst(X)", zsto,zout)
    127127
    128128         CALL histdef(nid_ins, "sols", "Solar rad. at surf.", "W/m2",
    129      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     129     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    130130     .                "inst(X)", zsto,zout)
    131131
    132132         CALL histdef(nid_ins, "soll", "IR rad. at surface", "W/m2",
    133      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     133     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    134134     .                "inst(X)", zsto,zout)
    135135
     
    137137
    138138         CALL histdef(nid_ins, "SWnet", "Net SW flux","W/m2",
    139      .                iim,jj_nb,nhori, klev,1,klev,nvert,
     139     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
    140140     .                32, "inst(X)", zsto,zout)
    141141
    142142         CALL histdef(nid_ins, "LWnet", "Net LW flux","W/m2",
    143      .                iim,jj_nb,nhori, klev,1,klev,nvert,
     143     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
    144144     .                32, "inst(X)", zsto,zout)
    145145
     
    149149           write(str2,'(i2.2)') k
    150150         CALL histdef(nid_ins,"thv"//str2,"Haze Opa Vis",
    151      .                "--",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     151     .                "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    152152     .                "ins(X)",zsto,zout)
    153153         ENDDO
     
    156156           write(str2,'(i2.2)') k
    157157         CALL histdef(nid_ins,"thi"//str2,"Haze Opa IR",
    158      .                "--",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     158     .                "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    159159     .                "ins(X)",zsto,zout)
    160160         ENDDO
     
    165165           write(str2,'(i2.2)') k
    166166         CALL histdef(nid_ins,"khv"//str2,"Haze ext Vis ",
    167      .                "m-1",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     167     .                "m-1",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    168168     .                "ins(X)",zsto,zout)
    169169         ENDDO
     
    172172           write(str2,'(i2.2)') k
    173173         CALL histdef(nid_ins,"khi"//str2,"Haze ext IR ",
    174      .                "m-1",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     174     .                "m-1",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    175175     .                "ins(X)",zsto,zout)
    176176         ENDDO
     
    181181           write(str2,'(i2.2)') k
    182182         CALL histdef(nid_ins,"tgv"//str2,"Haze Opa Vis",
    183      .                "--",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     183     .                "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    184184     .                "ins(X)",zsto,zout)
    185185         ENDDO
     
    188188           write(str2,'(i2.2)') k
    189189         CALL histdef(nid_ins,"tgi"//str2,"Haze Opa IR",
    190      .                "--",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     190     .                "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    191191     .                "ins(X)",zsto,zout)
    192192         ENDDO
     
    197197           write(str2,'(i2.2)') k
    198198         CALL histdef(nid_ins,"kgv"//str2,"Haze ext Vis ",
    199      .                "m-1",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     199     .                "m-1",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    200200     .                "ins(X)",zsto,zout)
    201201         ENDDO
     
    204204           write(str2,'(i2.2)') k
    205205         CALL histdef(nid_ins,"kgi"//str2,"Haze ext IR ",
    206      .                "m-1",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     206     .                "m-1",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    207207     .                "ins(X)",zsto,zout)
    208208         ENDDO
     
    214214
    215215         CALL histdef(nid_ins, "dtdyn", "Dynamics dT", "K/s",
    216      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     216     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    217217     .                "inst(X)", zsto,zout)
    218218
    219219         CALL histdef(nid_ins, "dtphy", "Physics dT", "K/s",
    220      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     220     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    221221     .                "inst(X)", zsto,zout)
    222222
    223223         CALL histdef(nid_ins, "dtvdf", "Boundary-layer dT", "K/s",
    224      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     224     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    225225     .                "inst(X)", zsto,zout)
    226226
    227227         CALL histdef(nid_ins, "dtajs", "Dry adjust. dT", "K/s",
    228      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     228     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    229229     .                "inst(X)", zsto,zout)
    230230
    231231         CALL histdef(nid_ins, "dtswr", "SW radiation dT", "K/s",
    232      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     232     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    233233     .                "inst(X)", zsto,zout)
    234234
    235235         CALL histdef(nid_ins, "dtlwr", "LW radiation dT", "K/s",
    236      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     236     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    237237     .                "inst(X)", zsto,zout)
    238238
    239239c        CALL histdef(nid_ins, "dtec", "Cinetic dissip dT", "K/s",
    240 c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     240c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    241241c    .                "inst(X)", zsto,zout)
    242242
    243243c        CALL histdef(nid_ins, "dvvdf", "Boundary-layer dV", "m/s2",
    244 c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     244c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    245245c    .                "inst(X)", zsto,zout)
    246246
     
    253253c        call histdef(nid_ins, "taux",
    254254c    $         "Zonal wind stress", "Pa", 
    255 c    $         iim,jj_nb,nhori, 1,1,1, nvert, 32,
     255c    $         nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    256256c    $         "inst(X)", zsto,zout)
    257257
    258258c        call histdef(nid_ins, "tauy",
    259259c    $         "Meridional xind stress", "Pa", 
    260 c    $         iim,jj_nb,nhori, 1,1,1, nvert, 32,
     260c    $         nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    261261c    $         "inst(X)", zsto,zout)
    262262
    263263c        CALL histdef(nid_ins, "cdrm", "Momentum drag coef.", "-",
    264 c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     264c    .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    265265c    .                "inst(X)", zsto,zout)
    266266
    267267c        CALL histdef(nid_ins, "cdrh", "Heat drag coef.", "-",
    268 c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     268c    .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    269269c    .                "inst(X)", zsto,zout)
    270270
  • trunk/LMDZ.TITAN/libf/phytitan/ini_histmth.h

    r1356 r1530  
    2222
    2323         CALL histdef(nid_mth, "phis", "Surface geop. height", "-",
    24      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     24     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    2525     .                "once",  zsto,zout)
    2626
    2727         CALL histdef(nid_mth, "aire", "Grid area", "-",
    28      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     28     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    2929     .                "once",  zsto,zout)
    3030
    3131ccccccc axe Ls
    3232         CALL histdef(nid_mth, "ls", "Solar longitude", "degrees",
    33      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     33     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    3434     .                "ave(X)", zsto,zout)
    3535
     
    3737
    3838         CALL histdef(nid_mth, "tsol", "Surface Temperature", "K",
    39      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     39     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    4040     .                "ave(X)", zsto,zout)
    4141
    4242         CALL histdef(nid_mth, "psol", "Surface Pressure", "Pa",
    43      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     43     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    4444     .                "ave(X)", zsto,zout)
    4545
    4646c        CALL histdef(nid_mth, "ue", "Zonal energy transport", "-",
    47 c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     47c    .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    4848c    .                "ave(X)", zsto,zout)
    4949
    5050c        CALL histdef(nid_mth, "ve", "Merid energy transport", "-",
    51 c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     51c    .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    5252c    .                "ave(X)", zsto,zout)
    5353
     
    6060
    6161         CALL histdef(nid_mth, "temp", "Air temperature", "K",
    62      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     62     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    6363     .                "ave(X)", zsto,zout)
    6464
    6565         CALL histdef(nid_mth, "pres", "Air pressure", "Pa",
    66      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     66     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    6767     .                "ave(X)", zsto,zout)
    6868
    6969         CALL histdef(nid_mth, "geop", "Geopotential height", "m",
    70      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     70     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    7171     .                "ave(X)", zsto,zout)
    7272
    7373         CALL histdef(nid_mth, "vitu", "Zonal wind", "m/s",
    74      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     74     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    7575     .                "ave(X)", zsto,zout)
    7676
    7777         CALL histdef(nid_mth, "vitv", "Meridional wind", "m/s",
    78      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     78     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    7979     .                "ave(X)", zsto,zout)
    8080
    8181         CALL histdef(nid_mth, "vitw", "Vertical wind", "Pa/s",
    82      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     82     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    8383     .                "ave(X)", zsto,zout)
    8484
    8585c        CALL histdef(nid_mth, "Kz", "vertical diffusion coef", "m2/s",
    86 c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     86c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    8787c    .                "ave(X)", zsto,zout)
    8888
    8989         CALL histdef(nid_mth, "tops", "Solar rad. at TOA", "W/m2",
    90      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     90     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    9191     .                "ave(X)", zsto1,zout)
    9292
    9393         CALL histdef(nid_mth, "duvdf", "Boundary-layer dU", "m/s2",
    94      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     94     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    9595     .                "ave(X)", zsto,zout)
    9696
    9797         CALL histdef(nid_mth, "dudyn", "Dynamics dU", "m/s2",
    98      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     98     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    9999     .                "ave(X)", zsto,zout)
    100100
     
    105105c           DO iq=1,nmicro
    106106c             CALL histdef(nid_mth, tname(iq), ttext(iq), "n/m2",
    107 c     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     107c     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    108108c     .                "ave(X)", zsto,zout)
    109109c           ENDDO
    110110             CALL histdef(nid_mth, "qaer","nb tot aer" , "n/m2",
    111      .                    iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     111     .                    nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    112112     .                    "ave(X)", zsto,zout)
    113113
    114114            if (clouds.eq.1) then
    115115             CALL histdef(nid_mth, "qnoy","nb tot noy" , "n/m2",
    116      .                    iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     116     .                    nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    117117     .                    "ave(X)", zsto,zout)
    118118             CALL histdef(nid_mth, "qgl1","V tot gl1" , "m3/m2",
    119      .                    iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     119     .                    nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    120120     .                    "ave(X)", zsto,zout)
    121121             CALL histdef(nid_mth, "qgl2","V tot gl2" , "m3/m2",
    122      .                    iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     122     .                    nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    123123     .                    "ave(X)", zsto,zout)
    124124             CALL histdef(nid_mth, "qgl3","V tot gl3" , "m3/m2",
    125      .                    iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     125     .                    nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    126126     .                    "ave(X)", zsto,zout)
    127127c--------------
    128128c ----- SATURATION ESP NUAGES
    129129               CALL histdef(nid_mth,"ch4sat", "saturation CH4", "--",
    130      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     130     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    131131     .                "ave(X)", zsto,zout)
    132132               CALL histdef(nid_mth,"c2h6sat", "saturation C2H6", "--",
    133      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     133     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    134134     .                "ave(X)", zsto,zout)
    135135               CALL histdef(nid_mth,"c2h2sat", "saturation C2H2", "--",
    136      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     136     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    137137     .                "ave(X)", zsto,zout)
    138138c --------------
    139139c ----- RESERVOIR DE SURFACE
    140140               CALL histdef(nid_mth, "reserv", "Reservoir surface","m",
    141      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     141     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    142142     .                "ave(X)", zsto,zout)
    143143c --------------
    144144c ----- ECHANGE GAZ SURF/ATM (evaporation)
    145145               CALL histdef(nid_mth, "evapch4", "Evaporation CH4","m",
    146      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     146     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    147147     .                "ave(X)", zsto,zout)
    148148c --------------
    149149c ----- PRECIPITATIONS (precipitations moyennes)
    150150               CALL histdef(nid_mth,"prech4","Precip CH4","um/s",
    151      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     151     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    152152     .                "ave(X)", zsto,zout)
    153153               CALL histdef(nid_mth,"prec2h6","Precip C2H6",
    154      .                "um/s",iim,jj_nb,nhori, 1,1,1, nvert, 32,
     154     .                "um/s",nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    155155     .                "ave(X)", zsto,zout)
    156156               CALL histdef(nid_mth,"prec2h2","Precip C2H2",
    157      .                "um/s",iim,jj_nb,nhori, 1,1,1, nvert, 32,
     157     .                "um/s",nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    158158     .                "ave(X)", zsto,zout)
    159159               CALL histdef(nid_mth,"prenoy","Precip NOY",
    160      .                "um/s",iim,jj_nb,nhori, 1,1,1, nvert, 32,
     160     .                "um/s",nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    161161     .                "ave(X)", zsto,zout)
    162162               CALL histdef(nid_mth,"preaer","Precip AER",
    163      .                "um/s",iim,jj_nb,nhori, 1,1,1, nvert, 32,
     163     .                "um/s",nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    164164     .                "ave(X)", zsto,zout)
    165165c --------------
    166166c ----- FLUX GLACE
    167167               CALL histdef(nid_mth,"flxgl1", "flux gl CH4",
    168      .              "kg/m2/s",iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     168     .              "kg/m2/s",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    169169     .              "ave(X)", zsto,zout)
    170170               CALL histdef(nid_mth,"flxgl2", "flux gl C2H6",
    171      .              "kg/m2/s",iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     171     .              "kg/m2/s",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    172172     .              "ave(X)", zsto,zout)
    173173               CALL histdef(nid_mth,"flxgl3", "flux gl C2H2",
    174      .              "kg/m2/s",iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     174     .              "kg/m2/s",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    175175     .              "ave(X)", zsto,zout)
    176176c --------------
    177177c ----- Source/puits GLACE
    178178               CALL histdef(nid_mth,"solch4", "dQ gl CH4",
    179      .              "m3/m3",iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     179     .              "m3/m3",nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    180180     .              "ave(X)", zsto,zout)
    181181               CALL histdef(nid_mth,"solc2h6", "dQ gl C2H6",
    182      .              "m3/m3",iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     182     .              "m3/m3",nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    183183     .              "ave(X)", zsto,zout)
    184184               CALL histdef(nid_mth,"solc2h2", "dQ gl C2H2",
    185      .              "m3/m3",iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     185     .              "m3/m3",nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    186186     .              "ave(X)", zsto,zout)
    187187c --------------
    188188c ----- RAYON DES GOUTTES
    189189               CALL histdef(nid_mth,"rcldbar", "rayon moyen goutte",
    190      .                "m",iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     190     .                "m",nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    191191     .                "ave(X)", zsto,zout)
    192192            endif
     
    197197           DO iq=nmicro+1,nqmax
    198198         CALL histdef(nid_mth, tname(iq), ttext(iq), "ppm",
    199      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     199     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    200200     .                "ave(X)", zsto,zout)
    201201           ENDDO
     
    203203c          DO iq=nmicro+1,nqmax
    204204c        CALL histdef(nid_mth, "c_"//tname(iq), "c_"//ttext(iq),
    205 c    .        "ppm/s",iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     205c    .        "ppm/s",nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    206206c    .                "ave(X)", zsto,zout)
    207207c          ENDDO
     
    219219
    220220         CALL histdef(nid_mth, "topl", "IR rad. at TOA", "W/m2",
    221      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     221     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    222222     .                "ave(X)", zsto1,zout)
    223223
    224224         CALL histdef(nid_mth, "sols", "Solar rad. at surf.", "W/m2",
    225      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     225     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    226226     .                "ave(X)", zsto1,zout)
    227227
    228228         CALL histdef(nid_mth, "soll", "IR rad. at surface", "W/m2",
    229      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     229     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    230230     .                "ave(X)", zsto1,zout)
    231231
     
    233233
    234234         CALL histdef(nid_mth, "SWnet", "Net SW flux","W/m2",
    235      .                iim,jj_nb,nhori, klev,1,klev,nvert,
     235     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
    236236     .                32, "ave(X)", zsto1,zout)
    237237
    238238c        CALL histdef(nid_mth, "SWup", "upward SW flux","W/m2",
    239 c    .                iim,jj_nb,nhori, klev,1,klev,nvert,
     239c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
    240240c    .                32, "ave(X)", zsto1,zout)
    241241
    242242c        CALL histdef(nid_mth, "SWdn", "downward SW flux","W/m2",
    243 c    .                iim,jj_nb,nhori, klev,1,klev,nvert,
     243c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
    244244c    .                32, "ave(X)", zsto1,zout)
    245245
    246246         CALL histdef(nid_mth, "LWnet", "Net LW flux","W/m2",
    247      .                iim,jj_nb,nhori, klev,1,klev,nvert,
     247     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
    248248     .                32, "ave(X)", zsto1,zout)
    249249
    250250c        CALL histdef(nid_mth, "LWup", "upward LW flux","W/m2",
    251 c    .                iim,jj_nb,nhori, klev,1,klev,nvert,
     251c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
    252252c    .                32, "ave(X)", zsto1,zout)
    253253
    254254c        CALL histdef(nid_mth, "LWdn", "downward LW flux","W/m2",
    255 c    .                iim,jj_nb,nhori, klev,1,klev,nvert,
     255c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
    256256c    .                32, "ave(X)", zsto1,zout)
    257257
    258258         CALL histdef(nid_mth, "fluxvdf", "PBL net flux","W/m2",
    259      .                iim,jj_nb,nhori, klev,1,klev,nvert,
     259     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
    260260     .                32, "ave(X)", zsto,zout)
    261261
    262262         CALL histdef(nid_mth, "fluxdyn", "Dyn. net flux","W/m2",
    263      .                iim,jj_nb,nhori, klev,1,klev,nvert,
     263     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
    264264     .                32, "ave(X)", zsto,zout)
    265265
    266266         CALL histdef(nid_mth, "fluxajs", "Dry adj. net flux","W/m2",
    267      .                iim,jj_nb,nhori, klev,1,klev,nvert,
     267     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
    268268     .                32, "ave(X)", zsto,zout)
    269269
    270270c        CALL histdef(nid_mth, "fluxec", "Cin. net flux","W/m2",
    271 c    .                iim,jj_nb,nhori, klev,1,klev,nvert,
     271c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
    272272c    .                32, "ave(X)", zsto,zout)
    273273
     
    277277           write(str2,'(i2.2)') k
    278278         CALL histdef(nid_mth,"thv"//str2,"Haze Opa Vis",
    279      .                "--",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     279     .                "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    280280     .                "ave(X)",zsto1,zout)
    281281         ENDDO
     
    284284           write(str2,'(i2.2)') k
    285285         CALL histdef(nid_mth,"thi"//str2,"Haze Opa IR",
    286      .                "--",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     286     .                "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    287287     .                "ave(X)",zsto1,zout)
    288288         ENDDO
     
    293293           write(str2,'(i2.2)') k
    294294         CALL histdef(nid_mth,"khv"//str2,"Haze ext Vis ",
    295      .                "m-1",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     295     .                "m-1",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    296296     .                "ave(X)",zsto1,zout)
    297297         ENDDO
     
    300300           write(str2,'(i2.2)') k
    301301         CALL histdef(nid_mth,"khi"//str2,"Haze ext IR ",
    302      .                "m-1",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     302     .                "m-1",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    303303     .                "ave(X)",zsto1,zout)
    304304         ENDDO
     
    309309           write(str2,'(i2.2)') k
    310310         CALL histdef(nid_mth,"tgv"//str2,"Gas Opa Vis",
    311      .                "--",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     311     .                "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    312312     .                "ave(X)",zsto1,zout)
    313313         ENDDO
     
    316316           write(str2,'(i2.2)') k
    317317         CALL histdef(nid_mth,"tgi"//str2,"Haze Opa IR",
    318      .                "--",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     318     .                "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    319319     .                "ave(X)",zsto1,zout)
    320320         ENDDO
     
    325325           write(str2,'(i2.2)') k
    326326         CALL histdef(nid_mth,"kgv"//str2,"Gas ext Vis ",
    327      .                "m-1",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     327     .                "m-1",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    328328     .                "ave(X)",zsto1,zout)
    329329         ENDDO
     
    332332           write(str2,'(i2.2)') k
    333333         CALL histdef(nid_mth,"kgi"//str2,"Gas ext IR ",
    334      .                "m-1",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     334     .                "m-1",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    335335     .                "ave(X)",zsto1,zout)
    336336         ENDDO
     
    340340         if (clouds.eq.1) then
    341341           CALL histdef(nid_mth,"tcld","Cld Opa proxy",
    342      .                "--",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     342     .                "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    343343     .                "ave(X)",zsto,zout)
    344344
     
    346346c ----- EXTINCTION NUAGES
    347347           CALL histdef(nid_mth,"kcld","Cld Ext proxy",
    348      .                "m-1",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     348     .                "m-1",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    349349     .                "ave(X)",zsto,zout)
    350350         endif
     
    355355             write(str2,'(i2.2)') k
    356356             CALL histdef(nid_mth,"occcld"//str2,"occ cld",
    357      .       "--",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     357     .       "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    358358     .       "ave(X)",zsto,zout)
    359359           enddo
     
    365365
    366366         CALL histdef(nid_mth, "dtdyn", "Dynamics dT", "K/s",
    367      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     367     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    368368     .                "ave(X)", zsto,zout)
    369369
    370370         CALL histdef(nid_mth, "dtphy", "Physics dT", "K/s",
    371      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     371     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    372372     .                "ave(X)", zsto,zout)
    373373
    374374         CALL histdef(nid_mth, "dtvdf", "Boundary-layer dT", "K/s",
    375      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     375     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    376376     .                "ave(X)", zsto,zout)
    377377
    378378         CALL histdef(nid_mth, "dtajs", "Dry adjust. dT", "K/s",
    379      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     379     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    380380     .                "ave(X)", zsto,zout)
    381381
    382382         CALL histdef(nid_mth, "dtswr", "SW radiation dT", "K/s",
    383      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     383     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    384384     .                "ave(X)", zsto,zout)
    385385
    386386         CALL histdef(nid_mth, "dtlwr", "LW radiation dT", "K/s",
    387      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     387     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    388388     .                "ave(X)", zsto,zout)
    389389
    390390c        CALL histdef(nid_mth, "dtec", "Cinetic dissip dT", "K/s",
    391 c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     391c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    392392c    .                "ave(X)", zsto,zout)
    393393
    394394c        CALL histdef(nid_mth, "dvvdf", "Boundary-layer dV", "m/s2",
    395 c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     395c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    396396c    .                "ave(X)", zsto,zout)
    397397
     
    404404c        call histdef(nid_mth, "taux",
    405405c    $         "Zonal wind stress", "Pa", 
    406 c    $         iim,jj_nb,nhori, 1,1,1, nvert, 32,
     406c    $         nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    407407c    $         "ave(X)", zsto,zout)
    408408
    409409c        call histdef(nid_mth, "tauy",
    410410c    $         "Meridional xind stress", "Pa", 
    411 c    $         iim,jj_nb,nhori, 1,1,1, nvert, 32,
     411c    $         nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    412412c    $         "ave(X)", zsto,zout)
    413413
    414414c        CALL histdef(nid_mth, "cdrm", "Momentum drag coef.", "-",
    415 c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     415c    .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    416416c    .                "ave(X)", zsto,zout)
    417417
    418418c        CALL histdef(nid_mth, "cdrh", "Heat drag coef.", "-",
    419 c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     419c    .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    420420c    .                "ave(X)", zsto,zout)
    421421
  • trunk/LMDZ.TITAN/libf/phytitan/inicondens.F

    r102 r1530  
    1010
    1111      use dimphy
     12      use mod_grid_phy_lmdz, only: nbp_lev
    1213      IMPLICIT NONE
    13 #include "dimensions.h"
    1414#include "YOMCST.h"
    1515
     
    1717c    -----------
    1818      INTEGER ny
    19       real press(llm),temp(llm)  ! pressure in mbar !
    20       REAL yc(llm,ny)
     19      real press(nbp_lev),temp(nbp_lev)  ! pressure in mbar !
     20      REAL yc(nbp_lev,ny)
    2121      character*10 nomy(ny+1)
    2222     
     
    2828      do ic=1,ny
    2929       print*, 'traceur CH(', ic, ')=', nomy(ic),'------------'
    30            do l=1,llm
     30           do l=1,nbp_lev
    3131
    3232c  Par defaut, yc est a 1 c'est a dire qu'on ne condense pas
  • trunk/LMDZ.TITAN/libf/phytitan/iophy.F90

    r1056 r1530  
    2424
    2525  subroutine init_iophy_new(rlat,rlon)
    26   USE dimphy
     26  USE dimphy, only: klon
    2727  USE mod_phys_lmdz_para
    28   USE mod_grid_phy_lmdz
     28  USE mod_grid_phy_lmdz, only: nbp_lon, nbp_lat, klon_glo
    2929  USE ioipsl
    3030  implicit none
    31   include 'dimensions.h'   
    3231    real,dimension(klon),intent(in) :: rlon
    3332    real,dimension(klon),intent(in) :: rlat
     
    5150   
    5251!$OMP MASTER 
    53     ALLOCATE(io_lat(jjm+1-1/(iim*jjm)))
    54     io_lat(1)=rlat_glo(1)
    55     io_lat(jjm+1-1/(iim*jjm))=rlat_glo(klon_glo)
    56     IF ((iim*jjm) > 1) then
    57       DO i=2,jjm
    58         io_lat(i)=rlat_glo(2+(i-2)*iim)
     52    ALLOCATE(io_lat(nbp_lat))
     53    IF (klon_glo == 1) THEN
     54      io_lat(1)=rlat_glo(1)
     55    ELSE
     56      io_lat(1)=rlat_glo(1)
     57      io_lat(nbp_lat)=rlat_glo(klon_glo)
     58      DO i=2,nbp_lat-1
     59        io_lat(i)=rlat_glo(2+(i-2)*nbp_lon)
    5960      ENDDO
    6061    ENDIF
    6162
    62     ALLOCATE(io_lon(iim))
    63     io_lon(:)=rlon_glo(2-1/(iim*jjm):iim+1-1/(iim*jjm))
     63    ALLOCATE(io_lon(nbp_lon))
     64    IF (klon_glo == 1) THEN
     65      io_lon(1)=rlon_glo(1)
     66    ELSE
     67      io_lon(1:nbp_lon)=rlon_glo(2:nbp_lon+1)
     68    ENDIF
    6469
    6570    ddid=(/ 1,2 /)
    66     dsg=(/ iim, jjm+1-1/(iim*jjm) /)
    67     dsl=(/ iim, jj_nb /)
     71    dsg=(/ nbp_lon, nbp_lat /)
     72    dsl=(/ nbp_lon, jj_nb /)
    6873    dpf=(/ 1,jj_begin /)
    69     dpl=(/ iim, jj_end /)
     74    dpl=(/ nbp_lon, jj_end /)
    7075    dhs=(/ ii_begin-1,0 /)
    71     if (mpi_rank==mpi_size-1) then
     76    IF (mpi_rank==mpi_size-1) THEN
    7277      dhe=(/0,0/)
    73     else
    74       dhe=(/ iim-ii_end,0 /) 
    75     endif
     78    ELSE
     79      dhe=(/ nbp_lon-ii_end,0 /) 
     80    ENDIF
    7681   
    7782    call flio_dom_set(mpi_size,mpi_rank,ddid,dsg,dsl,dpf,dpl,dhs,dhe, &
     
    8590  USE dimphy
    8691  USE mod_phys_lmdz_para
    87   use ioipsl
    88   implicit none
    89   include 'dimensions.h'   
    90     real,dimension(iim),intent(in) :: lon
    91     real,dimension(jjm+1-1/(iim*jjm)),intent(in) :: lat
     92  use ioipsl, only: flio_dom_set
     93  USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat
     94  implicit none
     95    real,dimension(nbp_lon),intent(in) :: lon
     96    real,dimension(nbp_lat),intent(in) :: lat
    9297
    9398    INTEGER,DIMENSION(2) :: ddid
     
    100105
    101106!$OMP MASTER 
    102     allocate(io_lat(jjm+1-1/(iim*jjm)))
     107    allocate(io_lat(nbp_lat))
    103108    io_lat(:)=lat(:)
    104     allocate(io_lon(iim))
     109    allocate(io_lon(nbp_lon))
    105110    io_lon(:)=lon(:)
    106111   
    107112    ddid=(/ 1,2 /)
    108     dsg=(/ iim, jjm+1-1/(iim*jjm) /)
    109     dsl=(/ iim, jj_nb /)
     113    dsg=(/ nbp_lon, nbp_lat /)
     114    dsl=(/ nbp_lon, jj_nb /)
    110115    dpf=(/ 1,jj_begin /)
    111     dpl=(/ iim, jj_end /)
     116    dpl=(/ nbp_lon, jj_end /)
    112117    dhs=(/ ii_begin-1,0 /)
    113118    if (mpi_rank==mpi_size-1) then
    114119      dhe=(/0,0/)
    115120    else
    116       dhe=(/ iim-ii_end,0 /) 
     121      dhe=(/ nbp_lon-ii_end,0 /) 
    117122    endif
    118123   
     
    125130 
    126131  subroutine histbeg_phy(name,itau0,zjulian,dtime,nhori,nid_day)
    127   USE dimphy
    128   USE mod_phys_lmdz_para
    129   use ioipsl
     132  USE mod_phys_lmdz_para, only: jj_begin, jj_end, jj_nb, is_sequential
     133  USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat
     134  use ioipsl, only: histbeg
    130135  use write_field
    131136  implicit none
    132   include 'dimensions.h'
    133137   
    134138    character*(*), intent(IN) :: name
     
    141145!$OMP MASTER   
    142146    if (is_sequential) then
    143       call histbeg(name,iim,io_lon, jj_nb,io_lat(jj_begin:jj_end), &
    144                    1,iim,1,jj_nb,itau0, zjulian, dtime, nhori, nid_day)
     147      call histbeg(name,nbp_lon,io_lon, jj_nb,io_lat(jj_begin:jj_end), &
     148                   1,nbp_lon,1,jj_nb,itau0, zjulian, dtime, nhori, nid_day)
    145149    else
    146       call histbeg(name,iim,io_lon, jj_nb,io_lat(jj_begin:jj_end), &
    147                    1,iim,1,jj_nb,itau0, zjulian, dtime, nhori, nid_day,phys_domain_id)
     150      call histbeg(name,nbp_lon,io_lon, jj_nb,io_lat(jj_begin:jj_end), &
     151                   1,nbp_lon,1,jj_nb,itau0, zjulian, dtime, nhori, nid_day,phys_domain_id)
    148152    endif
    149153!$OMP END MASTER
     
    154158             plon,plat,plon_bounds,plat_bounds, &
    155159             nname,itau0,zjulian,dtime,nnhori,nnid_day)
    156   USE dimphy
     160  USE dimphy, only: klon
    157161  USE mod_phys_lmdz_para
    158   USE mod_grid_phy_lmdz
    159   use ioipsl
    160   use write_field
    161   implicit none
    162   include 'dimensions.h'
     162  USE mod_grid_phy_lmdz, only: klon_glo, nbp_lon, nbp_lat
     163  use ioipsl, only: histbeg
     164
     165  implicit none
    163166
    164167    real,dimension(klon),intent(in) :: rlon
     
    186189    REAL, allocatable, dimension(:) :: npplat, npplon
    187190    REAL, allocatable, dimension(:,:) :: npplat_bounds, npplon_bounds
    188     INTEGER, PARAMETER :: jjmp1=jjm+1-1/jjm
    189     REAL, dimension(iim,jjmp1) :: zx_lon, zx_lat
     191    REAL, dimension(nbp_lon,nbp_lat) :: zx_lon, zx_lat
    190192
    191193    CALL gather(rlat,rlat_glo)
     
    212214     endif
    213215!
    214      IF ( tabij(i).LE.iim) THEN
     216     IF ( tabij(i).LE.nbp_lon) THEN
    215217      plat_bounds(i,1)=rlat_glo(tabij(i))
    216218     ELSE
    217       plat_bounds(i,1)=rlat_glo(tabij(i)-iim)
     219      plat_bounds(i,1)=rlat_glo(tabij(i)-nbp_lon)
    218220     ENDIF
    219      plat_bounds(i,2)=rlat_glo(tabij(i)+iim)
     221     plat_bounds(i,2)=rlat_glo(tabij(i)+nbp_lon)
    220222!
    221223!    print*,'CFMIP_iophy point i lon lon_bds',i,plon_bounds(i,1),rlon_glo(tabij(i)),plon_bounds(i,2)
     
    233235     ENDDO
    234236
    235        CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon_glo,zx_lon)
    236        if ((iim*jjm).gt.1) then
    237        DO i = 1, iim
     237       CALL gr_fi_ecrit(1,klon,nbp_lon,nbp_lat,rlon_glo,zx_lon)
     238       if ((nbp_lon*nbp_lat).gt.1) then
     239       DO i = 1, nbp_lon
    238240         zx_lon(i,1) = rlon_glo(i+1)
    239          zx_lon(i,jjmp1) = rlon_glo(i+1)
     241         zx_lon(i,nbp_lat) = rlon_glo(i+1)
    240242       ENDDO
    241243       endif
    242        CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat_glo,zx_lat)
     244       CALL gr_fi_ecrit(1,klon,nbp_lon,nbp_lat,rlat_glo,zx_lat)
    243245
    244246    DO i=1,pim
     
    249251
    250252     if (ipt(i).EQ.1) then
    251       plon_bounds(i,1)=zx_lon(iim,jpt(i))
     253      plon_bounds(i,1)=zx_lon(nbp_lon,jpt(i))
    252254      plon_bounds(i,2)=360.+zx_lon(ipt(i)+1,jpt(i))
    253255     endif
    254256 
    255      if (ipt(i).EQ.iim) then
     257     if (ipt(i).EQ.nbp_lon) then
    256258      plon_bounds(i,2)=360.+zx_lon(1,jpt(i))
    257259     endif
     
    265267     endif
    266268 
    267      if (jpt(i).EQ.jjmp1) then
    268       plat_bounds(i,1)=zx_lat(ipt(i),jjmp1)+0.001
    269       plat_bounds(i,2)=zx_lat(ipt(i),jjmp1)-0.001
     269     if (jpt(i).EQ.nbp_lat) then
     270      plat_bounds(i,1)=zx_lat(ipt(i),nbp_lat)+0.001
     271      plat_bounds(i,2)=zx_lat(ipt(i),nbp_lat)-0.001
    270272     endif
    271273!
     
    323325 
    324326  subroutine histwrite2d_phy(nid,lpoint,name,itau,field)
    325   USE dimphy
     327  USE dimphy, only: klon
    326328  USE mod_phys_lmdz_para
    327   USE ioipsl
    328   implicit none
    329   include 'dimensions.h'
     329  USE ioipsl, only: histwrite
     330  USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat
     331  implicit none
    330332   
    331333    integer,intent(in) :: nid
     
    336338    REAL,dimension(klon_mpi) :: buffer_omp
    337339    INTEGER, allocatable, dimension(:) :: index2d
    338     REAL :: Field2d(iim,jj_nb)
     340    REAL :: Field2d(nbp_lon,jj_nb)
    339341
    340342    integer :: ip
     
    347349    CALL grid1Dto2D_mpi(buffer_omp,Field2d)
    348350    if(.NOT.lpoint) THEN
    349      ALLOCATE(index2d(iim*jj_nb))
    350      ALLOCATE(fieldok(iim*jj_nb))
    351      CALL histwrite(nid,name,itau,Field2d,iim*jj_nb,index2d)
     351     ALLOCATE(index2d(nbp_lon*jj_nb))
     352     ALLOCATE(fieldok(nbp_lon*jj_nb))
     353     CALL histwrite(nid,name,itau,Field2d,nbp_lon*jj_nb,index2d)
    352354    else
    353355     ALLOCATE(fieldok(npstn))
     
    378380
    379381  subroutine histwrite3d_phy(nid,lpoint,name,itau,field)
    380   USE dimphy
     382  USE dimphy, only: klon
    381383  USE mod_phys_lmdz_para
    382 
    383   use ioipsl
    384   implicit none
    385   include 'dimensions.h'
     384  USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat
     385  use ioipsl, only: histwrite
     386  implicit none
    386387   
    387388    integer,intent(in) :: nid
     
    391392    real,dimension(:,:),intent(in) :: field  ! --> field(klon,:)
    392393    REAL,dimension(klon_mpi,size(field,2)) :: buffer_omp
    393     REAL :: Field3d(iim,jj_nb,size(field,2))
     394    REAL :: Field3d(nbp_lon,jj_nb,size(field,2))
    394395    INTEGER :: ip, n, nlev
    395396    INTEGER, ALLOCATABLE, dimension(:) :: index3d
     
    409410    CALL grid1Dto2D_mpi(buffer_omp,field3d)
    410411    if(.NOT.lpoint) THEN
    411      ALLOCATE(index3d(iim*jj_nb*nlev))
    412      ALLOCATE(fieldok(iim*jj_nb,nlev))
    413      CALL histwrite(nid,name,itau,Field3d,iim*jj_nb*nlev,index3d)
     412     ALLOCATE(index3d(nbp_lon*jj_nb*nlev))
     413     ALLOCATE(fieldok(nbp_lon*jj_nb,nlev))
     414     CALL histwrite(nid,name,itau,Field3d,nbp_lon*jj_nb*nlev,index3d)
    414415    else
    415416      nlev=size(field,2)
  • trunk/LMDZ.TITAN/libf/phytitan/orodrag.F

    r1056 r1530  
    1010      use dimphy
    1111      IMPLICIT NONE
    12 
    13 #include "dimensions.h"
    14 #include "paramet.h"
    1512
    1613c
  • trunk/LMDZ.TITAN/libf/phytitan/orosetup.F

    r1056 r1530  
    9797      use dimphy
    9898      implicit none
    99 
    100 #include "dimensions.h"
    101 #include "paramet.h"
    10299
    103100#include "YOMCST.h"
  • trunk/LMDZ.TITAN/libf/phytitan/phyetat0.F90

    r1525 r1530  
    2121! Objet: Lecture de l'etat initial pour la physique
    2222!======================================================================
    23 #include "dimensions.h"
    2423#include "netcdf.inc"
    2524#include "dimsoil.h"
  • trunk/LMDZ.TITAN/libf/phytitan/phyredem.F90

    r1525 r1530  
    2020! Objet: Lecture de l'etat initial pour la physique
    2121!======================================================================
    22 #include "dimensions.h"
    2322#include "netcdf.inc"
    2423#include "dimsoil.h"
  • trunk/LMDZ.TITAN/libf/phytitan/physiq.F

    r1525 r1530  
    7373      USE time_phylmdz_mod, only: itau_phy,day_ref,annee_ref,nday
    7474      USE logic_mod, only: iflag_trac,moyzon_ch,moyzon_mu
     75      USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat, nbp_lev
    7576      IMPLICIT none
    7677c======================================================================
  • trunk/LMDZ.TITAN/libf/phytitan/radlwsw.F

    r1356 r1530  
    4141      USE write_field_phy
    4242       IMPLICIT none
    43 #include "dimensions.h"
    4443#include "YOMCST.h"
    4544#include "clesphys.h"
  • trunk/LMDZ.TITAN/libf/phytitan/soil.F

    r102 r1530  
    4848      use dimphy
    4949      IMPLICIT NONE
    50 #include "dimensions.h"
    5150#include "YOMCST.h"
    5251#include "dimsoil.h"
  • trunk/LMDZ.TITAN/libf/phytitan/sources.F

    r808 r1530  
    2323         use dimphy
    2424         IMPLICIT NONE
    25 #include "dimensions.h"
    26 #include "microtab.h"
    2725#include "YOMCST.h"
    2826c
  • trunk/LMDZ.TITAN/libf/phytitan/sugwd.F

    r1056 r1530  
    7676      use dimphy
    7777      IMPLICIT NONE
    78 
    79 #include "dimensions.h"
    80 #include "paramet.h"
    8178
    8279#include "YOEGWD.h"
  • trunk/LMDZ.TITAN/libf/phytitan/tgmdat_mod.F90

    r1461 r1530  
    11MODULE  TGMDAT_MOD
    22      IMPLICIT NONE
    3       include "dimensions.h"
    43
    5       INTEGER,PARAMETER :: NLAYER=llm,NLEVEL=NLAYER+1
    64      REAL,SAVE :: PI=3.14159265358979323846
    75! RGAS IS THE UNIVERSAL GAS CONSTANT IN UNITS OF: M SEC-2 AMU K-1 KM
  • trunk/LMDZ.TITAN/libf/phytitan/ustarhb.F

    r495 r1530  
    1818      use dimphy
    1919      IMPLICIT none
    20 #include "dimensions.h"
    2120#include "YOMCST.h"
    2221c
  • trunk/LMDZ.TITAN/libf/phytitan/vdif_kcay.F

    r102 r1530  
    88      use dimphy
    99      IMPLICIT NONE
    10 #include "dimensions.h"
    1110c.......................................................................
    1211c
  • trunk/LMDZ.TITAN/libf/phytitan/yamada.F

    r102 r1530  
    88      use dimphy
    99      IMPLICIT NONE
    10 #include "dimensions.h"
    1110c.......................................................................
    1211c
  • trunk/LMDZ.TITAN/libf/phytitan/yamada4.F

    r814 r1530  
    99      use dimphy
    1010      IMPLICIT NONE
    11 #include "dimensions.h"
    1211#include "YOMCST.h"
    1312c.......................................................................
  • trunk/LMDZ.TITAN/libf/phytitan/zenang.F

    r119 r1530  
    2626      use dimphy
    2727      IMPLICIT none
    28 #include "dimensions.h"
    2928#include "YOMCST.h"
    3029#include "comorbit.h"
  • trunk/LMDZ.VENUS/libf/phyvenus/aaam_bud.F

    r101 r1530  
    77c
    88      use dimphy
     9      use mod_grid_phy_lmdz, only: nbp_lon, nbp_lat, klon_glo
    910      implicit none
    1011c======================================================================
     
    5253c ===================
    5354c
    54 c iim--common-I: Number of longitude intervals
    55 c jjm--common-I: Number of latitude intervals
     55c nbp_lon--common-I: Number of longitude intervals
     56c nbp_lat-1--common-I: Number of latitude intervals
    5657c klon-common-I: Number of points seen by the physics
    57 c                iim*(jjm-1)+2 for instance
     58c                nbp_lon*(nbp_lat-1-1)+2 for instance
    5859c klev-common-I: Number of vertical layers
    5960c======================================================================
     
    7879c======================================================================
    7980
    80 #include "dimensions.h"
    8181c
    8282c ARGUMENTS
     
    106106C  PUT AAM QUANTITIES AT ZERO:
    107107C
    108       if(iim+1.gt.801.or.jjm+1.gt.401)then
     108      if(nbp_lon+1.gt.801.or.nbp_lat.gt.401)then
    109109      print *,' Pb de dimension dans aaam_bud'
    110110      stop
     
    114114      hadley=1.e18
    115115      hadday=1.e18*1.e7
    116       dlat=xpi/float(jjm)
    117       dlon=2.*xpi/float(iim)
     116      IF (klon_glo.EQ.1) THEN
     117        dlat=xpi
     118      ELSE
     119        dlat=xpi/float(nbp_lat-1)
     120      ENDIF
     121      dlon=2.*xpi/float(nbp_lon)
    118122     
    119123      do iax=1,3
     
    140144          zlat(1)=plat(l)*xpi/180.
    141145
    142         do i=1,iim+1
     146        do i=1,nbp_lon+1
    143147
    144148          zs(i,1)=phis(l)/rg
     
    154158
    155159
    156       do j = 2,jjm
     160      do j = 2,nbp_lat-1
    157161
    158162C Values at Greenwich (Periodicity)
    159163
    160       zs(iim+1,j)=phis(l+1)/rg
    161       ps(iim+1,j)=p(l+1,1)
    162           ssou(iim+1,j)=dragu(l+1)+liftu(l+1)
    163           ssov(iim+1,j)=dragv(l+1)+liftv(l+1)
    164           blsu(iim+1,j)=clu(l+1)
    165           blsv(iim+1,j)=clv(l+1)
    166       zlon(iim+1)=-plon(l+1)*xpi/180.
     164      zs(nbp_lon+1,j)=phis(l+1)/rg
     165      ps(nbp_lon+1,j)=p(l+1,1)
     166          ssou(nbp_lon+1,j)=dragu(l+1)+liftu(l+1)
     167          ssov(nbp_lon+1,j)=dragv(l+1)+liftv(l+1)
     168          blsu(nbp_lon+1,j)=clu(l+1)
     169          blsv(nbp_lon+1,j)=clv(l+1)
     170      zlon(nbp_lon+1)=-plon(l+1)*xpi/180.
    167171      zlat(j)=plat(l+1)*xpi/180.
    168172
    169       ub(iim+1,j)=0.
    170       vb(iim+1,j)=0.
     173      ub(nbp_lon+1,j)=0.
     174      vb(nbp_lon+1,j)=0.
    171175         do k=1,nlev
    172          ub(iim+1,j)=ub(iim+1,j)+u(l+1,k)*(p(l+1,k)-p(l+1,k+1))/rg
    173          vb(iim+1,j)=vb(iim+1,j)+v(l+1,k)*(p(l+1,k)-p(l+1,k+1))/rg
     176         ub(nbp_lon+1,j)=ub(nbp_lon+1,j)+u(l+1,k)*
     177     &                                   (p(l+1,k)-p(l+1,k+1))/rg
     178         vb(nbp_lon+1,j)=vb(nbp_lon+1,j)+v(l+1,k)*
     179     &                                   (p(l+1,k)-p(l+1,k+1))/rg
    174180         enddo
    175181     
    176182
    177       do i=1,iim
     183      do i=1,nbp_lon
    178184
    179185      l=l+1
     
    201207
    202208      l=l+1
    203       ub(1,jjm+1)=0.
    204       vb(1,jjm+1)=0.
     209      ub(1,nbp_lat)=0.
     210      vb(1,nbp_lat)=0.
    205211      do k=1,nlev
    206          ub(1,jjm+1)=ub(1,jjm+1)+u(l,k)*(p(l,k)-p(l,k+1))/rg
    207          vb(1,jjm+1)=vb(1,jjm+1)+v(l,k)*(p(l,k)-p(l,k+1))/rg
    208       enddo
    209       zlat(jjm+1)=plat(l)*xpi/180.
    210 
    211       do i=1,iim+1
    212       zs(i,jjm+1)=phis(l)/rg
    213       ps(i,jjm+1)=p(l,1)
    214           ssou(i,jjm+1)=dragu(l)+liftu(l)
    215           ssov(i,jjm+1)=dragv(l)+liftv(l)
    216           blsu(i,jjm+1)=clu(l)
    217           blsv(i,jjm+1)=clv(l)
    218       ub(i,jjm+1)=ub(1,jjm+1)                               
    219       vb(i,jjm+1)=vb(1,jjm+1)                               
     212         ub(1,nbp_lat)=ub(1,nbp_lat)+u(l,k)*(p(l,k)-p(l,k+1))/rg
     213         vb(1,nbp_lat)=vb(1,nbp_lat)+v(l,k)*(p(l,k)-p(l,k+1))/rg
     214      enddo
     215      zlat(nbp_lat)=plat(l)*xpi/180.
     216
     217      do i=1,nbp_lon+1
     218      zs(i,nbp_lat)=phis(l)/rg
     219      ps(i,nbp_lat)=p(l,1)
     220          ssou(i,nbp_lat)=dragu(l)+liftu(l)
     221          ssov(i,nbp_lat)=dragv(l)+liftv(l)
     222          blsu(i,nbp_lat)=clu(l)
     223          blsv(i,nbp_lat)=clv(l)
     224      ub(i,nbp_lat)=ub(1,nbp_lat)                               
     225      vb(i,nbp_lat)=vb(1,nbp_lat)                               
    220226      enddo
    221227
     
    223229C  MOMENT ANGULAIRE
    224230C
    225         DO j=1,jjm   
    226         DO i=1,iim
     231        DO j=1,nbp_lat-1   
     232        DO i=1,nbp_lon
    227233
    228234           raam(1)=raam(1)-rea**3*dlon*dlat*0.5*
     
    261267C
    262268
    263         DO j=1,jjm
    264         DO i=1,iim
     269        DO j=1,nbp_lat-1
     270        DO i=1,nbp_lon
    265271           tmou(1)=tmou(1)-rea**2*dlon*0.5*sin(zlon(i))
    266272     c  *(zs(i,j)-zs(i,j+1))
     
    272278        ENDDO
    273279           
    274         DO j=2,jjm
    275         DO i=1,iim
     280        DO j=2,nbp_lat-1
     281        DO i=1,nbp_lon
    276282           tmou(1)=tmou(1)+rea**2*dlat*0.5*sin(zlat(j))
    277283     c  *(zs(i+1,j)-zs(i,j))
     
    288294C
    289295        l=1
    290         DO j=2,jjm
    291         DO i=1,iim
     296        DO j=2,nbp_lat-1
     297        DO i=1,nbp_lon
    292298        l=l+1
    293299           tsso(1)=tsso(1)-rea**3*cos(zlat(j))*dlon*dlat*
     
    335341100   format(F12.5,5(1x,F12.5))
    336342
    337       write(iam+1,*)((zs(i,j),i=1,iim),j=1,jjm+1)
    338       write(iam+1,*)((ps(i,j),i=1,iim),j=1,jjm+1)
    339       write(iam+1,*)((ub(i,j),i=1,iim),j=1,jjm+1)
    340       write(iam+1,*)((vb(i,j),i=1,iim),j=1,jjm+1)
    341       write(iam+1,*)((ssou(i,j),i=1,iim),j=1,jjm+1)
    342       write(iam+1,*)((ssov(i,j),i=1,iim),j=1,jjm+1)
    343       write(iam+1,*)((blsu(i,j),i=1,iim),j=1,jjm+1)
    344       write(iam+1,*)((blsv(i,j),i=1,iim),j=1,jjm+1)
     343      write(iam+1,*)((zs(i,j),i=1,nbp_lon),j=1,nbp_lat)
     344      write(iam+1,*)((ps(i,j),i=1,nbp_lon),j=1,nbp_lat)
     345      write(iam+1,*)((ub(i,j),i=1,nbp_lon),j=1,nbp_lat)
     346      write(iam+1,*)((vb(i,j),i=1,nbp_lon),j=1,nbp_lat)
     347      write(iam+1,*)((ssou(i,j),i=1,nbp_lon),j=1,nbp_lat)
     348      write(iam+1,*)((ssov(i,j),i=1,nbp_lon),j=1,nbp_lat)
     349      write(iam+1,*)((blsu(i,j),i=1,nbp_lon),j=1,nbp_lat)
     350      write(iam+1,*)((blsv(i,j),i=1,nbp_lon),j=1,nbp_lat)
    345351
    346352      RETURN
  • trunk/LMDZ.VENUS/libf/phyvenus/ajsec.F

    r1301 r1530  
    77
    88      use dimphy
     9      use mod_grid_phy_lmdz, only: nbp_lev
    910      use cpdet_mod, only: t2tpot, tpot2t
    1011      IMPLICIT none
     
    2728c d_qfi-----output-R-Incrementation des traceurs
    2829c======================================================================
    29 #include "dimensions.h"
    3030#include "YOMCST.h"
    3131      REAL paprs(klon,klev+1), pplay(klon,klev)
     
    6666c------------------------------------- passage en temperature potentielle
    6767! ADAPTATION GCM POUR CP(T)
    68       call t2tpot(klon*llm,tfi,zh,ppk)
     68      call t2tpot(klon*nbp_lev,tfi,zh,ppk)
    6969c
    7070      DO k = limbas, limhau
     
    168168c------------------------------------- et calcul du d_t
    169169! ADAPTATION GCM POUR CP(T)
    170       call tpot2t(klon*llm,zh,zt,ppk)
     170      call tpot2t(klon*nbp_lev,zh,zt,ppk)
    171171
    172172      DO k = limbas, limhau
  • trunk/LMDZ.VENUS/libf/phyvenus/ballon.F

    r1442 r1530  
    33
    44      use dimphy
     5      use mod_grid_phy_lmdz, only: nbp_lon, nbp_lat
    56      implicit none
    67
     
    9495c======================================================================
    9596
    96 #include "dimensions.h"
    9797#include "YOMCST.h"
    9898c
     
    110110      INTEGER jj,ii,ll
    111111
    112       REAL zlon(iim+1),zlat(jjm+1)
    113       save zlon,zlat
     112      REAL,SAVE,ALLOCATABLE :: zlon(:),zlat(:)
    114113
    115114      REAL time
     
    141140      print*,"BALLOONS ACTIVATED"
    142141
     142      allocate(zlon(nbp_lon+1))
     143      allocate(zlat(nbp_lat))
     144     
    143145C Latitudes:
    144146      zlat(1)=plat(1)*RPI/180.
    145       do j = 2,jjm
    146          k=(j-2)*iim+2
     147      do j = 2,nbp_lat-1
     148         k=(j-2)*nbp_lon+2
    147149         zlat(j)=plat(k)*RPI/180.
    148150      enddo
    149       zlat(jjm+1)=plat(klon)*RPI/180.
     151      zlat(nbp_lat)=plat(klon)*RPI/180.
    150152
    151153C Longitudes:
    152       do i = 1,iim
     154      do i = 1,nbp_lon
    153155         k=i+1
    154156         zlon(i)=plon(k)*RPI/180.
    155157      enddo
    156       zlon(iim+1)=zlon(1)+2.*RPI
     158      zlon(nbp_lon+1)=zlon(1)+2.*RPI
    157159
    158160c verif init     lat de 90 à -90, lon de -180 à 180
     
    258260
    259261      use dimphy
     262      use mod_grid_phy_lmdz, only: nbp_lon, nbp_lat
    260263      implicit none
    261264
     
    285288c======================================================================
    286289
    287 #include "dimensions.h"
    288290#include "YOMCST.h"
    289291c
     
    291293c
    292294      real map_u(klon),map_v(klon),map_a(klon)
    293       real latit(jjm+1),longit(iim)
     295      real latit(nbp_lat),longit(nbp_lon)
    294296      real phi,lam,ubal,vbal,abal
    295297c
     
    298300      INTEGER i,j,k
    299301      INTEGER jj,ii
    300       REAL    ujj(iim+1),vjj(iim+1),ajj(iim+1)
     302      REAL    ujj(nbp_lon+1),vjj(nbp_lon+1),ajj(nbp_lon+1)
    301303      REAL    factlat,factlon
    302304
     
    304306c-------------------------------------------------
    305307        jj=1  ! POLE NORD
    306         do j=2,jjm
     308        do j=2,nbp_lat-1
    307309          if (phi.lt.latit(j)) jj=j
    308310        enddo
     
    311313c pole nord
    312314        if (jj.eq.1) then
    313          do i=1,iim
     315         do i=1,nbp_lon
    314316          ujj(i) = map_u(i+1)*factlat + map_u(1)*(1-factlat)
    315317          vjj(i) = map_v(i+1)*factlat + map_v(1)*(1-factlat)
     
    317319         enddo
    318320c pole sud
    319         elseif (jj.eq.jjm) then
    320          do i=1,iim
    321           k = (jj-2)*iim+1+i
     321        elseif (jj.eq.nbp_lat-1) then
     322         do i=1,nbp_lon
     323          k = (jj-2)*nbp_lon+1+i
    322324          ujj(i) = map_u(klon)*factlat + map_u(k)*(1-factlat)
    323325          vjj(i) = map_v(klon)*factlat + map_v(k)*(1-factlat)
     
    326328c autres latitudes
    327329        else
    328          do i=1,iim
    329           k = (jj-2)*iim+1+i
    330           ujj(i) = map_u(k+iim)*factlat + map_u(k)*(1-factlat)
    331           vjj(i) = map_v(k+iim)*factlat + map_v(k)*(1-factlat)
    332           ajj(i) = map_a(k+iim)*factlat + map_a(k)*(1-factlat)
     330         do i=1,nbp_lon
     331          k = (jj-2)*nbp_lon+1+i
     332          ujj(i) = map_u(k+nbp_lon)*factlat + map_u(k)*(1-factlat)
     333          vjj(i) = map_v(k+nbp_lon)*factlat + map_v(k)*(1-factlat)
     334          ajj(i) = map_a(k+nbp_lon)*factlat + map_a(k)*(1-factlat)
    333335         enddo
    334336        endif
    335         ujj(iim+1)=ujj(1)
    336         vjj(iim+1)=vjj(1)
    337         ajj(iim+1)=ajj(1)
     337        ujj(nbp_lon+1)=ujj(1)
     338        vjj(nbp_lon+1)=vjj(1)
     339        ajj(nbp_lon+1)=ajj(1)
    338340
    339341c Interpolation in longitudes
    340342c-------------------------------------------------
    341343        ii=1  ! lon=-180
    342         do i=2,iim
     344        do i=2,nbp_lon
    343345          if (lam.gt.longit(i)) ii=i
    344346        enddo
     
    372374c======================================================================
    373375
    374 #include "dimensions.h"
    375376#include "YOMCST.h"
    376377c
  • trunk/LMDZ.VENUS/libf/phyvenus/blendrad.F

    r1310 r1530  
    1515      use dimphy
    1616      implicit none
    17 #include "dimensions.h"
    1817c#include "dimradmars.h"
    1918#include "nlteparams.h"
  • trunk/LMDZ.VENUS/libf/phyvenus/clmain.F

    r1443 r1530  
    3535      USE interface_surf
    3636      use dimphy
     37      use mod_grid_phy_lmdz, only: nbp_lev
    3738      use cpdet_mod, only: t2tpot
    3839      IMPLICIT none
     
    6667cAA la premiere couche
    6768c======================================================================
    68 #include "dimensions.h"
    6969c$$$ PB ajout pour soil
    7070#include "dimsoil.h"
     
    324324
    325325! ADAPTATION GCM POUR CP(T)
    326          call t2tpot(knon*llm,yt,yteta,ppk)
     326         call t2tpot(knon*nbp_lev,yt,yteta,ppk)
    327327
    328328         yzlev(1:knon,1)=0.
     
    479479      USE interface_surf
    480480      use dimphy
     481      use mod_grid_phy_lmdz, only: nbp_lon, nbp_lat, nbp_lev
    481482      use cpdet_mod, only: t2tpot,tpot2t,cpdet
    482483
     
    486487c Objet: diffusion verticale de "h"
    487488c======================================================================
    488 #include "dimensions.h"
    489489#include "YOMCST.h"
    490490#include "dimsoil.h"
     
    594594
    595595c passage en enthalpie potentielle
    596       call t2tpot(knon*llm,t,local_h,ppk)
     596      call t2tpot(knon*nbp_lev,t,local_h,ppk)
    597597c     print*,"tpot en entree de clqh=",local_h(klon/2,:)
    598598
     
    685685! ADAPTATION GCM POUR CP(T)
    686686      CALL interfsurf_hq(itime, dtime, rmu0,
    687      e klon, iim, jjm, knon,
     687     e klon, nbp_lon, nbp_lat-1, knon,
    688688     e rlon, rlat, cufi, cvfi,
    689689     e debut, lafin, soil_model, nsoilmx,tsoil,
     
    733733      ENDDO
    734734      ENDDO
    735       call tpot2t(knon*llm,local_h,d_t,ppk)
     735      call tpot2t(knon*nbp_lev,local_h,d_t,ppk)
    736736
    737737c     print*,"tpot en sortie de clqh=",local_h(klon/2,:)
     
    781781c flux_v---output-R- (diagnostic) flux du vent: (kg m/s)/(m**2 s)
    782782c======================================================================
    783 #include "dimensions.h"
    784783#include "iniprint.h"
    785784      INTEGER knon
     
    908907c pcfh-----output-R- coefficients a calculer (chaleur et humidite)
    909908c======================================================================
    910 #include "dimensions.h"
    911909#include "YOMCST.h"
    912910#include "iniprint.h"
     
    11721170
    11731171      use dimphy
     1172      use mod_grid_phy_lmdz, only: nbp_lev
    11741173      use cpdet_mod, only: cpdet
    11751174      IMPLICIT none
     
    11881187c pcfh-----output-R- coefficients a calculer (chaleur et humidite)
    11891188c======================================================================
    1190 #include "dimensions.h"
    11911189#include "YOMCST.h"
    11921190#include "iniprint.h"
  • trunk/LMDZ.VENUS/libf/phyvenus/cltrac.F

    r1442 r1530  
    2727c flux_tr--output-R- flux de tr
    2828c======================================================================
    29 #include "dimensions.h"
    3029      REAL dtime
    3130      REAL coef(klon,klev)
  • trunk/LMDZ.VENUS/libf/phyvenus/coefkzmin.F

    r101 r1530  
    99      IMPLICIT NONE
    1010
    11 #include "dimensions.h"
    1211#include "YOMCST.h"
    1312
  • trunk/LMDZ.VENUS/libf/phyvenus/concentrations2.F

    r1452 r1530  
    2222!     declarations
    2323 
    24 #include "dimensions.h"
    2524#include "YOMCST.h"
    2625#include "clesphys.h"
  • trunk/LMDZ.VENUS/libf/phyvenus/conduction.F

    r1310 r1530  
    1717c   declarations:
    1818c-----------------------------------------------------------------------
    19 
    20 !#include "dimensions.h"
    21 !#include "dimphys.h"
    22 !#include "comcstfi.h"
    23 !#include "surfdat.h"
    24 !#include "chimiedata.h"
    25 !#include "conc.h"
    2619
    2720c   arguments:
  • trunk/LMDZ.VENUS/libf/phyvenus/diagphy.F

    r1017 r1530  
    5151      implicit none
    5252
    53 #include "dimensions.h"
    5453#include "YOMCST.h"
    5554C
     
    212211      IMPLICIT NONE
    213212C
    214 #include "dimensions.h"
    215213#include "YOMCST.h"
    216214C
  • trunk/LMDZ.VENUS/libf/phyvenus/drag_noro.F

    r1301 r1530  
    1111      use dimphy
    1212      IMPLICIT none
    13 
    14 #include "dimensions.h"
    15 #include "paramet.h"
    1613
    1714c======================================================================
  • trunk/LMDZ.VENUS/libf/phyvenus/euvheat.F90

    r1442 r1530  
    3232!    ------------------
    3333!
    34 !#include "dimensions.h"
    3534#include "YOMCST.h"
    3635#include "clesphys.h"
  • trunk/LMDZ.VENUS/libf/phyvenus/flott_gwd_ran.F90

    r826 r1530  
    1313      use dimphy
    1414      implicit none
    15 
    16 #include "dimensions.h"
    17 #include "paramet.h"
    1815
    1916#include "YOEGWD.h"
  • trunk/LMDZ.VENUS/libf/phyvenus/grid_noro.F

    r1301 r1530  
    5252C=======================================================================
    5353
     54      use mod_grid_phy_lmdz, only: nbp_lon, nbp_lat
    5455      IMPLICIT none
    5556     
    56 #include "dimensions.h"
    5757#include "YOMCST.h"
    5858
     
    7979C INTERMEDIATE FIELDS  (CORRELATIONS OF OROGRAPHY GRADIENT)
    8080
    81       REAL ztz(iim+1,jjm+1),zxtzx(iim+1,jjm+1)
    82       REAL zytzy(iim+1,jjm+1),zxtzy(iim+1,jjm+1)
    83       REAL weight(iim+1,jjm+1)
     81      REAL ztz(nbp_lon+1,nbp_lat),zxtzx(nbp_lon+1,nbp_lat)
     82      REAL zytzy(nbp_lon+1,nbp_lat),zxtzy(nbp_lon+1,nbp_lat)
     83      REAL weight(nbp_lon+1,nbp_lat)
    8484
    8585C CORRELATIONS OF USN OROGRAPHY GRADIENTS
     
    107107c   
    108108c
    109       if(iim.ne.imar) STOP 'Problem dim. x'
    110       if(jjm.ne.jmar-1) STOP 'Problem dim. y'
     109      if(nbp_lon.ne.imar) STOP 'Problem dim. x'
     110      if(nbp_lat-1.ne.jmar-1) STOP 'Problem dim. y'
    111111      IF (imar.GT.2200 .OR. jmar.GT.1100) THEN
    112112         PRINT*, 'imar or jmar too big', imar, jmar
     
    114114      ENDIF
    115115
    116       IF(imar+1.ne.iim+1.or.jmar.ne.jjm+1)THEN
     116      IF(imar+1.ne.nbp_lon+1.or.jmar.ne.nbp_lat)THEN
    117117        print *,' imar or jmar bad dimensions:',imar,jmar
    118118        call abort
     
    297297C  FIRST FILTER, MOVING AVERAGE OVER 9 POINTS.
    298298
    299        CALL MVA9(zmea,iim+1,jjm+1)
    300        CALL MVA9(zstd,iim+1,jjm+1)
    301        CALL MVA9(zpic,iim+1,jjm+1)
    302        CALL MVA9(zval,iim+1,jjm+1)
    303        CALL MVA9(zxtzx,iim+1,jjm+1)
    304        CALL MVA9(zxtzy,iim+1,jjm+1)
    305        CALL MVA9(zytzy,iim+1,jjm+1)
     299       CALL MVA9(zmea,nbp_lon+1,nbp_lat)
     300       CALL MVA9(zstd,nbp_lon+1,nbp_lat)
     301       CALL MVA9(zpic,nbp_lon+1,nbp_lat)
     302       CALL MVA9(zval,nbp_lon+1,nbp_lat)
     303       CALL MVA9(zxtzx,nbp_lon+1,nbp_lat)
     304       CALL MVA9(zxtzy,nbp_lon+1,nbp_lat)
     305       CALL MVA9(zytzy,nbp_lon+1,nbp_lat)
    306306
    307307       DO ii = 1, imar
  • trunk/LMDZ.VENUS/libf/phyvenus/gwprofil.F

    r808 r1530  
    3838      use dimphy
    3939      IMPLICIT NONE
    40 
    41 #include "dimensions.h"
    42 #include "paramet.h"
    4340
    4441#include "YOMCST.h"
  • trunk/LMDZ.VENUS/libf/phyvenus/gwstress.F

    r101 r1530  
    5454      use dimphy
    5555      implicit none
    56 
    57 #include "dimensions.h"
    58 #include "paramet.h"
    5956
    6057#include "YOMCST.h"
  • trunk/LMDZ.VENUS/libf/phyvenus/hgardfou.F

    r101 r1530  
    99c Verifier la temperature
    1010c======================================================================
    11 #include "dimensions.h"
    1211#include "YOMCST.h"
    1312      REAL t(klon,klev), tsol(klon)
  • trunk/LMDZ.VENUS/libf/phyvenus/hrtherm.F

    r1442 r1530  
    1515
    1616
    17 #include "dimensions.h"
    1817#include "param.h"
    1918#include "param_v4.h"
  • trunk/LMDZ.VENUS/libf/phyvenus/ini_histday.h

    r902 r1530  
    2323c
    2424         CALL histdef(nid_day, "phis", "Surface geop. height", "-",
    25      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     25     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    2626     .                "once",  zsto,zout)
    2727c
    2828         CALL histdef(nid_day, "aire", "Grid area", "-",
    29      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     29     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    3030     .                "once",  zsto,zout)
    3131c
    3232         CALL histdef(nid_day, "tsol", "Surface Temperature", "K",
    33      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     33     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    3434     .                "ave(X)", zsto,zout)
    3535c
    3636         CALL histdef(nid_day, "psol", "Surface Pressure", "Pa",
    37      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     37     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    3838     .                "ave(X)", zsto,zout)
    3939c
    4040c        CALL histdef(nid_day, "ue", "Zonal energy transport", "-",
    41 c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     41c    .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    4242c    .                "ave(X)", zsto,zout)
    4343c
    4444c        CALL histdef(nid_day, "ve", "Merid energy transport", "-",
    45 c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     45c    .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    4646c    .                "ave(X)", zsto,zout)
    4747c
    4848c        CALL histdef(nid_day, "cdragh", "Drag coef on T", "-",
    49 c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     49c    .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    5050c    .                "ave(X)", zsto,zout)
    5151c
    5252c        CALL histdef(nid_day, "cdragm", "Drag coef on U", "-",
    53 c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     53c    .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    5454c    .                "ave(X)", zsto,zout)
    5555c
     
    6262c
    6363         CALL histdef(nid_day, "temp", "Air temperature", "K",
    64      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     64     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    6565     .                "ave(X)", zsto,zout)
    6666c
    6767         CALL histdef(nid_day, "pres", "Air pressure", "Pa",
    68      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     68     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    6969     .                "ave(X)", zsto,zout)
    7070c
    7171         CALL histdef(nid_day, "geop", "Geopotential height", "m",
    72      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     72     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    7373     .                "ave(X)", zsto,zout)
    7474c
    7575         CALL histdef(nid_day, "vitu", "Zonal wind", "m/s",
    76      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     76     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    7777     .                "ave(X)", zsto,zout)
    7878c
    7979         CALL histdef(nid_day, "vitv", "Meridional wind", "m/s",
    80      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     80     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    8181     .                "ave(X)", zsto,zout)
    8282c
    8383         CALL histdef(nid_day, "vitw", "Vertical wind", "Pa/s",
    84      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     84     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    8585     .                "ave(X)", zsto,zout)
    8686c
    8787         CALL histdef(nid_day, "dudyn", "Dynamics dU", "m/s2",
    88      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     88     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    8989     .                "ave(X)", zsto,zout)
    9090c
    9191         CALL histdef(nid_day, "duvdf", "Boundary-layer dU", "m/s2",
    92      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     92     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    9393     .                "ave(X)", zsto,zout)
    9494c
    9595c        CALL histdef(nid_day, "mang", "Angular momentum", "kg m2/s",
    96 c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     96c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    9797c    .                "ave(X)", zsto,zout)
    9898c
    9999c        CALL histdef(nid_day, "Kz", "vertical diffusion coef", "m2/s",
    100 c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     100c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    101101c    .                "ave(X)", zsto,zout)
    102102c
     
    107107          WRITE(str2,'(i2.2)') iq
    108108          CALL histdef(nid_day, tname(iq), ttext(iq), "ppm",
    109      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     109     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    110110     .                "ave(X)", zsto,zout)
    111111             ELSE
     
    117117c
    118118         CALL histdef(nid_day, "tops", "Solar rad. at TOA", "W/m2",
    119      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     119     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    120120     .                "ave(X)", zsto,zout)
    121121c
     
    130130c
    131131         CALL histdef(nid_day, "topl", "IR rad. at TOA", "W/m2",
    132      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     132     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    133133     .                "ave(X)", zsto,zout)
    134134c
    135135         CALL histdef(nid_day, "sols", "Solar rad. at surf.", "W/m2",
    136      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     136     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    137137     .                "ave(X)", zsto,zout)
    138138c
    139139         CALL histdef(nid_day, "soll", "IR rad. at surface", "W/m2",
    140      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     140     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    141141     .                "ave(X)", zsto,zout)
    142142c
     
    144144c
    145145         CALL histdef(nid_day, "SWnet", "Net SW flux","W/m2",
    146      .                iim,jj_nb,nhori, klev,1,klev,nvert,
     146     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
    147147     .                32, "ave(X)", zsto,zout)
    148148c
    149149         CALL histdef(nid_day, "LWnet", "Net LW flux","W/m2",
    150      .                iim,jj_nb,nhori, klev,1,klev,nvert,
     150     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
    151151     .                32, "ave(X)", zsto,zout)
    152152c
    153153         CALL histdef(nid_day, "fluxvdf", "PBL net flux","W/m2",
    154      .                iim,jj_nb,nhori, klev,1,klev,nvert,
     154     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
    155155     .                32, "ave(X)", zsto,zout)
    156156c
    157157         CALL histdef(nid_day, "fluxdyn", "Dyn. net flux","W/m2",
    158      .                iim,jj_nb,nhori, klev,1,klev,nvert,
     158     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
    159159     .                32, "ave(X)", zsto,zout)
    160160c
    161161         CALL histdef(nid_day, "fluxajs", "Dry adj. net flux","W/m2",
    162      .                iim,jj_nb,nhori, klev,1,klev,nvert,
     162     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
    163163     .                32, "ave(X)", zsto,zout)
    164164c
    165165c        CALL histdef(nid_day, "fluxec", "Cin. net flux","W/m2",
    166 c    .                iim,jj_nb,nhori, klev,1,klev,nvert,
     166c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
    167167c    .                32, "ave(X)", zsto,zout)
    168168c
     
    173173c
    174174         CALL histdef(nid_day, "dtdyn", "Dynamics dT", "K/s",
    175      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     175     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    176176     .                "ave(X)", zsto,zout)
    177177c
    178178c        CALL histdef(nid_day, "dtphy", "Physics dT", "K/s",
    179 c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     179c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    180180c    .                "ave(X)", zsto,zout)
    181181c
    182182         CALL histdef(nid_day, "dtvdf", "Boundary-layer dT", "K/s",
    183      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     183     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    184184     .                "ave(X)", zsto,zout)
    185185c
    186186         CALL histdef(nid_day, "dtajs", "Dry adjust. dT", "K/s",
    187      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     187     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    188188     .                "ave(X)", zsto,zout)
    189189c
    190190         CALL histdef(nid_day, "dtswr", "SW radiation dT", "K/s",
    191      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     191     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    192192     .                "ave(X)", zsto,zout)
    193193c
    194194         CALL histdef(nid_day, "dtlwr", "LW radiation dT", "K/s",
    195      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     195     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    196196     .                "ave(X)", zsto,zout)
    197197c
    198198c        CALL histdef(nid_day, "dtec", "Cinetic dissip dT", "K/s",
    199 c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     199c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    200200c    .                "ave(X)", zsto,zout)
    201201c
    202202         CALL histdef(nid_mth, "duajs", "Dry convection dU", "m/s2",
    203      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     203     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    204204     .                "ave(X)", zsto,zout)
    205205c
    206206         CALL histdef(nid_mth, "dugwo", "GW oro dU", "m/s2",
    207      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     207     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    208208     .                "ave(X)", zsto,zout)
    209209c
    210210         CALL histdef(nid_mth, "dugwno", "GW non-oro dU", "m/s2",
    211      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     211     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    212212     .                "ave(X)", zsto,zout)
    213213c
    214214c        CALL histdef(nid_mth, "dvvdf", "Boundary-layer dV", "m/s2",
    215 c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     215c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    216216c    .                "ave(X)", zsto,zout)
    217217c
     
    223223c        call histdef(nid_day, "taux",
    224224c    $         "Zonal wind stress", "Pa", 
    225 c    $         iim,jj_nb,nhori, 1,1,1, nvert, 32,
     225c    $         nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    226226c    $         "ave(X)", zsto,zout)
    227227c
    228228c        call histdef(nid_day, "tauy",
    229229c    $         "Meridional xind stress", "Pa", 
    230 c    $         iim,jj_nb,nhori, 1,1,1, nvert, 32,
     230c    $         nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    231231c    $         "ave(X)", zsto,zout)
    232232c
    233233c        CALL histdef(nid_day, "cdrm", "Momentum drag coef.", "-",
    234 c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     234c    .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    235235c    .                "ave(X)", zsto,zout)
    236236c
    237237c        CALL histdef(nid_day, "cdrh", "Heat drag coef.", "-",
    238 c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     238c    .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    239239c    .                "ave(X)", zsto,zout)
    240240c
  • trunk/LMDZ.VENUS/libf/phyvenus/ini_histins.h

    r1468 r1530  
    2323c
    2424         CALL histdef(nid_ins, "phis", "Surface geop. height", "-",
    25      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     25     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    2626     .                "once",  zsto,zout)
    2727c
    2828         CALL histdef(nid_ins, "aire", "Grid area", "-",
    29      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     29     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    3030     .                "once",  zsto,zout)
    3131c
    3232         CALL histdef(nid_ins, "tsol", "Surface Temperature", "K",
    33      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     33     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    3434     .                "ins(X)", zsto,zout)
    3535c
    3636         CALL histdef(nid_ins, "psol", "Surface Pressure", "Pa",
    37      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     37     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    3838     .                "ins(X)", zsto,zout)
    3939c
    4040c        CALL histdef(nid_ins, "ue", "Zonal energy transport", "-",
    41 c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     41c    .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    4242c    .                "ins(X)", zsto,zout)
    4343c
    4444c        CALL histdef(nid_ins, "ve", "Merid energy transport", "-",
    45 c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     45c    .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    4646c    .                "ins(X)", zsto,zout)
    4747c
    4848c        CALL histdef(nid_ins, "cdragh", "Drag coef on T", "-",
    49 c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     49c    .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    5050c    .                "ins(X)", zsto,zout)
    5151c
    5252c        CALL histdef(nid_ins, "cdragm", "Drag coef on U", "-",
    53 c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     53c    .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    5454c    .                "ins(X)", zsto,zout)
    5555c
     
    6262c
    6363         CALL histdef(nid_ins, "temp", "Air temperature", "K",
    64      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     64     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    6565     .                "ins(X)", zsto,zout)
    6666c
    6767         CALL histdef(nid_ins, "pres", "Air pressure", "Pa",
    68      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     68     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    6969     .                "ins(X)", zsto,zout)
    7070c
    7171         CALL histdef(nid_ins, "geop", "Geopotential height", "m",
    72      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     72     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    7373     .                "ins(X)", zsto,zout)
    7474c
    7575         CALL histdef(nid_ins, "vitu", "Zonal wind", "m/s",
    76      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     76     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    7777     .                "ins(X)", zsto,zout)
    7878c
    7979         CALL histdef(nid_ins, "vitv", "Meridional wind", "m/s",
    80      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     80     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    8181     .                "ins(X)", zsto,zout)
    8282c
    8383         CALL histdef(nid_ins, "vitw", "Vertical wind", "Pa/s",
    84      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     84     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    8585     .                "ins(X)", zsto,zout)
    8686c
    8787         CALL histdef(nid_ins, "dudyn", "Dynamics dU", "m/s2",
    88      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     88     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    8989     .                "ins(X)", zsto,zout)
    9090c
    9191         CALL histdef(nid_ins, "duvdf", "Boundary-layer dU", "m/s2",
    92      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     92     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    9393     .                "ins(X)", zsto,zout)
    9494c
    9595c        CALL histdef(nid_ins, "mang", "Angular momentum", "kg m2/s",
    96 c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     96c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    9797c    .                "ins(X)", zsto,zout)
    9898c
    9999c        CALL histdef(nid_ins, "Kz", "vertical diffusion coef", "m2/s",
    100 c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     100c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    101101c    .                "ins(X)", zsto,zout)
    102102c
    103103         CALL histdef(nid_ins, "mmean", "Mean molecular mass", "g/mol",
    104      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     104     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    105105     .                "ins(X)", zsto,zout)
    106106
    107107         CALL histdef(nid_ins, "rho", "Air density [mass/Vol]", "kg/m3",
    108      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     108     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    109109     .                "ins(X)", zsto,zout)
    110110
     
    116116          WRITE(str2,'(i2.2)') iq
    117117          CALL histdef(nid_ins, tname(iq), ttext(iq), "mol/mol",
    118      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     118     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    119119     .                "ins(X)", zsto,zout)
    120120             ELSE
     
    127127       IF (callthermos .and. ok_chem) THEN
    128128          CALL histdef(nid_ins, "d_qmoldif CO2", "Dif molec" , "kg/kg",
    129      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     129     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    130130     .                "ins(X)", zsto,zout)
    131131          CALL histdef(nid_ins, "d_qmoldif O3p", "Dif molec" , "kg/kg",
    132      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     132     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    133133     .                "ins(X)", zsto,zout)
    134134          CALL histdef(nid_ins, "d_qmoldif N2", "Dif molec" , "kg/kg",
    135      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     135     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    136136     .                "ins(X)", zsto,zout)
    137137       ENDIF
    138138c
    139139         CALL histdef(nid_ins, "tops", "Solar rad. at TOA", "W/m2",
    140      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     140     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    141141     .                "ins(X)", zsto,zout)
    142142c
     
    147147c
    148148         CALL histdef(nid_ins, "NBRTOTm1", "Nbr total droplet",
    149      .                "#/cm3", iim,jj_nb,nhori, klev,1,klev, nvert, 32,
     149     .                "#/cm3",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    150150     .                "ins(X)", zsto,zout)
    151151c
     
    155155c     .    for log normal distribution" ,
    156156c     .                "fraction",
    157 c     .                iim,jj_nb,nhori, klev,1,klev, nvert, 32,
     157c     .                nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32,
    158158c     .                "ins(X)", zsto,zout)
    159159c
     
    163163c     .    for log normal distribution",
    164164c     .                "fraction",
    165 c     .                iim,jj_nb,nhori, klev,1,klev, nvert, 32,
     165c     .                nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32,
    166166c     .                "ins(X)", zsto,zout)
    167167c
     
    171171c
    172172         CALL histdef(nid_ins, "NBRTOTm2", "Nbr total droplet",
    173      .                "#/cm3", iim,jj_nb,nhori, klev,1,klev, nvert, 32,
     173     .                "#/cm3",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    174174     .                "ins(X)", zsto,zout)
    175175c
     
    179179c     .    for log normal distribution" ,
    180180c     .                "fraction",
    181 c     .                iim,jj_nb,nhori, klev,1,klev, nvert, 32,
     181c     .                nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32,
    182182c     .                "ins(X)", zsto,zout)
    183183c
     
    187187c     .    for log normal distribution",
    188188c     .                "fraction",
    189 c     .                iim,jj_nb,nhori, klev,1,klev, nvert, 32,
     189c     .                nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32,
    190190c     .                "ins(X)", zsto,zout)
    191191c
     
    195195c
    196196         CALL histdef(nid_ins, "NBRTOTm3", "Nbr total droplet", "#/cm3",
    197      .                iim,jj_nb,nhori, klev,1,klev, nvert, 32,
     197     .                nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32,
    198198     .                "ins(X)", zsto,zout)
    199199c
     
    203203c     .   for log normal distribution" ,
    204204c     .                "fraction",
    205 c     .                iim,jj_nb,nhori, klev,1,klev, nvert, 32,
     205c     .                nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32,
    206206c     .                "ins(X)", zsto,zout)
    207207c
     
    211211c     .    for log normal distribution",
    212212c     .                "fraction",
    213 c     .                iim,jj_nb,nhori, klev,1,klev, nvert, 32,
     213c     .                nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32,
    214214c     .                "ins(X)", zsto,zout)
    215215c
     
    222222         CALL histdef(nid_ins, "WH2SO4", "Weight fraction H2SO4",
    223223     .                "fraction",
    224      .                iim,jj_nb,nhori, klev,1,klev, nvert, 32,
     224     .                nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32,
    225225     .                "ins(X)", zsto,zout)
    226226c
     
    229229         CALL histdef(nid_ins, "rho_droplet", "density cloud droplet",
    230230     .                "kg.m-3",
    231      .                iim,jj_nb,nhori, klev,1,klev, nvert, 32,
     231     .                nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32,
    232232     .                "ins(X)", zsto,zout)
    233233c
     
    239239         CALL histdef(nid_ins, "d_tr_sed_H2SO4", "var mmr from sedim",
    240240     .                "kg/kg",
    241      .                iim,jj_nb,nhori, klev,1,klev, nvert, 32,
     241     .                nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32,
    242242     .                "ins(X)", zsto,zout)
    243243c
     
    246246         CALL histdef(nid_ins, "d_tr_sed_H2O", "var mmr from sedim",
    247247     .                "kg/kg",
    248      .                iim,jj_nb,nhori, klev,1,klev, nvert, 32,
     248     .                nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32,
    249249     .                "ins(X)", zsto,zout)
    250250c
     
    253253         CALL histdef(nid_ins, "F_sedim", "tendency from sedim",
    254254     .                "kg.m-2.s-1",
    255      .                iim,jj_nb,nhori, klev,1,klev, nvert, 32,
     255     .                nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32,
    256256     .                "ins(X)", zsto,zout)
    257257c
     
    267267c
    268268         CALL histdef(nid_ins, "topl", "IR rad. at TOA", "W/m2",
    269      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     269     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    270270     .                "ins(X)", zsto,zout)
    271271c
    272272         CALL histdef(nid_ins, "sols", "Solar rad. at surf.", "W/m2",
    273      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     273     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    274274     .                "ins(X)", zsto,zout)
    275275c
    276276         CALL histdef(nid_ins, "soll", "IR rad. at surface", "W/m2",
    277      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     277     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    278278     .                "ins(X)", zsto,zout)
    279279c
     
    281281c
    282282         CALL histdef(nid_ins, "SWnet", "Net SW flux","W/m2",
    283      .                iim,jj_nb,nhori, klev,1,klev,nvert,
     283     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
    284284     .                32, "ins(X)", zsto,zout)
    285285c
    286286         CALL histdef(nid_ins, "LWnet", "Net LW flux","W/m2",
    287      .                iim,jj_nb,nhori, klev,1,klev,nvert,
     287     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
    288288     .                32, "ins(X)", zsto,zout)
    289289c
    290290         CALL histdef(nid_ins, "fluxvdf", "PBL net flux","W/m2",
    291      .                iim,jj_nb,nhori, klev,1,klev,nvert,
     291     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
    292292     .                32, "ins(X)", zsto,zout)
    293293c
    294294         CALL histdef(nid_ins, "fluxdyn", "Dyn. net flux","W/m2",
    295      .                iim,jj_nb,nhori, klev,1,klev,nvert,
     295     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
    296296     .                32, "ins(X)", zsto,zout)
    297297c
    298298         CALL histdef(nid_ins, "fluxajs", "Dry adj. net flux","W/m2",
    299      .                iim,jj_nb,nhori, klev,1,klev,nvert,
     299     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
    300300     .                32, "ins(X)", zsto,zout)
    301301c
    302302c        CALL histdef(nid_ins, "fluxec", "Cin. net flux","W/m2",
    303 c    .                iim,jj_nb,nhori, klev,1,klev,nvert,
     303c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
    304304c    .                32, "ins(X)", zsto,zout)
    305305c
     
    310310c
    311311         CALL histdef(nid_ins, "dtdyn", "Dynamics dT", "K/s",
    312      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     312     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    313313     .                "ins(X)", zsto,zout)
    314314c
    315315c        CALL histdef(nid_ins, "dtphy", "Physics dT", "K/s",
    316 c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     316c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    317317c    .                "ins(X)", zsto,zout)
    318318c
    319319         CALL histdef(nid_ins, "dtvdf", "Boundary-layer dT", "K/s",
    320      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     320     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    321321     .                "ins(X)", zsto,zout)
    322322c
    323323         CALL histdef(nid_ins, "dtajs", "Dry adjust. dT", "K/s",
    324      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     324     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    325325     .                "ins(X)", zsto,zout)
    326326c
    327327         CALL histdef(nid_ins, "dtswr", "SW radiation dT", "K/s",
    328      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     328     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    329329     .                "ins(X)", zsto,zout)
    330330c
    331331         CALL histdef(nid_ins, "dtlwr", "LW radiation dT", "K/s",
    332      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     332     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    333333     .                "ins(X)", zsto,zout)
    334334c
    335335c        CALL histdef(nid_ins, "dtec", "Cinetic dissip dT", "K/s",
    336 c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     336c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    337337c    .                "ins(X)", zsto,zout)
    338338c
    339339         CALL histdef(nid_ins, "duajs", "Dry convection dU", "m/s2",
    340      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     340     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    341341     .                "ins(X)", zsto,zout)
    342342c
    343343         CALL histdef(nid_ins, "dugwo", "GW oro dU", "m/s2",
    344      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     344     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    345345     .                "ins(X)", zsto,zout)
    346346c
    347347         CALL histdef(nid_ins, "dugwno", "GW non-oro dU", "m/s2",
    348      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     348     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    349349     .                "ins(X)", zsto,zout)
    350350c
    351351c        CALL histdef(nid_ins, "dvvdf", "Boundary-layer dV", "m/s2",
    352 c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     352c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    353353c    .                "ins(X)", zsto,zout)
    354354c
     
    360360c        call histdef(nid_ins, "taux",
    361361c    $         "Zonal wind stress", "Pa", 
    362 c    $         iim,jj_nb,nhori, 1,1,1, nvert, 32,
     362c    $         nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    363363c    $         "ins(X)", zsto,zout)
    364364c
    365365c        call histdef(nid_ins, "tauy",
    366366c    $         "Meridional xind stress", "Pa", 
    367 c    $         iim,jj_nb,nhori, 1,1,1, nvert, 32,
     367c    $         nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    368368c    $         "ins(X)", zsto,zout)
    369369c
    370370c        CALL histdef(nid_ins, "cdrm", "Momentum drag coef.", "-",
    371 c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     371c    .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    372372c    .                "ins(X)", zsto,zout)
    373373c
    374374c        CALL histdef(nid_ins, "cdrh", "Heat drag coef.", "-",
    375 c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     375c    .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    376376c    .                "ins(X)", zsto,zout)
    377377c
  • trunk/LMDZ.VENUS/libf/phyvenus/ini_histmth.h

    r1518 r1530  
    2323c
    2424         CALL histdef(nid_mth, "phis", "Surface geop. height", "-",
    25      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     25     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    2626     .                "once",  zsto,zout)
    2727c
    2828         CALL histdef(nid_mth, "aire", "Grid area", "-",
    29      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     29     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    3030     .                "once",  zsto,zout)
    3131c
    3232         CALL histdef(nid_mth, "tsol", "Surface Temperature", "K",
    33      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     33     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    3434     .                "ave(X)", zsto,zout)
    3535c
    3636         CALL histdef(nid_mth, "psol", "Surface Pressure", "Pa",
    37      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     37     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    3838     .                "ave(X)", zsto,zout)
    3939c
    4040c        CALL histdef(nid_mth, "ue", "Zonal energy transport", "-",
    41 c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     41c    .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    4242c    .                "ave(X)", zsto,zout)
    4343c
    4444c        CALL histdef(nid_mth, "ve", "Merid energy transport", "-",
    45 c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     45c    .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    4646c    .                "ave(X)", zsto,zout)
    4747c
    4848c        CALL histdef(nid_mth, "cdragh", "Drag coef on T", "-",
    49 c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     49c    .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    5050c    .                "ave(X)", zsto,zout)
    5151c
    5252c        CALL histdef(nid_mth, "cdragm", "Drag coef on U", "-",
    53 c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     53c    .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    5454c    .                "ave(X)", zsto,zout)
    5555c
     
    6262c
    6363         CALL histdef(nid_mth, "temp", "Air temperature", "K",
    64      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     64     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    6565     .                "ave(X)", zsto,zout)
    6666c
    6767         CALL histdef(nid_mth, "pres", "Air pressure", "Pa",
    68      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     68     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    6969     .                "ave(X)", zsto,zout)
    7070c
    7171         CALL histdef(nid_mth, "geop", "Geopotential height", "m",
    72      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     72     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    7373     .                "ave(X)", zsto,zout)
    7474c
    7575         CALL histdef(nid_mth, "vitu", "Zonal wind", "m/s",
    76      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     76     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    7777     .                "ave(X)", zsto,zout)
    7878c
    7979         CALL histdef(nid_mth, "vitv", "Meridional wind", "m/s",
    80      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     80     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    8181     .                "ave(X)", zsto,zout)
    8282c
    8383         CALL histdef(nid_mth, "vitw", "Vertical wind", "Pa/s",
    84      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     84     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    8585     .                "ave(X)", zsto,zout)
    8686c
    8787         CALL histdef(nid_mth, "dudyn", "Dynamics dU", "m/s2",
    88      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     88     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    8989     .                "ave(X)", zsto,zout)
    9090c
    9191         CALL histdef(nid_mth, "duvdf", "Boundary-layer dU", "m/s2",
    92      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     92     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    9393     .                "ave(X)", zsto,zout)
    9494c
    9595c        CALL histdef(nid_mth, "mang", "Angular momentum", "kg m2/s",
    96 c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     96c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    9797c    .                "ave(X)", zsto,zout)
    9898c
    9999c        CALL histdef(nid_mth, "Kz", "vertical diffusion coef", "m2/s",
    100 c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     100c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    101101c    .                "ave(X)", zsto,zout)
    102102c
    103103         CALL histdef(nid_mth, "mmean", "Mean molecular mass", "g/mol",
    104      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     104     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    105105     .                "ave(X)", zsto,zout)
    106106
    107107c        CALL histdef(nid_mth, "rho", "Air density [mass/Vol]", "kg/m3",
    108 c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     108c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    109109c    .                "ave(X)", zsto,zout)
    110110
     
    116116          WRITE(str2,'(i2.2)') iq
    117117          CALL histdef(nid_mth, tname(iq), ttext(iq), "mol/mol",
    118      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     118     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    119119     .                "ave(X)", zsto,zout)
    120120             ELSE
     
    127127       IF (callthermos .and. ok_chem) THEN
    128128          CALL histdef(nid_mth, "d_qmoldifCO2", "Dif molec" , "kg/kg",
    129      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     129     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    130130     .                "ave(X)", zsto,zout)
    131131          CALL histdef(nid_mth, "d_qmoldifO3p", "Dif molec" , "kg/kg",
    132      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     132     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    133133     .                "ave(X)", zsto,zout)
    134134          CALL histdef(nid_mth, "d_qmoldifN2", "Dif molec" , "kg/kg",
    135      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     135     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    136136     .                "ave(X)", zsto,zout)
    137137       ENDIF
    138138c
    139139         CALL histdef(nid_mth, "tops", "Solar rad. at TOA", "W/m2",
    140      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     140     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    141141     .                "ave(X)", zsto,zout)
    142142c
     
    151151c
    152152         CALL histdef(nid_mth, "topl", "IR rad. at TOA", "W/m2",
    153      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     153     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    154154     .                "ave(X)", zsto,zout)
    155155c
    156156         CALL histdef(nid_mth, "sols", "Solar rad. at surf.", "W/m2",
    157      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     157     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    158158     .                "ave(X)", zsto,zout)
    159159c
    160160         CALL histdef(nid_mth, "soll", "IR rad. at surface", "W/m2",
    161      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     161     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    162162     .                "ave(X)", zsto,zout)
    163163c
     
    165165c
    166166         CALL histdef(nid_mth, "SWnet", "Net SW flux","W/m2",
    167      .                iim,jj_nb,nhori, klev,1,klev,nvert,
     167     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
    168168     .                32, "ave(X)", zsto,zout)
    169169c
    170170         CALL histdef(nid_mth, "LWnet", "Net LW flux","W/m2",
    171      .                iim,jj_nb,nhori, klev,1,klev,nvert,
     171     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
    172172     .                32, "ave(X)", zsto,zout)
    173173c
    174174c        CALL histdef(nid_mth, "fluxvdf", "PBL net flux","W/m2",
    175 c    .                iim,jj_nb,nhori, klev,1,klev,nvert,
     175c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
    176176c    .                32, "ave(X)", zsto,zout)
    177177c
    178178c        CALL histdef(nid_mth, "fluxdyn", "Dyn. net flux","W/m2",
    179 c    .                iim,jj_nb,nhori, klev,1,klev,nvert,
     179c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
    180180c    .                32, "ave(X)", zsto,zout)
    181181c
    182182c        CALL histdef(nid_mth, "fluxajs", "Dry adj. net flux","W/m2",
    183 c    .                iim,jj_nb,nhori, klev,1,klev,nvert,
     183c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
    184184c    .                32, "ave(X)", zsto,zout)
    185185c
    186186c        CALL histdef(nid_mth, "fluxec", "Cin. net flux","W/m2",
    187 c    .                iim,jj_nb,nhori, klev,1,klev,nvert,
     187c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
    188188c    .                32, "ave(X)", zsto,zout)
    189189c
     
    194194c
    195195         CALL histdef(nid_mth, "dtdyn", "Dynamics dT", "K/s",
    196      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     196     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    197197     .                "ave(X)", zsto,zout)
    198198c
    199199c        CALL histdef(nid_mth, "dtphy", "Physics dT", "K/s",
    200 c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     200c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    201201c    .                "ave(X)", zsto,zout)
    202202c
    203203         CALL histdef(nid_mth, "dtvdf", "Boundary-layer dT", "K/s",
    204      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     204     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    205205     .                "ave(X)", zsto,zout)
    206206c
    207207         CALL histdef(nid_mth, "dtajs", "Dry adjust. dT", "K/s",
    208      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     208     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    209209     .                "ave(X)", zsto,zout)
    210210c
    211211         CALL histdef(nid_mth, "dtswr", "SW radiation dT", "K/s",
    212      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     212     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    213213     .                "ave(X)", zsto,zout)
    214214c
    215215         CALL histdef(nid_mth,"dtswrNLTE", "SWNLTE radiation dT",
    216      .                "K/s", iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     216     .                "K/s",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    217217     .                "ave(X)", zsto,zout)
    218218c
    219219         CALL histdef(nid_mth,"dtswrDCrisp","SWDCrisp radiation dT",
    220      .                "K/s", iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     220     .                "K/s",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    221221     .                "ave(X)", zsto,zout)
    222222c
    223223         CALL histdef(nid_mth, "dtlwr", "LW radiation dT", "K/s",
    224      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     224     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    225225     .                "ave(X)", zsto,zout)
    226226c
    227227         CALL histdef(nid_mth,"dtlwrNLTE", "LWNLTE radiation dT",
    228      .                "K/s", iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     228     .                "K/s",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    229229     .                "ave(X)", zsto,zout)
    230230c
    231231         CALL histdef(nid_mth,"dtlwrLTE","LW_LTE radiation dT",
    232      .                "K/s", iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     232     .                "K/s",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    233233     .                "ave(X)", zsto,zout)
    234234c
    235235         CALL histdef(nid_mth, "dteuv", "UV radiation dT", "K/s",
    236      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     236     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    237237     .                "ave(X)", zsto,zout)
    238238         CALL histdef(nid_mth, "dtcond", "Therm conduction", "K/s",
    239      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     239     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    240240     .                "ave(X)", zsto,zout)
    241241         CALL histdef(nid_mth, "dumolvis", "molec viscosity (u)"
    242      .                ,"m/s2", iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     242     .                ,"m/s2",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    243243     .                "ave(X)", zsto,zout)
    244244         CALL histdef(nid_mth, "dvmolvis", "molec viscosity (v)",
    245      .                "m/s2", iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     245     .                "m/s2",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    246246     .                "ave(X)", zsto,zout)
    247247
    248248c
    249249c        CALL histdef(nid_mth, "dtec", "Cinetic dissip dT", "K/s",
    250 c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     250c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    251251c    .                "ave(X)", zsto,zout)
    252252c
    253253         CALL histdef(nid_mth, "duajs", "Dry convection dU", "m/s2",
    254      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     254     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    255255     .                "ave(X)", zsto,zout)
    256256c
    257257         CALL histdef(nid_mth, "dugwo", "GW oro dU", "m/s2",
    258      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     258     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    259259     .                "ave(X)", zsto,zout)
    260260c
    261261         CALL histdef(nid_mth, "dugwno", "GW non-oro dU", "m/s2",
    262      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     262     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    263263     .                "ave(X)", zsto,zout)
    264264
    265265c         CALL histdef(nid_mth, "co2_vmr", "volume mixture ratio co2",
    266 c     .         "cm-3",iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     266c     .         "cm-3",nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    267267c     .                "ave(X)", zsto,zout)
    268268
    269269c         CALL histdef(nid_mth, "o_vmr", "volume mixture ratio o",
    270 c     .                "cm-3", iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     270c     .                "cm-3",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    271271c     .                "ave(X)", zsto,zout)
    272272c         CALL histdef(nid_mth, "co_vmr", "volume mixture ratio co",
    273 c     .         "cm-3",iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     273c     .         "cm-3",nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    274274c     .                "ave(X)", zsto,zout)
    275275
    276276c         CALL histdef(nid_mth, "n2_vmr", "volume mixture ratio n2",
    277 c     .                "cm-3", iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     277c     .                "cm-3",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    278278c     .                "ave(X)", zsto,zout)
    279279c
    280280c        CALL histdef(nid_mth, "dvvdf", "Boundary-layer dV", "m/s2",
    281 c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     281c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    282282c    .                "ave(X)", zsto,zout)
    283283c
     
    289289c        call histdef(nid_mth, "taux",
    290290c    $         "Zonal wind stress", "Pa", 
    291 c    $         iim,jj_nb,nhori, 1,1,1, nvert, 32,
     291c    $         nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    292292c    $         "ave(X)", zsto,zout)
    293293c
    294294c        call histdef(nid_mth, "tauy",
    295295c    $         "Meridional xind stress", "Pa", 
    296 c    $         iim,jj_nb,nhori, 1,1,1, nvert, 32,
     296c    $         nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    297297c    $         "ave(X)", zsto,zout)
    298298c
    299299c        CALL histdef(nid_mth, "cdrm", "Momentum drag coef.", "-",
    300 c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     300c    .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    301301c    .                "ave(X)", zsto,zout)
    302302c
    303303c        CALL histdef(nid_mth, "cdrh", "Heat drag coef.", "-",
    304 c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     304c    .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    305305c    .                "ave(X)", zsto,zout)
    306306c
  • trunk/LMDZ.VENUS/libf/phyvenus/iophy.F90

    r892 r1530  
    2424
    2525  subroutine init_iophy_new(rlat,rlon)
    26   USE dimphy
     26  USE dimphy, only: klon
    2727  USE mod_phys_lmdz_para
    28   USE mod_grid_phy_lmdz
     28  USE mod_grid_phy_lmdz, only: nbp_lon, nbp_lat, klon_glo
    2929  USE ioipsl
    3030  implicit none
    31   include 'dimensions.h'   
    3231    real,dimension(klon),intent(in) :: rlon
    3332    real,dimension(klon),intent(in) :: rlat
     
    5150   
    5251!$OMP MASTER 
    53     ALLOCATE(io_lat(jjm+1-1/(iim*jjm)))
    54     io_lat(1)=rlat_glo(1)
    55     io_lat(jjm+1-1/(iim*jjm))=rlat_glo(klon_glo)
    56     IF ((iim*jjm) > 1) then
    57       DO i=2,jjm
    58         io_lat(i)=rlat_glo(2+(i-2)*iim)
     52    ALLOCATE(io_lat(nbp_lat))
     53    IF (klon_glo == 1) THEN
     54      io_lat(1)=rlat_glo(1)
     55    ELSE
     56      io_lat(1)=rlat_glo(1)
     57      io_lat(nbp_lat)=rlat_glo(klon_glo)
     58      DO i=2,nbp_lat-1
     59        io_lat(i)=rlat_glo(2+(i-2)*nbp_lon)
    5960      ENDDO
    6061    ENDIF
    6162
    62     ALLOCATE(io_lon(iim))
    63     io_lon(:)=rlon_glo(2-1/(iim*jjm):iim+1-1/(iim*jjm))
     63    ALLOCATE(io_lon(nbp_lon))
     64    IF (klon_glo == 1) THEN
     65      io_lon(1)=rlon_glo(1)
     66    ELSE
     67      io_lon(1:nbp_lon)=rlon_glo(2:nbp_lon+1)
     68    ENDIF
    6469
    6570! POUR VENUS, LA CARTE EST A L ENVERS !!
     
    6873
    6974    ddid=(/ 1,2 /)
    70     dsg=(/ iim, jjm+1-1/(iim*jjm) /)
    71     dsl=(/ iim, jj_nb /)
     75    dsg=(/ nbp_lon, nbp_lat /)
     76    dsl=(/ nbp_lon, jj_nb /)
    7277    dpf=(/ 1,jj_begin /)
    73     dpl=(/ iim, jj_end /)
     78    dpl=(/ nbp_lon, jj_end /)
    7479    dhs=(/ ii_begin-1,0 /)
    75     if (mpi_rank==mpi_size-1) then
     80    IF (mpi_rank==mpi_size-1) THEN
    7681      dhe=(/0,0/)
    77     else
    78       dhe=(/ iim-ii_end,0 /) 
    79     endif
     82    ELSE
     83      dhe=(/ nbp_lon-ii_end,0 /) 
     84    ENDIF
    8085   
    8186    call flio_dom_set(mpi_size,mpi_rank,ddid,dsg,dsl,dpf,dpl,dhs,dhe, &
     
    8994  USE dimphy
    9095  USE mod_phys_lmdz_para
    91   use ioipsl
    92   implicit none
    93   include 'dimensions.h'   
    94     real,dimension(iim),intent(in) :: lon
    95     real,dimension(jjm+1-1/(iim*jjm)),intent(in) :: lat
     96  use ioipsl, only: flio_dom_set
     97  USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat
     98  implicit none
     99    real,dimension(nbp_lon),intent(in) :: lon
     100    real,dimension(nbp_lat),intent(in) :: lat
    96101
    97102    INTEGER,DIMENSION(2) :: ddid
     
    104109
    105110!$OMP MASTER 
    106     allocate(io_lat(jjm+1-1/(iim*jjm)))
     111    allocate(io_lat(nbp_lat))
    107112    io_lat(:)=lat(:)
    108     allocate(io_lon(iim))
     113    allocate(io_lon(nbp_lon))
    109114    io_lon(:)=lon(:)
    110115   
    111116    ddid=(/ 1,2 /)
    112     dsg=(/ iim, jjm+1-1/(iim*jjm) /)
    113     dsl=(/ iim, jj_nb /)
     117    dsg=(/ nbp_lon, nbp_lat /)
     118    dsl=(/ nbp_lon, jj_nb /)
    114119    dpf=(/ 1,jj_begin /)
    115     dpl=(/ iim, jj_end /)
     120    dpl=(/ nbp_lon, jj_end /)
    116121    dhs=(/ ii_begin-1,0 /)
    117122    if (mpi_rank==mpi_size-1) then
    118123      dhe=(/0,0/)
    119124    else
    120       dhe=(/ iim-ii_end,0 /) 
     125      dhe=(/ nbp_lon-ii_end,0 /) 
    121126    endif
    122127   
     
    129134 
    130135  subroutine histbeg_phy(name,itau0,zjulian,dtime,nhori,nid_day)
    131   USE dimphy
    132   USE mod_phys_lmdz_para
    133   use ioipsl
     136  USE mod_phys_lmdz_para, only: jj_begin, jj_end, jj_nb, is_sequential
     137  USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat
     138  use ioipsl, only: histbeg
    134139  use write_field
    135140  implicit none
    136   include 'dimensions.h'
    137141   
    138142    character*(*), intent(IN) :: name
     
    145149!$OMP MASTER   
    146150    if (is_sequential) then
    147       call histbeg(name,iim,io_lon, jj_nb,io_lat(jj_begin:jj_end), &
    148                    1,iim,1,jj_nb,itau0, zjulian, dtime, nhori, nid_day)
     151      call histbeg(name,nbp_lon,io_lon, jj_nb,io_lat(jj_begin:jj_end), &
     152                   1,nbp_lon,1,jj_nb,itau0, zjulian, dtime, nhori, nid_day)
    149153    else
    150       call histbeg(name,iim,io_lon, jj_nb,io_lat(jj_begin:jj_end), &
    151                    1,iim,1,jj_nb,itau0, zjulian, dtime, nhori, nid_day,phys_domain_id)
     154      call histbeg(name,nbp_lon,io_lon, jj_nb,io_lat(jj_begin:jj_end), &
     155                   1,nbp_lon,1,jj_nb,itau0, zjulian, dtime, nhori, nid_day,phys_domain_id)
    152156    endif
    153157!$OMP END MASTER
     
    158162             plon,plat,plon_bounds,plat_bounds, &
    159163             nname,itau0,zjulian,dtime,nnhori,nnid_day)
    160   USE dimphy
     164  USE dimphy, only: klon
    161165  USE mod_phys_lmdz_para
    162   USE mod_grid_phy_lmdz
    163   use ioipsl
    164   use write_field
    165   implicit none
    166   include 'dimensions.h'
     166  USE mod_grid_phy_lmdz, only: klon_glo, nbp_lon, nbp_lat
     167  use ioipsl, only: histbeg
     168
     169  implicit none
    167170
    168171    real,dimension(klon),intent(in) :: rlon
     
    190193    REAL, allocatable, dimension(:) :: npplat, npplon
    191194    REAL, allocatable, dimension(:,:) :: npplat_bounds, npplon_bounds
    192     INTEGER, PARAMETER :: jjmp1=jjm+1-1/jjm
    193     REAL, dimension(iim,jjmp1) :: zx_lon, zx_lat
     195    REAL, dimension(nbp_lon,nbp_lat) :: zx_lon, zx_lat
    194196
    195197    CALL gather(rlat,rlat_glo)
     
    216218     endif
    217219!
    218      IF ( tabij(i).LE.iim) THEN
     220     IF ( tabij(i).LE.nbp_lon) THEN
    219221      plat_bounds(i,1)=rlat_glo(tabij(i))
    220222     ELSE
    221       plat_bounds(i,1)=rlat_glo(tabij(i)-iim)
     223      plat_bounds(i,1)=rlat_glo(tabij(i)-nbp_lon)
    222224     ENDIF
    223      plat_bounds(i,2)=rlat_glo(tabij(i)+iim)
     225     plat_bounds(i,2)=rlat_glo(tabij(i)+nbp_lon)
    224226!
    225227!    print*,'CFMIP_iophy point i lon lon_bds',i,plon_bounds(i,1),rlon_glo(tabij(i)),plon_bounds(i,2)
     
    237239     ENDDO
    238240
    239        CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon_glo,zx_lon)
    240        if ((iim*jjm).gt.1) then
    241        DO i = 1, iim
     241       CALL gr_fi_ecrit(1,klon,nbp_lon,nbp_lat,rlon_glo,zx_lon)
     242       if ((nbp_lon*nbp_lat).gt.1) then
     243       DO i = 1, nbp_lon
    242244         zx_lon(i,1) = rlon_glo(i+1)
    243          zx_lon(i,jjmp1) = rlon_glo(i+1)
     245         zx_lon(i,nbp_lat) = rlon_glo(i+1)
    244246       ENDDO
    245247       endif
    246        CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat_glo,zx_lat)
     248       CALL gr_fi_ecrit(1,klon,nbp_lon,nbp_lat,rlat_glo,zx_lat)
    247249
    248250    DO i=1,pim
     
    253255
    254256     if (ipt(i).EQ.1) then
    255       plon_bounds(i,1)=zx_lon(iim,jpt(i))
     257      plon_bounds(i,1)=zx_lon(nbp_lon,jpt(i))
    256258      plon_bounds(i,2)=360.+zx_lon(ipt(i)+1,jpt(i))
    257259     endif
    258260 
    259      if (ipt(i).EQ.iim) then
     261     if (ipt(i).EQ.nbp_lon) then
    260262      plon_bounds(i,2)=360.+zx_lon(1,jpt(i))
    261263     endif
     
    269271     endif
    270272 
    271      if (jpt(i).EQ.jjmp1) then
    272       plat_bounds(i,1)=zx_lat(ipt(i),jjmp1)+0.001
    273       plat_bounds(i,2)=zx_lat(ipt(i),jjmp1)-0.001
     273     if (jpt(i).EQ.nbp_lat) then
     274      plat_bounds(i,1)=zx_lat(ipt(i),nbp_lat)+0.001
     275      plat_bounds(i,2)=zx_lat(ipt(i),nbp_lat)-0.001
    274276     endif
    275277!
     
    327329 
    328330  subroutine histwrite2d_phy(nid,lpoint,name,itau,field)
    329   USE dimphy
     331  USE dimphy, only: klon
    330332  USE mod_phys_lmdz_para
    331   USE ioipsl
    332   implicit none
    333   include 'dimensions.h'
     333  USE ioipsl, only: histwrite
     334  USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat
     335  implicit none
    334336   
    335337    integer,intent(in) :: nid
     
    340342    REAL,dimension(klon_mpi) :: buffer_omp
    341343    INTEGER, allocatable, dimension(:) :: index2d
    342     REAL :: Field2d(iim,jj_nb)
     344    REAL :: Field2d(nbp_lon,jj_nb)
    343345
    344346    integer :: ip
     
    351353    CALL grid1Dto2D_mpi(buffer_omp,Field2d)
    352354    if(.NOT.lpoint) THEN
    353      ALLOCATE(index2d(iim*jj_nb))
    354      ALLOCATE(fieldok(iim*jj_nb))
    355      CALL histwrite(nid,name,itau,Field2d,iim*jj_nb,index2d)
     355     ALLOCATE(index2d(nbp_lon*jj_nb))
     356     ALLOCATE(fieldok(nbp_lon*jj_nb))
     357     CALL histwrite(nid,name,itau,Field2d,nbp_lon*jj_nb,index2d)
    356358    else
    357359     ALLOCATE(fieldok(npstn))
     
    382384
    383385  subroutine histwrite3d_phy(nid,lpoint,name,itau,field)
    384   USE dimphy
     386  USE dimphy, only: klon
    385387  USE mod_phys_lmdz_para
    386 
    387   use ioipsl
    388   implicit none
    389   include 'dimensions.h'
     388  USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat
     389  use ioipsl, only: histwrite
     390  implicit none
    390391   
    391392    integer,intent(in) :: nid
     
    395396    real,dimension(:,:),intent(in) :: field  ! --> field(klon,:)
    396397    REAL,dimension(klon_mpi,size(field,2)) :: buffer_omp
    397     REAL :: Field3d(iim,jj_nb,size(field,2))
     398    REAL :: Field3d(nbp_lon,jj_nb,size(field,2))
    398399    INTEGER :: ip, n, nlev
    399400    INTEGER, ALLOCATABLE, dimension(:) :: index3d
     
    413414    CALL grid1Dto2D_mpi(buffer_omp,field3d)
    414415    if(.NOT.lpoint) THEN
    415      ALLOCATE(index3d(iim*jj_nb*nlev))
    416      ALLOCATE(fieldok(iim*jj_nb,nlev))
    417      CALL histwrite(nid,name,itau,Field3d,iim*jj_nb*nlev,index3d)
     416     ALLOCATE(index3d(nbp_lon*jj_nb*nlev))
     417     ALLOCATE(fieldok(nbp_lon*jj_nb,nlev))
     418     CALL histwrite(nid,name,itau,Field3d,nbp_lon*jj_nb*nlev,index3d)
    418419    else
    419420      nlev=size(field,2)
  • trunk/LMDZ.VENUS/libf/phyvenus/jthermcalc.F

    r1442 r1530  
    1818
    1919c     common variables and constants
    20       include "dimensions.h"
    2120      include "param.h"
    2221      include "param_v4.h"
     
    10621061
    10631062c     common variables and constants
    1064 #include "dimensions.h"
    10651063c#include "tracer.h"
    10661064#include "param.h"
     
    14081406
    14091407c     common variables and constants
    1410 #include "dimensions.h"
    14111408#include "param.h"
    14121409#include "param_v4.h"
     
    17021699
    17031700!     common variables and constants
    1704       include "dimensions.h"
    17051701      include "param.h"
    17061702      include 'param_v4.h'
  • trunk/LMDZ.VENUS/libf/phyvenus/jthermcalc_e107.F

    r1310 r1530  
    1616
    1717c     common variables and constants
    18 #include "dimensions.h"
    1918#include "param.h"
    2019#include "param_v4.h"
  • trunk/LMDZ.VENUS/libf/phyvenus/load_ksi.F

    r1301 r1530  
    44      IMPLICIT none
    55
    6 #include "dimensions.h"
    76#include "YOMCST.h"
    87#include "comcstVE.h"
  • trunk/LMDZ.VENUS/libf/phyvenus/lw_venus_ve.F

    r1301 r1530  
    99      IMPLICIT none
    1010
    11 #include "dimensions.h"
    1211#include "YOMCST.h"
    1312C
  • trunk/LMDZ.VENUS/libf/phyvenus/lwi.F

    r1017 r1530  
    66
    77
    8 #include "dimensions.h"
    98#include "YOMCST.h"
    109#include "timerad.h"
  • trunk/LMDZ.VENUS/libf/phyvenus/moldiff_red.F90

    r1442 r1530  
    88implicit none
    99
    10 !#include "dimensions.h"
    1110!#include "dimphys.h"
    1211#include "comcstfi.h"
     
    965964        use infotrac
    966965        IMPLICIT NONE
    967         !#include "dimensions.h"
    968966
    969967        INTEGER :: nl,nq,l, nn
     
    10511049        use infotrac
    10521050        IMPLICIT NONE
    1053 !#include "dimensions.h"
    10541051       
    10551052        INTEGER :: nl,nq,il,l,i,iq,nlx,iz,ig
     
    13941391        use infotrac
    13951392        IMPLICIT NONE
    1396 !#include "dimensions.h"
    13971393        INTEGER :: nl,nq,nlx,il,nn,iP,ig,compteur
    13981394        INTEGER,DIMENSION(1) :: indP
     
    14971493!        use infotrac
    14981494        IMPLICIT NONE
    1499 !#include "dimensions.h"
    15001495        INTEGER :: nl,nq,nlx,il,nn,iP,ig,compteur
    15011496        INTEGER,DIMENSION(1) :: indP
  • trunk/LMDZ.VENUS/libf/phyvenus/moldiffcoeff_red.F

    r1442 r1530  
    1616c
    1717c=======================================================================
    18 !#include "dimensions.h"
    19 !#include "dimphys.h"
    20 !#include "callkeys.h"
    21 !#include "comdiurn.h"
    22 !#include "chimiedata.h"
    23 !#include "tracer.h"
    24 !#include "conc.h"
    2518#include "diffusion.h"
    2619
  • trunk/LMDZ.VENUS/libf/phyvenus/molvis.F

    r1310 r1530  
    2020c   declarations:
    2121c-----------------------------------------------------------------------
    22 
    23 !#include "dimensions.h"
    24 !#include "dimphys.h"
    25 !#include "comcstfi.h"
    26 !#include "surfdat.h"
    27 !#include "chimiedata.h"
    28 !#include "conc.h"
    2922
    3023c   arguments:
  • trunk/LMDZ.VENUS/libf/phyvenus/new_cloud_sedim.F

    r1525 r1530  
    1919c   -------------
    2020#include "YOMCST.h"     
    21 c#include "dimensions.h"
    2221c#include "dimphys.h"
    2322c#include "comcstfi.h"
  • trunk/LMDZ.VENUS/libf/phyvenus/nirco2abs.F

    r1442 r1530  
    5454c
    5555
    56 !#include "dimensions.h"
    5756#include "YOMCST.h"
    5857#include "clesphys.h"
  • trunk/LMDZ.VENUS/libf/phyvenus/nlte_tcool.F

    r1442 r1530  
    3333      implicit none
    3434
    35       include "dimensions.h"
    3635      include "nlte_paramdef.h"
    3736      include "nlte_commons.h"
  • trunk/LMDZ.VENUS/libf/phyvenus/nltecool.F

    r1310 r1530  
    3939
    4040#include "nltedata.h" ! (Equivalent to the reading of the "nlte_escape.dat" file)
    41 #include "dimensions.h"
    42 
    43 !#include "dimphys.h"
    44 !#include "chimiedata.h"
    45 !#include "conc.h" !Added to have "dynamic composition" in the scheme
    46 !#include "tracer.h" !"
    47 !#include "callkeys.h"
    4841
    4942c Input and output variables
  • trunk/LMDZ.VENUS/libf/phyvenus/nlthermeq.F

    r1310 r1530  
    1010      use dimphy
    1111      implicit none
    12 #include "dimensions.h"
    1312c#include "dimradmars.h"
    1413#include "nlteparams.h"
  • trunk/LMDZ.VENUS/libf/phyvenus/orodrag.F

    r101 r1530  
    1010      use dimphy
    1111      IMPLICIT NONE
    12 
    13 #include "dimensions.h"
    14 #include "paramet.h"
    1512
    1613c
  • trunk/LMDZ.VENUS/libf/phyvenus/orosetup.F

    r101 r1530  
    9797      use dimphy
    9898      implicit none
    99 
    100 #include "dimensions.h"
    101 #include "paramet.h"
    10299
    103100#include "YOMCST.h"
  • trunk/LMDZ.VENUS/libf/phyvenus/param_read.F

    r1452 r1530  
    77
    88c     common variables and constants
    9       include "dimensions.h"
    109      include 'param.h'
    1110      include 'param_v4.h'
  • trunk/LMDZ.VENUS/libf/phyvenus/param_read_e107.F

    r1310 r1530  
    77
    88c     common variables and constants
    9       include "dimensions.h"
    109      include 'param.h'
    1110      include 'param_v4.h'
  • trunk/LMDZ.VENUS/libf/phyvenus/phyetat0.F90

    r1525 r1530  
    2121! Objet: Lecture de l'etat initial pour la physique
    2222!======================================================================
    23 #include "dimensions.h"
    2423#include "netcdf.inc"
    2524#include "dimsoil.h"
  • trunk/LMDZ.VENUS/libf/phyvenus/phyredem.F90

    r1525 r1530  
    2020! Objet: Lecture de l'etat initial pour la physique
    2121!======================================================================
    22 #include "dimensions.h"
    2322#include "netcdf.inc"
    2423#include "dimsoil.h"
  • trunk/LMDZ.VENUS/libf/phyvenus/physiq.F

    r1525 r1530  
    7070!      use ieee_arithmetic
    7171      use time_phylmdz_mod, only: annee_ref, day_ref, itau_phy
     72      use mod_grid_phy_lmdz, only: nbp_lon
    7273      use logic_mod, only: iflag_trac
    7374      IMPLICIT none
     
    8081#define histins
    8182c======================================================================
    82 #include "dimensions.h"
    83       integer jjmp1
    84       parameter (jjmp1=jjm+1-1/jjm)
    8583#include "dimsoil.h"
    8684#include "clesphys.h"
  • trunk/LMDZ.VENUS/libf/phyvenus/phytrac_chimie.F

    r1525 r1530  
    2828      IMPLICIT none
    2929     
    30 c#include "dimensions.h"
    3130#include "clesphys.h"
    3231c#include "temps.h"
  • trunk/LMDZ.VENUS/libf/phyvenus/phytrac_emiss.F

    r1525 r1530  
    3232      USE comgeomphy
    3333      USE chemparam_mod,only:M_tr
     34      USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat
    3435      IMPLICIT none
    3536#include "YOMCST.h"
    36 #include "dimensions.h"
    3737#include "clesphys.h"
    38 #include "paramet.h"
    3938c======================================================================
    4039
     
    124123         
    125124         ig_volcan(ilat,ilon)= 0
    126          deltalat = 180./jjm
    127          deltalon = 360./iim
     125         if ((nbp_lon*nbp_lat)==1) then ! running a 1D simulation
     126           deltalat=180.
     127           deltalon=360.
     128         else
     129           deltalat = 180./(nbp_lat-1)
     130           deltalon = 360./nbp_lon
     131         endif
    128132
    129133         do i=1,nlon
  • trunk/LMDZ.VENUS/libf/phyvenus/phytrac_relax.F

    r1525 r1530  
    3131      IMPLICIT none
    3232#include "YOMCST.h"
    33 #include "dimensions.h"
    3433#include "clesphys.h"
    35 #include "paramet.h"
    3634c======================================================================
    3735
  • trunk/LMDZ.VENUS/libf/phyvenus/radlwsw.F

    r1442 r1530  
    3232      use write_field_phy
    3333      IMPLICIT none
    34 #include "dimensions.h"
    3534#include "YOMCST.h"
    3635#include "clesphys.h"
  • trunk/LMDZ.VENUS/libf/phyvenus/radlwsw.NewtonCool

    r1310 r1530  
    2929      USE write_field_phy
    3030      IMPLICIT none
    31 #include "dimensions.h"
    3231#include "YOMCST.h"
    3332#include "clesphys.h"
  • trunk/LMDZ.VENUS/libf/phyvenus/soil.F

    r101 r1530  
    4848      use dimphy
    4949      IMPLICIT NONE
    50 #include "dimensions.h"
    5150#include "YOMCST.h"
    5251#include "dimsoil.h"
  • trunk/LMDZ.VENUS/libf/phyvenus/sugwd.F

    r101 r1530  
    7676      use dimphy
    7777      IMPLICIT NONE
    78 
    79 #include "dimensions.h"
    80 #include "paramet.h"
    8178
    8279#include "YOEGWD.h"
  • trunk/LMDZ.VENUS/libf/phyvenus/sw_venus_cl.F

    r1442 r1530  
    88      IMPLICIT none
    99
    10 #include "dimensions.h"
    1110#include "YOMCST.h"
    1211C
  • trunk/LMDZ.VENUS/libf/phyvenus/sw_venus_cl_1Dglobave.F

    r1442 r1530  
    88      IMPLICIT none
    99
    10 #include "dimensions.h"
    1110#include "YOMCST.h"
    1211C
  • trunk/LMDZ.VENUS/libf/phyvenus/sw_venus_dc.F

    r1442 r1530  
    88      IMPLICIT none
    99
    10 #include "dimensions.h"
    1110#include "YOMCST.h"
    1211C
  • trunk/LMDZ.VENUS/libf/phyvenus/sw_venus_dc_1Dglobave.F

    r1442 r1530  
    88      IMPLICIT none
    99
    10 #include "dimensions.h"
    1110#include "YOMCST.h"
    1211C
  • trunk/LMDZ.VENUS/libf/phyvenus/sw_venus_ve.F

    r1442 r1530  
    88      IMPLICIT none
    99
    10 #include "dimensions.h"
    1110#include "YOMCST.h"
    1211C
  • trunk/LMDZ.VENUS/libf/phyvenus/sw_venus_ve_1Dglobave.F

    r1442 r1530  
    88      IMPLICIT none
    99
    10 #include "dimensions.h"
    1110#include "YOMCST.h"
    1211C
  • trunk/LMDZ.VENUS/libf/phyvenus/transp.F

    r101 r1530  
    1313      use dimphy
    1414      IMPLICIT none
    15 #include "dimensions.h"
    1615#include "YOMCST.h"
    1716c
  • trunk/LMDZ.VENUS/libf/phyvenus/ustarhb.F

    r892 r1530  
    1818      use dimphy
    1919      IMPLICIT none
    20 #include "dimensions.h"
    2120#include "YOMCST.h"
    2221c
  • trunk/LMDZ.VENUS/libf/phyvenus/vdif_kcay.F

    r101 r1530  
    88      use dimphy
    99      IMPLICIT NONE
    10 #include "dimensions.h"
    1110c.......................................................................
    1211c
  • trunk/LMDZ.VENUS/libf/phyvenus/yamada.F

    r101 r1530  
    88      use dimphy
    99      IMPLICIT NONE
    10 #include "dimensions.h"
    1110c.......................................................................
    1211c
  • trunk/LMDZ.VENUS/libf/phyvenus/yamada4.F

    r101 r1530  
    99      use dimphy
    1010      IMPLICIT NONE
    11 #include "dimensions.h"
    1211c.......................................................................
    1312c
  • trunk/LMDZ.VENUS/libf/phyvenus/zenang.F

    r1301 r1530  
    2626      use dimphy
    2727      IMPLICIT none
    28 #include "dimensions.h"
    2928#include "YOMCST.h"
    3029c================================================================
Note: See TracChangeset for help on using the changeset viewer.