Changeset 15650
- Timestamp:
- 2022-01-18T09:50:49+01:00 (3 years ago)
- Location:
- branches/UKMO/dev_r5518_obs_oper_WISA_osse/NEMOGCM/NEMO/OPA_SRC/OBS
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_obs_oper_WISA_osse/NEMOGCM/NEMO/OPA_SRC/OBS/diaobs.F90
r14591 r15650 959 959 & ltype_clim = .TRUE. 960 960 961 IF ( (TRIM(cobstypessurf(jtype)) == 'sla') .OR. & 962 & (TRIM(cobstypessurf(jtype)) == 'sit') .OR. & 961 IF ( TRIM(cobstypessurf(jtype)) == 'sla') THEN 962 ltype_clim = .TRUE. 963 nvarssurf(jtype) = 1 964 nextrsurf(jtype) = 2 965 ELSE IF ( (TRIM(cobstypessurf(jtype)) == 'sit') .OR. & 963 966 & (TRIM(cobstypessurf(jtype)) == 'fbd') ) THEN 964 967 nvarssurf(jtype) = 1 … … 1467 1470 CASE('sla') 1468 1471 zsurfvar(:,:,1) = sshn(:,:) 1472 IF ( kstp == 0 ) THEN 1473 IF ( lwp .AND. surfdataqc(jtype)%nsstpmpp(1) > 0 ) THEN 1474 CALL ctl_warn( 'Sea-ice not initialised on zeroth '// & 1475 & 'time-step but some obs are valid then.' ) 1476 WRITE(numout,*)surfdataqc(jtype)%nsstpmpp(1), & 1477 & ' SLA obs will not be checked against ice coverage' 1478 ENDIF 1479 zsurfclim(:,:,1) = 0.0 1480 ELSE 1481 ! Pass sea-ice concentration in through SLA climatology 1482 ! for use in rejecting SLA under ice. 1483 zsurfclim(:,:,1) = fr_i(:,:) 1484 ENDIF 1485 1486 1469 1487 CASE('sss') 1470 1488 zsurfvar(:,:,1) = tsn(:,:,1,jp_sal) -
branches/UKMO/dev_r5518_obs_oper_WISA_osse/NEMOGCM/NEMO/OPA_SRC/OBS/obs_oper.F90
r13393 r15650 954 954 IF ( surfdataqc%lclim ) surfdataqc%rclm(jobs,1) = zclm(1) 955 955 956 IF ( surfdataqc%lclim .AND. TRIM(surfdataqc%cvars(1)) == 'SLA' ) THEN 957 IF (surfdataqc%rclm(jobs,1) > 0.01) THEN 958 surfdataqc%nqc(jobs) = 4 959 END IF 960 ENDIF 961 956 962 IF ( zext(1) == obfillflt ) THEN 957 963 ! If the observation value is a fill value, set QC flag to bad -
branches/UKMO/dev_r5518_obs_oper_WISA_osse/NEMOGCM/NEMO/OPA_SRC/OBS/obs_write.F90
r14591 r15650 494 494 495 495 CALL alloc_obfbdata( fbdata, 1, surfdata%nsurf, 1, & 496 & 1 + iadd_mdt + iadd_std + iadd , &496 & 1 + iadd_mdt + iadd_std + iadd_clm + iadd, & 497 497 & 1 + iext, .TRUE. ) 498 498 … … 516 516 fbdata%cgrid(1) = 'T' 517 517 DO ja = 1, iadd 518 fbdata%caddname(1+iadd_mdt+iadd_std+ ja) = padd%cdname(ja)519 fbdata%caddlong(1+iadd_mdt+iadd_std+ ja,1) = padd%cdlong(ja,1)520 fbdata%caddunit(1+iadd_mdt+iadd_std+ ja,1) = padd%cdunit(ja,1)518 fbdata%caddname(1+iadd_mdt+iadd_std+iadd_clm+ja) = padd%cdname(ja) 519 fbdata%caddlong(1+iadd_mdt+iadd_std+iadd_clm+ja,1) = padd%cdlong(ja,1) 520 fbdata%caddunit(1+iadd_mdt+iadd_std+iadd_clm+ja,1) = padd%cdunit(ja,1) 521 521 END DO 522 522
Note: See TracChangeset
for help on using the changeset viewer.