Changeset 398
- Timestamp:
- Jul 22, 2002, 5:50:38 PM (22 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ.3.3/branches/rel-LF/libf/phylmd/interface_surf.F90
r397 r398 43 43 real, allocatable, dimension(:),save :: coastalflow, riverflow 44 44 !!$PB 45 !!REAL, ALLOCATABLE, DIMENSION(:,:), SAVE :: tmp_rriv, tmp_rcoa45 REAL, ALLOCATABLE, DIMENSION(:,:), SAVE :: tmp_rriv, tmp_rcoa 46 46 !! pour simuler la fonte des glaciers antarctiques 47 47 REAL, ALLOCATABLE, DIMENSION(:,:), SAVE :: coeff_iceberg … … 276 276 endif 277 277 !!$PB 278 !!ALLOCATE (tmp_rriv(iim,jjm+1), stat=error)279 ! !if (error /= 0) then280 !!abort_message='Pb allocation tmp_rriv'281 !!call abort_gcm(modname,abort_message,1)282 !!endif283 !!ALLOCATE (tmp_rcoa(iim,jjm+1), stat=error)284 !!if (error /= 0) then285 !!abort_message='Pb allocation tmp_rcoa'286 !!call abort_gcm(modname,abort_message,1)287 !!endif278 ALLOCATE (tmp_rriv(iim,jjm+1), stat=error) 279 if (error /= 0) then 280 abort_message='Pb allocation tmp_rriv' 281 call abort_gcm(modname,abort_message,1) 282 endif 283 ALLOCATE (tmp_rcoa(iim,jjm+1), stat=error) 284 if (error /= 0) then 285 abort_message='Pb allocation tmp_rcoa' 286 call abort_gcm(modname,abort_message,1) 287 endif 288 288 !!$ 289 289 else if (size(coastalflow) /= knon) then … … 1022 1022 bidule=0. 1023 1023 bidule(1:knon)=riverflow(1:knon) 1024 ! call gath2cpl(bidule, tmp_rriv, klon, knon,iim,jjm,ktindex)1024 call gath2cpl(bidule, tmp_rriv, klon, knon,iim,jjm,knindex) 1025 1025 bidule=0. 1026 1026 bidule(1:knon)=coastalflow(1:knon) 1027 ! call gath2cpl(bidule, tmp_rcoa, klon, knon,iim,jjm,ktindex)1027 call gath2cpl(bidule, tmp_rcoa, klon, knon,iim,jjm,knindex) 1028 1028 alb_new(1:knon) = albedo_out(1:knon,1) 1029 1029 alblw(1:knon) = albedo_out(1:knon,2) … … 1144 1144 real, allocatable, dimension(:,:,:),save :: tmp_snow, tmp_evap, tmp_tsol 1145 1145 real, allocatable, dimension(:,:,:),save :: tmp_fder, tmp_albe, tmp_taux 1146 real, allocatable, dimension(:,:,:),save :: tmp_tauy, tmp_rriv, tmp_rcoa1147 !!REAL, ALLOCATABLE, DIMENSION(:,:,:),SAVE :: tmp_tauy1146 !!$ real, allocatable, dimension(:,:,:),save :: tmp_tauy, tmp_rriv, tmp_rcoa 1147 REAL, ALLOCATABLE, DIMENSION(:,:,:),SAVE :: tmp_tauy 1148 1148 ! variables a passer au coupleur 1149 1149 real, dimension(iim, jjm+1) :: wri_sol_ice, wri_sol_sea, wri_nsol_ice … … 1336 1336 cpl_tauy(ig,cpl_index) = cpl_tauy(ig,cpl_index) & 1337 1337 & + tauy(ig) / FLOAT(nexca) 1338 cpl_rriv(ig,cpl_index) = cpl_rriv(ig,cpl_index) &1339 & + riverflow(ig) / FLOAT(nexca)/dtime1340 cpl_rcoa(ig,cpl_index) = cpl_rcoa(ig,cpl_index) &1341 & + coastalflow(ig) / FLOAT(nexca)/dtime1338 !!$ cpl_rriv(ig,cpl_index) = cpl_rriv(ig,cpl_index) & 1339 !!$ & + riverflow(ig) / FLOAT(nexca)/dtime 1340 !!$ cpl_rcoa(ig,cpl_index) = cpl_rcoa(ig,cpl_index) & 1341 !!$ & + coastalflow(ig) / FLOAT(nexca)/dtime 1342 1342 enddo 1343 !IF (cpl_index .EQ. 1) THEN1344 !cpl_rriv(:,:) = cpl_rriv(:,:) + tmp_rriv(:,:) / FLOAT(nexca)1345 !cpl_rcoa(:,:) = cpl_rcoa(:,:) + tmp_rcoa(:,:) / FLOAT(nexca)1346 !ENDIF1343 IF (cpl_index .EQ. 1) THEN 1344 cpl_rriv(:,:) = cpl_rriv(:,:) + tmp_rriv(:,:) / FLOAT(nexca) 1345 cpl_rcoa(:,:) = cpl_rcoa(:,:) + tmp_rcoa(:,:) / FLOAT(nexca) 1346 ENDIF 1347 1347 endif 1348 1348 … … 1442 1442 allocate(tmp_taux(iim,jjm+1,2), stat=error); sum_error = sum_error + error 1443 1443 allocate(tmp_tauy(iim,jjm+1,2), stat=error); sum_error = sum_error + error 1444 allocate(tmp_rriv(iim,jjm+1,2), stat=error); sum_error = sum_error + error1445 allocate(tmp_rcoa(iim,jjm+1,2), stat=error); sum_error = sum_error + error1444 !!$ allocate(tmp_rriv(iim,jjm+1,2), stat=error); sum_error = sum_error + error 1445 !!$ allocate(tmp_rcoa(iim,jjm+1,2), stat=error); sum_error = sum_error + error 1446 1446 if (sum_error /= 0) then 1447 1447 abort_message='Pb allocation variables couplees pour l''ecriture' … … 1465 1465 call gath2cpl(cpl_taux(1,cpl_index), tmp_taux(1,1,cpl_index), klon, knon,iim,jjm, knindex) 1466 1466 call gath2cpl(cpl_tauy(1,cpl_index), tmp_tauy(1,1,cpl_index), klon, knon,iim,jjm, knindex) 1467 call gath2cpl(cpl_rriv(1,cpl_index), tmp_rriv(1,1,cpl_index), klon, knon,iim,jjm, knindex)1468 call gath2cpl(cpl_rcoa(1,cpl_index), tmp_rcoa(1,1,cpl_index), klon, knon,iim,jjm, knindex)1467 !!$ call gath2cpl(cpl_rriv(1,cpl_index), tmp_rriv(1,1,cpl_index), klon, knon,iim,jjm, knindex) 1468 !!$ call gath2cpl(cpl_rcoa(1,cpl_index), tmp_rcoa(1,1,cpl_index), klon, knon,iim,jjm, knindex) 1469 1469 1470 1470 ! … … 1495 1495 & tmp_snow(:,:,2) * tamp_srf(:,:,2) / deno 1496 1496 !!$PB 1497 wri_rriv = tmp_rriv(:,:,1) * tamp_srf(:,:,1) / deno + &1498 & tmp_rriv(:,:,2) * tamp_srf(:,:,2) / deno1499 wri_rcoa = tmp_rcoa(:,:,1) * tamp_srf(:,:,1) / deno + &1500 & tmp_rcoa(:,:,2) * tamp_srf(:,:,2) / deno1497 !!$ wri_rriv = tmp_rriv(:,:,1) * tamp_srf(:,:,1) / deno + & 1498 !!$ & tmp_rriv(:,:,2) * tamp_srf(:,:,2) / deno 1499 !!$ wri_rcoa = tmp_rcoa(:,:,1) * tamp_srf(:,:,1) / deno + & 1500 !!$ & tmp_rcoa(:,:,2) * tamp_srf(:,:,2) / deno 1501 1501 wri_taux = tmp_taux(:,:,1) * tamp_srf(:,:,1) / deno + & 1502 1502 & tmp_taux(:,:,2) * tamp_srf(:,:,2) / deno … … 1583 1583 deallocate(tmp_tauy, stat=error); sum_error = sum_error + error 1584 1584 !!$PB 1585 deallocate(tmp_rriv, stat=error); sum_error = sum_error + error1586 deallocate(tmp_rcoa, stat=error); sum_error = sum_error + error1585 !!$ deallocate(tmp_rriv, stat=error); sum_error = sum_error + error 1586 !!$ deallocate(tmp_rcoa, stat=error); sum_error = sum_error + error 1587 1587 if (sum_error /= 0) then 1588 1588 abort_message='Pb deallocation variables couplees'
Note: See TracChangeset
for help on using the changeset viewer.