Changeset 1403 for LMDZ4/trunk/libf/phylmd/calwake.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/calwake.F (modified) (6 diffs)

LMDZ4/trunk/libf/phylmd/calwake.F

@@ +1 @@
+!
+! $Id$
+!
       SUBROUTINE CALWAKE(paprs,pplay,dtime
      :             ,t,q,omgb
@@ -189 +192 @@
      $                ,Cstar,d_deltat_gw
      $                ,d_deltatw,d_deltaqw)
-
-      DO i=1,klon 
-       IF (ktopw(i) .GT. 0) THEN
-         DO l=1,klev
+c
+      DO l=1,klev
+       DO i=1,klon 
+        IF (ktopw(i) .GT. 0) THEN
            wake_deltat(i,l)= dtw(i,l)
            wake_deltaq(i,l)= dqw(i,l)
@@ -212 +215 @@
            wake_ddeltat(i,l) = d_deltatw(i,l)
            wake_ddeltaq(i,l) = d_deltaqw(i,l)
-         ENDDO
-       ELSE
-         DO l = 1,klev
+        ELSE
            wake_deltat(i,l)= 0.
            wake_deltaq(i,l)= 0.
@@ -222 +223 @@
            wake_dtKE(i,l) = 0.
            wake_dqKE(i,l) = 0.
+           wake_dtPBL(i,l) = 0.
+           wake_dqPBL(i,l) = 0.
            wake_omg(i,l) = 0.
            wake_dp_deltomg(i,l) = 0.
@@ -230 +233 @@
            undi_t(i,l)=te(i,l)
            undi_q(i,l)=qe(i,l)
-         ENDDO
-       ENDIF
-
+           wake_ddeltat(i,l) = 0.
+           wake_ddeltaq(i,l) = 0.
+        ENDIF
+       ENDDO
+      ENDDO
+c
+      DO i=1,klon 
        wake_h(i)= hw(i)
        wake_s(i)= sigmaw(i)
@@ -241 +248 @@
        wake_Cstar(i) = Cstar(i)
        wake_dens(i) = wdens(i)
-c
-cIM 290108 999  CONTINUE
-c
-      ENDDO
+      ENDDO
+c
       RETURN
       END
+
       SUBROUTINE CALWAKE_scal(paprs,pplay,dtime
      :             ,t,q,omgb