source: LMDZ4/trunk/libf/phylmd/ini_wake.F @ 894

Last change on this file since 894 was 880, checked in by Laurent Fairhead, 17 years ago

Correction petit oubli
LF

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 2.1 KB
Line 
1      SUBROUTINE INI_WAKE(wape,fip,it_wape_prescr,
2     :     wape_prescr, fip_prescr, alp_bl_prescr, ale_bl_prescr)
3***************************************************************
4*                                                             *
5*        INI_WAKE : variables d'initialisation de la poche    *
6*                   froide, necessaires au declenchement      *
7*                   de la convection.                         *
8*                                                             *
9*                                                             *
10***************************************************************
11c Arguments
12c =========
13c Input
14c -----
15c   wape           : valeur de l'energie potentielle de la poche (WAPE)
16c                    dans l'etat initial
17c   fip            : valeur de la puissance incidente sur le front (FIP)
18c                    dans l'etat initial
19c Output
20c ------
21c   it_wape_prescr : nombre de pas de temps pendant lesquels la WAPE
22c            doit etre imposee.
23c   wape_prescr    : valeur prescrite de la WAPE.
24c   fip_prescr     : valeur prescrite de la FIP.
25c
26c Variables internes
27c ==================
28c   it = nbre de pas de temps lu
29c   w  = WAPE lue
30c   f  = FIP lue
31c
32cdeclarations
33      real ale_bl_prescr
34      real alp_bl_prescr
35      real it
36cCR: on rajoute ale et alp de la PBL precrits
37      open (99,file='wake.data',form='formatted')
38      read (99,*) it
39      read (99,*) w
40      read (99,*) f
41      read (99,*) u
42      read (99,*) p
43      close (99)
44c
45      print *,' it,w ',it,w
46      it_wape_prescr = it
47      if (w .lt. 0) then
48         wape_prescr = wape
49         fip_prescr = fip
50      else
51         wape_prescr = w
52         fip_prescr = f
53      endif
54c
55      print *,' u,p ',u,p
56      alp_bl_prescr=u
57      ale_bl_prescr=p
58      print *,'Initialisation de la poche : WAPE, FIP imposees ='
59     $               ,wape_prescr, fip_prescr
60      print *, '                   pendant ',it_wape_prescr,' steps'
61c
62      print *,'Initialisation de la BL: ALP, ALE imposees ='
63     $               ,alp_bl_prescr, ale_bl_prescr
64      return
65      end
Note: See TracBrowser for help on using the repository browser.