source: LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/ini_wake.F @ 1373

Last change on this file since 1373 was 1336, checked in by idelkadi, 15 years ago
  • Remise de la valeur par default de alp_offset a 0.
  • Ajout de la possiblite de lecture des parametres de ini_wake.F dans le fichier ini_wake_param.data
  • Ajout de la possiblite de lecture des parametres de wake.F dans le fichier wake_param.data
  • Correction dans la partie convection (nouvelle physique)
  • 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 1336 2010-04-01 12:27:52Z musat $
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.