Changeset 11863 for branches/UKMO/dev_r5518_obs_oper_update_SlaAssim/NEMOGCM/NEMO/OPA_SRC/OBS/obs_read_prof.F90
- Timestamp:
- 2019-11-05T18:22:07+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_obs_oper_update_SlaAssim/NEMOGCM/NEMO/OPA_SRC/OBS/obs_read_prof.F90
r11468 r11863 46 46 & kvars, kextr, kstp, ddobsini, ddobsend, & 47 47 & ldvar, ldignmis, ldsatt, & 48 & ldmod, ldclim, kdailyavtypes )48 & ldmod, ldclim, cdvars, kdailyavtypes ) 49 49 !!--------------------------------------------------------------------- 50 50 !! … … 81 81 REAL(dp), INTENT(IN) :: ddobsini ! Obs. ini time in YYYYMMDD.HHMMSS 82 82 REAL(dp), INTENT(IN) :: ddobsend ! Obs. end time in YYYYMMDD.HHMMSS 83 CHARACTER(len=8), DIMENSION(kvars), INTENT(IN) :: cdvars 83 84 INTEGER, DIMENSION(imaxavtypes), OPTIONAL :: & 84 85 & kdailyavtypes ! Types of daily average observations … … 87 88 CHARACTER(LEN=15), PARAMETER :: cpname='obs_rea_prof' 88 89 CHARACTER(len=8) :: clrefdate 89 CHARACTER(len=8), DIMENSION(:), ALLOCATABLE :: clvars 90 CHARACTER(len=8), DIMENSION(:), ALLOCATABLE :: clvarsin 90 91 INTEGER :: jvar 91 92 INTEGER :: ji … … 222 223 223 224 IF ( jj == 1 ) THEN 224 ALLOCATE( clvars ( inpfiles(jj)%nvar ) )225 ALLOCATE( clvarsin( inpfiles(jj)%nvar ) ) 225 226 DO ji = 1, inpfiles(jj)%nvar 226 clvars(ji) = inpfiles(jj)%cname(ji) 227 clvarsin(ji) = inpfiles(jj)%cname(ji) 228 IF ( clvarsin(ji) /= cdvars(ji) ) THEN 229 CALL ctl_stop( 'Feedback file variables do not match', & 230 & ' expected variable names for this type' ) 231 ENDIF 227 232 END DO 228 233 ELSE 229 234 DO ji = 1, inpfiles(jj)%nvar 230 IF ( inpfiles(jj)%cname(ji) /= clvars (ji) ) THEN235 IF ( inpfiles(jj)%cname(ji) /= clvarsin(ji) ) THEN 231 236 CALL ctl_stop( 'Feedback file variables not consistent', & 232 237 & ' with previous files for this type' ) … … 507 512 profdata%nprof = 0 508 513 profdata%nvprot(:) = 0 509 profdata%cvars(:) = clvars (:)514 profdata%cvars(:) = clvarsin(:) 510 515 iprof = 0 511 516 … … 809 814 ! Deallocate temporary data 810 815 !----------------------------------------------------------------------- 811 DEALLOCATE( ifileidx, iprofidx, zdat, clvars )816 DEALLOCATE( ifileidx, iprofidx, zdat, clvarsin ) 812 817 813 818 !-----------------------------------------------------------------------
Note: See TracChangeset
for help on using the changeset viewer.