- Timestamp:
- 2013-09-20T18:21:40+02:00 (11 years ago)
- Location:
- branches/2013/dev_r3987_UKMO4_OBS/NEMOGCM/NEMO/OOO_SRC
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2013/dev_r3987_UKMO4_OBS/NEMOGCM/NEMO/OOO_SRC/off_netcdf.F90
r4030 r4031 119 119 WRITE(cf_end, FMT="(I4.4,'.nc')") knproc 120 120 cfname = trim(cf_stem)//trim(cf_end) 121 122 ! ##################### DEBUG ############################123 WRITE(*,*) 'Debug :', trim(cdate), trim(cmodel), trim(cname)124 WRITE(*,*) 'Debug :', trim(cf_end)125 WRITE(*,*) 'Debug :', trim(cfname)126 ! ##################### DEBUG ############################127 121 128 122 ! QC attribute string settings -
branches/2013/dev_r3987_UKMO4_OBS/NEMOGCM/NEMO/OOO_SRC/off_oper.F90
r4030 r4031 24 24 !! Initialise NEMO 25 25 CALL nemo_init 26 WRITE(*,*) "DEBUG: ln_cl4 = ", ln_cl427 26 !! Initialise Offline obs_oper data 28 27 CALL off_data_init( ln_cl4 ) … … 46 45 CALL dia_obs_init 47 46 END IF 48 WRITE(*,*) "DEBUG: Read model counterpart"49 47 !! Read next model counterpart 50 48 CALL off_rea_dri 51 WRITE(*,*) "DEBUG: Interpolate matchup"52 49 !! Interpolate to observation space 53 50 CALL off_oper_interp 54 WRITE(*,*) "DEBUG: Write matchup"55 51 !! Pipe to output files 56 52 CALL dia_obs_wri -
branches/2013/dev_r3987_UKMO4_OBS/NEMOGCM/NEMO/OOO_SRC/off_read.F90
r4030 r4031 46 46 ELSE IF (TRIM(cmatchname) == '') THEN 47 47 ! Feedback file settings 48 WRITE(*,*) "DEBUG: Feedback file read"49 48 CALL off_read_dailymean(TRIM(cdfilename), kindex) 50 49 CALL off_read_juld(TRIM(cdfilename), kindex, cl4_modjuld) … … 171 170 WRITE(*,*) "Opening :", TRIM(filename) 172 171 ! Open Netcdf file to find dimension id 173 istat = nf90_open( trim(filename),nf90_noclobber,ncid)172 istat = nf90_open(path=TRIM(filename), mode=nf90_nowrite, ncid=ncid) 174 173 IF ( istat /= nf90_noerr ) THEN 175 174 WRITE(*,*) "WARNING: Could not open ", trim(filename) 175 WRITE(*,*) "ERROR: ", nf90_strerror(istat) 176 176 ENDIF 177 177 istat = nf90_inq_dimid(ncid,'x',xdim) … … 194 194 start_n = (/ nimpp, njmpp, 1, ifcst /) 195 195 count_n = (/ nlci, nlcj, jpk, 1 /) 196 WRITE(*,*) "DEBUG: start_n ", start_n197 WRITE(*,*) "DEBUG: count_n ", count_n198 196 start_s = (/ nimpp, njmpp, ifcst /) 199 197 count_s = (/ nlci, nlcj, 1 /) … … 204 202 istat = nf90_get_att(ncid, varid, '_FillValue', fill_val) 205 203 istat = nf90_get_var(ncid, varid, temp_tn, start_n, count_n) 206 ! Diagnose reading problems207 WRITE(*,*) "DEBUG: Max. temp_tn = ", MAXVAL(temp_tn)208 WRITE(*,*) "DEBUG: Ave. temp_tn = ", SUM(temp_tn, temp_tn>0)/(MAX(1,COUNT(temp_tn>0)))209 WRITE(*,*) "DEBUG: Min. temp_tn = ", MINVAL(temp_tn)210 204 WHERE(temp_tn(:,:,:) == fill_val) temp_tn(:,:,:) = fbrmdi 211 205 … … 249 243 sshn(:,nlcj+1:) = 0 250 244 END IF 251 252 ! Diagnose reading problems253 WRITE(*,*) "DEBUG: Max. tsn = ", MAXVAL(tsn)254 WRITE(*,*) "DEBUG: Ave. tsn = ", SUM(tsn, tsn>0)/(MAX(1,COUNT(tsn>0)))255 WRITE(*,*) "DEBUG: Min. tsn = ", MINVAL(tsn)256 245 257 246 ! Deallocate arrays -
branches/2013/dev_r3987_UKMO4_OBS/NEMOGCM/NEMO/OOO_SRC/off_write.F90
r4030 r4031 52 52 ! Open netcdf file 53 53 CALL chkerr(nf90_open(trim(cdfilename), nf90_write, ncid), cpname, __LINE__ ) 54 WRITE(*, *) 'DEBUG: off_wri_extra_3d: ', TRIM(cdfilename)55 54 ! Write data 56 55 CALL chkerr(nf90_inq_varid(ncid,TRIM(cdvarname), varid),cpname, __LINE__ ) … … 88 87 ! Open netcdf file 89 88 CALL chkerr(nf90_open(trim(cdfilename), nf90_write, ncid), cpname, __LINE__ ) 90 WRITE(*, *) 'DEBUG: off_wri_extra_4d_index: ', TRIM(cdfilename), TRIM(cdvarname)91 89 ! Write data 92 90 CALL chkerr(nf90_inq_varid(ncid,TRIM(cdvarname), varid),cpname, __LINE__ ) … … 121 119 ! Open netcdf file 122 120 CALL chkerr(nf90_open(trim(cdfilename), nf90_write, ncid), cpname, __LINE__ ) 123 WRITE(*, *) 'DEBUG: off_wri_extra_4d: ', TRIM(cdfilename), TRIM(cdvarname)124 121 ! Write data 125 122 CALL chkerr(nf90_inq_varid(ncid,TRIM(cdvarname), varid),cpname, __LINE__ )
Note: See TracChangeset
for help on using the changeset viewer.