Changeset 3401 for LMDZ6/trunk


Ignore:
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:
2 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/phylmd/iostart.F90

    r3070 r3401  
    283283  SUBROUTINE open_restartphy(filename)
    284284  USE netcdf
    285   USE mod_phys_lmdz_para
    286   USE mod_grid_phy_lmdz
    287   USE dimphy
     285  USE mod_phys_lmdz_para, ONLY: is_master
     286  USE mod_grid_phy_lmdz, ONLY: klon_glo
     287  USE dimphy, ONLY: klev, klevp1
     288  USE print_control_mod, ONLY: lunout
    288289  IMPLICIT NONE
    289290    CHARACTER(LEN=*),INTENT(IN) :: filename
    290291    INTEGER                     :: ierr
    291292   
    292     IF (is_mpi_root .AND. is_omp_root) THEN
    293       ierr = NF90_CREATE(filename, NF90_CLOBBER, nid_restart)
     293    IF (is_master) THEN
     294      ierr = NF90_CREATE(filename, IOR(NF90_CLOBBER,NF90_64BIT_OFFSET), &
     295                          nid_restart)
    294296      IF (ierr/=NF90_NOERR) THEN
    295         write(6,*)' Pb d''ouverture du fichier '//filename
    296         write(6,*)' ierr = ', ierr
    297         CALL abort_physic("", "", 1)
     297        write(lunout,*)'open_restartphy: problem creating file '//trim(filename)
     298        write(lunout,*)trim(nf90_strerror(ierr))
     299        CALL abort_physic("open_restartphy", trim(nf90_strerror(ierr)), 1)
    298300      ENDIF
    299301
  • LMDZ6/trunk/libf/phylmd/phyaqua_mod.F90

    r2979 r3401  
    526526      PRINT *, 'writelim: Ecriture du fichier limit'
    527527
    528       ierr = nf_create('limit.nc', nf_clobber, nid)
     528      ierr = nf_create('limit.nc', IOR(NF_CLOBBER,NF_64BIT_OFFSET), nid)
    529529
    530530      ierr = nf_put_att_text(nid, nf_global, 'title', 30, &
Note: See TracChangeset for help on using the changeset viewer.