Changeset 5160 for LMDZ6/branches/Amaury_dev/libf/phylmd/alpale_wk.F90
- Timestamp:
- Aug 3, 2024, 2:56:58 PM (7 weeks ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/alpale_wk.F90
r5144 r5160 49 49 ALLOCATE (cellrad(klon)) 50 50 ! Compute pseudo grid-cell radius cellrad, such that pi*cellrad^2=cell_area 51 print*, 'alpale_wk: cell_area(1) ', cell_area(1)51 PRINT *, 'alpale_wk: cell_area(1) ', cell_area(1) 52 52 cellrad(:) = sqrt(cell_area(:) / rpi) 53 53 first = .FALSE. … … 55 55 56 56 ! Compute wake radius 57 !! print*,'alpale_wk: sigmaw(1), wdens(1) ', sigmaw(1), wdens(1)57 !! PRINT *,'alpale_wk: sigmaw(1), wdens(1) ', sigmaw(1), wdens(1) 58 58 DO i = 1, klon 59 59 IF (zoccur(i) >= 1) THEN … … 65 65 66 66 ! Compute probability that the grid-cell is intersected by a gust front 67 !! print*,'alpale_wk: wkrad(1), cellrad(1) ', wkrad(1), cellrad(1)67 !! PRINT *,'alpale_wk: wkrad(1), cellrad(1) ', wkrad(1), cellrad(1) 68 68 !! proba_gf(:) = exp(-wdens(:)*rpi*max(wkrad(:)-cellrad(:),0.)**2) - & ! Formules 69 69 !! exp(-wdens(:)*rpi*(wkrad(:)+cellrad(:))**2) ! fausses ! … … 74 74 ! Compute Fip conditionned on the presence of some gust front within the 75 75 ! grid-cell 76 !! print*,'alpale_wk: proba_gf(1), fip(1), ', proba_gf(1), fip(1)76 !! PRINT *,'alpale_wk: proba_gf(1), fip(1), ', proba_gf(1), fip(1) 77 77 fip_cond(:) = fip(:) / proba_gf(:) 78 !! print*,'alpale_wk: wkrad(1), cellrad(1), proba_gf(1), fip(1), fip_cond(1) ', &78 !! PRINT *,'alpale_wk: wkrad(1), cellrad(1), proba_gf(1), fip(1), fip_cond(1) ', & 79 79 !! wkrad(1), cellrad(1), proba_gf(1), fip(1), fip_cond(1) 80 80
Note: See TracChangeset
for help on using the changeset viewer.