Changeset 1403 for LMDZ4/trunk/libf/phylmd/ini_wake.F

Timestamp:

Jul 1, 2010, 11:02:53 AM (14 years ago)

Author:

Laurent Fairhead

Message:

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

Location:

LMDZ4/trunk

Files:

• . (modified) (1 prop)
• libf/phylmd/ini_wake.F (modified) (5 diffs)

LMDZ4/trunk/libf/phylmd/ini_wake.F

@@ +1 @@
+!
+! $Id$
+!
       SUBROUTINE INI_WAKE(wape,fip,it_wape_prescr,
      :     wape_prescr, fip_prescr, alp_bl_prescr, ale_bl_prescr)
@@ -23 +26 @@
 c   wape_prescr    : valeur prescrite de la WAPE.
 c   fip_prescr     : valeur prescrite de la FIP.
+c   ale_bl_prescr  : valeur prescrite de la Ale de PBL.
+c   alp_bl_prescr  : valeur prescrite de la Alp de PBL.
 c
 c Variables internes
@@ -29 +34 @@
 c   w  = WAPE lue
 c   f  = FIP lue
+c   alebl  = Ale de PBL lue
+c   alpbl  = Alp de PBL lue
 c
+      include 'iniprint.h'
 cdeclarations
       real ale_bl_prescr
       real alp_bl_prescr
       real it
-cCR: on rajoute ale et alp de la PBL precrits
-c     open (99,file='wake.data',form='formatted')
-c     read (99,*) it
-c     read (99,*) w
-c     read (99,*) f
-c     read (99,*) u
-c     read (99,*) p
-c     close (99)
 
 ! FH A mettre si besoin dans physiq.def
@@ -48 +48 @@
       w=4.
       f=0.1
-      u=0.1
-      p=4.
+      alebl=4.
+      alpbl=0.1
 c
-      print *,' it,w ',it,w
+cCR: on rajoute ale et alp de la PBL precrits
+      open (99,file='ini_wake_param.data',form='formatted',
+     s      status='old',err=902)
+      read (99,*) it
+      read (99,*) w
+      read (99,*) f
+      read (99,*,end=901) alebl
+      read (99,*,end=901) alpbl
+901   close (99)
+902   continue
+c
+      write(lunout,*)' it,wape ',it,wape
       it_wape_prescr = it
       if (w .lt. 0) then
@@ -61 +72 @@
       endif
 c
-      print *,' u,p ',u,p
-      alp_bl_prescr=u
-      ale_bl_prescr=p
+      write(lunout,*)' alebl, alpbl ',alebl,alpbl
+      ale_bl_prescr=alebl
+      alp_bl_prescr=alpbl
       print *,'Initialisation de la poche : WAPE, FIP imposees ='
      $               ,wape_prescr, fip_prescr