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 |