- Timestamp:
- 2017-07-25T19:44:54+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/limmp.F90
r8369 r8373 249 249 !!------------------------------------------------------------------- 250 250 251 INTEGER, POINTER, DIMENSION(:) :: indxi ! compressed indices for cells with ice melting 252 INTEGER, POINTER, DIMENSION(:) :: indxj ! 253 254 REAL(wp), POINTER, DIMENSION(:,:) :: zwfx_mlw ! available meltwater for melt ponding 255 REAL(wp), POINTER, DIMENSION(:,:,:) :: zrfrac ! fraction of available meltwater retained for melt ponding 256 257 REAL(wp), PARAMETER :: zrmin = 0.15_wp ! minimum fraction of available meltwater retained for melt ponding 258 REAL(wp), PARAMETER :: zrmax = 0.70_wp ! maximum '' '' '' '' '' 259 REAL(wp), PARAMETER :: zrexp = 0.01_wp ! rate constant to refreeze melt ponds 260 REAL(wp), PARAMETER :: zpnd_aspect = 0.8_wp ! pond aspect ratio 261 262 REAL(wp) :: zhi ! dummy ice thickness 263 REAL(wp) :: zhs ! dummy snow depth 264 REAL(wp) :: zTp ! reference temperature 265 REAL(wp) :: zdTs ! dummy temperature difference 266 REAL(wp) :: z1_rhofw ! inverse freshwater density 267 REAL(wp) :: z1_zpnd_aspect ! inverse pond aspect ratio 268 REAL(wp) :: zvpold ! dummy pond volume 269 270 INTEGER :: ji, jj, jl, ij ! loop indices 271 INTEGER :: icells ! size of dummy array 272 251 INTEGER, DIMENSION(jpij) :: indxi ! compressed indices for cells with ice melting 252 INTEGER, DIMENSION(jpij) :: indxj ! 253 254 REAL(wp), DIMENSION(jpi,jpj) :: zwfx_mlw ! available meltwater for melt ponding 255 REAL(wp), DIMENSION(jpi,jpj,jpl) :: zrfrac ! fraction of available meltwater retained for melt ponding 256 257 REAL(wp), PARAMETER :: zrmin = 0.15_wp ! minimum fraction of available meltwater retained for melt ponding 258 REAL(wp), PARAMETER :: zrmax = 0.70_wp ! maximum '' '' '' '' '' 259 REAL(wp), PARAMETER :: zrexp = 0.01_wp ! rate constant to refreeze melt ponds 260 REAL(wp), PARAMETER :: zpnd_aspect = 0.8_wp ! pond aspect ratio 261 262 REAL(wp) :: zhi ! dummy ice thickness 263 REAL(wp) :: zhs ! dummy snow depth 264 REAL(wp) :: zTp ! reference temperature 265 REAL(wp) :: zdTs ! dummy temperature difference 266 REAL(wp) :: z1_rhofw ! inverse freshwater density 267 REAL(wp) :: z1_zpnd_aspect ! inverse pond aspect ratio 268 REAL(wp) :: zvpold ! dummy pond volume 269 270 INTEGER :: ji, jj, jl, ij ! loop indices 271 INTEGER :: icells ! size of dummy array 273 272 !!------------------------------------------------------------------- 274 275 CALL wrk_alloc( jpi*jpj, indxi, indxj)276 CALL wrk_alloc( jpi,jpj, zwfx_mlw )277 CALL wrk_alloc( jpi,jpj,jpl, zrfrac )278 279 273 z1_rhofw = 1. / rhofw 280 274 z1_zpnd_aspect = 1. / zpnd_aspect … … 390 384 391 385 ENDIF 392 393 CALL wrk_dealloc( jpi*jpj, indxi, indxj)394 CALL wrk_dealloc( jpi,jpj, zwfx_mlw )395 CALL wrk_dealloc( jpi,jpj,jpl, zrfrac )396 386 397 387 END SUBROUTINE lim_mp_cesm
Note: See TracChangeset
for help on using the changeset viewer.