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

Last change on this file since 1280 was 970, checked in by Laurent Fairhead, 16 years ago

Corrections petits bugs divers
LF

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 2.2 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
37c     open (99,file='wake.data',form='formatted')
38c     read (99,*) it
39c     read (99,*) w
40c     read (99,*) f
41c     read (99,*) u
42c     read (99,*) p
43c     close (99)
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      u=0.1
51      p=4.
52c
53      print *,' it,w ',it,w
54      it_wape_prescr = it
55      if (w .lt. 0) then
56         wape_prescr = wape
57         fip_prescr = fip
58      else
59         wape_prescr = w
60         fip_prescr = f
61      endif
62c
63      print *,' u,p ',u,p
64      alp_bl_prescr=u
65      ale_bl_prescr=p
66      print *,'Initialisation de la poche : WAPE, FIP imposees ='
67     $               ,wape_prescr, fip_prescr
68      print *, '                   pendant ',it_wape_prescr,' steps'
69c
70      print *,'Initialisation de la BL: ALP, ALE imposees ='
71     $               ,alp_bl_prescr, ale_bl_prescr
72      return
73      end
Note: See TracBrowser for help on using the repository browser.