Changeset 1217

Timestamp:

Apr 3, 2014, 12:52:48 PM (11 years ago)

Author:

milmd

Message:

Bug fix for when there is zero tracer.
MI

File:

• trunk/LMDZ.GENERIC/libf/phystd/iostart.F90 (modified) (1 diff)

trunk/LMDZ.GENERIC/libf/phystd/iostart.F90

@@ -528 +528 @@
       ENDIF
       
-      ierr=NF90_DEF_DIM(nid_restart,"number_of_advected_fields",nqtot,idim5)
-      IF (ierr/=NF90_NOERR) THEN
-        write(*,*)'open_restartphy: problem defining number_of_advected_fields dimension '
-        write(*,*)trim(nf90_strerror(ierr))
-        CALL ABORT
-      ENDIF
+      if (nqtot>0) then
+        ! only define a tracer dimension if there are tracers
+        ierr=NF90_DEF_DIM(nid_restart,"number_of_advected_fields",nqtot,idim5)
+        IF (ierr/=NF90_NOERR) THEN
+          write(*,*)'open_restartphy: problem defining number_of_advected_fields dimension '
+          write(*,*)trim(nf90_strerror(ierr))
+          CALL ABORT
+        ENDIF
+      endif
 
       ierr=NF90_DEF_DIM(nid_restart,"nlayer",klev,idim6)