Changeset 3401 for LMDZ6/trunk

Timestamp:

Oct 16, 2018, 10:29:07 AM (6 years ago)

Author:

Ehouarn Millour

Message:

Make produced "limit.nc" and "restartphy.nc" files be of "64-bit offset" type rather than "classic" NetCDF format to enable creation of large files (>2Gb).
EM

Location:

LMDZ6/trunk/libf/phylmd

Files:

• iostart.F90 (modified) (1 diff)
• phyaqua_mod.F90 (modified) (1 diff)

LMDZ6/trunk/libf/phylmd/iostart.F90

@@ -283 +283 @@
   SUBROUTINE open_restartphy(filename)
   USE netcdf
-  USE mod_phys_lmdz_para
-  USE mod_grid_phy_lmdz
-  USE dimphy
+  USE mod_phys_lmdz_para, ONLY: is_master
+  USE mod_grid_phy_lmdz, ONLY: klon_glo
+  USE dimphy, ONLY: klev, klevp1
+  USE print_control_mod, ONLY: lunout
   IMPLICIT NONE
     CHARACTER(LEN=*),INTENT(IN) :: filename
     INTEGER                     :: ierr
     
-    IF (is_mpi_root .AND. is_omp_root) THEN
-      ierr = NF90_CREATE(filename, NF90_CLOBBER, nid_restart)
+    IF (is_master) THEN
+      ierr = NF90_CREATE(filename, IOR(NF90_CLOBBER,NF90_64BIT_OFFSET), &
+                          nid_restart)
       IF (ierr/=NF90_NOERR) THEN
-        write(6,*)' Pb d''ouverture du fichier '//filename
-        write(6,*)' ierr = ', ierr
-        CALL abort_physic("", "", 1)
+        write(lunout,*)'open_restartphy: problem creating file '//trim(filename)
+        write(lunout,*)trim(nf90_strerror(ierr))
+        CALL abort_physic("open_restartphy", trim(nf90_strerror(ierr)), 1)
       ENDIF
 

LMDZ6/trunk/libf/phylmd/phyaqua_mod.F90

@@ -526 +526 @@
       PRINT *, 'writelim: Ecriture du fichier limit'
 
-      ierr = nf_create('limit.nc', nf_clobber, nid)
+      ierr = nf_create('limit.nc', IOR(NF_CLOBBER,NF_64BIT_OFFSET), nid)
 
       ierr = nf_put_att_text(nid, nf_global, 'title', 30, &