source: LMDZ4/trunk/libf/phylmd/moy_undefSTD.F @ 1664

Last change on this file since 1664 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: 4.0 KB
Line 
1!
2! $Id: moy_undefSTD.F 1403 2010-07-01 09:02:53Z fairhead $
3!
4      SUBROUTINE moy_undefSTD(itap,freq_outNMC,freq_moyNMC)
5      USE netcdf
6      USE dimphy
7      USE phys_state_var_mod ! Variables sauvegardees de la physique
8      IMPLICIT none
9c
10c====================================================================
11c
12c I. Musat : 09.2004
13c
14c Moyenne - a des frequences differentes - des valeurs bien definies
15c         (.NE.missing_val) des variables interpolees a un niveau de
16c         pression.
17c 1) les variables de type "day" (nout=1) ou "mth" (nout=2) sont sommees
18c    tous les pas de temps de la physique
19c
20c 2) les variables de type "NMC" (nout=3) sont calculees a partir
21c    des valeurs instantannees toutes les 6 heures
22c
23c
24c NB: mettre "inst(X)" dans le write_hist*NMC.h !
25c====================================================================
26cym#include "dimensions.h"
27cym      integer jjmp1
28cym      parameter (jjmp1=jjm+1-1/jjm)
29cym#include "dimphy.h"
30c
31c
32c variables Input
33c     INTEGER nlevSTD, klevSTD, itap
34c     PARAMETER(klevSTD=17)
35      INTEGER itap
36c
37c variables locales
38c     INTEGER i, k, nout, n
39c     PARAMETER(nout=3) !nout=1 day/nout=2 mth/nout=3 NMC
40      INTEGER i, k, n
41c     REAL dtime, freq_outNMC(nout), freq_moyNMC(nout)
42      REAL freq_outNMC(nout), freq_moyNMC(nout)
43c
44c variables Output
45c     REAL tnondef(klon,klevSTD,nout)
46c     REAL tsumSTD(klon,klevSTD,nout)
47c
48      REAL missing_val
49c
50      missing_val=nf90_fill_real
51c
52      DO n=1, nout
53c
54c calcul 1 fois par jour
55c
56       IF(MOD(itap,NINT(freq_outNMC(n)/dtime)).EQ.0) THEN
57c
58        DO k=1, nlevSTD
59         DO i=1, klon
60          IF(tnondef(i,k,n).NE.(freq_moyNMC(n))) THEN
61           tsumSTD(i,k,n)=tsumSTD(i,k,n)/
62     $     (freq_moyNMC(n)-tnondef(i,k,n))
63cIM BEG
64          usumSTD(i,k,n)=usumSTD(i,k,n)/
65     $     (freq_moyNMC(n)-tnondef(i,k,n))
66          vsumSTD(i,k,n)=vsumSTD(i,k,n)/
67     $     (freq_moyNMC(n)-tnondef(i,k,n))
68          wsumSTD(i,k,n)=wsumSTD(i,k,n)/
69     $     (freq_moyNMC(n)-tnondef(i,k,n))
70          phisumSTD(i,k,n)=phisumSTD(i,k,n)/
71     $     (freq_moyNMC(n)-tnondef(i,k,n))
72          qsumSTD(i,k,n)=qsumSTD(i,k,n)/
73     $     (freq_moyNMC(n)-tnondef(i,k,n))
74          rhsumSTD(i,k,n)=rhsumSTD(i,k,n)/
75     $     (freq_moyNMC(n)-tnondef(i,k,n))
76          uvsumSTD(i,k,n)=uvsumSTD(i,k,n)/
77     $     (freq_moyNMC(n)-tnondef(i,k,n))
78          vqsumSTD(i,k,n)=vqsumSTD(i,k,n)/
79     $     (freq_moyNMC(n)-tnondef(i,k,n))
80          vTsumSTD(i,k,n)=vTsumSTD(i,k,n)/
81     $     (freq_moyNMC(n)-tnondef(i,k,n))
82          wqsumSTD(i,k,n)=wqsumSTD(i,k,n)/
83     $     (freq_moyNMC(n)-tnondef(i,k,n))
84          vphisumSTD(i,k,n)=vphisumSTD(i,k,n)/
85     $     (freq_moyNMC(n)-tnondef(i,k,n))
86          wTsumSTD(i,k,n)=wTsumSTD(i,k,n)/
87     $     (freq_moyNMC(n)-tnondef(i,k,n))
88          u2sumSTD(i,k,n)=u2sumSTD(i,k,n)/
89     $     (freq_moyNMC(n)-tnondef(i,k,n))
90          v2sumSTD(i,k,n)=v2sumSTD(i,k,n)/
91     $     (freq_moyNMC(n)-tnondef(i,k,n))
92          T2sumSTD(i,k,n)=T2sumSTD(i,k,n)/
93     $     (freq_moyNMC(n)-tnondef(i,k,n))
94          O3sumSTD(i,k,n)=O3sumSTD(i,k,n)/
95     $     (freq_moyNMC(n)-tnondef(i,k,n))
96          O3daysumSTD(i,k,n)=O3daysumSTD(i,k,n)/
97     $     (freq_moyNMC(n)-tnondef(i,k,n))
98cIM END
99          ELSE
100           tsumSTD(i,k,n)=missing_val
101           usumSTD(i,k,n)=missing_val
102           vsumSTD(i,k,n)=missing_val
103           wsumSTD(i,k,n)=missing_val
104           phisumSTD(i,k,n)=missing_val
105           qsumSTD(i,k,n)=missing_val
106           rhsumSTD(i,k,n)=missing_val
107           uvsumSTD(i,k,n)=missing_val
108           vqsumSTD(i,k,n)=missing_val
109           vTsumSTD(i,k,n)=missing_val
110           wqsumSTD(i,k,n)=missing_val
111           vphisumSTD(i,k,n)=missing_val
112           wTsumSTD(i,k,n)=missing_val
113           u2sumSTD(i,k,n)=missing_val
114           v2sumSTD(i,k,n)=missing_val
115           T2sumSTD(i,k,n)=missing_val
116           O3sumSTD(i,k,n)=missing_val
117           O3daysumSTD(i,k,n)=missing_val
118          ENDIF !tnondef(i,k,n).NE.(freq_moyNMC(n))
119         ENDDO !i
120        ENDDO !k
121       ENDIF !MOD(itap,NINT(freq_outNMC(n)/dtime)).EQ.0
122c
123      ENDDO !n
124c
125      RETURN
126      END 
Note: See TracBrowser for help on using the repository browser.