Ignore:
Timestamp:
Aug 3, 2024, 2:56:58 PM (7 weeks ago)
Author:
abarral
Message:

Put .h into modules

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/branches/Amaury_dev/libf/phylmd/alpale_wk.F90

    r5144 r5160  
    4949    ALLOCATE (cellrad(klon))
    5050    !  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)
    5252    cellrad(:) = sqrt(cell_area(:) / rpi)
    5353    first = .FALSE.
     
    5555
    5656  !  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)
    5858  DO i = 1, klon
    5959    IF (zoccur(i) >= 1) THEN
     
    6565
    6666  !  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)
    6868  !!  proba_gf(:) = exp(-wdens(:)*rpi*max(wkrad(:)-cellrad(:),0.)**2) - &   ! Formules
    6969  !!                exp(-wdens(:)*rpi*(wkrad(:)+cellrad(:))**2)             ! fausses !
     
    7474  !  Compute Fip conditionned on the presence of some gust front within the
    7575  !  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)
    7777  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) ', &
    7979  !!                        wkrad(1), cellrad(1), proba_gf(1), fip(1), fip_cond(1)
    8080
Note: See TracChangeset for help on using the changeset viewer.