- Timestamp:
- 2017-04-13T16:21:08+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r7881_no_wrk_alloc/NEMOGCM/NEMO/LIM_SRC_2/limthd_zdf_2.F90
r5407 r7910 21 21 USE in_out_manager 22 22 USE lib_mpp ! MPP library 23 USE wrk_nemo ! work arrays24 23 USE lib_fortran ! Fortran utilities (allows no signed zero when 'key_nosignedzero' defined) 25 24 … … 73 72 !! 74 73 INTEGER :: ji ! dummy loop indices 75 REAL(wp), POINTER, DIMENSION(:) :: zqcmlts ! energy due to surface melting76 REAL(wp), POINTER, DIMENSION(:) :: zqcmltb ! energy due to bottom melting77 REAL(wp), POINTER, DIMENSION(:) :: ztsmlt ! snow/ice surface melting temperature78 REAL(wp), POINTER, DIMENSION(:) :: ztbif ! int. temp. at the mid-point of the 1st layer of the snow/ice sys.79 REAL(wp), POINTER, DIMENSION(:) :: zksn ! effective conductivity of snow80 REAL(wp), POINTER, DIMENSION(:) :: zkic ! effective conductivity of ice81 REAL(wp), POINTER, DIMENSION(:) :: zksndh ! thermal cond. at the mid-point of the 1st layer of the snow/ice sys.82 REAL(wp), POINTER, DIMENSION(:) :: zfcsu ! conductive heat flux at the surface of the snow/ice system83 REAL(wp), POINTER, DIMENSION(:) :: zfcsudt ! = zfcsu * dt84 REAL(wp), POINTER, DIMENSION(:) :: zi0 ! frac. of the net SW rad. which is not absorbed at the surface85 REAL(wp), POINTER, DIMENSION(:) :: z1mi0 ! fraction of the net SW radiation absorbed at the surface86 REAL(wp), POINTER, DIMENSION(:) :: zqmax ! maximum energy stored in brine pockets87 REAL(wp), POINTER, DIMENSION(:) :: zrcpdt ! h_su*rho_su*cp_su/dt(h_su being the thick. of surf. layer)88 REAL(wp), POINTER, DIMENSION(:) :: zts_old ! previous surface temperature89 REAL(wp), POINTER, DIMENSION(:) :: zidsn , z1midsn , zidsnic ! temporary variables90 REAL(wp), POINTER, DIMENSION(:) :: zfnet ! net heat flux at the top surface( incl. conductive heat flux)91 REAL(wp), POINTER, DIMENSION(:) :: zsprecip ! snow accumulation92 REAL(wp), POINTER, DIMENSION(:) :: zhsnw_old ! previous snow thickness93 REAL(wp), POINTER, DIMENSION(:) :: zdhictop ! change in ice thickness due to top surf ablation/accretion94 REAL(wp), POINTER, DIMENSION(:) :: zdhicbot ! change in ice thickness due to bottom surf abl/acc95 REAL(wp), POINTER, DIMENSION(:) :: zqsup ! energy transmitted to ocean (coming from top surf abl/acc)96 REAL(wp), POINTER, DIMENSION(:) :: zqocea ! energy transmitted to ocean (coming from bottom sur abl/acc)97 REAL(wp), POINTER, DIMENSION(:) :: zfrl_old ! previous sea/ice fraction98 REAL(wp), POINTER, DIMENSION(:) :: zfrld_1d ! new sea/ice fraction99 REAL(wp), POINTER, DIMENSION(:) :: zep ! internal temperature of the 2nd layer of the snow/ice system74 REAL(wp), DIMENSION(jpij) :: zqcmlts ! energy due to surface melting 75 REAL(wp), DIMENSION(jpij) :: zqcmltb ! energy due to bottom melting 76 REAL(wp), DIMENSION(jpij) :: ztsmlt ! snow/ice surface melting temperature 77 REAL(wp), DIMENSION(jpij) :: ztbif ! int. temp. at the mid-point of the 1st layer of the snow/ice sys. 78 REAL(wp), DIMENSION(jpij) :: zksn ! effective conductivity of snow 79 REAL(wp), DIMENSION(jpij) :: zkic ! effective conductivity of ice 80 REAL(wp), DIMENSION(jpij) :: zksndh ! thermal cond. at the mid-point of the 1st layer of the snow/ice sys. 81 REAL(wp), DIMENSION(jpij) :: zfcsu ! conductive heat flux at the surface of the snow/ice system 82 REAL(wp), DIMENSION(jpij) :: zfcsudt ! = zfcsu * dt 83 REAL(wp), DIMENSION(jpij) :: zi0 ! frac. of the net SW rad. which is not absorbed at the surface 84 REAL(wp), DIMENSION(jpij) :: z1mi0 ! fraction of the net SW radiation absorbed at the surface 85 REAL(wp), DIMENSION(jpij) :: zqmax ! maximum energy stored in brine pockets 86 REAL(wp), DIMENSION(jpij) :: zrcpdt ! h_su*rho_su*cp_su/dt(h_su being the thick. of surf. layer) 87 REAL(wp), DIMENSION(jpij) :: zts_old ! previous surface temperature 88 REAL(wp), DIMENSION(jpij) :: zidsn , z1midsn , zidsnic ! temporary variables 89 REAL(wp), DIMENSION(jpij) :: zfnet ! net heat flux at the top surface( incl. conductive heat flux) 90 REAL(wp), DIMENSION(jpij) :: zsprecip ! snow accumulation 91 REAL(wp), DIMENSION(jpij) :: zhsnw_old ! previous snow thickness 92 REAL(wp), DIMENSION(jpij) :: zdhictop ! change in ice thickness due to top surf ablation/accretion 93 REAL(wp), DIMENSION(jpij) :: zdhicbot ! change in ice thickness due to bottom surf abl/acc 94 REAL(wp), DIMENSION(jpij) :: zqsup ! energy transmitted to ocean (coming from top surf abl/acc) 95 REAL(wp), DIMENSION(jpij) :: zqocea ! energy transmitted to ocean (coming from bottom sur abl/acc) 96 REAL(wp), DIMENSION(jpij) :: zfrl_old ! previous sea/ice fraction 97 REAL(wp), DIMENSION(jpij) :: zfrld_1d ! new sea/ice fraction 98 REAL(wp), DIMENSION(jpij) :: zep ! internal temperature of the 2nd layer of the snow/ice system 100 99 REAL(wp), DIMENSION(3) :: & 101 100 zplediag & ! principle diagonal, subdiag. and supdiag. of the … … 164 163 , zibmlt, ziqr, zihgnew, zind, ztmp ! temporary scalars 165 164 !!---------------------------------------------------------------------- 166 CALL wrk_alloc( jpij, ztsmlt, ztbif , zksn , zkic , zksndh , zfcsu , zfcsudt , zi0 , z1mi0 , zqmax )167 CALL wrk_alloc( jpij, zrcpdt, zts_old, zidsn , z1midsn , zidsnic, zfnet , zsprecip, zhsnw_old, zdhictop, zdhicbot )168 CALL wrk_alloc( jpij, zqsup , zqocea , zfrl_old, zfrld_1d, zep , zqcmlts, zqcmltb )169 165 170 166 !----------------------------------------------------------------------- … … 809 805 END DO 810 806 ! 811 CALL wrk_dealloc( jpij, ztsmlt, ztbif , zksn , zkic , zksndh , zfcsu , zfcsudt , zi0 , z1mi0 , zqmax )812 CALL wrk_dealloc( jpij, zrcpdt, zts_old, zidsn , z1midsn , zidsnic, zfnet , zsprecip, zhsnw_old, zdhictop, zdhicbot )813 CALL wrk_dealloc( jpij, zqsup , zqocea , zfrl_old, zfrld_1d, zep , zqcmlts, zqcmltb )814 807 ! 815 808 END SUBROUTINE lim_thd_zdf_2
Note: See TracChangeset
for help on using the changeset viewer.