source: LMDZ5/branches/LMDZ5_SPLA/libf/phylmd/ini_wake.F90 @ 2404

Last change on this file since 2404 was 1992, checked in by lguez, 11 years ago

Converted to free source form files in libf/phylmd which were still in
fixed source form. The conversion was done using the polish mode of
the NAG Fortran Compiler.

In addition to converting to free source form, the processing of the
files also:

-- indented the code (including comments);

-- set Fortran keywords to uppercase, and set all other identifiers
to lower case;

-- added qualifiers to end statements (for example "end subroutine
conflx", instead of "end");

-- changed the terminating statements of all DO loops so that each
loop ends with an ENDDO statement (instead of a labeled continue).

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