Changeset 3940 for LMDZ6/trunk/libf/phylmdiso/fonte_neige_mod.F90
- Timestamp:
- Jun 15, 2021, 1:18:14 PM (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmdiso/fonte_neige_mod.F90
r3927 r3940 28 28 REAL, PRIVATE :: tau_calv 29 29 !$OMP THREADPRIVATE(tau_calv) 30 REAL, ALLOCATABLE, DIMENSION(:,:) , PRIVATE:: ffonte_global30 REAL, ALLOCATABLE, DIMENSION(:,:) :: ffonte_global 31 31 !$OMP THREADPRIVATE(ffonte_global) 32 REAL, ALLOCATABLE, DIMENSION(:,:) , PRIVATE:: fqfonte_global32 REAL, ALLOCATABLE, DIMENSION(:,:) :: fqfonte_global 33 33 !$OMP THREADPRIVATE(fqfonte_global) 34 REAL, ALLOCATABLE, DIMENSION(:,:) , PRIVATE:: fqcalving_global34 REAL, ALLOCATABLE, DIMENSION(:,:) :: fqcalving_global 35 35 !$OMP THREADPRIVATE(fqcalving_global) 36 REAL, ALLOCATABLE, DIMENSION(:) , PRIVATE:: runofflic_global36 REAL, ALLOCATABLE, DIMENSION(:) :: runofflic_global 37 37 !$OMP THREADPRIVATE(runofflic_global) 38 38 #ifdef ISO … … 60 60 ! The variable run_off_lic_0 is initialized to the field read from 61 61 ! restart file. The other variables are initialized to zero. 62 63 62 ! 64 63 !**************************************************************************************** … … 84 83 run_off_lic_0(:) = restart_runoff(:) 85 84 86 87 85 !**************************************************************************************** 88 86 ! Allocate other variables and initilize to zero … … 130 128 ENDIF 131 129 runofflic_global(:) = 0.0 132 133 134 130 135 131 !**************************************************************************************** … … 348 344 349 345 snow_evap = 0. 350 346 351 347 #ifdef ISOVERIF 352 348 write(*,*) 'klon,snow_evap(413)=',klon,snow_evap(413) … … 406 402 ! Y'a-t-il fonte de neige? 407 403 neige_fond = (snow(i)>epsfra .OR. nisurf==is_sic .OR. nisurf==is_lic) .AND. tsurf_new(i)>=RTT 408 IF (neige_fond) THEN 404 IF (neige_fond) THEN 409 405 fq_fonte = MIN( MAX((tsurf_new(i)-RTT )/chasno,0.0),snow(i)) 410 406 ffonte(i) = fq_fonte * RLMLT/dtime … … 412 408 snow(i) = MAX(0., snow(i) - fq_fonte) 413 409 bil_eau_s(i) = bil_eau_s(i) + fq_fonte 414 tsurf_new(i) = tsurf_new(i) - fq_fonte * chasno 410 tsurf_new(i) = tsurf_new(i) - fq_fonte * chasno 415 411 #ifdef ISO 416 412 fq_fonte_diag(i)=fq_fonte
Note: See TracChangeset
for help on using the changeset viewer.