Changeset 3296 for LMDZ6


Ignore:
Timestamp:
Mar 23, 2018, 10:13:38 PM (7 years ago)
Author:
oboucher
Message:

a USE statement was missing in mo_simple_plumes.F90
z array has been changed to reflect Stephanie's recommendation
that z is the altitude of the middle of the layers

Location:
LMDZ6/trunk/libf/phylmd
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/phylmd/macv2sp.F90

    r3287 r3296  
    3636  REAL,DIMENSION(klon,klev) :: aod_prof, ssa_prof, asy_prof
    3737  REAL,DIMENSION(klon,klev) :: z, dz
    38   REAL,DIMENSION(klon)      :: oro, zrho
     38  REAL,DIMENSION(klon)      :: oro, zrho, zt
    3939  !
    4040  INTEGER, PARAMETER :: nmon = 12
     
    6363  ! define the height and dheight arrays
    6464  !
    65   z(:,1)  = 0.                                      ! altitude of surface taken as 0
    66   DO k = 1, klev-1
    67     oro(:) = pphis(:)/RG                            ! surface height in m
    68     zrho(:) = pplay(:, k)/t_seri(:, k)/RD           ! air density in kg/m3
    69     dz(:,k) = (paprs(:,k)-paprs(:,k+1))/zrho(:)/RG  ! layer thickness in m
    70     z(:,k+1) = z(:,k) + dz(:,k)                     ! height of interfaces in m, starting from 0 at surface
     65  oro(:)  = pphis(:)/RG                             ! surface height in m
     66  !
     67  DO k = 1, klev
     68    zrho(:) = pplay(:,k)/t_seri(:,k)/RD                         ! air density in kg/m3
     69    dz(:,k) = (paprs(:,k)-paprs(:,k+1))/zrho(:)/RG              ! layer thickness in m
     70    IF (k==1) THEN
     71       z(:,1) = oro(:) + (paprs(:,1)-pplay(:,1))/zrho(:)/RG     ! altitude middle of first layer in m
     72       zt(:)  = oro(:) + dz(:,1)                                ! altitude top of first layer in m
     73    ELSE
     74      z(:,k) = zt(:) + (paprs(:,k)-pplay(:,k))/zrho(:)/RG       ! altitude middle of layer k in m
     75      zt(:)  = zt(:) + dz(:,k)                                  ! altitude top of layer k in m
     76    ENDIF
    7177  ENDDO
    7278  !
  • LMDZ6/trunk/libf/phylmd/mo_simple_plumes.F90

    r3295 r3296  
    7373    USE mod_phys_lmdz_mpi_data, ONLY: is_mpi_root
    7474    USE mod_phys_lmdz_omp_data, ONLY: is_omp_root
     75    USE mod_phys_lmdz_para
    7576    !
    7677    ! ----------
Note: See TracChangeset for help on using the changeset viewer.