Changeset 1403 for LMDZ4/trunk/libf/bibio/writedynav.F

Timestamp:

Jul 1, 2010, 11:02:53 AM (15 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/bibio/writedynav.F (modified) (9 diffs)

LMDZ4/trunk/libf/bibio/writedynav.F

@@ -2 +2 @@
 ! $Id$
 !
-      subroutine writedynav( histid, time, vcov, 
-     ,                          ucov,teta,ppk,phi,q,masse,ps,phis)
+      subroutine writedynav(time, vcov, 
+     ,                ucov,teta,ppk,phi,q,masse,ps,phis)
 
 #ifdef CPP_IOIPSL
@@ -9 +9 @@
 #endif
       USE infotrac, ONLY : nqtot, ttext
+      use com_io_dyn_mod, only : histaveid,histvaveid,histuaveid
       implicit none
 
@@ -17 +18 @@
 C
 C   Entree:
-C      histid: ID du fichier histoire
 C      time: temps de l'ecriture
 C      vcov: vents v covariants
@@ -29 +29 @@
 C      
 C
-C   Sortie:
-C      fileid: ID du fichier netcdf cree
 C
 C   L. Fairhead, LMD, 03/99
@@ -53 +51 @@
 C
 
-      INTEGER histid
       REAL vcov(ip1jm,llm),ucov(ip1jmp1,llm) 
-      REAL teta(ip1jmp1*llm),phi(ip1jmp1,llm),ppk(ip1jmp1*llm)                  
+      REAL teta(ip1jmp1*llm),phi(ip1jmp1,llm),ppk(ip1jmp1*llm)     
       REAL ps(ip1jmp1),masse(ip1jmp1,llm)                   
       REAL phis(ip1jmp1)                  
@@ -66 +63 @@
 C   Variables locales
 C
-      integer ndex2d(iip1*jjp1),ndex3d(iip1*jjp1*llm),iq, ii, ll
-      real us(ip1jmp1*llm), vs(ip1jmp1*llm)
+      integer ndex2d(ip1jmp1),ndexu(ip1jmp1*llm),ndexv(ip1jm*llm)
+      INTEGER iq, ii, ll
       real tm(ip1jmp1*llm)
       REAL vnat(ip1jm,llm),unat(ip1jmp1,llm) 
@@ -75 +72 @@
 C  Initialisations
 C
-      ndex3d = 0
+      ndexu = 0
+      ndexv = 0
       ndex2d = 0
       ok_sync = .TRUE.
-      us = 999.999
-      vs = 999.999
       tm = 999.999
       vnat = 999.999
@@ -91 +87 @@
 C  Appels a histwrite pour l'ecriture des variables a sauvegarder
 C
-C  Vents U scalaire
+C  Vents U 
 C
-      call gr_u_scal(llm, unat, us)
-      call histwrite(histid, 'u', itau_w, us, 
-     .               iip1*jjp1*llm, ndex3d)
+      call histwrite(histuaveid, 'u', itau_w, unat, 
+     .               iip1*jjp1*llm, ndexu)
 C
-C  Vents V scalaire
+C  Vents V
 C
-      call gr_v_scal(llm, vnat, vs)
-      call histwrite(histid, 'v', itau_w, vs, 
-     .               iip1*jjp1*llm, ndex3d)
+      call histwrite(histvaveid, 'v', itau_w, vnat, 
+     .               iip1*jjm*llm, ndexv)
 C
 C  Temperature potentielle moyennee
 C
-      call histwrite(histid, 'theta', itau_w, teta, 
-     .                iip1*jjp1*llm, ndex3d)
+      call histwrite(histaveid, 'theta', itau_w, teta, 
+     .                iip1*jjp1*llm, ndexu)
 C
 C  Temperature moyennee
@@ -113 +107 @@
         tm(ii) = teta(ii) * ppk(ii)/cpp
       enddo
-      call histwrite(histid, 'temp', itau_w, tm, 
-     .                iip1*jjp1*llm, ndex3d)
+      call histwrite(histaveid, 'temp', itau_w, tm, 
+     .                iip1*jjp1*llm, ndexu)
 C
 C  Geopotentiel
 C
-      call histwrite(histid, 'phi', itau_w, phi, 
-     .                iip1*jjp1*llm, ndex3d)
+      call histwrite(histaveid, 'phi', itau_w, phi, 
+     .                iip1*jjp1*llm, ndexu)
 C
 C  Traceurs
 C
-        DO iq=1,nqtot
-          call histwrite(histid, ttext(iq), itau_w, q(:,:,iq), 
-     .                   iip1*jjp1*llm, ndex3d)
-        enddo
+!        DO iq=1,nqtot
+!          call histwrite(histaveid, ttext(iq), itau_w, q(:,:,iq), 
+!     .                   iip1*jjp1*llm, ndexu)
+!        enddo
 C
 C  Masse
 C
-       call histwrite(histid, 'masse', itau_w, masse, iip1*jjp1, ndex2d)
+       call histwrite(histaveid, 'masse', itau_w, masse, 
+     $                   iip1*jjp1*llm, ndexu)
 C
 C  Pression au sol
 C
-       call histwrite(histid, 'ps', itau_w, ps, iip1*jjp1, ndex2d)
+       call histwrite(histaveid, 'ps', itau_w, ps, iip1*jjp1, ndex2d)
 C
 C  Geopotentiel au sol
 C
-       call histwrite(histid, 'phis', itau_w, phis, iip1*jjp1, ndex2d)
+!       call histwrite(histaveid,'phis',itau_w, phis,iip1*jjp1, ndex2d)
 C
 C  Fin
 C
-      if (ok_sync) call histsync(histid)
+      if (ok_sync) then
+          call histsync(histaveid)
+          call histsync(histvaveid)
+          call histsync(histuaveid)
+      ENDIF
 
 #else