source: LMDZ4/trunk/libf/dyn3dpar/interpre.F @ 5385

Last change on this file since 5385 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: 3.6 KB
RevLine 
[630]1!
[1403]2! $Id: interpre.F 1403 2010-07-01 09:02:53Z abarral $
[630]3!
4       subroutine interpre(q,qppm,w,fluxwppm,masse,
5     s            apppm,bpppm,massebx,masseby,pbaru,pbarv,
6     s            unatppm,vnatppm,psppm)
7
[1403]8      USE control_mod
9      implicit none
[630]10
11#include "dimensions.h"
12c#include "paramr2.h"
13#include "paramet.h"
14#include "comconst.h"
15#include "comdissip.h"
16#include "comvert.h"
17#include "comgeom2.h"
18#include "logic.h"
19#include "temps.h"
20#include "ener.h"
21#include "description.h"
22
23c---------------------------------------------------
24c Arguments     
25      real   apppm(llm+1),bpppm(llm+1)
26      real   q(iip1,jjp1,llm),qppm(iim,jjp1,llm)
27c---------------------------------------------------
28      real   masse(iip1,jjp1,llm)
29      real   massebx(iip1,jjp1,llm),masseby(iip1,jjm,llm)     
30      real   w(iip1,jjp1,llm+1)
31      real   fluxwppm(iim,jjp1,llm)
32      real   pbaru(iip1,jjp1,llm )
33      real   pbarv(iip1,jjm,llm)
34      real   unatppm(iim,jjp1,llm)
35      real   vnatppm(iim,jjp1,llm)
36      real   psppm(iim,jjp1)
37c---------------------------------------------------
38c Local
39      real   vnat(iip1,jjp1,llm)
40      real   unat(iip1,jjp1,llm)
41      real   fluxw(iip1,jjp1,llm)
42      real   smass(iip1,jjp1)
43c----------------------------------------------------
44      integer l,ij,i,j
45
46c       CALCUL DE LA PRESSION DE SURFACE
47c       Les coefficients ap et bp sont passés en common
48c       Calcul de la pression au sol en mb optimisée pour
49c       la vectorialisation
50                   
51         do j=1,jjp1
52             do i=1,iip1
53                smass(i,j)=0.
54             enddo
55         enddo
56
57         do l=1,llm
58             do j=1,jjp1
59                 do i=1,iip1
60                    smass(i,j)=smass(i,j)+masse(i,j,l)
61                 enddo
62             enddo
63         enddo
64     
65         do j=1,jjp1
66             do i=1,iim
67                 psppm(i,j)=smass(i,j)/aire(i,j)*g*0.01
68             end do
69         end do                       
70       
71c RECONSTRUCTION DES CHAMPS CONTRAVARIANTS
72c Le programme ppm3d travaille avec les composantes
73c de vitesse et pas les flux, on doit donc passer de l'un à l'autre
74c Dans le même temps, on fait le changement d'orientation du vent en v
75      do l=1,llm
76          do j=1,jjm
77              do i=1,iip1
78                  vnat(i,j,l)=-pbarv(i,j,l)/masseby(i,j,l)*cv(i,j)             
79              enddo
80          enddo
81          do  i=1,iim
82          vnat(i,jjp1,l)=0.
83          enddo
84          do j=1,jjp1
85              do i=1,iip1
86                  unat(i,j,l)=pbaru(i,j,l)/massebx(i,j,l)*cu(i,j)
87              enddo
88          enddo
89      enddo
90             
91c CALCUL DU FLUX MASSIQUE VERTICAL
92c Flux en l=1 (sol) nul
93      fluxw=0.       
94      do l=1,llm
95           do j=1,jjp1
96              do i=1,iip1             
97               fluxw(i,j,l)=w(i,j,l)*g*0.01/aire(i,j)
98C               print*,i,j,l,'fluxw(i,j,l)=',fluxw(i,j,l),
99C     c                      'w(i,j,l)=',w(i,j,l)
100              enddo
101           enddo
102      enddo
103     
104c INVERSION DES NIVEAUX
105c le programme ppm3d travaille avec une 3ème coordonnée inversée par rapport
106c de celle du LMDZ: z=1<=>niveau max, z=llm+1<=>surface
107c On passe donc des niveaux du LMDZ à ceux de Lin
108     
109      do l=1,llm+1
110          apppm(l)=ap(llm+2-l)
111          bpppm(l)=bp(llm+2-l)         
112      enddo
113     
114      do l=1,llm
115          do j=1,jjp1
116             do i=1,iim     
117                 unatppm(i,j,l)=unat(i,j,llm-l+1)
118                 vnatppm(i,j,l)=vnat(i,j,llm-l+1)
119                 fluxwppm(i,j,l)=fluxw(i,j,llm-l+1)
120                 qppm(i,j,l)=q(i,j,llm-l+1)                             
121             enddo
122          enddo                               
123      enddo
124   
125      return
126      end
127
128
129
130
131
132
Note: See TracBrowser for help on using the repository browser.