- Timestamp:
- 2015-04-13T12:13:10+02:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5204_CNRS_PISCES_dcy/NEMOGCM/NEMO/LIM_SRC_2/limthd_2.F90
r4990 r5206 114 114 CALL wrk_alloc( jpi, jpj, jpk, zmsk ) 115 115 116 IF( kt == nit000 ) CALL lim_thd_init_2 ! Initialization (first time-step only) 116 IF( kt == nit000 ) THEN 117 CALL lim_thd_init_2 ! Initialization (first time-step only) 118 IF( l_trcdm2dc ) ALLOCATE( fstric_mean(jpi,jpj), fstbif_mean_1d(jpij), qsr_ice_mean_1d(jpij) ) 119 ENDIF 117 120 118 121 !-------------------------------------------! … … 137 140 rdq_ice(:,:) = 0.e0 ! heat content associated with rdm_ice 138 141 zmsk (:,:,:) = 0.e0 139 IF( ltrcdm2dc_ice ) fstric_daymean (:,:) = 0.e0 ! part of solar radiation absorbing inside the ice 142 ! 143 IF( l_trcdm2dc ) fstric_mean(:,:) = 0.e0 ! part of solar radiation absorbing inside the ice 140 144 141 145 ! set to zero snow thickness smaller than epsi04 … … 285 289 CALL tab_2d_1d_2( nbpb, tbif_1d (1:nbpb , 3 ), tbif(:,:,3) , jpi, jpj, npb(1:nbpb) ) 286 290 CALL tab_2d_1d_2( nbpb, qsr_ice_1d (1:nbpb) , qsr_ice(:,:,1) , jpi, jpj, npb(1:nbpb) ) 287 IF( l trcdm2dc_ice) &291 IF( l_trcdm2dc ) & 288 292 & CALL tab_2d_1d_2( nbpb, qsr_ice_mean_1d (1:nbpb), qsr_ice_mean(:,:,1), jpi, jpj, npb(1:nbpb) ) 289 293 CALL tab_2d_1d_2( nbpb, fr1_i0_1d (1:nbpb) , fr1_i0 , jpi, jpj, npb(1:nbpb) ) … … 336 340 CALL tab_1d_2d_2( nbpb, qsr_ice(:,:,1), npb, qsr_ice_1d(1:nbpb) , jpi, jpj ) 337 341 CALL tab_1d_2d_2( nbpb, qns_ice(:,:,1), npb, qns_ice_1d(1:nbpb) , jpi, jpj ) 338 IF( l trcdm2dc_ice )THEN339 CALL tab_1d_2d_2( nbpb, fstric_ daymean , npb, fstbif_daymean_1d (1:nbpb), jpi, jpj )340 CALL tab_1d_2d_2( nbpb, qsr_ice_mean(:,:,1), npb, qsr_ice_mean_1d(1:nbpb), jpi, jpj )342 IF( l_trcdm2dc ) & 343 CALL tab_1d_2d_2( nbpb, fstric_mean , npb, fstbif_mean_1d (1:nbpb), jpi, jpj ) 344 CALL tab_1d_2d_2( nbpb, qsr_ice_mean(:,:,1), npb, qsr_ice_mean_1d(1:nbpb), jpi, jpj ) 341 345 ENDIF 342 IF( .NOT. lk_cpl ) CALL tab_1d_2d_2( nbpb, qla_ice(:,:,1), npb, qla_ice_1d(1:nbpb) 346 IF( .NOT. lk_cpl ) CALL tab_1d_2d_2( nbpb, qla_ice(:,:,1), npb, qla_ice_1d(1:nbpb), jpi, jpj ) 343 347 ! 344 348 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.