Changeset 4990 for trunk/NEMOGCM/NEMO/OPA_SRC/OBS/obs_read_sla.F90
- Timestamp:
- 2014-12-15T17:42:49+01:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/OBS/obs_read_sla.F90
r3651 r4990 111 111 & zphi, & 112 112 & zlam 113 REAL(dp), DIMENSION(:), ALLOCATABLE :: &113 real(wp), DIMENSION(:), ALLOCATABLE :: & 114 114 & zdat 115 115 LOGICAL :: llvalprof … … 117 117 TYPE(obfbdata), POINTER, DIMENSION(:) :: & 118 118 & inpfiles 119 REAL(dp), DIMENSION(knumfiles) :: &119 real(wp), DIMENSION(knumfiles) :: & 120 120 & djulini, & 121 121 & djulend … … 125 125 INTEGER :: iobs 126 126 INTEGER :: iobstot 127 INTEGER :: ios 128 INTEGER :: ioserrcount 127 129 CHARACTER(len=8) :: cl_refdate 128 130 … … 391 393 & iindx ) 392 394 393 CALL obs_surf_alloc( sladata, iobs, kvars, kextr, & 394 & jpi, jpj, kstp ) 395 CALL obs_surf_alloc( sladata, iobs, kvars, kextr, kstp, jpi, jpj ) 395 396 396 397 ! * Read obs/positions, QC, all variable and assign to sladata … … 400 401 ityp (:) = 0 401 402 itypmpp(:) = 0 403 404 ioserrcount = 0 402 405 403 406 DO jk = 1, iobstot … … 451 454 452 455 ! Instrument type 453 READ( inpfiles(jj)%cdtyp(ji), '(I4)' ) itype 456 READ( inpfiles(jj)%cdtyp(ji), '(I4)', IOSTAT = ios, ERR = 901 ) itype 457 901 IF ( ios /= 0 ) THEN 458 IF (ioserrcount == 0) CALL ctl_warn ( 'Problem converting an instrument type to integer. Setting type to zero' ) 459 ioserrcount = ioserrcount + 1 460 itype = 0 461 ENDIF 454 462 sladata%ntyp(iobs) = itype 455 463 ityp(itype+1) = ityp(itype+1) + 1
Note: See TracChangeset
for help on using the changeset viewer.