Changeset 533 for LMDZ4/trunk/libf

Timestamp:

Jun 1, 2004, 4:29:52 PM (21 years ago)

Author:

lmdzadmin

Message:

Problème de phasage avec LMDZ.3.3 qui buggait create_etat0_limit
LF

File:

• LMDZ4/trunk/libf/dyn3d/startvar.F (modified) (3 diffs)

LMDZ4/trunk/libf/dyn3d/startvar.F

@@ -558 +558 @@
     !  LOCAL
     !
-      REAL :: lev(1), date, dt
+!ac     REAL :: lev(1), date, dt
+      REAL :: date, dt
+      REAL, DIMENSION(:), ALLOCATABLE :: levphys_ini
+!ac
       INTEGER :: itau(1)
       INTEGER ::  llm_tmp, ttm_tmp
       INTEGER :: i, j
-      INTEGER :: iret
     !
       CHARACTER*25 title
@@ -579 +581 @@
      . ttm_tmp, fid_phys)
     !
-    !
       ALLOCATE (lat_phys(iml_phys,jml_phys))
       ALLOCATE (lon_phys(iml_phys,jml_phys))
+!ac
+      ALLOCATE (levphys_ini(llm_tmp))
+    !
+!      CALL flinopen(physfname, .FALSE., iml_phys, jml_phys, 
+!     . llm_tmp, lon_phys, lat_phys, lev, ttm_tmp, 
+!     . itau, date, dt, fid_phys)
     !
       CALL flinopen(physfname, .FALSE., iml_phys, jml_phys, 
-     . llm_tmp, lon_phys, lat_phys, lev, ttm_tmp, 
+     . llm_tmp, lon_phys, lat_phys, levphys_ini, ttm_tmp, 
      . itau, date, dt, fid_phys)
     !
+      DEALLOCATE (levphys_ini)
+!ac
+    !
     ! Allocate the space we will need to get the data out of this file
     !
@@ -593 +603 @@
     !   In case we have a file which is in degrees we do the transformation
     !
-      if (allocated(lon_rad)) DEALLOCATE(lon_rad, stat=iret)
-      if (allocated(lon_ini)) DEALLOCATE(lon_ini, stat=iret)
-
       ALLOCATE(lon_rad(iml_phys))
       ALLOCATE(lon_ini(iml_phys))