GCC Code Coverage Report
Directory: ./ Exec Total Coverage
File: phylmd/ini_wake.F90 Lines: 21 29 72.4 %
Date: 2023-06-30 12:56:34 Branches: 2 6 33.3 %

Line Branch Exec Source
1
2
! $Id: ini_wake.F90 2311 2015-06-25 07:45:24Z emillour $
3
4
1
SUBROUTINE ini_wake(wape, fip, it_wape_prescr, wape_prescr, fip_prescr, &
5
    alp_bl_prescr, ale_bl_prescr)
6
  USE print_control_mod, ONLY: lunout
7
  IMPLICIT NONE
8
  ! **************************************************************
9
  ! *
10
  ! INI_WAKE : variables d'initialisation de la poche    *
11
  ! froide, necessaires au declenchement      *
12
  ! de la convection.                         *
13
  ! *
14
  ! *
15
  ! **************************************************************
16
  ! Arguments
17
  ! =========
18
  ! Input
19
  ! -----
20
  ! wape           : valeur de l'energie potentielle de la poche (WAPE)
21
  ! dans l'etat initial
22
  ! fip            : valeur de la puissance incidente sur le front (FIP)
23
  ! dans l'etat initial
24
  ! Output
25
  ! ------
26
  ! it_wape_prescr : nombre de pas de temps pendant lesquels la WAPE
27
  ! doit etre imposee.
28
  ! wape_prescr    : valeur prescrite de la WAPE.
29
  ! fip_prescr     : valeur prescrite de la FIP.
30
  ! ale_bl_prescr  : valeur prescrite de la Ale de PBL.
31
  ! alp_bl_prescr  : valeur prescrite de la Alp de PBL.
32
33
  ! Variables internes
34
  ! ==================
35
  ! it = nbre de pas de temps lu
36
  ! w  = WAPE lue
37
  ! f  = FIP lue
38
  ! alebl  = Ale de PBL lue
39
  ! alpbl  = Alp de PBL lue
40
41
  ! declarations
42
  REAL wape, fip, wape_prescr, fip_prescr
43
  INTEGER it_wape_prescr
44
  REAL ale_bl_prescr
45
  REAL alp_bl_prescr
46
  REAL it
47
  REAL w,f,alebl,alpbl
48
49
  ! FH A mettre si besoin dans physiq.def
50
  ! FH : voir avec JYG
51
1
  it = 0.
52
1
  w = 4.
53
1
  f = 0.1
54
1
  alebl = 4.
55
1
  alpbl = 0.1
56
57
  ! CR: on rajoute ale et alp de la PBL precrits
58
  OPEN (99, FILE='ini_wake_param.data', FORM='formatted', STATUS='old', &
59
1
    ERR=902)
60
  READ (99, *) it
61
  READ (99, *) w
62
  READ (99, *) f
63
  READ (99, *, END=901) alebl
64
  READ (99, *, END=901) alpbl
65
901 CLOSE (99)
66
902 CONTINUE
67
68
1
  WRITE (lunout, *) ' it,wape ', it, wape
69
1
  it_wape_prescr = it
70
1
  IF (w<0) THEN
71
    wape_prescr = wape
72
    fip_prescr = fip
73
  ELSE
74
1
    wape_prescr = w
75
1
    fip_prescr = f
76
  END IF
77
78
1
  WRITE (lunout, *) ' alebl, alpbl ', alebl, alpbl
79
1
  ale_bl_prescr = alebl
80
1
  alp_bl_prescr = alpbl
81
1
  PRINT *, 'Initialisation de la poche : WAPE, FIP imposees =', wape_prescr, &
82
2
    fip_prescr
83
1
  PRINT *, '                   pendant ', it_wape_prescr, ' steps'
84
85
1
  PRINT *, 'Initialisation de la BL: ALP, ALE imposees =', alp_bl_prescr, &
86
2
    ale_bl_prescr
87
1
  RETURN
88
END SUBROUTINE ini_wake