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/diaobs.F90

    r10181 r10276  
    2222   USE obs_read_surf            ! Reading and allocation of surface obs 
    2323   USE obs_readmdt              ! Reading and allocation of MDT for SLA. 
     24   USE obs_readsnowdepth        ! Get model snow depth for conversion of freeboard to ice thickness. 
    2425   USE obs_prep                 ! Preparation of obs. (grid search etc). 
    2526   USE obs_oper                 ! Observation operators 
     
    9899      & surfdata, &          !: Initial surface data 
    99100      & surfdataqc           !: Surface data after quality control 
     101       
    100102   TYPE(obs_prof), PUBLIC, POINTER, DIMENSION(:) :: & 
    101103      & profdata, &          !: Initial profile data 
     
    134136      !!        !  15-02  (M. Martin) Simplification of namelist and code 
    135137      !!---------------------------------------------------------------------- 
    136  
     138#if defined key_cice 
     139      USE sbc_oce, ONLY : &        ! CICE variables 
     140         & thick_s                 ! snow depth for freeboard conversion        
     141#endif          
    137142      IMPLICIT NONE 
    138143 
     
    831836            nextrsurf(jtype) = 0 
    832837            IF ( TRIM(cobstypessurf(jtype)) == 'sla' ) nextrsurf(jtype) = 2 
     838            IF ( TRIM(cobstypessurf(jtype)) == 'sit' ) nextrsurf(jtype) = 2             
    833839 
    834840            !Read in surface obs types 
     
    844850               IF ( ln_altbias ) & 
    845851                  & CALL obs_rea_altbias ( surfdataqc(jtype), n2dintsurf(jtype), cn_altbiasfile ) 
     852            ENDIF 
     853             
     854            IF ( TRIM(cobstypessurf(jtype)) == 'sit' ) THEN 
     855               CALL obs_rea_snowdepth( surfdataqc(jtype), n2dintsurf(jtype), thick_s(:,:) ) 
    846856            ENDIF 
    847857 
     
    12511261               ELSE             
    12521262#if defined key_cice 
    1253                ! Will need to convert freeboard to thickness here 
    12541263                  zsurfvar(:,:) = thick_i(:,:) 
    12551264#elif defined key_lim2 || defined key_lim3 
     
    15331542               &               nit000, idaystp, zsurfvar, zsurfmask,    & 
    15341543               &               n2dintsurf(jtype), llnightav(jtype),     & 
    1535                &               ravglamscl(jtype), ravgphiscl(jtype),     & 
     1544               &               ravglamscl(jtype), ravgphiscl(jtype),    & 
    15361545               &               lfpindegs(jtype) ) 
    15371546 
     1547            ! Change label of data from FBD ("freeboard") to SIT ("Sea Ice Thickness") 
     1548            IF ( TRIM(surfdataqc(jtype)%cvars(1)) == 'FBD' ) THEN 
     1549              surfdata(jtype)%cvars(1) = 'SIT'  
     1550            ENDIF 
     1551             
    15381552         END DO 
    15391553 
Note: See TracChangeset for help on using the changeset viewer.