Changeset 4768 for branches/2014
- Timestamp:
- 2014-09-17T11:05:03+02:00 (10 years ago)
- Location:
- branches/2014/dev_r4650_UKMO14.5_SST_BIAS_CORRECTION/NEMOGCM
- Files:
-
- 1 added
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r4650_UKMO14.5_SST_BIAS_CORRECTION/NEMOGCM/CONFIG/SHARED/namelist_ref
r4384 r4768 1112 1112 ! sstfiles GHRSST input observation file name 1113 1113 ! ! sstfbfiles: Feedback SST input observation file name 1114 sstfbfiles = 'sst_01.nc' 'sst_02.nc' 'sst_03.nc' 'sst_04.nc' 'sst_05.nc'1114 sstfbfiles = 'sst_01.nc' 1115 1115 ! seaicefiles Sea Ice input observation file name 1116 1116 ! velavcurfiles Vel. cur. daily av. input file name … … 1132 1132 ln_grid_global = .true. 1133 1133 ln_grid_search_lookup = .false. 1134 ln_sstbias = .false. 1135 sstbias_files = 'sstbias.nc' 1134 1136 / 1135 1137 !----------------------------------------------------------------------- -
branches/2014/dev_r4650_UKMO14.5_SST_BIAS_CORRECTION/NEMOGCM/NEMO/OPA_SRC/OBS/diaobs.F90
r4624 r4768 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 89 90 LOGICAL, PUBLIC :: ln_sstbias !: Logical switch for bias corection of SST 91 90 92 REAL(KIND=dp), PUBLIC :: dobsini !: Observation window start date YYYYMMDD.HHMMSS 91 93 REAL(KIND=dp), PUBLIC :: dobsend !: Observation window end date YYYYMMDD.HHMMSS … … 102 104 ! !: If so use endailyavtypes 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? … … 135 139 !! ! 06-10 (A. Weaver) Cleaning and add controls 136 140 !! ! 07-03 (K. Mogensen) General handling of profiles 141 !! ! 14-08 (J.While) Incorporated SST bias correction 137 142 !!---------------------------------------------------------------------- 138 143 … … 144 149 CHARACTER(len=128) :: profbfiles(MaxNumFiles) 145 150 CHARACTER(len=128) :: sstfiles(MaxNumFiles) 146 CHARACTER(len=128) :: sstfbfiles(MaxNumFiles) 151 CHARACTER(len=128) :: sstfbfiles(MaxNumFiles) 152 CHARACTER(len=128) :: sstbias_files(MaxNumFiles) 147 153 CHARACTER(len=128) :: slafilesact(MaxNumFiles) 148 154 CHARACTER(len=128) :: slafilespas(MaxNumFiles) … … 180 186 & ln_velhradcp, velhradcpfiles, & 181 187 & ln_velfb, velfbfiles, ln_velfb_av, & 182 & ln_profb_enatim, ln_ignmis, ln_cl4 188 & ln_profb_enatim, ln_ignmis, ln_cl4, & 189 & ln_sstbias, sstbias_files 183 190 184 191 INTEGER :: jprofset … … 193 200 INTEGER :: jnumsst 194 201 INTEGER :: jnumsstfb 202 INTEGER :: jnumsstbias 195 203 INTEGER :: jnumseaice 196 204 INTEGER :: jnumvelavcur … … 208 216 !----------------------------------------------------------------------- 209 217 210 velcurfiles(:) = '' 218 !Initalise all values in namelist arrays 219 enactfiles(:) = '' 220 coriofiles(:) = '' 221 profbfiles(:) = '' 222 slafilesact(:) = '' 223 slafilespas(:) = '' 224 slafbfiles(:) = '' 225 sstfiles(:) = '' 226 sstfbfiles(:) = '' 227 seaicefiles(:) = '' 228 velcurfiles(:) = '' 229 veladcpfiles(:) = '' 230 velavcurfiles(:) = '' 231 velhrcurfiles(:) = '' 232 velavadcpfiles(:) = '' 233 velhradcpfiles(:) = '' 234 velfbfiles(:) = '' 235 velcurfiles(:) = '' 211 236 veladcpfiles(:) = '' 237 sstbias_files(:) = '' 238 endailyavtypes(:) = -1 239 endailyavtypes(1) = 820 240 ln_profb_ena(:) = .FALSE. 241 ln_profb_enatim(:) = .TRUE. 242 ln_velfb_av(:) = .FALSE. 243 ln_ignmis = .FALSE. 212 244 CALL ini_date( dobsini ) 213 245 CALL fin_date( dobsend ) … … 264 296 lmask(:) = .FALSE. 265 297 ENDIF 298 IF (ln_sstbias) THEN 299 lmask(:) = .FALSE. 300 WHERE (sstbias_files(:) /= '') lmask(:) = .TRUE. 301 jnumsstbias = COUNT(lmask) 302 lmask(:) = .FALSE. 303 ENDIF 266 304 IF (ln_seaice) THEN 267 305 lmask(:) = .FALSE. … … 315 353 WRITE(numout,*) ' Logical switch for GHRSST observations ln_ghrsst = ', ln_ghrsst 316 354 WRITE(numout,*) ' Logical switch for feedback SST data ln_sstfb = ', ln_sstfb 355 WRITE(numout,*) ' Logical switch for SST bias correction ln_sstbias = ', ln_sstbias 317 356 WRITE(numout,*) ' Logical switch for night-time SST obs ln_sstnight = ', ln_sstnight 318 357 WRITE(numout,*) ' Logical switch for SSS observations ln_sss = ', ln_sss … … 754 793 755 794 ENDIF 756 795 796 !Read in bias field and correct SST. 797 IF ( ln_sstbias ) THEN 798 IF ( jnumsstbias == 0 ) CALL ctl_stop("ln_sstbias set,"// & 799 " but no bias"// & 800 " files to read in") 801 CALL obs_app_sstbias( nsstsets, sstdatqc, n2dint, & 802 jnumsstbias, & 803 sstbias_files(1:jnumsstbias) ) 804 ENDIF 805 757 806 ENDIF 758 807
Note: See TracChangeset
for help on using the changeset viewer.