Changeset 11510


Ignore:
Timestamp:
2019-09-06T18:08:57+02:00 (13 months ago)
Author:
clem
Message:

add ponds for bdy

Location:
NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/cfgs/SHARED/namelist_ref

    r11507 r11510  
    639639   bn_s_i      = 'NOT USED'              ,         24.       , 'iicesalt',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
    640640   ! if bn_t_i, bn_tsu and bn_s_i are "not used", then define arbitrary temperatures and salinity 
    641    rn_ice_tem    = 270.       !  arbitrary temperature of incoming sea ice 
    642    rn_ice_sal    = 10.        !       --   salinity           -- 
    643    rn_ice_age    = 30.        !       --   age                -- 
     641   rn_ice_tem  = 270.         !  arbitrary temperature of incoming sea ice 
     642   rn_ice_sal  = 10.          !       --   salinity           -- 
     643   rn_ice_age  = 30.          !       --   age                -- 
     644   rn_ice_apnd = 0.2          !       --   pond conc.         -- 
     645   rn_ice_hpnd = 0.05         !       --   pond depth         -- 
    644646/ 
    645647!----------------------------------------------------------------------- 
  • NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/src/OCE/BDY/bdy_oce.F90

    r11507 r11510  
    111111   REAL(wp), DIMENSION(jp_bdy) ::   rice_sal                !: salinity    of incoming sea ice 
    112112   REAL(wp), DIMENSION(jp_bdy) ::   rice_age                !: age         of incoming sea ice 
     113   REAL(wp), DIMENSION(jp_bdy) ::   rice_apnd               !: pond conc.  of incoming sea ice 
     114   REAL(wp), DIMENSION(jp_bdy) ::   rice_hpnd               !: pond thick. of incoming sea ice 
    113115   ! 
    114116   !!---------------------------------------------------------------------- 
  • NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/src/OCE/BDY/bdydta.F90

    r11507 r11510  
    362362      !                                                         ! =F => baroclinic velocities in 3D boundary data 
    363363      LOGICAL                                ::   ln_zinterp    ! =T => requires a vertical interpolation of the bdydta 
    364       REAL(wp)                               ::   rn_ice_tem, rn_ice_sal, rn_ice_age 
     364      REAL(wp)                               ::   rn_ice_tem, rn_ice_sal, rn_ice_age, rn_ice_apnd, rn_ice_hpnd  
    365365      INTEGER                                ::   ipk,ipl       ! 
    366366      INTEGER                                ::   idvar         ! variable ID 
     
    380380      NAMELIST/nambdy_dta/ cn_dir, bn_tem, bn_sal, bn_u3d, bn_v3d, bn_ssh, bn_u2d, bn_v2d  
    381381      NAMELIST/nambdy_dta/ bn_a_i, bn_h_i, bn_h_s, bn_t_i, bn_t_s, bn_tsu, bn_s_i 
    382       NAMELIST/nambdy_dta/ rn_ice_tem, rn_ice_sal, rn_ice_age 
     382      NAMELIST/nambdy_dta/ rn_ice_tem, rn_ice_sal, rn_ice_age, rn_ice_apnd, rn_ice_hpnd 
    383383      NAMELIST/nambdy_dta/ ln_full_vel, ln_zinterp 
    384384      !!--------------------------------------------------------------------------- 
     
    434434         ENDIF 
    435435 
    436          ! temp, sal and age of incoming ice 
    437          rice_tem(jbdy) = rn_ice_tem 
    438          rice_sal(jbdy) = rn_ice_sal 
    439          rice_age(jbdy) = rn_ice_age 
     436         ! temp, salt, age and ponds of incoming ice 
     437         rice_tem (jbdy) = rn_ice_tem 
     438         rice_sal (jbdy) = rn_ice_sal 
     439         rice_age (jbdy) = rn_ice_age 
     440         rice_apnd(jbdy) = rn_ice_apnd 
     441         rice_hpnd(jbdy) = rn_ice_hpnd 
    440442          
    441443          
  • NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/src/OCE/BDY/bdyice.F90

    r11507 r11510  
    225225                  ! 
    226226               CASE( 0 )   ! velocity is inward 
    227                   oa_i(ji,jj,jl) = rice_age(jbdy) * a_i(ji,jj,jl)   ! age 
    228                   a_ip(ji,jj,jl) = 0._wp                            ! pond concentration 
    229                   v_ip(ji,jj,jl) = 0._wp                            ! pond volume 
     227                  oa_i(ji,jj,jl) = rice_age(jbdy) * a_i(ji,jj,jl)     ! age 
     228                  IF( ln_pnd_CST .OR. ln_pnd_H12 ) THEN 
     229                     a_ip(ji,jj,jl) = rice_apnd(jbdy)                 ! pond concentration 
     230                     v_ip(ji,jj,jl) = rice_apnd(jbdy)*rice_hpnd(jbdy) ! pond volume 
     231                  ENDIF 
    230232               CASE( 1 )   ! velocity is outward 
    231                   oa_i(ji,jj,jl) = oa_i(ib,jb,jl) ! age 
    232                   a_ip(ji,jj,jl) = a_ip(ib,jb,jl) ! pond concentration 
    233                   v_ip(ji,jj,jl) = v_ip(ib,jb,jl) ! pond volume 
     233                  oa_i(ji,jj,jl) = oa_i(ib,jb,jl)                     ! age 
     234                  IF( ln_pnd_CST .OR. ln_pnd_H12 ) THEN 
     235                     a_ip(ji,jj,jl) = a_ip(ib,jb,jl)                  ! pond concentration 
     236                     v_ip(ji,jj,jl) = v_ip(ib,jb,jl)                  ! pond volume 
     237                  ENDIF 
    234238               END SELECT 
    235239               ! 
Note: See TracChangeset for help on using the changeset viewer.