Changeset 1451 for LMDZ4

Timestamp:

Nov 12, 2010, 4:13:39 PM (13 years ago)

Author:

jghattas

Message:

Bug corrections :
dyn3dpar/initfluxsto_p.F, leapfrog_p.F and fluxstokenc_p.F : for offline option
phylmd/calcul_divers.h : corrected initialization, for case ecrit_mth=dtime

Location:

LMDZ4/branches/LMDZ4_AR5/libf

Files:

• dyn3dpar/fluxstokenc_p.F (modified) (1 diff)
• dyn3dpar/initfluxsto_p.F (modified) (2 diffs)
• dyn3dpar/leapfrog_p.F (modified) (1 diff)
• phylmd/calcul_divers.h (modified) (3 diffs)

LMDZ4/branches/LMDZ4_AR5/libf/dyn3dpar/fluxstokenc_p.F

@@ -57 +57 @@
 
 c AC initialisations
-cym      pbarug(:,:)   = 0.
+      pbarug(:,:)   = 0.
 cym      pbarvg(:,:,:) = 0.
 cym      wg(:,:)       = 0.

LMDZ4/branches/LMDZ4_AR5/libf/dyn3dpar/initfluxsto_p.F

@@ -203 +203 @@
      .              llm, nivsigs, zvertiid)
 c pour le fichier def
-      nivd(1) = 1
-      call histvert(filedid, 'sig_s', 'Niveaux sigma',
-     .  'sigma_level',
-     .              1, nivd, dvertiid)
-
+      if (mpi_rank==0) then 
+         nivd(1) = 1
+         call histvert(filedid, 'sig_s', 'Niveaux sigma',
+     .        'sigma_level',
+     .        1, nivd, dvertiid)
+      endif
 C
 C  Appels a histdef pour la definition des variables a sauvegarder
@@ -282 +283 @@
       call histend(fileid)
       call histend(filevid)
-      call histend(filedid)
+      if (mpi_rank==0) call histend(filedid)
       if (ok_sync) then
         call histsync(fileid)
         call histsync(filevid)
-        call histsync(filedid)
+        if (mpi_rank==0) call histsync(filedid)
       endif
         

LMDZ4/branches/LMDZ4_AR5/libf/dyn3dpar/leapfrog_p.F

@@ -588 +588 @@
      .        flxw,pk, iapptrac)
 
-       IF (offline) THEN
+       IF (offline .AND. .NOT. adjust) THEN
 Cmaf stokage du flux de masse pour traceurs OFF-LINE
 

LMDZ4/branches/LMDZ4_AR5/libf/phylmd/calcul_divers.h

@@ -3 +3 @@
 c
 c
-c initialisations diverses au "debut" du mois
+c initialisations diverses au "debut" du simulation
 c
-      IF(MOD(itap,INT(ecrit_mth/dtime)).EQ.1) THEN
+      IF(debut) THEN
          DO i=1, klon
           nday_rain(i)=0.
@@ -17 +17 @@
        ENDDO
 c
-      ENDIF !MOD(itap,INT(ecrit_mth)).EQ.1
+      ENDIF
 c
-      IF(MOD(itap,INT(ecrit_day/dtime)).EQ.0) THEN
+      IF(MOD(itap,INT(un_jour/dtime)).EQ.0) THEN
 c
-cIM calcul total_rain, nday_rain
+cIM calcul une fois par jour : total_rain, nday_rain
 c
        DO i = 1, klon
@@ -27 +27 @@
         IF(total_rain(i).GT.0.) nday_rain(i)=nday_rain(i)+1.
        ENDDO
-      ENDIF !itap.EQ.ecrit_mth
+      ENDIF