source: LMDZ6/branches/LMDZ_ECRad/libf/phylmd/ini_wake.F90 @ 5306

Last change on this file since 5306 was 2311, checked in by Ehouarn Millour, 9 years ago

Further modifications to enforce physics/dynamics separation:

  • moved iniprint.h and misc_mod back to dyn3d_common, as these should only be used by dynamics.
  • created print_control_mod in the physics to store flags prt_level, lunout, debug to be local to physics (should be used rather than iniprint.h)
  • created abort_physic.F90 , which does the same job as abort_gcm() did, but should be used instead when in physics.
  • reactivated inifis (turned it into a module, inifis_mod.F90) to initialize physical constants and print_control_mod flags.

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