New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 7333 – NEMO

Changeset 7333


Ignore:
Timestamp:
2016-11-25T10:05:40+01:00 (7 years ago)
Author:
deazer
Message:

Correct Precision issue with I/O

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/CO6_Restartable_Tidal_Analysis/NEMOGCM/NEMO/OPA_SRC/DIA/diaharmana.F90

    r6885 r7333  
    553553         write (clfinal,'(a,''_'',i4.4)') trim(clpath)//trim(clname),nproc 
    554554         open(66,file=TRIM(clfinal)) 
    555          write(66,'(1e15.5)') cc 
     555         write(66,'(1e20.9)') cc 
    556556         !These Three which contain the most data can be moved to a regular 
    557557         !restart file 
     
    571571         enddo 
    572572 
    573          write(66,'(1e15.5)') anau 
    574          write(66,'(1e15.5)') anav 
    575          write(66,'(1e15.5)') anaf 
     573         write(66,'(1e20.9)') anau 
     574         write(66,'(1e20.9)') anav 
     575         write(66,'(1e20.9)') anaf 
    576576         write(66,'(1e25.20)') fjulday_startharm 
    577577         close(66) 
     
    614614         write (clfinal,'(a,''_'',i4.4)') trim(clpath)//trim(clname),nproc 
    615615         open(66,file=TRIM(clfinal),status='old') 
    616          read(66,'(1e15.5)') cc 
     616         read(66,'(1e20.9)') cc 
    617617!2D regular fields can be read from normal restart this saves space and handy to 
    618618!view in netcdf format also. 
     
    632632            CALL iom_get(   numror,jpdom_autoglo, TRIM(Wave(ntide(ih))%cname_tide)//'bvbar_sin'   , bvbar(ih*2+1, : , : )    ) 
    633633         enddo 
    634          read(66,'(1e15.5)') anau 
    635          read(66,'(1e15.5)') anav 
    636          read(66,'(1e15.5)') anaf 
     634         read(66,'(1e20.9)') anau 
     635         read(66,'(1e20.9)') anav 
     636         read(66,'(1e20.9)') anaf 
    637637         read(66,'(1e25.20)') fjulday_startharm 
    638638      WRITE(numout,*) 'Checking anaf is correct' 
Note: See TracChangeset for help on using the changeset viewer.