source: LMDZ5/trunk/libf/phylmd/ini_wake.F90 @ 2201

Last change on this file since 2201 was 2197, checked in by Ehouarn Millour, 10 years ago

Added 'implicit none' statements and proper variable definitions where they were missing.
EM

  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 2.3 KB
Line 
1
2! $Id: ini_wake.F90 2197 2015-02-09 07:13:05Z crio $
3
4SUBROUTINE ini_wake(wape, fip, it_wape_prescr, wape_prescr, fip_prescr, &
5    alp_bl_prescr, ale_bl_prescr)
6  IMPLICIT NONE
7  ! **************************************************************
8  ! *
9  ! INI_WAKE : variables d'initialisation de la poche    *
10  ! froide, necessaires au declenchement      *
11  ! de la convection.                         *
12  ! *
13  ! *
14  ! **************************************************************
15  ! Arguments
16  ! =========
17  ! Input
18  ! -----
19  ! wape           : valeur de l'energie potentielle de la poche (WAPE)
20  ! dans l'etat initial
21  ! fip            : valeur de la puissance incidente sur le front (FIP)
22  ! dans l'etat initial
23  ! Output
24  ! ------
25  ! it_wape_prescr : nombre de pas de temps pendant lesquels la WAPE
26  ! doit etre imposee.
27  ! wape_prescr    : valeur prescrite de la WAPE.
28  ! fip_prescr     : valeur prescrite de la FIP.
29  ! ale_bl_prescr  : valeur prescrite de la Ale de PBL.
30  ! alp_bl_prescr  : valeur prescrite de la Alp de PBL.
31
32  ! Variables internes
33  ! ==================
34  ! it = nbre de pas de temps lu
35  ! w  = WAPE lue
36  ! f  = FIP lue
37  ! alebl  = Ale de PBL lue
38  ! alpbl  = Alp de PBL lue
39
40  include 'iniprint.h'
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.