source: LMDZ5/branches/LMDZ5V1.0-dev/libf/phylmd/ini_wake.F @ 1523

Last change on this file since 1523 was 1403, checked in by Laurent Fairhead, 14 years ago

Merged LMDZ4V5.0-dev branch changes r1292:r1399 to trunk.

Validation:
Validation consisted in compiling the HEAD revision of the trunk,
LMDZ4V5.0-dev branch and the merged sources and running different
configurations on local and SX8 machines comparing results.

Local machine: bench configuration, 32x24x11, gfortran

  • IPSLCM5A configuration (comparison between trunk and merged sources):
    • numerical convergence on dynamical fields over 3 days
    • start files are equivalent (except for RN and PB fields)
    • daily history files equivalent
  • MH07 configuration, new physics package (comparison between LMDZ4V5.0-dev branch and merged sources):
    • numerical convergence on dynamical fields over 3 days
    • start files are equivalent (except for RN and PB fields)
    • daily history files equivalent

SX8 machine (brodie), 96x95x39 on 4 processors:

  • IPSLCM5A configuration:
    • start files are equivalent (except for RN and PB fields)
    • monthly history files equivalent
  • MH07 configuration:
    • start files are equivalent (except for RN and PB fields)
    • monthly history files equivalent

Changes to the makegcm and create_make_gcm scripts to take into account
main programs in F90 files


Fusion de la branche LMDZ4V5.0-dev (r1292:r1399) au tronc principal

Validation:
La validation a consisté à compiler la HEAD de le trunk et de la banche
LMDZ4V5.0-dev et les sources fusionnées et de faire tourner le modéle selon
différentes configurations en local et sur SX8 et de comparer les résultats

En local: 32x24x11, config bench/gfortran

  • pour une config IPSLCM5A (comparaison tronc/fusion):
    • convergence numérique sur les champs dynamiques après 3 jours
    • restart et restartphy égaux (à part sur RN et Pb)
    • fichiers histoire égaux
  • pour une config nlle physique (MH07) (comparaison LMDZ4v5.0-dev/fusion):
    • convergence numérique sur les champs dynamiques après 3 jours
    • restart et restartphy égaux
    • fichiers histoire équivalents

Sur brodie, 96x95x39 sur 4 proc:

  • pour une config IPSLCM5A:
    • restart et restartphy égaux (à part sur RN et PB)
    • pas de différence dans les fichiers histmth.nc
  • pour une config MH07
    • restart et restartphy égaux (à part sur RN et PB)
    • pas de différence dans les fichiers histmth.nc

Changement sur makegcm et create_make-gcm pour pouvoir prendre en compte des
programmes principaux en *F90

  • 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 1403 2010-07-01 09:02:53Z idelkadi $
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.