- Timestamp:
- Jul 20, 2024, 12:00:23 PM (2 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/fonte_neige_mod.F90
r5022 r5087 178 178 IF (iso_eau > 0) THEN 179 179 CALL iso_verif_egalite_vect1D( & 180 &xtrun_off_lic_0,run_off_lic_0,'fonte_neige 100', &181 &niso,klon)180 xtrun_off_lic_0,run_off_lic_0,'fonte_neige 100', & 181 niso,klon) 182 182 ENDIF !IF (iso_eau > 0) THEN 183 183 #endif … … 233 233 snow, qsol, tsurf_new, evap & 234 234 #ifdef ISO 235 &,fq_fonte_diag,fqfonte_diag,snow_evap_diag,fqcalving_diag &236 &,max_eau_sol_diag,runoff_diag,run_off_lic_diag,coeff_rel_diag &237 #endif 238 &)235 ,fq_fonte_diag,fqfonte_diag,snow_evap_diag,fqcalving_diag & 236 ,max_eau_sol_diag,runoff_diag,run_off_lic_diag,coeff_rel_diag & 237 #endif 238 ) 239 239 240 240 USE indice_sol_mod … … 471 471 SUBROUTINE fonte_neige_final(restart_runoff & 472 472 #ifdef ISO 473 &,xtrestart_runoff &473 ,xtrestart_runoff & 474 474 #endif 475 &)475 ) 476 476 ! 477 477 ! This subroutine returns run_off_lic_0 for later writing to restart file. … … 505 505 DO i=1,klon 506 506 IF (iso_verif_egalite_nostop(run_off_lic_0(i) & 507 &,xtrun_off_lic_0(iso_eau,i) &508 &,'fonte_neige 413') &509 &== 1) then507 ,xtrun_off_lic_0(iso_eau,i) & 508 ,'fonte_neige 413') & 509 == 1) then 510 510 WRITE(*,*) 'i=',i 511 511 STOP … … 546 546 fqfonte_out, ffonte_out, run_off_lic_out & 547 547 #ifdef ISO 548 &,fxtcalving_out, fxtfonte_out,xtrun_off_lic_out &548 ,fxtcalving_out, fxtfonte_out,xtrun_off_lic_out & 549 549 #endif 550 &)550 ) 551 551 552 552 … … 626 626 #ifdef ISO 627 627 SUBROUTINE gestion_neige_besoin_varglob_fonte_neige(klon,knon, & 628 &xtprecip_snow,xtprecip_rain, &629 &fxtfonte_neige,fxtcalving, &630 &knindex,nisurf,run_off_lic_diag,coeff_rel_diag)628 xtprecip_snow,xtprecip_rain, & 629 fxtfonte_neige,fxtcalving, & 630 knindex,nisurf,run_off_lic_diag,coeff_rel_diag) 631 631 632 632 ! dans cette routine, on a besoin des variables globales de … … 661 661 j = knindex(i) 662 662 CALL iso_verif_egalite(xtrun_off_lic_0(iso_eau,j), & 663 &run_off_lic_0(j),'gestion_neige_besoin_varglob_fonte_neige 625')663 run_off_lic_0(j),'gestion_neige_besoin_varglob_fonte_neige 625') 664 664 ENDDO 665 665 ENDIF … … 676 676 DO ixt = 1, niso 677 677 xtrun_off_lic(ixt,i) = (coeff_rel_diag * fxtcalving(ixt,i)) & 678 &+(1. - coeff_rel_diag) * xtrun_off_lic_0(ixt,j)678 +(1. - coeff_rel_diag) * xtrun_off_lic_0(ixt,j) 679 679 xtrun_off_lic_0(ixt,j) = xtrun_off_lic(ixt,i) 680 680 xtrun_off_lic(ixt,i) = xtrun_off_lic(ixt,i) + fxtfonte_neige(ixt,i) + xtprecip_rain(ixt,i) … … 683 683 IF (iso_eau > 0) THEN 684 684 IF (iso_verif_egalite_choix_nostop(xtrun_off_lic(iso_eau,i), & 685 &run_off_lic_diag(i),'gestion_neige_besoin_varglob_fonte_neige 1201a', &686 & errmax,errmaxrel) == 1) THEN685 run_off_lic_diag(i),'gestion_neige_besoin_varglob_fonte_neige 1201a', & 686 errmax,errmaxrel) == 1) THEN 687 687 WRITE(*,*) 'i,j=',i,j 688 688 WRITE(*,*) 'coeff_rel_diag=',coeff_rel_diag
Note: See TracChangeset
for help on using the changeset viewer.