Changeset 6326
- Timestamp:
- 2016-02-18T16:10:11+01:00 (8 years ago)
- Location:
- branches/UKMO/dev_r4650_general_vert_coord_obsoper/NEMOGCM
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r4650_general_vert_coord_obsoper/NEMOGCM/CONFIG/SHARED/namelist_ref
r5837 r6326 1253 1253 ln_grid_global = .true. 1254 1254 ln_grid_search_lookup = .false. 1255 ln_sstbias = .false. 1256 sstbias_files = 'sstbias.nc' 1255 1257 / 1256 1258 !----------------------------------------------------------------------- -
branches/UKMO/dev_r4650_general_vert_coord_obsoper/NEMOGCM/NEMO/OPA_SRC/OBS/diaobs.F90
r6016 r6326 25 25 USE obs_read_sla ! Reading and allocation of SLA observations 26 26 USE obs_read_sst ! Reading and allocation of SST observations 27 USE obs_sstbias ! Bias correction routine for SST 27 28 USE obs_readmdt ! Reading and allocation of MDT for SLA. 28 29 USE obs_read_seaice ! Reading and allocation of Sea Ice observations … … 87 88 LOGICAL, PUBLIC :: ln_ignmis !: Logical switch for ignoring missing files 88 89 LOGICAL, PUBLIC :: ln_s_at_t !: Logical switch to compute model S at T observations 90 LOGICAL, PUBLIC :: ln_sstbias !: Logical switch for bias corection of SST 89 91 90 92 REAL(KIND=dp), PUBLIC :: dobsini !: Observation window start date YYYYMMDD.HHMMSS … … 103 105 & ln_profb_enatim !: Change tim for 820 enact data set. 104 106 107 INTEGER, DIMENSION(MaxNumFiles), PUBLIC :: sstbias_type !SST bias type 108 105 109 LOGICAL, DIMENSION(MaxNumFiles) :: & 106 110 & ln_velfb_av !: Is the velocity feedback files daily average? … … 137 141 !! ! 06-10 (A. Weaver) Cleaning and add controls 138 142 !! ! 07-03 (K. Mogensen) General handling of profiles 143 !! ! 14-08 (J.While) Incorporated SST bias correction 139 144 !!---------------------------------------------------------------------- 140 145 … … 146 151 CHARACTER(len=128) :: profbfiles(MaxNumFiles) 147 152 CHARACTER(len=128) :: sstfiles(MaxNumFiles) 148 CHARACTER(len=128) :: sstfbfiles(MaxNumFiles) 153 CHARACTER(len=128) :: sstfbfiles(MaxNumFiles) 154 CHARACTER(len=128) :: sstbias_files(MaxNumFiles) 149 155 CHARACTER(len=128) :: slafilesact(MaxNumFiles) 150 156 CHARACTER(len=128) :: slafilespas(MaxNumFiles) … … 182 188 & ln_velhradcp, velhradcpfiles, & 183 189 & ln_velfb, velfbfiles, ln_velfb_av, & 184 & ln_profb_enatim, ln_ignmis, ln_cl4 190 & ln_profb_enatim, ln_ignmis, ln_cl4, & 191 & ln_sstbias, sstbias_files 185 192 186 193 INTEGER :: jprofset … … 195 202 INTEGER :: jnumsst 196 203 INTEGER :: jnumsstfb 204 INTEGER :: jnumsstbias 197 205 INTEGER :: jnumseaice 198 206 INTEGER :: jnumvelavcur … … 229 237 velcurfiles(:) = '' 230 238 veladcpfiles(:) = '' 239 sstbias_files(:) = '' 231 240 endailyavtypes(:) = -1 232 241 endailyavtypes(1) = 820 … … 290 299 lmask(:) = .FALSE. 291 300 ENDIF 301 IF (ln_sstbias) THEN 302 lmask(:) = .FALSE. 303 WHERE (sstbias_files(:) /= '') lmask(:) = .TRUE. 304 jnumsstbias = COUNT(lmask) 305 lmask(:) = .FALSE. 306 ENDIF 292 307 IF (ln_seaice) THEN 293 308 lmask(:) = .FALSE. … … 341 356 WRITE(numout,*) ' Logical switch for GHRSST observations ln_ghrsst = ', ln_ghrsst 342 357 WRITE(numout,*) ' Logical switch for feedback SST data ln_sstfb = ', ln_sstfb 358 WRITE(numout,*) ' Logical switch for SST bias correction ln_sstbias = ', ln_sstbias 343 359 WRITE(numout,*) ' Logical switch for night-time SST obs ln_sstnight = ', ln_sstnight 344 360 WRITE(numout,*) ' Logical switch for SSS observations ln_sss = ', ln_sss … … 776 792 777 793 ENDIF 794 795 !Read in bias field and correct SST. 796 IF ( ln_sstbias ) THEN 797 IF ( jnumsstbias == 0 ) CALL ctl_stop("ln_sstbias set,"// & 798 " but no bias"// & 799 " files to read in") 800 CALL obs_app_sstbias( nsstsets, sstdatqc, n2dint, & 801 jnumsstbias, & 802 sstbias_files(1:jnumsstbias) ) 803 ENDIF 778 804 779 805 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.