Changeset 1463 for trunk/NEMO/LIM_SRC_2
- Timestamp:
- 2009-06-09T16:45:31+02:00 (15 years ago)
- Location:
- trunk/NEMO/LIM_SRC_2
- Files:
-
- 6 edited
-
ice_2.F90 (modified) (1 diff)
-
iceini_2.F90 (modified) (1 diff)
-
limsbc_2.F90 (modified) (4 diffs)
-
limthd_2.F90 (modified) (2 diffs)
-
limwri_2.F90 (modified) (1 diff)
-
limwri_dimg_2.h90 (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/LIM_SRC_2/ice_2.F90
r1229 r1463 66 66 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: rdvonif !: Lateral Variation of ice volume (only used for outputs) 67 67 68 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: sist !: Sea-Ice Surface Temperature (Kelvin ??? degree ??? I don't know)68 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: sist !: Sea-Ice Surface Temperature (Kelvin) 69 69 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: tfu !: Freezing/Melting point temperature of sea water at SSS 70 70 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: hicif !: Ice thickness -
trunk/NEMO/LIM_SRC_2/iceini_2.F90
r1229 r1463 68 68 ENDIF 69 69 70 tn_ice(:,: ) = sist(:,:) ! initialisation of ice temperature70 tn_ice(:,:,1) = sist(:,:) ! initialisation of ice temperature 71 71 fr_i (:,:) = 1.0 - frld(:,:) ! initialisation of sea-ice fraction 72 72 # if defined key_coupled 73 alb_ice(:,: ) = albege(:,:) ! sea-ice albedo73 alb_ice(:,:,1) = albege(:,:) ! sea-ice albedo 74 74 # endif 75 75 ! -
trunk/NEMO/LIM_SRC_2/limsbc_2.F90
r1370 r1463 173 173 ! computation the solar flux at ocean surface 174 174 #if defined key_coupled 175 zqsr = qsr_tot(ji,jj) + ( fstric(ji,jj) - qsr_ice(ji,jj ) ) * ( 1.0 - pfrld(ji,jj) )175 zqsr = qsr_tot(ji,jj) + ( fstric(ji,jj) - qsr_ice(ji,jj,1) ) * ( 1.0 - pfrld(ji,jj) ) 176 176 #else 177 177 zqsr = pfrld(ji,jj) * qsr(ji,jj) + ( 1. - pfrld(ji,jj) ) * fstric(ji,jj) … … 286 286 IF ( lk_cpl ) THEN 287 287 ! Ice surface temperature 288 tn_ice(:,: ) = sist(:,:) ! sea-ice surface temperature288 tn_ice(:,:,1) = sist(:,:) ! sea-ice surface temperature 289 289 ! Computation of snow/ice and ocean albedo 290 290 ! INTERFACE 3D versus 2D … … 292 292 zhicif(:,:,1) = hicif(:,:) ; zhsnif(:,:,1) = hsnif(:,:) 293 293 CALL albedo_ice( zsist, zhicif, zhsnif, zalbp, zalb ) 294 alb_ice(:,: ) = 0.5 * ( zalbp(:,:,1) + zalb (:,:,1) ) ! Ice albedo (mean clear and overcast skys)294 alb_ice(:,:,1) = 0.5 * ( zalbp(:,:,1) + zalb (:,:,1) ) ! Ice albedo (mean clear and overcast skys) 295 295 ENDIF 296 296 … … 300 300 CALL prt_ctl(tab2d_1=utau , clinfo1=' lim_sbc: utau : ', mask1=umask, & 301 301 & tab2d_2=vtau , clinfo2=' vtau : ' , mask2=vmask ) 302 CALL prt_ctl(tab2d_1=fr_i , clinfo1=' lim_sbc: fr_i : ', tab2d_2=tn_ice , clinfo2=' tn_ice : ')302 CALL prt_ctl(tab2d_1=fr_i , clinfo1=' lim_sbc: fr_i : ', tab2d_2=tn_ice(:,:,1), clinfo2=' tn_ice : ') 303 303 ENDIF 304 304 -
trunk/NEMO/LIM_SRC_2/limthd_2.F90
r1228 r1463 195 195 zfi = 1.0 - pfrld(ji,jj) 196 196 qldif(ji,jj) = tms(ji,jj) * rdt_ice & 197 & * ( ( qsr_tot(ji,jj) - qsr_ice(ji,jj ) * zfi ) * ( 1.0 - thcm(ji,jj) ) &198 & + ( qns_tot(ji,jj) - qns_ice(ji,jj ) * zfi ) &197 & * ( ( qsr_tot(ji,jj) - qsr_ice(ji,jj,1) * zfi ) * ( 1.0 - thcm(ji,jj) ) & 198 & + ( qns_tot(ji,jj) - qns_ice(ji,jj,1) * zfi ) & 199 199 & + frld(ji,jj) * ( fdtcn(ji,jj) + ( 1.0 - zindb ) * fsbbq(ji,jj) ) ) 200 200 #else … … 261 261 ! 262 262 ! put the variable in a 1-D array for thermodynamics process 263 CALL tab_2d_1d_2( nbpb, frld_1d (1:nbpb) , frld , jpi, jpj, npb(1:nbpb) ) 264 CALL tab_2d_1d_2( nbpb, h_ice_1d (1:nbpb) , hicif , jpi, jpj, npb(1:nbpb) ) 265 CALL tab_2d_1d_2( nbpb, h_snow_1d (1:nbpb) , hsnif , jpi, jpj, npb(1:nbpb) ) 266 CALL tab_2d_1d_2( nbpb, sist_1d (1:nbpb) , sist , jpi, jpj, npb(1:nbpb) ) 267 CALL tab_2d_1d_2( nbpb, tbif_1d (1:nbpb , 1 ), tbif(:,:,1), jpi, jpj, npb(1:nbpb) ) 268 CALL tab_2d_1d_2( nbpb, tbif_1d (1:nbpb , 2 ), tbif(:,:,2), jpi, jpj, npb(1:nbpb) ) 269 CALL tab_2d_1d_2( nbpb, tbif_1d (1:nbpb , 3 ), tbif(:,:,3), jpi, jpj, npb(1:nbpb) ) 270 CALL tab_2d_1d_2( nbpb, qsr_ice_1d (1:nbpb) , qsr_ice , jpi, jpj, npb(1:nbpb) ) 271 CALL tab_2d_1d_2( nbpb, fr1_i0_1d (1:nbpb) , fr1_i0 , jpi, jpj, npb(1:nbpb) ) 272 CALL tab_2d_1d_2( nbpb, fr2_i0_1d (1:nbpb) , fr2_i0 , jpi, jpj, npb(1:nbpb) ) 273 CALL tab_2d_1d_2( nbpb, qns_ice_1d (1:nbpb) , qns_ice , jpi, jpj, npb(1:nbpb) ) 263 CALL tab_2d_1d_2( nbpb, frld_1d (1:nbpb) , frld , jpi, jpj, npb(1:nbpb) ) 264 CALL tab_2d_1d_2( nbpb, h_ice_1d (1:nbpb) , hicif , jpi, jpj, npb(1:nbpb) ) 265 CALL tab_2d_1d_2( nbpb, h_snow_1d (1:nbpb) , hsnif , jpi, jpj, npb(1:nbpb) ) 266 CALL tab_2d_1d_2( nbpb, sist_1d (1:nbpb) , sist , jpi, jpj, npb(1:nbpb) ) 267 CALL tab_2d_1d_2( nbpb, tbif_1d (1:nbpb , 1 ), tbif(:,:,1) , jpi, jpj, npb(1:nbpb) ) 268 CALL tab_2d_1d_2( nbpb, tbif_1d (1:nbpb , 2 ), tbif(:,:,2) , jpi, jpj, npb(1:nbpb) ) 269 CALL tab_2d_1d_2( nbpb, tbif_1d (1:nbpb , 3 ), tbif(:,:,3) , jpi, jpj, npb(1:nbpb) ) 270 CALL tab_2d_1d_2( nbpb, qsr_ice_1d (1:nbpb) , qsr_ice(:,:,1) , jpi, jpj, npb(1:nbpb) ) 271 CALL tab_2d_1d_2( nbpb, fr1_i0_1d (1:nbpb) , fr1_i0 , jpi, jpj, npb(1:nbpb) ) 272 CALL tab_2d_1d_2( nbpb, fr2_i0_1d (1:nbpb) , fr2_i0 , jpi, jpj, npb(1:nbpb) ) 273 CALL tab_2d_1d_2( nbpb, qns_ice_1d(1:nbpb) , qns_ice(:,:,1), jpi, jpj, npb(1:nbpb) ) 274 CALL tab_2d_1d_2( nbpb, dqns_ice_1d(1:nbpb) , dqns_ice(:,:,1), jpi, jpj, npb(1:nbpb) ) 274 275 IF( .NOT. lk_cpl ) THEN 275 CALL tab_2d_1d_2( nbpb, qla_ice_1d (1:nbpb) , qla_ice, jpi, jpj, npb(1:nbpb) )276 CALL tab_2d_1d_2( nbpb, dqla_ice_1d(1:nbpb) , dqla_ice 276 CALL tab_2d_1d_2( nbpb, qla_ice_1d (1:nbpb) , qla_ice(:,:,1), jpi, jpj, npb(1:nbpb) ) 277 CALL tab_2d_1d_2( nbpb, dqla_ice_1d(1:nbpb) , dqla_ice(:,:,1), jpi, jpj, npb(1:nbpb) ) 277 278 ENDIF 278 CALL tab_2d_1d_2( nbpb, dqns_ice_1d(1:nbpb) , dqns_ice , jpi, jpj, npb(1:nbpb) )279 279 CALL tab_2d_1d_2( nbpb, tfu_1d (1:nbpb) , tfu , jpi, jpj, npb(1:nbpb) ) 280 280 CALL tab_2d_1d_2( nbpb, sprecip_1d (1:nbpb) , sprecip , jpi, jpj, npb(1:nbpb) ) -
trunk/NEMO/LIM_SRC_2/limwri_2.F90
r1359 r1463 158 158 zcmo(ji,jj,15) = utaui_ice(ji,jj) 159 159 zcmo(ji,jj,16) = vtaui_ice(ji,jj) 160 zcmo(ji,jj,17) = qsr_ice(ji,jj )161 zcmo(ji,jj,18) = qns_ice(ji,jj )160 zcmo(ji,jj,17) = qsr_ice(ji,jj,1) 161 zcmo(ji,jj,18) = qns_ice(ji,jj,1) 162 162 zcmo(ji,jj,19) = sprecip(ji,jj) 163 163 END DO -
trunk/NEMO/LIM_SRC_2/limwri_dimg_2.h90
r1156 r1463 123 123 zcmo(ji,jj,15) = utaui_ice(ji,jj) 124 124 zcmo(ji,jj,16) = vtaui_ice(ji,jj) 125 zcmo(ji,jj,17) = qsr_ice(ji,jj )126 zcmo(ji,jj,18) = qns_ice(ji,jj )125 zcmo(ji,jj,17) = qsr_ice(ji,jj,1) 126 zcmo(ji,jj,18) = qns_ice(ji,jj,1) 127 127 zcmo(ji,jj,19) = sprecip(ji,jj) 128 128 END DO … … 166 166 rcmoy(ji,jj,15) = utaui_ice(ji,jj) 167 167 rcmoy(ji,jj,16) = vtaui_ice(ji,jj) 168 rcmoy(ji,jj,17) = qsr_ice(ji,jj )169 rcmoy(ji,jj,18) = qns_ice(ji,jj )168 rcmoy(ji,jj,17) = qsr_ice(ji,jj,1) 169 rcmoy(ji,jj,18) = qns_ice(ji,jj,1) 170 170 rcmoy(ji,jj,19) = sprecip(ji,jj) 171 171 END DO
Note: See TracChangeset
for help on using the changeset viewer.
