bedload_main.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\sisyphe\bedload_main.f
00002 !
00074                      SUBROUTINE BEDLOAD_MAIN
00075 !                    ***********************
00076 !
00077      &(ACLADM,KSP,KSR, V2DPAR,UNSV2D,CF,EBOR,FW,HN,LIQBOR,
00078      & MASK, MASKEL, MASKPT, Q, QBOR, U2D,
00079      & V2D, S,UNLADM,UW,THETAW,MU,TOB,TOBW,TW,ZF,
00080      & DEBUG, HIDFAC, ICF, IELMT, ISOUS, KDDL, KDIR,
00081      & KENT, KINC, KLOG, KNEU, KSORT, LOADMETH, LT,
00082      & NPOIN, NPTFR, NSICLA, OPTBAN, LS0, BETA, FD90, FDM,
00083      & GRAV, HIDI, HMIN, VCE, CSF_SABLE, XMVE, XMVS, XWC,
00084      & PI, KARMAN, ZERO, KARIM_HOLLY_YANG,MSK, SUSP, VF,
00085      & ENTET, CONST_ALAYER, LCONDIS, LGRAFED, MESH,
00086      & ELAY, LIEBOR, LIMTEC, MASKTR,
00087      & IT1, T1, T2, T3, T4, T5, T6, T7, T8, T9, T10, T11,
00088      & T12,T13,UNORM,AC, AT0, DTS, ELAY0, FRACSED_GF,
00089      & AVAIL, BREACH, CALFA, COEFPN,
00090      & DZF_GF, HIDING, QSCL_C, QSCL_S, QS_C,
00091      & QSCLXC, QSXC, QSCLYC, QSYC, SALFA, ZF_C, ZFCL_C, NSOUS,
00092      & ENTETS, SECCURRENT, SLOPEFF,
00093      & PHISED, DEVIA, BETA2, BIJK,SEDCO,HOULE,
00094      & U3D,V3D,CODE,FLBCLA,MAXADV)
00095 !
00096 !***********************************************************************
00097 ! SISYPHE   V6P1                                   21/07/2011
00098 !***********************************************************************
00099 !
00100 !
00101 !
00102 !
00103 !
00104 !
00105 !
00106 !
00107 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00108 !| AC             |<->| CRITICAL SHIELDS PARAMETER
00109 !| ACLADM         |-->| MEAN DIAMETER OF SEDIMENT
00110 !| AT0            |<->| TIME IN S
00111 !| AVAIL          |<->| VOLUME PERCENT OF EACH CLASS
00112 !| BETA           |-->| COEFFICIENT FOR SLOPING BED EFFECT ( KOCH AND FLOKSTRA)
00113 !| BETA2          |-->| COEFFICIENT FOR THE DEVIATION  (TALMON ET AL.)
00114 !| BIJK           |-->| COEFFICIENT OF THE BIJKER FORMULA
00115 !| BREACH         |<->| INDICATOR FOR NON ERODIBLE BED (FINITE VOLUMES SHEMES)
00116 !| CALFA          |<->| COSINUS OF THE ANGLE BETWEEN MEAN FLOW AND TRANSPORT
00117 !| CF             |-->| QUADRATIC FRICTION COEFFICIENT
00118 !| CODE           |-->| HYDRODYNAMIC CODE IN CASE OF COUPLING
00119 !| COEFPN         |<->| CORRECTION OF TRANSORT FOR SLOPING BED EFFECT
00120 !| CONST_ALAYER   |-->| CONSTANT ACTIVE LAYER THICKNESS OR NOT
00121 !| DEBUG          |-->| FLAG FOR DEBUGGING
00122 !| DEVIA          |-->| SLOPE EFFECT FORMULA FOR DEVIATION
00123 !| DTS            |<->| TIME STEP FOR SUSPENSION
00124 !| DZF_GF         |<->| (A SUPPRIMER)
00125 !| EBOR           |<->| IMPOSED BOUNDARY CONDITION FOR BED EVOLUTION (DIRICHLET)
00126 !| ELAY           |<->| THICKNESS OF SURFACE LAYER
00127 !| ELAY0          |<->| ACTIVE LAYER THICKNESS
00128 !| ENTET          |-->| LOGICAL, IF YES INFORMATION IS GIVEN ON MASS CONSERVATION
00129 !| ENTETS         |<->| LOGICAL, IF YES INFORMATION IS GIVEN ON MASS CONSERVATION FOR SUSPENSION
00130 !| FD90           |-->| DIAMETER D90
00131 !| FDM            |-->| DIAMETER DM FOR EACH CLASS
00132 !| FLBCLA         |-->| BLOCK OF FLUXES AT BOUNDARY FOR EACH CLASS
00133 !| FRACSED_GF     |<->| (A SUPPRIMER)
00134 !| FW             |-->| WAVE FRICTION FACTOR
00135 !| GRAV           |-->| ACCELERATION OF GRAVITY
00136 !| HIDFAC         |-->| HIDING FACTOR FORMULAS
00137 !| HIDI           |-->| HIDING FACTOR FOR PARTICULAR SIZE CLASS (HIDFAC =0)
00138 !| HIDING         |-->| HIDING FACTOR CORRECTION
00139 !| HMIN           |-->| MINIMUM VALUE OF WATER DEPTH
00140 !| HN             |-->| WATER DEPTH
00141 !| HOULE          |-->| LOGICAL, FOR WAVE EFFECTS
00142 !| ICF            |-->| BED-LOAD OR TOTAL LOAD TRANSPORT FORMULAS
00143 !| IELMT          |-->| NUMBER OF ELEMENTS
00144 !| ISOUS          |-->| SUB-ITERATIONS
00145 !| IT1            |<->| INTEGER WORK ARRAY IN A BIEF_OBJ STRUCTURE
00146 !| KARMAN         |-->| VON KARMAN CONSTANT
00147 !| KDDL           |-->| CONVENTION FOR DEGREE OF FREEDOM
00148 !| KDIR           |-->| CONVENTION FOR DIRICHLET POINT
00149 !| KENT           |-->| CONVENTION FOR LIQUID INPUT WITH PRESCRIBED VALUE
00150 !| KINC           |-->| CONVENTION FOR INCIDENT WAVE BOUNDARY CONDITION
00151 !| KLOG           |-->| CONVENTION FOR SOLID BOUNDARY
00152 !| KNEU           |-->| CONVENTION FOR NEUMANN CONDITION
00153 !| KSORT          |-->| CONVENTION FOR FREE OUTPUT
00154 !| KSP            |-->| BED SKIN ROUGHNESS
00155 !| KSR            |-->| RIPPLE BED ROUGHNESS
00156 !| LCONDIS        |-->| LOGICAL, CONSTANT FLOW DISCHARGE
00157 !| LGRAFED        |-->| (A SUPPRIMER)
00158 !| LIEBOR         |<->| TYPE OF BOUNDARY CONDITIONS FOR BED EVOLUTION
00159 !| LIMTEC         |<->| TECHNICAL BOUNDARY CONDITION (NEUMAN...)
00160 !| LIQBOR         |-->| TYPE OF BOUNDARY CONDITION FOR QS
00161 !| LOADMETH       |-->| (A SUPPRIMER)
00162 !| LS0            |-->| (A SUPPRIMER)
00163 !| LT             |-->| ITERATION
00164 !| MASK           |-->| BLOCK OF MASKS, EVERY ONE FOR A TYPE OF BOUNDARY
00165 !|                |   | SEE DIFFIN.F IN LIBRARY BIEF.
00166 !| MASKEL         |-->| MASKING OF ELEMENTS
00167 !| MASKPT         |-->| MASKING PER POINT
00168 !| MASKTR         |<->| MASKING FOR TRACERS, PER POINT
00169 !| MESH           |<->| MESH STRUCTURE
00170 !| MSK            |-->| IF YES, THERE IS MASKED ELEMENTS
00171 !| MU             |<->| CORRECTION FACTOR FOR BED ROUGHNESS
00172 !| NPOIN          |-->| NUMBER OF POINTS
00173 !| NPTFR          |-->| NUMBER OF BOUNDARY POINTS
00174 !| NSICLA         |-->| NUMBER OF SIZE CLASSES FOR BED MATERIALS
00175 !| NSOUS          |<->| NUMBER OF SUB-ITERATIONS
00176 !| OPTBAN         |-->| OPTION FOR TIDAL FLATS
00177 !| PHISED         |-->| ANGLE OF REPOSE OF THE SEDIMENT
00178 !| PI             |-->| PI
00179 !| Q              |-->| FLOW DISCHARGE
00180 !| QBOR           |-->| BOUNDARY CONDITION FOR TRANSPORT RATE
00181 !| QSCLXC         |<->| TRANSPORT RATE FOR EACH CLASS X-DIRECTION
00182 !| QSCLYC         |<->| TRANSPORT RATE FOR EACH CLASS Y-DIRECTION
00183 !| QSCL_C         |<->| BEDLOAD TRANSPORT RATE
00184 !| QSCL_S         |<->| SUSPENDED LOAD TRANSPORT RATE
00185 !| QSXC           |<->| BEDLOAD TRANSPORT RATE X-DIRECTION
00186 !| QSYC           |<->| BEDLOAD TRANSPORT RATE Y-DIRECTION
00187 !| QS_C           |<->| BEDLOAD TRANSPORT RATE
00188 !| S              |-->| VOID STRUCTURE
00189 !| SALFA          |<->| SINUS OF THE ANGLE BETWEEN TRANSPORT RATE AND CURRENT
00190 !| SECCURRENT     |-->| LOGICAL, PARAMETRISATION FOR SECONDARY CURRENTS
00191 !| SEDCO          |-->| LOGICAL, SEDIMENT COHESIVE OR NOT
00192 !| SLOPEFF        |-->| LOGICAL, SLOPING BED EFFECT OR NOT
00193 !| SUSP           |-->| LOGICAL, SUSPENSION
00194 !| T1             |<->| WORK BIEF_OBJ STRUCTURE
00195 !| T10            |<->| WORK BIEF_OBJ STRUCTURE
00196 !| T11            |<->| WORK BIEF_OBJ STRUCTURE
00197 !| T12            |<->| WORK BIEF_OBJ STRUCTURE
00198 !| T13            |<->| WORK BIEF_OBJ STRUCTURE
00199 !| T2             |<->| WORK BIEF_OBJ STRUCTURE
00200 !| T3             |<->| WORK BIEF_OBJ STRUCTURE
00201 !| T4             |<->| WORK BIEF_OBJ STRUCTURE
00202 !| T5             |<->| WORK BIEF_OBJ STRUCTURE
00203 !| T6             |<->| WORK BIEF_OBJ STRUCTURE
00204 !| T7             |<->| WORK BIEF_OBJ STRUCTURE
00205 !| T8             |<->| WORK BIEF_OBJ STRUCTURE
00206 !| T9             |<->| WORK BIEF_OBJ STRUCTURE
00207 !| THETAW         |-->| ANGLE BETWEEN WAVE AND CURRENT
00208 !| TOB            |<->| BED SHEAR STRESS (TOTAL FRICTION)
00209 !| TOBW           |-->| WAVE INDUCED SHEAR STRESS
00210 !| TW             |-->| WAVE PERIOD
00211 !| U2D            |<->| MEAN FLOW VELOCITY X-DIRECTION
00212 !| U3D            |-->| THREE-DIMENSIONAL VELOCITY X-DIRECTION
00213 !| UNLADM         |-->| MEAN DIAMETER OF ACTIVE STRATUM LAYER
00214 !| UNORM          |<->| NORM OF THE MEAN FLOW VELOCITY
00215 !| UNSV2D         |-->| INVERSE OF INTEGRALS OF TEST FUNCTIONS
00216 !| UW             |-->| ORBITAL WAVE VELOCITY
00217 !| V2D            |<->| MEAN FLOW VELOCITY Y-DIRECTION
00218 !| V2DPAR         |-->| INTEGRAL OF TEST FUNCTIONS, ASSEMBLED IN PARALLEL
00219 !| V3D            |-->| THREE-DIMENSIONAL VELOCITY Y-DIRECTION
00220 !| VCE            |-->| WATER VISCOSITY
00221 !| VF             |-->| LOGICAL, FINITE VOLUMES OR NOT
00222 !| CSF_SABLE      |-->| BED VOLUME CONCENTRATION CSF = (1-POROSITY)
00223 !| XMVE           |-->| FLUID DENSITY
00224 !| XMVS           |-->| SEDIMENT DENSITY
00225 !| XWC            |-->| SETTLING VELOCITY
00226 !| ZERO           |-->| ZERO
00227 !| ZF             |-->| ELEVATION OF BOTTOM
00228 !| ZFCL_C         |<->| BEDLOAD EVOLUTION FOR EACH SEDIMENT CLASS
00229 !| ZF_C           |<->| BEDLOAD EVOLUTION
00230 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00231 !
00232       USE BIEF
00233       USE INTERFACE_SISYPHE, EX_BEDLOAD_MAIN => BEDLOAD_MAIN
00234       USE DECLARATIONS_SISYPHE, ONLY : DREDGESIM,NOMBLAY
00235       IMPLICIT NONE
00236       INTEGER LNG,LU
00237       COMMON/INFO/LNG,LU
00238 !
00239 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00240 !
00241       TYPE(BIEF_OBJ),   INTENT(IN)    :: ACLADM, KSR,V2DPAR,UNSV2D
00242       TYPE(BIEF_OBJ),   INTENT(IN)    :: CF,FW,KSP,HN,LIQBOR
00243       TYPE(BIEF_OBJ),   INTENT(IN)    :: MASK, MASKEL, MASKPT
00244       TYPE(BIEF_OBJ),   INTENT(IN)    :: Q, QBOR
00245       TYPE(BIEF_OBJ),   INTENT(INOUT) :: U2D,V2D,TOB,MU,UNORM,EBOR
00246       TYPE(BIEF_OBJ),   INTENT(IN)    :: S,UNLADM
00247       TYPE(BIEF_OBJ),   INTENT(IN)    :: UW, THETAW,  TOBW, TW
00248       TYPE(BIEF_OBJ),   INTENT(IN)    :: ZF
00249       INTEGER,          INTENT(IN)    :: DEBUG, HIDFAC, ICF,MAXADV
00250       INTEGER,          INTENT(IN)    :: IELMT, ISOUS, KDDL, KDIR, KENT
00251       INTEGER,          INTENT(IN)    :: KINC, KLOG, KNEU, KSORT
00252       INTEGER,          INTENT(IN)    :: LOADMETH, LT,NPOIN, NPTFR
00253       INTEGER,          INTENT(IN)    :: NSICLA, OPTBAN
00254       DOUBLE PRECISION, INTENT(IN)    :: LS0, BETA, FD90(NSICLA)
00255       DOUBLE PRECISION, INTENT(IN)    :: FDM(NSICLA),GRAV
00256       DOUBLE PRECISION, INTENT(IN)    :: HIDI(NSICLA),HMIN,VCE
00257       DOUBLE PRECISION, INTENT(IN)    :: CSF_SABLE,XMVE,XMVS,XWC(NSICLA)
00258       DOUBLE PRECISION, INTENT(IN)    :: PI,KARMAN,ZERO
00259       DOUBLE PRECISION, INTENT(IN)    :: KARIM_HOLLY_YANG
00260       LOGICAL,          INTENT(IN)    :: MSK, SUSP, VF
00261       LOGICAL,          INTENT(IN)    :: ENTET, CONST_ALAYER
00262       LOGICAL,          INTENT(IN)    :: LCONDIS, LGRAFED,SECCURRENT
00263       LOGICAL,          INTENT(IN)    :: SEDCO(NSICLA),HOULE
00264       TYPE(BIEF_MESH),  INTENT(INOUT) :: MESH
00265       TYPE(BIEF_OBJ),   INTENT(INOUT) :: ELAY,FLBCLA
00266       TYPE(BIEF_OBJ),   INTENT(INOUT) :: LIEBOR, LIMTEC, MASKTR
00267       TYPE(BIEF_OBJ),   INTENT(INOUT) :: IT1,T1,T2,T3,T4,T5,T6,T7
00268       TYPE(BIEF_OBJ),   INTENT(INOUT) :: T8,T9,T10,T11,T12,T13
00269       DOUBLE PRECISION, INTENT(INOUT) :: AC(NSICLA), AT0, DTS, ELAY0
00270       DOUBLE PRECISION, INTENT(INOUT) :: FRACSED_GF(NSICLA)
00271       DOUBLE PRECISION, INTENT(INOUT) :: AVAIL(NPOIN,NOMBLAY,NSICLA)
00272       TYPE(BIEF_OBJ),   INTENT(INOUT) :: BREACH, CALFA, COEFPN
00273       TYPE(BIEF_OBJ),   INTENT(INOUT) :: DZF_GF
00274       TYPE(BIEF_OBJ),   INTENT(INOUT) :: HIDING
00275       TYPE(BIEF_OBJ),   INTENT(INOUT) :: QSCL_C,QSCL_S
00276       TYPE(BIEF_OBJ),   INTENT(INOUT) :: QS_C, QSCLXC, QSXC, QSCLYC
00277       TYPE(BIEF_OBJ),   INTENT(INOUT) :: QSYC, SALFA, ZF_C, ZFCL_C
00278       INTEGER,          INTENT(INOUT) :: NSOUS
00279       LOGICAL,          INTENT(INOUT) :: ENTETS
00280       DOUBLE PRECISION,   INTENT(IN)  :: BETA2, PHISED
00281       INTEGER, INTENT (IN)            :: SLOPEFF, DEVIA
00282       DOUBLE PRECISION, INTENT(IN)    :: BIJK
00283       TYPE(BIEF_OBJ),    INTENT(IN)    :: U3D,V3D
00284       CHARACTER(LEN=24), INTENT(IN)    :: CODE
00285 !
00286 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00287 !
00288       INTEGER I
00289 !
00290 !
00291 !======================================================================!
00292 !======================================================================!
00293 !                               PROGRAM                                !
00294 !======================================================================!
00295 !======================================================================!
00296 !
00297 !
00298 !     INITIALISES TECHNICAL BOUNDARY CONDITIONS
00299 !
00300       IF (DEBUG > 0) WRITE(LU,*) 'BEDLOAD_DIFFIN'
00301       CALL BEDLOAD_DIFFIN
00302      &        (U2D, V2D, MESH%NBOR, MESH%XNEBOR, MESH%YNEBOR,
00303      &         MASKEL, MESH%NELBOR, NPTFR, KENT, KSORT, KLOG,
00304      &         KDIR, KDDL, KNEU, MSK, IT1, LIEBOR, MASKTR,LIMTEC,
00305      &         MESH%IKLBOR%I,MESH%NELEB,MESH%NELEBX)
00306       IF (DEBUG > 0) WRITE(LU,*) 'END_BEDLOAD_DIFFIN'
00307 !
00308       DO I = 1, NSICLA
00309 !
00310 !       FOR SAND
00311         IF(.NOT.SEDCO(I)) THEN
00312           IF (DEBUG > 0) WRITE(LU,*)
00313      &      'BEDLOAD_SOLIDISCHARGE : ',I,'/',NSICLA
00314           CALL BEDLOAD_SOLIDISCHARGE
00315      &       (MESH, U2D, V2D, UNORM,HN, TW, UW, MU,TOB,CF,
00316      &         TOBW,FW,THETAW,AVAIL(1:NPOIN,1,I),
00317      &         MASKPT, MASKEL, ACLADM,
00318      &         UNLADM,KSP,KSR, LIQBOR, QBOR%ADR(I)%P, DEBUG, NPOIN,
00319      &         NPTFR, IELMT, ICF, KENT, OPTBAN, HIDFAC, GRAV,
00320      &         FDM(I), FD90(I), XWC(I), XMVE, XMVS, VCE, HMIN,
00321      &         HIDI(I),KARMAN,ZERO,PI,
00322      &         KARIM_HOLLY_YANG,SUSP,MSK,T1,T2,
00323      &         T3, T4, T5, T6, T7, T8, T9, T10, T11,T12, AC(I),
00324      &         HIDING,QSCL_C%ADR(I)%P,QSCL_S%ADR(I)%P,
00325      &         SLOPEFF,COEFPN,PHISED,CALFA,SALFA,BETA,ZF,S,
00326      &         DEVIA, BETA2 , SECCURRENT, BIJK,HOULE,UNSV2D,
00327      &         U3D,V3D,CODE)
00328           IF(DEBUG > 0) WRITE(LU,*) 'END_BEDLOAD_SOLIDISCHARGE'
00329         ELSE
00330 !         FOR COHESIVE SEDIMENT: ZERO BEDLOAD TRANSPORT RATE
00331 !         JMH: IS THIS USEFUL ???
00332           CALL OS('X=0     ',X=QSCL_C%ADR(I)%P)
00333           CALL OS('X=0     ',X=QSCLXC%ADR(I)%P)
00334           CALL OS('X=0     ',X=QSCLYC%ADR(I)%P)
00335         ENDIF
00336 !
00337       ENDDO
00338 !
00339 !     COMPUTES THE EVOLUTION FOR EACH CLASS
00340 !
00341       DO I = 1, NSICLA
00342 !
00343         IF(.NOT.SEDCO(I)) THEN
00344 !
00345           IF (DEBUG > 0) WRITE(LU,*) 'BEDLOAD_EVOL : ',I,'/',NSICLA
00346           CALL BEDLOAD_EVOL(HN,Q,S,ELAY,ACLADM,AVAIL(1:NPOIN,1,I),
00347      &                      COEFPN,CALFA,SALFA,LIMTEC,
00348      &                      EBOR%ADR(I)%P,MASKEL,MASK,
00349      &                      V2DPAR,UNSV2D,DEBUG,NPOIN,NPTFR,IELMT,
00350      &                      KENT,KDIR,KDDL,LOADMETH,
00351      &                      DTS,FDM(I),FD90(I),HMIN,LS0,GRAV,XMVS,XMVE,
00352      &                      VCE,VF,ENTETS,MSK,CONST_ALAYER,
00353      &                      LCONDIS,MESH,QSCL_C%ADR(I)%P,
00354      &                      T1,T2,T3,T4,T5,T6,T7,T8,T9,T10,T11,T12,
00355      &                      T13,CSF_SABLE,BREACH,QSCLXC%ADR(I)%P,
00356      &                      QSCLYC%ADR(I)%P,ZFCL_C%ADR(I)%P,SLOPEFF,
00357      &                      I,FLBCLA,LIQBOR,QBOR%ADR(I)%P,MAXADV)
00358           IF(DEBUG.GT.0) WRITE(LU,*) 'END_BEDLOAD_EVOL'
00359 !
00360 !         NOW DIVIDING BY CSF_SABLE TO GET THE EVOLUTION OF BED
00361 !         INCLUDING VOIDS
00362           CALL OS('X=CX    ',X= ZFCL_C%ADR(I)%P,C=1.D0/CSF_SABLE)
00363 !
00364         ELSE
00365 !
00366 !         NO EVOLUTION FOR COHESIVE SEDIMENT
00367           CALL OS('X=0     ',X=ZFCL_C%ADR(I)%P)
00368 !
00369         ENDIF
00370 !
00371       ENDDO
00372 !
00373 !     CALLS DREDGESIM
00374 !
00375       IF(DREDGESIM) CALL DREDGESIM_INTERFACE(2)
00376       ! *********************************************** !
00377       ! II - EVOLUTIONS AND QS FOR EACH CLASS ARE ADDED !
00378       ! *********************************************** !
00379       ! II.1 - INITIALISES
00380       ! ---------------------
00381       CALL OS('X=0     ', X=QS_C)
00382       CALL OS('X=0     ', X=ZF_C)
00383       ! II.2 - ADDS THE CLASSES
00384       ! ----------------------
00385       !
00386       DO I=1,NSICLA
00387         IF(.NOT.SEDCO(I)) THEN
00388           CALL OS('X=X+Y   ', X=QS_C, Y=QSCL_C%ADR(I)%P)
00389           CALL OS('X=X+Y   ', X=ZF_C, Y=ZFCL_C%ADR(I)%P)
00390         ENDIF
00391       ENDDO
00392 !
00393 !     TIDAL FLATS WITH MASKING     JMH ON 27/07/2006
00394 !
00395       IF(OPTBAN.EQ.2) CALL OS('X=XY    ',X=ZF_C,Y=MASKPT)
00396 !
00397       ! II.3 - SLOPE EFFECT FOR THE SUM OF THE QS
00398       ! -----------------------------------------
00399       ! QS : COEFPN ALREADY ADDED IN QSCL_C
00400       CALL OS('X=YZ    ', X=QSXC, Y=QS_C, Z=CALFA)
00401       CALL OS('X=YZ    ', X=QSYC, Y=QS_C, Z=SALFA)
00402 !
00403 !======================================================================!
00404 !======================================================================!
00405 !
00406       RETURN
00407       END

Generated on Fri Aug 31 2013 18:12:58 by S.E.Bourban (HRW) using doxygen 1.7.0