Changeset 297 for codes

Timestamp:

11/27/14 15:01:16 (10 years ago)

Author:

millour

Message:

Make restart file format be netcdf4, so that output files of size greater than 2GB may be written.
EM

File:

• codes/icosagcm/branches/SATURN_DYNAMICO/ICOSAGCM/src/restart.f90 (modified) (4 diffs)

codes/icosagcm/branches/SATURN_DYNAMICO/ICOSAGCM/src/restart.f90

@@ -74 +74 @@
 
     IF (is_mpi_root) THEN
-      status = NF90_CREATE(TRIM(ADJUSTL(restart_file_name))//'.nc', NF90_CLOBBER, ncid)
+      status = NF90_CREATE(TRIM(ADJUSTL(restart_file_name))//'.nc', IOR(NF90_CLOBBER,NF90_NETCDF4), ncid)
+      IF (status/=NF90_NOERR) THEN
+        WRITE(*,*) "write_restart error: failed to create ",TRIM(ADJUSTL(restart_file_name))//'.nc'
+        WRITE(*,*) TRIM(NF90_STRERROR(status))
+      ELSE
+        WRITE(*,*) "write_restart: created ",TRIM(ADJUSTL(restart_file_name))//'.nc'
+      ENDIF
       status = NF90_DEF_DIM(ncid,'cell',ncell_glo,cellId)
       status = NF90_DEF_DIM(ncid,'edge',3*ncell_glo,edgeId)
@@ -119 +125 @@
           ENDIF
         ENDIF
+        IF (status/=NF90_NOERR) THEN
+          WRITE(*,*) "write_restart error: for ",TRIM(ADJUSTL(field(1)%name))
+          WRITE(*,*) TRIM(NF90_STRERROR(status))
+        ENDIF
       ENDDO
           
       
       status = NF90_ENDDEF(ncid)
+      IF (status/=NF90_NOERR) THEN
+         WRITE(*,*) "write_restart error: failed to switch out of netcdf define mode"
+         WRITE(*,*) TRIM(NF90_STRERROR(status))
+         STOP
+      ENDIF
       
       ALLOCATE(lon(ncell_glo),lat(ncell_glo),bounds_lon(0:nvert-1,ncell_glo),bounds_lat(0:nvert-1,ncell_glo))
@@ -147 +162 @@
       status=NF90_PUT_VAR(ncid,bounds_latId,REAL(bounds_lat,r8),start=(/ 1,1 /),count=(/ nvert,ncell_glo /))
       status=NF90_PUT_VAR(ncid,levAxisId,REAL(presnivs,r8),start=(/ 1 /),count=(/ llm /))
-    ENDIF
+    ENDIF ! of IF (is_mpi_root)
 
     DO nf=1,nfield
@@ -163 +178 @@
     IF (is_mpi_root) THEN
       status = NF90_CLOSE(ncid)      
+      IF (status/=NF90_NOERR) THEN
+         WRITE(*,*) "write_restart error: failed to close ",TRIM(ADJUSTL(restart_file_name))//'.nc'
+         WRITE(*,*) TRIM(NF90_STRERROR(status))
+         STOP
+      ELSE
+         WRITE(*,*) "write_restart: finished writing ",TRIM(ADJUSTL(restart_file_name))//'.nc'
+      ENDIF
     ENDIF
 !$OMP END MASTER