Changeset 11935


Ignore:
Timestamp:
2019-11-20T11:37:14+01:00 (6 months ago)
Author:
dcarneir
Message:

Changing SBC scripts (sbc_oce and sbcice_cice) to include sea ice thickness

Location:
branches/UKMO/dev_r5518_obs_oper_update_sit/NEMOGCM/NEMO/OPA_SRC/SBC
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_obs_oper_update_sit/NEMOGCM/NEMO/OPA_SRC/SBC/sbc_oce.F90

    r7960 r11935  
    121121   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   sprecip           !: solid precipitation                          [Kg/m2/s] 
    122122   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   fr_i              !: ice fraction = 1 - lead fraction      (between 0 to 1) 
     123   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   thick_i           !: ice thickness [m] 
     124   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   thick_s           !: snow depth [m] 
     125   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   vol_i             !: ice volume [m3] 
    123126#if defined key_cpl_carbon_cycle 
    124127   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   atm_co2           !: atmospheric pCO2                             [ppm] 
     
    167170         ! 
    168171      ALLOCATE( tprecip(jpi,jpj) , sprecip(jpi,jpj) , fr_i(jpi,jpj) ,     & 
     172         &      thick_i(jpi,jpj) , thick_s(jpi,jpj) , vol_i(jpi,jpj),     & 
    169173#if defined key_cpl_carbon_cycle 
    170174         &      atm_co2(jpi,jpj) ,                                        & 
  • branches/UKMO/dev_r5518_obs_oper_update_sit/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90

    r7960 r11935  
    229229      CALL lbc_lnk ( fr_iu , 'U', 1. ) 
    230230      CALL lbc_lnk ( fr_iv , 'V', 1. ) 
     231 
     232! Snow and ice thickness 
     233! vice, vsno is volume per unit area of grid cell = thickness 
     234      CALL cice2nemo(vice,thick_i,'T', 1. ) 
     235      CALL cice2nemo(vsno,thick_s,'T', 1. ) 
     236       
     237      IF ( (ksbc == jp_flx) .OR. (ksbc == jp_purecpl) ) THEN 
     238         DO jl = 1,ncat 
     239            CALL cice2nemo(vsnon(:,:,jl,:),ht_s(:,:,jl),'T', 1.) 
     240            CALL cice2nemo(vicen(:,:,jl,:),ht_i(:,:,jl),'T', 1.) 
     241         ENDDO 
     242      ENDIF 
     243       
     244! T point to U point 
     245! T point to V point 
     246 
     247! Sea ice thickness 
     248      thick_iu(:,:)=0.0 
     249      thick_iv(:,:)=0.0 
     250      DO jj=1,jpjm1 
     251         DO ji=1,jpim1 
     252            thick_iu(ji,jj)=0.5*(thick_i(ji,jj)+thick_i(ji+1,jj))*umask(ji,jj,1) 
     253            thick_iv(ji,jj)=0.5*(thick_i(ji,jj)+thick_i(ji,jj+1))*vmask(ji,jj,1) 
     254         ENDDO 
     255      ENDDO 
     256 
     257      CALL lbc_lnk ( thick_iu , 'U', 1. ) 
     258      CALL lbc_lnk ( thick_iv , 'V', 1. ) 
     259      
     260! Snow depth 
     261      thick_su(:,:)=0.0 
     262      thick_sv(:,:)=0.0 
     263      DO jj=1,jpjm1 
     264         DO ji=1,jpim1 
     265            thick_su(ji,jj)=0.5*(thick_s(ji,jj)+thick_s(ji+1,jj))*umask(ji,jj,1) 
     266            thick_sv(ji,jj)=0.5*(thick_s(ji,jj)+thick_s(ji,jj+1))*vmask(ji,jj,1) 
     267         ENDDO 
     268      ENDDO 
     269         
     270      CALL lbc_lnk ( thick_su , 'U', 1. ) 
     271      CALL lbc_lnk ( thick_sv , 'V', 1. )      
    231272 
    232273      !                                      ! embedded sea ice 
Note: See TracChangeset for help on using the changeset viewer.