!copyright (C) 2001  MSC-RPN COMM  %%%RPNPHY%%%
*** S/P SHALLCONV3
#include "phy_macros_f.h"

      subroutine shallconv3 (d, dsiz, f, fsiz, v, vsiz, kount,trnch, 1,1
     1                       cdt1, ni, nk              )
*
#include "impnone.cdk"
*
      INTEGER ni, nk, kount, trnch
      INTEGER dsiz,fsiz,vsiz
      REAL d(dsiz), f(fsiz), v(vsiz)
      REAL cdt1
*
*Author
*          A-M Leduc. (March 2002)
*
*Revisions
* 001      S-Belair, A-M.Leduc (nov 2002) add convective counter 
*                                      kshal for ktrsnt.
* 002      A-M.Leduc (nov 2002)  - add switch ishlcvt(2) and remove
*                             delt as argument--->shallconv2
*Object
*          This is the interface routine for shallow convection schemes
*
*Arguments
*
*          - Input -
* D        "dynamic" bus (dynamics input fields)
*
*          - Input/Output -
* F        "permanent" bus (field of permanent physics variables)
* V        "volatile" bus 
*
*          - Input -
* DSIZ     dimension of d
* FSIZ     dimension of f
* VSIZ     dimension of v
* TRNCH    row number
* KOUNT    timestep number
* CDT1     timestep * factdt
*          see common block "options"
* NI       horizontal running length
* NK       vertical dimension
*
**
*
#include "indx_sfc.cdk"
#include "phybus.cdk"
#include "options.cdk"

      EXTERNAL ktrsnt3
*
      INTEGER i,k,ik
      INTEGER zilab(ni,nk)
      REAL    zdbdt(ni)
*
**
      IF(ISHLCVT(2).EQ.1)THEN


        DO k=1,nk
        DO i=1,ni
          zilab(i,k) = 0
        END DO
        END DO
*
        DO i=1,ni
         zdbdt(i) = 0.
        END DO
*
        CALL ktrsnt3( v(tshal), v(hushal), zilab, v(ckt), v(qktl),
     1               v(qkti), zdbdt,f(tlcs),f(tscs),v(qcz),
     1               d(tplus), d(tmoins), d(huplus), d(humoins),
     1               d(gzmoins6), v(qdifv), d(pplus), d(pmoins), 
     1               d(sigm), cdt1, v(kshal),ni, nk-1  )
*
*       TRANSFORMATION EN HAUTEUR D'EAU (M) - DIVISER PAR DENSITE EAU
*
        DO i=1,ni
           f(tscs+i-1) = f(tscs+i-1) * 1.e-03
           f(tlcs+i-1) = f(tlcs+i-1) * 1.e-03
        END DO
      ENDIF
*
*
      RETURN
      END