Changeset 1531 for trunk/LMDZ.GENERIC/libf/phystd/inistats.F

Timestamp:

Apr 7, 2016, 9:17:30 AM (9 years ago)

Author:

mturbet

Message:

Generic GCM:

• Fix buggy ouputs in 1D introduced by previous code modifications.

EM+MT

File:

• trunk/LMDZ.GENERIC/libf/phystd/inistats.F (modified) (3 diffs)

trunk/LMDZ.GENERIC/libf/phystd/inistats.F

@@ -16 +16 @@
       integer :: l,nsteppd
       real, dimension(nbp_lev) ::  sig_s
-      real :: lon_reg_ext(nbp_lon+1) ! extended longitudes
+      real,allocatable :: lon_reg_ext(:) ! extended longitudes
       integer :: idim_lat,idim_lon,idim_llm,idim_llmp1,idim_time
       real, dimension(istime) :: lt
       integer :: nvarid
 
+
+      IF (nbp_lon*nbp_lat==1) THEN
+        ! 1D model
+        ALLOCATE(lon_reg_ext(1))
+      ELSE
+        ! 3D model
+        ALLOCATE(lon_reg_ext(nbp_lon+1))
+      ENDIF
+      
       write (*,*) 
       write (*,*) '                        || STATS ||'
@@ -46 +55 @@
       
       lon_reg_ext(1:nbp_lon)=lon_reg(1:nbp_lon)
-      !add extra redundant point (180 degrees, since lon_reg starts at -180
-      lon_reg_ext(nbp_lon+1)=-lon_reg_ext(1)
+      IF (nbp_lon*nbp_lat/=1) THEN
+        ! In 3D, add extra redundant point (180 degrees,
+        ! since lon_reg starts at -180)
+        lon_reg_ext(nbp_lon+1)=-lon_reg_ext(1)
+      ENDIF
 
       if (is_master) then
@@ -59 +71 @@
 
       ierr = NF_DEF_DIM (nid, "latitude", nbp_lat, idim_lat)
-      ierr = NF_DEF_DIM (nid, "longitude", nbp_lon+1, idim_lon)
+      IF (nbp_lon*nbp_lat==1) THEN
+        ierr = NF_DEF_DIM (nid, "longitude", 1, idim_lon)
+      ELSE
+        ierr = NF_DEF_DIM (nid, "longitude", nbp_lon+1, idim_lon)
+      ENDIF
       ierr = NF_DEF_DIM (nid, "altitude", nbp_lev, idim_llm)
       ierr = NF_DEF_DIM (nid, "llmp1", nbp_lev+1, idim_llmp1)