Directory: | ./ |
---|---|
File: | phys/ini_wake.f90 |
Date: | 2022-01-11 19:19:34 |
Exec | Total | Coverage | |
---|---|---|---|
Lines: | 21 | 29 | 72.4% |
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/2✗ Branch 1 not taken.
✓ Branch 2 taken 1 times.
|
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/2✗ Branch 0 not taken.
✓ Branch 1 taken 1 times.
|
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 | ||
89 |