!copyright (C) 2001 MSC-RPN COMM %%%RPNPHY%%% ***S/P LIN_INICHAMP1_TL *SUBROUTINE LIN_INICHAMP1_TL(E, ESIZ, F, FSIZ, 1 $ V, VSIZ, D, DSIZ, $ TRAV2D, SELOC, KOUNT, TRNCH, $ DT, CDT1, NI, NK) * #include "impnone.cdk"
* INTEGER I, IK, K, ESIZ, FSIZ, VSIZ, DSIZ, KOUNT, NI, NK, TRNCH REAL E(ESIZ), F(FSIZ), D(DSIZ), V(VSIZ), DT, CDT1 REAL SELOC(NI,NK) REAL TRAV2D (NI,NK) * *Author * S. Laroche (Dec 2002) Initialization for TLM * *Revision * * *Object * To initialize arrays. * * Arguments * * - Input - * F field for permanent physics variables * FSIZ dimension of F * V volatile bus * VSIZ dimension of V * D dynamics bus * DSIZ dimension of D * TRAV2D work field * SELOC intermediate (staggered) sigma levels (2D) * KOUNT timestep number * TRNCH row number * DT length of timestep * CDT1 = DT for 2-time level models * = 2*DT for 3-time level models * NI horizontal dimension * NK vertical dimension * ** * #include "indx_sfc.cdk"
#include "options.cdk"
#include "phy_macros_f.h"
#include "phybus.cdk"
#include "consphy.cdk"
INTEGER NIK * EXTERNAL INISURF * NIK = NI*NK * * ************************************************************************ * INITIALISATIONS FAITES A KOUNT = 1 SEULEMENT * * -------------------------------------------- * ************************************************************************ * IF (KOUNT.EQ.1) THEN * DO I=0,NI-1 F(ASC + I) = 0.0 F(ALC + I) = 0.0 F(ASS + I) = 0.0 F(ALS + I) = 0.0 END DO * DO I=1,VSIZ V(I) = 0.0 END DO * ENDIF * * ************************************************************************ * INITIALISATIONS FAITES A KOUNT > 1 SEULEMENT * * -------------------------------------------- * ************************************************************************ * IF (KOUNT.GT.1) THEN * DO I=1,VSIZ V(I) = 0.0 END DO * ENDIF * RETURN END