source: LMDZ5/trunk/libf/phylmd/ini_wake.F @ 1907

Last change on this file since 1907 was 1907, checked in by lguez, 10 years ago

Added a copyright property to every file of the distribution, except
for the fcm files (which have their own copyright). Use svn propget on
a file to see the copyright. For instance:

$ svn propget copyright libf/phylmd/physiq.F90
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

Also added the files defining the CeCILL version 2 license, in French
and English, at the top of the LMDZ tree.
  • 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.5 KB
Line 
1!
2! $Id: ini_wake.F 1907 2013-11-26 13:10:46Z lguez $
3!
4      SUBROUTINE INI_WAKE(wape,fip,it_wape_prescr,
5     :     wape_prescr, fip_prescr, 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***************************************************************
14c Arguments
15c =========
16c Input
17c -----
18c   wape           : valeur de l'energie potentielle de la poche (WAPE)
19c                    dans l'etat initial
20c   fip            : valeur de la puissance incidente sur le front (FIP)
21c                    dans l'etat initial
22c Output
23c ------
24c   it_wape_prescr : nombre de pas de temps pendant lesquels la WAPE
25c            doit etre imposee.
26c   wape_prescr    : valeur prescrite de la WAPE.
27c   fip_prescr     : valeur prescrite de la FIP.
28c   ale_bl_prescr  : valeur prescrite de la Ale de PBL.
29c   alp_bl_prescr  : valeur prescrite de la Alp de PBL.
30c
31c Variables internes
32c ==================
33c   it = nbre de pas de temps lu
34c   w  = WAPE lue
35c   f  = FIP lue
36c   alebl  = Ale de PBL lue
37c   alpbl  = Alp de PBL lue
38c
39      include 'iniprint.h'
40cdeclarations
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
52c
53cCR: on rajoute ale et alp de la PBL precrits
54      open (99,file='ini_wake_param.data',form='formatted',
55     s      status='old',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
63c
64      write(lunout,*)' it,wape ',it,wape
65      it_wape_prescr = it
66      if (w .lt. 0) then
67         wape_prescr = wape
68         fip_prescr = fip
69      else
70         wape_prescr = w
71         fip_prescr = f
72      endif
73c
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 ='
78     $               ,wape_prescr, fip_prescr
79      print *, '                   pendant ',it_wape_prescr,' steps'
80c
81      print *,'Initialisation de la BL: ALP, ALE imposees ='
82     $               ,alp_bl_prescr, ale_bl_prescr
83      return
84      end
Note: See TracBrowser for help on using the repository browser.