Changeset 147 for LMDZ.3.3/branches


Ignore:
Timestamp:
Oct 24, 2000, 9:55:53 AM (24 years ago)
Author:
lmdzadmin
Message:

Dimensionnement de tous les tableaux a klon plutot qu'a knon pour cause
de compilateur NEC deficient (?)
LF

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ.3.3/branches/rel-LF/libf/phylmd/interface_surf.F90

    r143 r147  
    136136  integer, intent(IN) :: nisurf
    137137  integer, intent(IN) :: knon
    138   integer, dimension(knon), intent(in) :: knindex
     138  integer, dimension(klon), intent(in) :: knindex
    139139  real, dimension(klon,nbsrf), intent(IN) :: pctsrf
    140140  logical, intent(IN) :: debut, lafin, ok_veget
    141141  real, dimension(klon), intent(IN) :: rlon, rlat
    142   real, dimension(knon), intent(IN) :: zlev
    143   real, dimension(knon), intent(IN) :: u1_lay, v1_lay
    144   real, dimension(knon), intent(IN) :: temp_air, spechum
    145   real, dimension(knon), intent(IN) :: hum_air, ccanopy
    146   real, dimension(knon), intent(IN) :: tq_cdrag, petAcoef, peqAcoef
    147   real, dimension(knon), intent(IN) :: petBcoef, peqBcoef
    148   real, dimension(knon), intent(IN) :: precip_rain, precip_snow
    149   real, dimension(knon), intent(IN) :: lwdown, swnet, swdown, ps, albedo
    150   real, dimension(knon), intent(IN) :: tsurf, p1lay
    151   real, dimension(knon), intent(IN) :: radsol
     142  real, dimension(klon), intent(IN) :: zlev
     143  real, dimension(klon), intent(IN) :: u1_lay, v1_lay
     144  real, dimension(klon), intent(IN) :: temp_air, spechum
     145  real, dimension(klon), intent(IN) :: hum_air, ccanopy
     146  real, dimension(klon), intent(IN) :: tq_cdrag, petAcoef, peqAcoef
     147  real, dimension(klon), intent(IN) :: petBcoef, peqBcoef
     148  real, dimension(klon), intent(IN) :: precip_rain, precip_snow
     149  real, dimension(klon), intent(IN) :: lwdown, swnet, swdown, ps, albedo
     150  real, dimension(klon), intent(IN) :: tsurf, p1lay
     151  real, dimension(klon), intent(IN) :: radsol
    152152  real, dimension(klon), intent(IN) :: zmasq
    153153  real, dimension(klon), intent(IN) :: fder, taux, tauy
    154154  character (len = 6)  :: ocean
    155155  integer              :: npas, nexca ! nombre et pas de temps couplage
    156   real, dimension(knon), intent(INOUT) :: evap, snow, qsol
     156  real, dimension(klon), intent(INOUT) :: evap, snow, qsol
    157157
    158158! Parametres de sortie
    159   real, dimension(knon), intent(OUT):: fluxsens, fluxlat
    160   real, dimension(knon), intent(OUT):: tsol_rad, tsurf_new, alb_new
    161   real, dimension(knon), intent(OUT):: emis_new, z0_new
    162   real, dimension(knon), intent(OUT):: dflux_l, dflux_s
     159  real, dimension(klon), intent(OUT):: fluxsens, fluxlat
     160  real, dimension(klon), intent(OUT):: tsol_rad, tsurf_new, alb_new
     161  real, dimension(klon), intent(OUT):: emis_new, z0_new
     162  real, dimension(klon), intent(OUT):: dflux_l, dflux_s
    163163  real, dimension(klon,nbsrf), intent(OUT) :: pctsrf_new
    164164  real, dimension(klon), intent(INOUT):: agesno
    … …  
    170170  INTEGER              :: error, ii
    171171  logical              :: check = .true.
    172   real, dimension(knon):: cal, beta, dif_grnd, capsol
     172  real, dimension(klon):: cal, beta, dif_grnd, capsol
    173173  real, parameter      :: calice=1.0/(5.1444e+06*0.15), tau_gl=1./86400.*5.
    174174  real, parameter      :: calsno=1./(2.3867e+06*.15)
    175   real, dimension(knon):: alb_ice
    176   real, dimension(knon):: tsurf_temp
     175  real, dimension(klon):: alb_ice
     176  real, dimension(klon):: tsurf_temp
    177177  real, dimension(klon):: alb_neig_grid, alb_eau
    178   real, dimension(knon):: alb_neig
    179   REAL, DIMENSION(knon):: lmt_rug, lmt_alb
    180   real, DIMENSION(knon):: zfra
     178  real, dimension(klon):: alb_neig
     179  real, DIMENSION(klon):: zfra
    181180
    182181  if (check) write(*,*) 'Entree ', modname
    … …  
    247246      call calbeta(dtime, nisurf, knon, snow, qsol, beta, capsol, dif_grnd)
    248247      cal = RCPD * capsol
    249       call calcul_fluxs( knon, nisurf, dtime, &
     248      call calcul_fluxs( klon, knon, nisurf, dtime, &
    250249     &   tsurf, p1lay, cal, beta, tq_cdrag, ps, &
    251250     &   precip_rain, precip_snow, snow, qsol,  &
    … …  
    261260     & alb_new, z0_new)
    262261!
    263 ! Pb compilo sun
    264 !       alb_neig = alb_neig_grid(knindex)
    265 !      alb_new = alb_neig*zfra + lmt_alb(knindex)*(1.0-zfra)
    266 !      z0_new = lmt_rug(knindex)
    267262!
    268263       DO ii = 1, knon
    … …  
    327322    dif_grnd = 0.
    328323
    329     call calcul_fluxs( knon, nisurf, dtime, &
     324    call calcul_fluxs( klon, knon, nisurf, dtime, &
    330325     &   tsurf_temp, p1lay, cal, beta, tq_cdrag, ps, &
    331326     &   precip_rain, precip_snow, snow, qsol,  &
    … …  
    383378    where (snow > 0.0) cal = calsno
    384379
    385     call calcul_fluxs( knon, nisurf, dtime, &
     380    call calcul_fluxs( klon, knon, nisurf, dtime, &
    386381     &   tsurf_temp, p1lay, cal, beta, tq_cdrag, ps, &
    387382     &   precip_rain, precip_snow, snow, qsol,  &
    … …  
    413408    dif_grnd = 0.0
    414409
    415     call calcul_fluxs( knon, nisurf, dtime, &
     410    call calcul_fluxs( klon, knon, nisurf, dtime, &
    416411     &   tsurf, p1lay, cal, beta, tq_cdrag, ps, &
    417412     &   precip_rain, precip_snow, snow, qsol,  &
    … …  
    540535  integer, intent(IN) :: nisurf
    541536  integer, intent(IN) :: knon
    542   integer, dimension(knon), intent(IN) :: knindex
     537  integer, dimension(klon), intent(IN) :: knindex
    543538  logical, intent(IN) :: debut, lafin, ok_veget
    544539  real, dimension(klon), intent(IN) :: rlon, rlat
    545   real, dimension(knon), intent(IN) :: zlev
    546   real, dimension(knon), intent(IN) :: u1_lay, v1_lay
    547   real, dimension(knon), intent(IN) :: temp_air, spechum
    548   real, dimension(knon), intent(IN) :: hum_air, ccanopy
    549   real, dimension(knon), intent(IN) :: tq_cdrag, petAcoef, peqAcoef
    550   real, dimension(knon), intent(IN) :: petBcoef, peqBcoef
    551   real, dimension(knon), intent(IN) :: precip_rain, precip_snow
    552   real, dimension(knon), intent(IN) :: lwdown, swnet, swdown, ps
    553   real, dimension(knon), intent(IN) :: tsurf, p1lay
    554   real, dimension(knon), intent(IN) :: radsol
     540  real, dimension(klon), intent(IN) :: zlev
     541  real, dimension(klon), intent(IN) :: u1_lay, v1_lay
     542  real, dimension(klon), intent(IN) :: temp_air, spechum
     543  real, dimension(klon), intent(IN) :: hum_air, ccanopy
     544  real, dimension(klon), intent(IN) :: tq_cdrag, petAcoef, peqAcoef
     545  real, dimension(klon), intent(IN) :: petBcoef, peqBcoef
     546  real, dimension(klon), intent(IN) :: precip_rain, precip_snow
     547  real, dimension(klon), intent(IN) :: lwdown, swnet, swdown, ps
     548  real, dimension(klon), intent(IN) :: tsurf, p1lay
     549  real, dimension(klon), intent(IN) :: radsol
    555550! Parametres de sortie
    556   real, dimension(knon), intent(OUT):: evap, fluxsens, fluxlat
    557   real, dimension(knon), intent(OUT):: tsol_rad, tsurf_new, alb_new
    558   real, dimension(knon), intent(OUT):: emis_new, z0_new
    559   real, dimension(knon), intent(OUT):: dflux_s, dflux_l
     551  real, dimension(klon), intent(OUT):: evap, fluxsens, fluxlat
     552  real, dimension(klon), intent(OUT):: tsol_rad, tsurf_new, alb_new
     553  real, dimension(klon), intent(OUT):: emis_new, z0_new
     554  real, dimension(klon), intent(OUT):: dflux_s, dflux_l
    560555
    561556! Local
    … …  
    566561  character (len = 80) :: abort_message
    567562  logical              :: check = .true.
    568   real, dimension(knon) :: cal, beta, dif_grnd, capsol
     563  real, dimension(klon) :: cal, beta, dif_grnd, capsol
    569564! type de couplage dans sechiba
    570565!  character (len=10)   :: coupling = 'implicit'
    … …  
    581576  integer, save          :: rest_id_stom, hist_id_stom
    582577
    583   real, dimension(knon):: snow, qsol
     578  real, dimension(klon):: snow, qsol
    584579
    585580  if (check) write(*,*)'Entree ', modname
    … …  
    721716  integer, intent(IN) :: knon
    722717  real, dimension(klon,nbsrf), intent(IN) :: pctsrf
    723   integer, dimension(knon), intent(in) :: knindex
     718  integer, dimension(klon), intent(in) :: knindex
    724719  logical, intent(IN) :: debut, lafin
    725720  real, dimension(klon), intent(IN) :: rlon, rlat
    726721  character (len = 6)  :: ocean
    727   real, dimension(knon), intent(IN) :: lwdown, swdown
    728   real, dimension(knon), intent(IN) :: precip_rain, precip_snow
    729   real, dimension(knon), intent(IN) :: tsurf, fder, albsol, taux, tauy
     722  real, dimension(klon), intent(IN) :: lwdown, swdown
     723  real, dimension(klon), intent(IN) :: precip_rain, precip_snow
     724  real, dimension(klon), intent(IN) :: tsurf, fder, albsol, taux, tauy
    730725  INTEGER              :: nexca, npas, kstep
    731726  real, dimension(klon), intent(IN) :: zmasq
    732727
    733   real, dimension(knon), intent(INOUT) :: evap
     728  real, dimension(klon), intent(INOUT) :: evap
    734729
    735730! Parametres de sortie
    736   real, dimension(knon), intent(OUT):: tsurf_new, alb_new, alb_ice
     731  real, dimension(klon), intent(OUT):: tsurf_new, alb_new, alb_ice
    737732  real, dimension(klon,nbsrf), intent(OUT) :: pctsrf_new
    738733
    … …  
    803798    enddo
    804799    sum_error = 0
    805     allocate(cpl_sols(knon,2), stat = error); sum_error = sum_error + error
    806     allocate(cpl_nsol(knon,2), stat = error); sum_error = sum_error + error
    807     allocate(cpl_rain(knon,2), stat = error); sum_error = sum_error + error
    808     allocate(cpl_snow(knon,2), stat = error); sum_error = sum_error + error
    809     allocate(cpl_evap(knon,2), stat = error); sum_error = sum_error + error
    810     allocate(cpl_tsol(knon,2), stat = error); sum_error = sum_error + error
    811     allocate(cpl_fder(knon,2), stat = error); sum_error = sum_error + error
    812     allocate(cpl_albe(knon,2), stat = error); sum_error = sum_error + error
    813     allocate(cpl_taux(knon,2), stat = error); sum_error = sum_error + error
    814     allocate(cpl_tauy(knon,2), stat = error); sum_error = sum_error + error
    815     allocate(cpl_rcoa(knon,2), stat = error); sum_error = sum_error + error
    816     allocate(cpl_rriv(knon,2), stat = error); sum_error = sum_error + error
     800    allocate(cpl_sols(klon,2), stat = error); sum_error = sum_error + error
     801    allocate(cpl_nsol(klon,2), stat = error); sum_error = sum_error + error
     802    allocate(cpl_rain(klon,2), stat = error); sum_error = sum_error + error
     803    allocate(cpl_snow(klon,2), stat = error); sum_error = sum_error + error
     804    allocate(cpl_evap(klon,2), stat = error); sum_error = sum_error + error
     805    allocate(cpl_tsol(klon,2), stat = error); sum_error = sum_error + error
     806    allocate(cpl_fder(klon,2), stat = error); sum_error = sum_error + error
     807    allocate(cpl_albe(klon,2), stat = error); sum_error = sum_error + error
     808    allocate(cpl_taux(klon,2), stat = error); sum_error = sum_error + error
     809    allocate(cpl_tauy(klon,2), stat = error); sum_error = sum_error + error
     810    allocate(cpl_rcoa(klon,2), stat = error); sum_error = sum_error + error
     811    allocate(cpl_rriv(klon,2), stat = error); sum_error = sum_error + error
    817812    allocate(read_sst(iim, jjm+1), stat = error); sum_error = sum_error + error
    818813    allocate(read_sic(iim, jjm+1), stat = error); sum_error = sum_error + error
    … …  
    884879!
    885880    if (nisurf == is_oce) then
     881      if (check) write(*,*)'rentree fromcpl, itime-1 = ',itime-1
    886882      call fromcpl(itime-1,(jjm+1)*iim,                                  &
    887883     &        read_sst, read_sic, read_sit, read_alb_sic)
    … …  
    940936!
    941937! allocation memoire
    942     sum_error = 0
    943     allocate(tmp_sols(iim,jjm+1,2), stat=error); sum_error = sum_error + error
    944     allocate(tmp_nsol(iim,jjm+1,2), stat=error); sum_error = sum_error + error
    945     allocate(tmp_rain(iim,jjm+1,2), stat=error); sum_error = sum_error + error
    946     allocate(tmp_snow(iim,jjm+1,2), stat=error); sum_error = sum_error + error
    947     allocate(tmp_evap(iim,jjm+1,2), stat=error); sum_error = sum_error + error
    948     allocate(tmp_tsol(iim,jjm+1,2), stat=error); sum_error = sum_error + error
    949     allocate(tmp_fder(iim,jjm+1,2), stat=error); sum_error = sum_error + error
    950     allocate(tmp_albe(iim,jjm+1,2), stat=error); sum_error = sum_error + error
    951     allocate(tmp_taux(iim,jjm+1,2), stat=error); sum_error = sum_error + error
    952     allocate(tmp_tauy(iim,jjm+1,2), stat=error); sum_error = sum_error + error
    953     allocate(tmp_rriv(iim,jjm+1,2), stat=error); sum_error = sum_error + error
    954     allocate(tmp_rcoa(iim,jjm+1,2), stat=error); sum_error = sum_error + error
    955     if (sum_error /= 0) then
    956       abort_message='Pb allocation variables couplees'
    957       call abort_gcm(modname,abort_message,1)
     938    if (nisurf == is_oce) then
     939      sum_error = 0
     940      allocate(tmp_sols(iim,jjm+1,2), stat=error); sum_error = sum_error + error
     941      allocate(tmp_nsol(iim,jjm+1,2), stat=error); sum_error = sum_error + error
     942      allocate(tmp_rain(iim,jjm+1,2), stat=error); sum_error = sum_error + error
     943      allocate(tmp_snow(iim,jjm+1,2), stat=error); sum_error = sum_error + error
     944      allocate(tmp_evap(iim,jjm+1,2), stat=error); sum_error = sum_error + error
     945      allocate(tmp_tsol(iim,jjm+1,2), stat=error); sum_error = sum_error + error
     946      allocate(tmp_fder(iim,jjm+1,2), stat=error); sum_error = sum_error + error
     947      allocate(tmp_albe(iim,jjm+1,2), stat=error); sum_error = sum_error + error
     948      allocate(tmp_taux(iim,jjm+1,2), stat=error); sum_error = sum_error + error
     949      allocate(tmp_tauy(iim,jjm+1,2), stat=error); sum_error = sum_error + error
     950      allocate(tmp_rriv(iim,jjm+1,2), stat=error); sum_error = sum_error + error
     951      allocate(tmp_rcoa(iim,jjm+1,2), stat=error); sum_error = sum_error + error
     952      if (sum_error /= 0) then
     953        abort_message='Pb allocation variables couplees pour l''ecriture'
     954        call abort_gcm(modname,abort_message,1)
     955      endif
    958956    endif
    959957
    … …  
    961959! Mise sur la bonne grille des champs a passer au coupleur
    962960!
    963     call gath2cpl(cpl_sols(1,nisurf), tmp_sols(1,1,nisurf), klon, knon,iim,jjm,                  knindex)
    964     call gath2cpl(cpl_nsol(1,nisurf), tmp_nsol(1,1,nisurf), klon, knon,iim,jjm,                  knindex)
    965     call gath2cpl(cpl_rain(1,nisurf), tmp_rain(1,1,nisurf), klon, knon,iim,jjm,                  knindex)
    966     call gath2cpl(cpl_snow(1,nisurf), tmp_snow(1,1,nisurf), klon, knon,iim,jjm,                  knindex)
    967     call gath2cpl(cpl_evap(1,nisurf), tmp_evap(1,1,nisurf), klon, knon,iim,jjm,                  knindex)
    968     call gath2cpl(cpl_tsol(1,nisurf), tmp_tsol(1,1,nisurf), klon, knon,iim,jjm,                  knindex)
    969     call gath2cpl(cpl_fder(1,nisurf), tmp_fder(1,1,nisurf), klon, knon,iim,jjm,                  knindex)
    970     call gath2cpl(cpl_albe(1,nisurf), tmp_albe(1,1,nisurf), klon, knon,iim,jjm,                  knindex)
    971     call gath2cpl(cpl_taux(1,nisurf), tmp_taux(1,1,nisurf), klon, knon,iim,jjm,                  knindex)
    972     call gath2cpl(cpl_tauy(1,nisurf), tmp_tauy(1,1,nisurf), klon, knon,iim,jjm,                  knindex)
    973     call gath2cpl(cpl_rriv(1,nisurf), tmp_rriv(1,1,nisurf), klon, knon,iim,jjm,                  knindex)
    974     call gath2cpl(cpl_rcoa(1,nisurf), tmp_rcoa(1,1,nisurf), klon, knon,iim,jjm,                  knindex)
     961    cpl_index = 1
     962    if (nisurf == is_sic) cpl_index = 2
     963    call gath2cpl(cpl_sols(1,cpl_index), tmp_sols(1,1,cpl_index), klon, knon,iim,jjm,                  knindex)
     964    call gath2cpl(cpl_nsol(1,cpl_index), tmp_nsol(1,1,cpl_index), klon, knon,iim,jjm,                  knindex)
     965    call gath2cpl(cpl_rain(1,cpl_index), tmp_rain(1,1,cpl_index), klon, knon,iim,jjm,                  knindex)
     966    call gath2cpl(cpl_snow(1,cpl_index), tmp_snow(1,1,cpl_index), klon, knon,iim,jjm,                  knindex)
     967    call gath2cpl(cpl_evap(1,cpl_index), tmp_evap(1,1,cpl_index), klon, knon,iim,jjm,                  knindex)
     968    call gath2cpl(cpl_tsol(1,cpl_index), tmp_tsol(1,1,cpl_index), klon, knon,iim,jjm,                  knindex)
     969    call gath2cpl(cpl_fder(1,cpl_index), tmp_fder(1,1,cpl_index), klon, knon,iim,jjm,                  knindex)
     970    call gath2cpl(cpl_albe(1,cpl_index), tmp_albe(1,1,cpl_index), klon, knon,iim,jjm,                  knindex)
     971    call gath2cpl(cpl_taux(1,cpl_index), tmp_taux(1,1,cpl_index), klon, knon,iim,jjm,                  knindex)
     972    call gath2cpl(cpl_tauy(1,cpl_index), tmp_tauy(1,1,cpl_index), klon, knon,iim,jjm,                  knindex)
     973    call gath2cpl(cpl_rriv(1,cpl_index), tmp_rriv(1,1,cpl_index), klon, knon,iim,jjm,                  knindex)
     974    call gath2cpl(cpl_rcoa(1,cpl_index), tmp_rcoa(1,1,cpl_index), klon, knon,iim,jjm,                  knindex)
    975975
    976976!
    … …  
    11081108  integer, intent(IN) :: knon
    11091109  integer, intent(IN) :: klon
    1110   integer, dimension(knon), intent(in) :: knindex
     1110  integer, dimension(klon), intent(in) :: knindex
    11111111  logical, intent(IN) :: debut
    11121112
    11131113! Parametres de sortie
    1114   real, intent(out), dimension(knon) :: lmt_sst
     1114  real, intent(out), dimension(klon) :: lmt_sst
    11151115  real, intent(out), dimension(klon,nbsrf) :: pctsrf_new
    11161116
    … …  
    11291129  logical     :: check = .true.
    11301130! Champs lus dans le fichier de CL
    1131   real, allocatable , save, dimension(:) :: sst_lu, alb_lu, rug_lu, nat_lu
     1131  real, allocatable , save, dimension(:) :: sst_lu, rug_lu, nat_lu
    11321132  real, allocatable , save, dimension(:,:) :: pct_tmp
    11331133!
    … …  
    13511351  integer, intent(IN) :: knon
    13521352  integer, intent(IN) :: klon
    1353   integer, dimension(knon), intent(in) :: knindex
     1353  integer, dimension(klon), intent(in) :: knindex
    13541354  logical, intent(IN) :: debut
    13551355
    13561356! Parametres de sortie
    1357   real, intent(out), dimension(knon) :: lmt_alb
    1358   real, intent(out), dimension(knon) :: lmt_rug
     1357  real, intent(out), dimension(klon) :: lmt_alb
     1358  real, intent(out), dimension(klon) :: lmt_rug
    13591359
    13601360! Variables locales
    … …  
    14701470!
    14711471
    1472   SUBROUTINE calcul_fluxs( knon, nisurf, dtime, &
     1472  SUBROUTINE calcul_fluxs( klon, knon, nisurf, dtime, &
    14731473     & tsurf, p1lay, cal, beta, coef1lay, ps, &
    14741474     & precip_rain, precip_snow, snow, qsol, &
    … …  
    15151515
    15161516! Parametres d'entree
    1517   integer, intent(IN) :: knon, nisurf
     1517  integer, intent(IN) :: knon, nisurf, klon
    15181518  real   , intent(IN) :: dtime
    1519   real, dimension(knon), intent(IN) :: petAcoef, peqAcoef
    1520   real, dimension(knon), intent(IN) :: petBcoef, peqBcoef
    1521   real, dimension(knon), intent(IN) :: ps, q1lay
    1522   real, dimension(knon), intent(IN) :: tsurf, p1lay, cal, beta, coef1lay
    1523   real, dimension(knon), intent(IN) :: precip_rain, precip_snow
    1524   real, dimension(knon), intent(IN) :: radsol, dif_grnd
    1525   real, dimension(knon), intent(IN) :: t1lay, u1lay, v1lay
    1526   real, dimension(knon), intent(INOUT) :: snow, qsol
     1519  real, dimension(klon), intent(IN) :: petAcoef, peqAcoef
     1520  real, dimension(klon), intent(IN) :: petBcoef, peqBcoef
     1521  real, dimension(klon), intent(IN) :: ps, q1lay
     1522  real, dimension(klon), intent(IN) :: tsurf, p1lay, cal, beta, coef1lay
     1523  real, dimension(klon), intent(IN) :: precip_rain, precip_snow
     1524  real, dimension(klon), intent(IN) :: radsol, dif_grnd
     1525  real, dimension(klon), intent(IN) :: t1lay, u1lay, v1lay
     1526  real, dimension(klon), intent(INOUT) :: snow, qsol
    15271527
    15281528! Parametres sorties
    1529   real, dimension(knon), intent(OUT):: tsurf_new, evap, fluxsens, fluxlat
    1530   real, dimension(knon), intent(OUT):: dflux_s, dflux_l
     1529  real, dimension(klon), intent(OUT):: tsurf_new, evap, fluxsens, fluxlat
     1530  real, dimension(klon), intent(OUT):: dflux_s, dflux_l
    15311531
    15321532! Variables locales
    15331533  integer :: i
    1534   real, dimension(knon) :: zx_mh, zx_nh, zx_oh
    1535   real, dimension(knon) :: zx_mq, zx_nq, zx_oq
    1536   real, dimension(knon) :: zx_pkh, zx_dq_s_dt, zx_qsat, zx_coef
    1537   real, dimension(knon) :: zx_sl, zx_k1,  zx_dq,  zx_cq,  zx_dh, zx_ch
    1538   real, dimension(knon) :: zx_h_ts, zx_q_0 , d_ts
     1534  real, dimension(klon) :: zx_mh, zx_nh, zx_oh
     1535  real, dimension(klon) :: zx_mq, zx_nq, zx_oq
     1536  real, dimension(klon) :: zx_pkh, zx_dq_s_dt, zx_qsat, zx_coef
     1537  real, dimension(klon) :: zx_sl, zx_k1,  zx_dq,  zx_cq,  zx_dh, zx_ch
     1538  real, dimension(klon) :: zx_h_ts, zx_q_0 , d_ts
    15391539  real                  :: zdelta, zcvm5, zx_qs, zcor, zx_dq_s_dh
    15401540  real                  :: bilan_f, fq_fonte
    … …  
    16951695! input
    16961696  integer                   :: klon, knon, iim, jjm
    1697   real, dimension(knon)     :: champ_in
    1698   integer, dimension(knon)  :: knindex
     1697  real, dimension(klon)     :: champ_in
     1698  integer, dimension(klon)  :: knindex
    16991699! output
    17001700  real, dimension(iim,jjm+1)  :: champ_out
    … …  
    17421742  integer                   :: klon, knon, iim, jjm
    17431743  real, dimension(iim,jjm+1)     :: champ_in
    1744   integer, dimension(knon)  :: knindex
     1744  integer, dimension(klon)  :: knindex
    17451745! output
    1746   real, dimension(knon)  :: champ_out
     1746  real, dimension(klon)  :: champ_out
    17471747! local
    17481748  integer                   :: i, ig, j
Note: See TracChangeset for help on using the changeset viewer.