Changeset 2139

Timestamp:

Nov 4, 2014, 12:21:52 AM (10 years ago)

Author:

fhourdin

Message:

Correction dimensions tableau en 1D pour la lecture des restart.
Bug fixing for 1D restarts.

File:

• LMDZ5/trunk/libf/phylmd/1DUTILS.h (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/1DUTILS.h

@@ -362 +362 @@
       CHARACTER*(*) fichnom
 !Al1 plev tronque pour .nc mais plev(klev+1):=0
-      real :: plev(klon,klev),play (klon,klev),phi(klon,klev)
+      real :: plev(klon,klev+1),play (klon,klev),phi(klon,klev)
       real :: presnivs(klon,klev)
       real :: ucov(klon,klev),vcov(klon,klev),temp(klon,klev)
@@ -445 +445 @@
 !  Lecture des champs
 !
-      plev(1,klev+1)=0.
-      CALL get_field("plev",plev,found)
-      IF (.NOT. found) PRINT*, modname//'Le champ <Plev> est absent'
       CALL get_field("play",play,found)
       IF (.NOT. found) PRINT*, modname//'Le champ <Play> est absent'
@@ -464 +461 @@
       CALL get_field("omega2",omega2,found)
       IF (.NOT. found) PRINT*, modname//'Le champ <omega2> est absent'
+      plev(1,klev+1)=0.
+      CALL get_field("plev",plev(:,1:klev),found)
+      IF (.NOT. found) PRINT*, modname//'Le champ <Plev> est absent'
 
       Do iq=1,nqtot