Changeset 319 for codes/icosagcm/branches

Timestamp:

01/16/15 15:18:53 (9 years ago)

Author:

ymipsl

Message:

Iteration number is appened to the restart file name.

YM

File:

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

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

@@ -51 +51 @@
   REAL(rstd),ALLOCATABLE :: lon(:),lat(:),bounds_lon(:,:),bounds_lat(:,:)
   LOGICAL,SAVE :: first=.TRUE.
+  CHARACTER(LEN=16) :: str_it
 !$OMP THREADPRIVATE(first)
    
@@ -83 +84 @@
       
 
-    IF (is_mpi_root .AND. first) THEN
-      status = NF90_CREATE(TRIM(ADJUSTL(restart_file_name))//'.nc', IOR(NF90_CLOBBER,NF90_NETCDF4), ncid)
+    IF (is_mpi_root) THEN
+      WRITE(str_it,*) it
+      status = NF90_CREATE(TRIM(ADJUSTL(restart_file_name))//"_"//TRIM(ADJUSTL(str_it))//'.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(*,*) "write_restart error: failed to create ",TRIM(ADJUSTL(restart_file_name))//"_"//TRIM(ADJUSTL(str_it))//'.nc'
         WRITE(*,*) TRIM(NF90_STRERROR(status))
       ELSE
-        WRITE(*,*) "write_restart: created ",TRIM(ADJUSTL(restart_file_name))//'.nc'
+        WRITE(*,*) "write_restart: created ",TRIM(ADJUSTL(restart_file_name))//"_"//TRIM(ADJUSTL(str_it))//'.nc'
       ENDIF
       status = NF90_DEF_DIM(ncid,'cell',ncell_glo,cellId)
@@ -191 +193 @@
 
     IF (is_mpi_root) THEN
-      status = NF90_SYNC(ncid)      
+      status = NF90_CLOSE(ncid)      
       IF (status/=NF90_NOERR) THEN
-         WRITE(*,*) "write_restart error: failed to close ",TRIM(ADJUSTL(restart_file_name))//'.nc'
+         WRITE(*,*) "write_restart error: failed to close ",TRIM(ADJUSTL(restart_file_name))//"_"//TRIM(ADJUSTL(str_it))//'.nc'
          WRITE(*,*) TRIM(NF90_STRERROR(status))
          STOP
       ELSE
-         WRITE(*,*) "write_restart: finished writing ",TRIM(ADJUSTL(restart_file_name))//'.nc'
+         WRITE(*,*) "write_restart: finished writing ",TRIM(ADJUSTL(restart_file_name))//"_"//TRIM(ADJUSTL(str_it))//'.nc'
       ENDIF
     ENDIF