Ignore:
Timestamp:
2018-11-02T18:15:27+01:00 (3 years ago)
Author:
emmafiedler
Message:

Freeboard assimilation updates

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_obs_oper_update_icethick/NEMOGCM/NEMO/OPA_SRC/OBS/obs_oper.F90

    r9306 r10276  
    3737   USE obs_grid,      ONLY : &  
    3838      & obs_level_search      
    39  
     39#if defined key_cice 
     40   USE ice_constants, ONLY : &    ! For conversion from sea ice freeboard to thickness 
     41      & rhos, rhoi, rhow 
     42#endif  
    4043   IMPLICIT NONE 
    4144 
     
    798801         ENDIF 
    799802          
     803         IF ( TRIM(surfdataqc%cvars(1)) == 'FBD' ) THEN 
     804            ! Convert radar freeboard to true freeboard (add 1/4 snow depth; 1/4 based on ratio of speed of light in vacuum compared to snow (3.0e8 vs 2.4e8 m/s)) 
     805            surfdataqc%rext(jobs,1) = surfdataqc%robs(jobs,1)             
     806            surfdataqc%robs(jobs,1) = surfdataqc%rext(jobs,1) + 0.25*surfdataqc%rext(jobs,2) 
     807            ! Convert corrected freeboard to ice thickness following Tilling et al. (2016) 
     808            surfdataqc%robs(jobs,1) = (surfdataqc%robs(jobs,1)*rhow + surfdataqc%rext(jobs,2)*rhos)/ & 
     809                                      (rhow - rhoi)             
     810         ENDIF 
     811          
    800812         IF ( zext(1) == obfillflt ) THEN 
    801813            ! If the observation value is a fill value, set QC flag to bad 
Note: See TracChangeset for help on using the changeset viewer.