- Timestamp:
- 2021-11-09T16:27:42+01:00 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/UKMO/NEMO_4.0.4_generic_obs/src/OCE/OBS/obs_bias.F90
r15187 r15487 115 115 ENDIF 116 116 117 IF (lwp)WRITE(numout,*)118 IF (lwp)WRITE(numout,*) 'obs_app_bias : '119 IF (lwp)WRITE(numout,*) '----------------- '117 IF (lwp) WRITE(numout,*) 118 IF (lwp) WRITE(numout,*) 'obs_app_bias : ' 119 IF (lwp) WRITE(numout,*) '----------------- ' 120 120 IF ( ll_extvar ) THEN 121 IF (lwp)WRITE(numout,*) 'Read observation bias for ', TRIM(obsdata%cextvars(kvar))121 IF (lwp) WRITE(numout,*) 'Read observation bias for ', TRIM(obsdata%cextvars(kvar)) 122 122 ELSE 123 IF (lwp)WRITE(numout,*) 'Read observation bias for ', TRIM(obsdata%cvars(kvar))123 IF (lwp) WRITE(numout,*) 'Read observation bias for ', TRIM(obsdata%cvars(kvar)) 124 124 ENDIF 125 125 126 126 ! Open and read the files 127 z_obsbias(:,:,:) =0.0_wp127 z_obsbias(:,:,:) = 0.0_wp 128 128 DO jtype = 1, knumtypes 129 129 130 numobsbias =0131 IF (lwp)WRITE(numout,*) 'Opening ',cl_bias_files(jtype)130 numobsbias = 0 131 IF (lwp) WRITE(numout,*) 'Opening ', cl_bias_files(jtype) 132 132 CALL iom_open( cl_bias_files(jtype), numobsbias, ldstop=.FALSE. ) 133 133 IF (numobsbias > 0) THEN … … 138 138 !so have to use NETCDF 139 139 !routines directly - should be upgraded in the future 140 iret =NF90_OPEN(TRIM(cl_bias_files(jtype)), NF90_NOWRITE, incfile)140 iret = NF90_OPEN(TRIM(cl_bias_files(jtype)), NF90_NOWRITE, incfile) 141 141 IF ( .NOT. ll_extvar ) THEN 142 142 iret=NF90_GET_ATT( incfile, NF90_GLOBAL, TRIM(obsdata%cvars(kvar))//"_source", & … … 144 144 ibiastypes(jtype) = ifile_source 145 145 ENDIF 146 iret =NF90_CLOSE(incfile)146 iret = NF90_CLOSE(incfile) 147 147 IF ( iret /= 0 ) CALL ctl_stop( & 148 148 'obs_app_bias : Cannot read bias type from file '// & … … 196 196 & zmask_tmp(2,2,inumtype), & 197 197 & zbias( 2,2,inumtype ) ) 198 jt =1198 jt = 1 199 199 DO jobs = 1, obsdata%nsurf 200 200 IF ( obsdata%ntyp(jobs) == ibiastypes(jtype) ) THEN … … 211 211 & igrdi_tmp(:,:,:), igrdj_tmp(:,:,:), & 212 212 & z_obsbias(:,:,jtype), zbias(:,:,:) ) 213 jt =1213 jt = 1 214 214 DO jobs = 1, obsdata%nsurf 215 215 IF ( ( obsdata%ntyp(jobs) == ibiastypes(jtype) ) .OR. & … … 230 230 obsdata%robs(jobs,kvar) = obsdata%robs(jobs,kvar) - zext(1) 231 231 ENDIF 232 jt =jt+1232 jt = jt + 1 233 233 ENDIF 234 234 END DO … … 254 254 WRITE(numout,*) "Bias correction applied successfully" 255 255 IF ( .NOT. ll_extvar ) THEN 256 WRITE(numout,*) "Obs types: ", ibiastypes(:), &257 256 WRITE(numout,*) "Obs types: ", ibiastypes(:), & 257 & " Have all been bias corrected" 258 258 ENDIF 259 259 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.