Changeset 1471
- Timestamp:
- 2009-06-11T16:27:34+02:00 (15 years ago)
- Location:
- trunk/NEMO
- Files:
-
- 10 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/LIM_SRC_2/ice_2.F90
r1470 r1471 91 91 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: dmgwi !: Variation of the mass of snow ice 92 92 93 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: albege !: Albedo of the snow or ice (only for outputs)94 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: albecn !: Albedo of the ocean (only for outputs)95 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: tauc !: Cloud optical depth96 97 93 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: u_ice, v_ice !: two components of the ice velocity at I-point (m/s) 98 94 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: u_oce, v_oce !: two components of the ocean velocity at I-point (m/s) 99 95 100 REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpsmax) :: scal0 !: ???101 96 REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jplayersp1) :: tbif !: Temperature inside the ice/snow layer 102 103 !! REAL(wp), DIMENSION(jpi,jpj,0:jpkmax+1) :: reslum !: Relative absorption of solar radiation in each ocean level104 97 105 98 !!* moment used in the advection scheme -
trunk/NEMO/LIM_SRC_2/iceini_2.F90
r1465 r1471 68 68 tn_ice(:,:,1) = sist(:,:) ! initialisation of ice temperature 69 69 fr_i (:,:) = 1.0 - frld(:,:) ! initialisation of sea-ice fraction 70 # if defined key_coupled71 alb_ice(:,:,1) = albege(:,:) ! sea-ice albedo72 # endif73 70 ! 74 71 END SUBROUTINE ice_init_2 -
trunk/NEMO/LIM_SRC_2/limistate_2.F90
r1470 r1471 104 104 u_ice (:,:) = 0.e0 105 105 v_ice (:,:) = 0.e0 106 # if defined key_coupled107 albege(:,:) = 0.8 * tms(:,:)108 # endif109 106 110 107 !--- Moments for advection. -
trunk/NEMO/LIM_SRC_2/limrst_2.F90
r1470 r1471 112 112 CALL iom_rstput( iter, nitrst, numriw, 'frld' , frld (:,:) ) 113 113 CALL iom_rstput( iter, nitrst, numriw, 'sist' , sist (:,:) ) 114 # if defined key_coupled115 CALL iom_rstput( iter, nitrst, numriw, 'albege', albege(:,:) )116 # endif117 114 CALL iom_rstput( iter, nitrst, numriw, 'tbif1' , tbif (:,:,1) ) 118 115 CALL iom_rstput( iter, nitrst, numriw, 'tbif2' , tbif (:,:,2) ) … … 199 196 CALL iom_get( numrir, jpdom_autoglo, 'frld' , frld ) 200 197 CALL iom_get( numrir, jpdom_autoglo, 'sist' , sist ) 201 # if defined key_coupled202 CALL iom_get( numrir, jpdom_autoglo, 'albege', albege )203 # endif204 198 CALL iom_get( numrir, jpdom_autoglo, 'tbif1' , tbif(:,:,1) ) 205 199 CALL iom_get( numrir, jpdom_autoglo, 'tbif2' , tbif(:,:,2) ) -
trunk/NEMO/LIM_SRC_2/par_ice_2.F90
r1234 r1471 14 14 15 15 INTEGER, PUBLIC, PARAMETER :: jpl = 1 !: number of ice categories (only 1 in LIM-2) 16 INTEGER, PUBLIC, PARAMETER :: jpkmax = 1 !: ???17 INTEGER, PUBLIC, PARAMETER :: jpsmax = 2 !: ???18 16 19 17 INTEGER, PUBLIC, PARAMETER :: jplayers = 2 !: number of vertical ice layers -
trunk/NEMO/LIM_SRC_3/ice.F90
r1465 r1471 332 332 fhmec !: Heat flux due to snow loss during compression 333 333 334 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: & !:335 albege , & !: Albedo of the snow or ice (only for outputs)336 albecn , & !: Albedo of the ocean (only for outputs)337 tauc !: Cloud optical depth338 339 334 ! temporary arrays for dummy version of the code 340 335 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: & !: … … 382 377 vt_i_typ !: total volume contained in each ice type 383 378 384 REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpsmax) :: & !:385 scal0 !: ???386 387 379 REAL(wp), PUBLIC, DIMENSION(jpi,jpj,nlay_s,jpl) :: & !: 388 380 t_s, & !: Snow temperatures (K) … … 396 388 e_i, & !: Ice thermal contents [ Joules*10^9 ] 397 389 s_i !: Ice salinities 398 399 REAL(wp), DIMENSION(jpi,jpj,0:jpkmax+1) :: & !:400 reslum !: Relative absorption of solar radiation in each ocean level401 390 402 391 !!-------------------------------------------------------------------------- -
trunk/NEMO/LIM_SRC_3/iceini.F90
r1465 r1471 91 91 92 92 fr_i(:,:) = at_i(:,:) ! initialisation of sea-ice fraction 93 # if defined key_coupled94 Must be adpated to LIM395 alb_ice(:,:,:) = albege(:,:) ! sea-ice albedo96 # endif97 93 98 94 nstart = numit + nn_fsbc -
trunk/NEMO/LIM_SRC_3/limistate.F90
r1465 r1471 466 466 stress2_i(:,:) = 0.0 467 467 stress12_i(:,:) = 0.0 468 469 # if defined key_coupled470 albege(:,:) = 0.8 * tms(:,:)471 # endif472 468 473 469 !-------------------------------------------------------------------- -
trunk/NEMO/LIM_SRC_3/limrst.F90
r1469 r1471 137 137 CALL iom_rstput( iter, nitrst, numriw, znam , z2d ) 138 138 END DO 139 # if defined key_coupled 140 CALL iom_rstput( iter, nitrst, numriw, 'albege', albege(:,:) ) 141 # endif 139 142 140 DO jl = 1, jpl 143 141 WRITE(zchar,'(I1)') jl … … 495 493 ENDIF 496 494 497 # if defined key_coupled498 CALL iom_get( numrir, jpdom_autoglo, 'albege' , albege )499 # endif500 495 DO jl = 1, jpl 501 496 WRITE(zchar,'(I1)') jl -
trunk/NEMO/LIM_SRC_3/par_ice.F90
r1156 r1471 15 15 PUBLIC ! allows par_oce and par_kind to be known in ice modules 16 16 17 INTEGER, PARAMETER :: & !:18 jpkmax = 1 , & !: ???19 jpsmax = 2 !: ???20 21 17 !ICE THERMODYNAMICS 22 18 INTEGER , PARAMETER :: & !:
Note: See TracChangeset
for help on using the changeset viewer.