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/SBC/sbcice_cice.F90

    r10181 r10276  
    251251      CALL lbc_lnk ( fr_iv , 'V', 1. ) 
    252252       
    253 !!*****JUST ADDED BUT NEEDS CONFIRMING!!**** 
     253 
    254254! Snow and ice thickness 
    255       CALL cice2nemo(vice,thick_i,'T', 1. )  !!Do this for snow depth too? set thick_s. Duplicated below for embedded ice? 
    256  
    257 ! vice is volume per unit area of grid cell = thickness 
     255! vice, vsno is volume per unit area of grid cell = thickness 
     256      CALL cice2nemo(vice,thick_i,'T', 1. ) 
     257      CALL cice2nemo(vsno,thick_s,'T', 1. ) 
    258258       
    259       IF ( (ksbc == jp_flx) .OR. (ksbc == jp_purecpl) ) THEN ! Confirm what this if loop is for 
     259      IF ( (ksbc == jp_flx) .OR. (ksbc == jp_purecpl) ) THEN 
    260260         DO jl = 1,ncat 
    261261            CALL cice2nemo(vsnon(:,:,jl,:),ht_s(:,:,jl),'T', 1. ) 
     
    263263         ENDDO 
    264264      ENDIF 
    265        
     265 
     266 
    266267! T point to U point 
    267268! T point to V point 
     269 
     270! Sea ice thickness 
    268271      thick_iu(:,:)=0.0 
    269272      thick_iv(:,:)=0.0 
     
    278281      CALL lbc_lnk ( thick_iv , 'V', 1. ) 
    279282 
    280 !!**********************************************************      
     283! Snow depth 
     284      thick_su(:,:)=0.0 
     285      thick_sv(:,:)=0.0 
     286      DO jj=1,jpjm1 
     287         DO ji=1,jpim1 
     288            thick_su(ji,jj)=0.5*(thick_s(ji,jj)+thick_s(ji+1,jj))*umask(ji,jj,1) 
     289            thick_sv(ji,jj)=0.5*(thick_s(ji,jj)+thick_s(ji,jj+1))*vmask(ji,jj,1) 
     290         ENDDO 
     291      ENDDO 
     292 
     293      CALL lbc_lnk ( thick_su , 'U', 1. ) 
     294      CALL lbc_lnk ( thick_sv , 'V', 1. )       
    281295       
    282296 
Note: See TracChangeset for help on using the changeset viewer.