Ignore:
Timestamp:
2019-11-20T17:32:11+01:00 (10 months ago)
Author:
dcarneir
Message:

Changing GO6 package to include sea ice thickness DA

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_GO6_package_FOAMv14_sit/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90

    r10302 r11939  
    5757                flatn_f,fsurfn_f,fcondtopn_f,                    & 
    5858#ifdef key_asminc 
    59                 daice_da,fresh_da,fsalt_da,                    & 
     59                daice_da,dhi_da,fresh_da,fsalt_da,               & 
    6060#endif 
    6161                uatm,vatm,wind,fsw,flw,Tair,potT,Qa,rhoa,zlvl,   & 
     
    250250      CALL lbc_lnk ( fr_iu , 'U', 1. ) 
    251251      CALL lbc_lnk ( fr_iv , 'V', 1. ) 
     252 
     253! Snow and ice thickness 
     254! vice, vsno is volume per unit area of grid cell = thickness 
     255      CALL cice2nemo(vice,thick_i,'T', 1. ) 
     256      CALL cice2nemo(vsno,thick_s,'T', 1. ) 
     257      
     258      IF ( (ksbc == jp_flx) .OR. (ksbc == jp_purecpl) ) THEN 
     259         DO jl = 1,ncat 
     260            CALL cice2nemo(vsnon(:,:,jl,:),ht_s(:,:,jl),'T', 1. ) 
     261            CALL cice2nemo(vicen(:,:,jl,:),ht_i(:,:,jl),'T', 1. ) 
     262         ENDDO 
     263      ENDIF 
     264 
     265! T point to U point 
     266! T point to V point 
     267 
     268! Sea ice thickness 
     269      thick_iu(:,:)=0.0 
     270      thick_iv(:,:)=0.0 
     271      DO jj=1,jpjm1 
     272         DO ji=1,jpim1 
     273            thick_iu(ji,jj)=0.5*(thick_i(ji,jj)+thick_i(ji+1,jj))*umask(ji,jj,1) 
     274            thick_iv(ji,jj)=0.5*(thick_i(ji,jj)+thick_i(ji,jj+1))*vmask(ji,jj,1) 
     275         ENDDO 
     276      ENDDO 
     277 
     278      CALL lbc_lnk ( thick_iu , 'U', 1. ) 
     279      CALL lbc_lnk ( thick_iv , 'V', 1. ) 
     280 
     281! Snow depth 
     282      thick_su(:,:)=0.0 
     283      thick_sv(:,:)=0.0 
     284      DO jj=1,jpjm1 
     285         DO ji=1,jpim1 
     286            thick_su(ji,jj)=0.5*(thick_s(ji,jj)+thick_s(ji+1,jj))*umask(ji,jj,1) 
     287            thick_sv(ji,jj)=0.5*(thick_s(ji,jj)+thick_s(ji,jj+1))*vmask(ji,jj,1) 
     288         ENDDO 
     289      ENDDO 
     290 
     291      CALL lbc_lnk ( thick_su , 'U', 1. ) 
     292      CALL lbc_lnk ( thick_sv , 'V', 1. )      
    252293 
    253294      !                                      ! embedded sea ice 
     
    311352      nfsalt_da(:,:) = 0.0    
    312353      ndaice_da(:,:) = 0.0          
     354      ndsit_da(:,:)  = 0.0          
    313355#endif 
    314356      ! 
     
    469511      ztmp(:,:)=ndaice_da(:,:)*tmask(:,:,1) 
    470512      Call nemo2cice(ztmp,daice_da,'T', 1. ) 
     513!Ice thickness change (from assimilation) 
     514      ztmp(:,:)=ndsit_da(:,:)*tmask(:,:,1) 
     515      Call nemo2cice(ztmp,dhi_da,'T', 1. ) 
    471516#endif  
    472517 
Note: See TracChangeset for help on using the changeset viewer.