Changeset 15288 for NEMO/branches/UKMO/NEMO_4.0.4_generic_obs
- Timestamp:
- 2021-09-24T17:11:43+02:00 (3 years ago)
- Location:
- NEMO/branches/UKMO/NEMO_4.0.4_generic_obs/src/OCE/OBS
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/UKMO/NEMO_4.0.4_generic_obs/src/OCE/OBS/obs_read_surf.F90
r15240 r15288 606 606 607 607 ! QC flags 608 ! WHY IS THIS NOT A FUNCTION OF NUM VARS?609 608 surfdata%nqc(iobs) = inpfiles(jj)%ivqc(ji,jvar) 610 609 611 610 ! Observed value 612 611 surfdata%robs(iobs,jvar) = inpfiles(jj)%pob(1,ji,jvar) 613 614 ! THIS NEEDS SORTING615 ! ! Model and MDT is set to fbrmdi unless read from file616 ! IF ( ldmod ) THEN617 ! surfdata%rmod(iobs,jvar) = inpfiles(jj)%padd(1,ji,1,1)618 ! IF ( TRIM(surfdata%cvars(1)) == 'SLA' ) THEN619 ! surfdata%rext(iobs,1) = inpfiles(jj)%padd(1,ji,2,1)620 ! surfdata%rext(iobs,2) = inpfiles(jj)%pext(1,ji,1)621 ! ENDIF622 ! ELSE623 ! surfdata%rmod(iobs,jvar) = fbrmdi624 ! IF ( TRIM(surfdata%cvars(1)) == 'SLA' ) surfdata%rext(iobs,:) = fbrmdi625 ! ENDIF626 612 627 613 ! Additional variables -
NEMO/branches/UKMO/NEMO_4.0.4_generic_obs/src/OCE/OBS/obs_write.F90
r15285 r15288 364 364 ENDIF 365 365 END DO 366 !fbdata%cname(1) = surfdata%cvars(1)367 !fbdata%coblong(1) = 'Sea level anomaly'368 !fbdata%cobunit(1) = 'Metres'369 !fbdata%cextname(1) = 'MDT'370 !fbdata%cextlong(1) = 'Mean dynamic topography'371 !fbdata%cextunit(1) = 'Metres'372 !fbdata%caddname(2) = 'SSH'373 !fbdata%caddlong(2,1) = 'Model Sea surface height'374 !fbdata%caddunit(2,1) = 'Metres'375 !IF ( TRIM(surfdata%cvars(1)) == 'ICECONC' ) THEN376 ! fbdata%caddlong(1,1) = 'Model interpolated ICE'377 !ELSE378 ! fbdata%caddlong(1,1) = 'Model interpolated ' // TRIM(surfdata%cvars(1))379 !ENDIF380 366 381 367 WRITE(clfname, FMT="(A,'fb_fdbk_',I4.4,'.nc')") TRIM(clfiletype), nproc … … 434 420 DO jvar = 1, surfdata%nvar 435 421 fbdata%padd(1,jo,1,jvar) = surfdata%rmod(jo,jvar) 436 !IF ( TRIM(surfdata%cvars(1)) == 'SLA' ) fbdata%padd(1,jo,2,1) = surfdata%rext(jo,1)437 422 fbdata%pob(1,jo,jvar) = surfdata%robs(jo,jvar) 438 423 IF ( surfdata%nqc(jo) > 255 ) THEN … … 449 434 ENDIF 450 435 fbdata%iobsk(1,jo,jvar) = 0 451 !IF ( TRIM(surfdata%cvars(1)) == 'SLA' ) fbdata%pext(1,jo,1) = surfdata%rext(jo,2)452 436 IF (iadd > 0) THEN 453 437 DO ja = 1, iadd
Note: See TracChangeset
for help on using the changeset viewer.