source: dynamico_lmdz/aquaplanet/LMDZ5/libf/phylmd/ini_wake.F90 @ 3817

Last change on this file since 3817 was 3817, checked in by millour, 10 years ago

Further cleanup and removal of references to iniprint.h.
Also added bench testcase 48x36x19.
EM

File size: 2.3 KB
Line 
1
2! $Id: ini_wake.F90 2197 2015-02-09 07:13:05Z emillour $
3
4SUBROUTINE ini_wake(wape, fip, it_wape_prescr, wape_prescr, fip_prescr, &
5    alp_bl_prescr, ale_bl_prescr)
6  USE inifis_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  it = 0.
52  w = 4.
53  f = 0.1
54  alebl = 4.
55  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    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
65901 CLOSE (99)
66902 CONTINUE
67
68  WRITE (lunout, *) ' it,wape ', it, wape
69  it_wape_prescr = it
70  IF (w<0) THEN
71    wape_prescr = wape
72    fip_prescr = fip
73  ELSE
74    wape_prescr = w
75    fip_prescr = f
76  END IF
77
78  WRITE (lunout, *) ' alebl, alpbl ', alebl, alpbl
79  ale_bl_prescr = alebl
80  alp_bl_prescr = alpbl
81  PRINT *, 'Initialisation de la poche : WAPE, FIP imposees =', wape_prescr, &
82    fip_prescr
83  PRINT *, '                   pendant ', it_wape_prescr, ' steps'
84
85  PRINT *, 'Initialisation de la BL: ALP, ALE imposees =', alp_bl_prescr, &
86    ale_bl_prescr
87  RETURN
88END SUBROUTINE ini_wake
Note: See TracBrowser for help on using the repository browser.